Bottomonium suppression in the quark-gluon plasma
Michael Strickland
Kent State University Kent, OH USA
Sharif University of Technology July 14, 2020
Bottomonium suppression in the quark-gluon plasma Michael - - PowerPoint PPT Presentation
Bottomonium suppression in the quark-gluon plasma Michael Strickland Kent State University Kent, OH USA Sharif University of Technology July 14, 2020 Quarks are normally confined inside hadrons HADRONS Baryons Mesons Gluons
Kent State University Kent, OH USA
Sharif University of Technology July 14, 2020
Baryons Mesons HADRONS
Gluons hold the quarks together
… …
2
Name Mass [GeV/c2] Electric Charge Up u 0.0024 +2/3 Down d 0.0048
Strange s 0.104
Charm c 1.27 2/3 Bottom b 4.2
Top t 171.2 2/3
have electric charge and “color charge”
electric charge and “anti-color charge”
proton which is ~ 1 GeV
Diagram Symbol
quark antiquark flow of time q q
Blue Red Green Anti-Red Anti-Blue Anti-Green “Colorless”
3
4 Color ionized plasma
5
Crossover First Order Second Order
200 MeV à 2 x 1012 K
Tc ~ 154 MeV
T [MeV]
Color Superconductor
100 200 300 1000 Hadron Gas
Nuclear Matter
Quark Gluon Plasma
Quark Matter
Vacuum
μ [MeV]
B
300 600 900 1200 1500 Quark-Gluon Plasma
1 loop αs ;
MS
176 MeV μB 0 MeV
ideal
Wuppertal Budapest NNLO HTLpt
6
No fit parameters!
Andersen, Leganger, Su, and MS 1009.4644, 1103.2528
Hadron Resonance Gas Quark Gluon Plasma
RHIC 200 GeV LHC 5.023 TeV LHC 2.76 TeV
Resummed perturbation theory (Strickland et al) Lattice QCD (Wuppertal-Budapest group)
7
Crossover First Order Second Order
200 MeV à 2 x 1012 K
Tc ~ 154 MeV
T [MeV]
Color Superconductor
100 200 300 1000 Hadron Gas
Nuclear Matter
Quark Gluon Plasma
Quark Matter
Vacuum
μ [MeV]
B
300 600 900 1200 1500
neutron stars supernovae
LHC 5.023 TeV à T0 ~ 700 MeV LHC 2.76 TeV à T0 ~ 600 MeV RHIC 200 GeV à T0 ~ 400 MeV
fraction of the mass of observed hadronic matter.
from the strong force.
8
Bashir et al, Commun. Theor. Phys. 58 (2012) 79-134.
Logarithmic scale
Figure courtesy B. Muller
9
10
200 MeV à 2 x 1012 K
~ 10-14 m Entire event lasts ~ 10 fm/c which is ~ 3 x 10-23 s !!!
T = 200 MeV à 2 x 1012 K
~ 10-14 m Entire event lasts ~ 10 fm/c which is ~ 3 x 10-23 s !!!
11
200 MeV à 2 x 1012 K
~ 10-14 m Entire event lasts ~ 10 fm/c which is ~ 3 x 10-23 s !!!
T = 200 MeV à 2 x 1012 K
~ 10-14 m Entire event lasts ~ 10 fm/c which is ~ 3 x 10-23 s !!!
and allows us to extract transport coefficients like the shear viscosity
propagation; provides information about momentum diffusion and energy loss of partons in the QGP
bound state survival in the QGP
information about initial conditions
final particle distributions are thermalized
provides information about the size of the QGP and collective flow profiles
12
and allows us to extract transport coefficients like the shear viscosity
propagation; provides information about momentum diffusion and energy loss of partons in the QGP
bound state survival in the QGP
information about initial conditions
final particle distributions are thermalized
provides information about the size of the QGP and collective flow profiles
13
14
15
E288 exp @ Fermilab, 1977
16
Vacuum decay lifetime = 3654 fm/c ~ 10-20 s Vacuum bindng energy of Y(1s) is ~ 1 GeV
17 Color ionized plasma
18
and MS, 1302.2180
states one needs a fully relativistic treatment
sufficiently high then one can take the “heavy quark limit”
having to deal with a non- relativistic terms plus relativistic corrections
19
Name Mass [GeV/c2] Electric Charge Up u 0.0024 +2/3 Down d 0.0048
Strange s 0.104
Charm c 1.27 2/3 Bottom b 4.2
Top t 171.2 2/3
20
shows, it works quite well
masses of the bottomonium sates is 0.22%
21
Screening of electric interaction with screening length rD = 1/mD
A test charge polarizes the particles of the plasma and they “screen” its charge
Peter Debye 1884 - 1966
electric plasma occurs in the QGP
strong interactions of quarks and gluons
22
Coloumb Debye-screened
VColoumb(r) = −α r − → VDebye(r) = −α r e−mDr
23
Coloumb Debye-screened
VColoumb(r) = −α r − → VDebye(r) = −α r e−mDr
24
Coloumb Debye-screened
VColoumb(r) = −α r − → VDebye(r) = −α r e−mDr
25
Anisotropic potential calculation: Dumitru, Guo, and MS, 0711.4722 and 0903.4703 Gluon propagator in an anisotropic plasma: Romatschke and MS, hep-ph/0304092
Using the real-time formalism one can express the potential in terms of the static advanced, retarded, and Feynman propagators Real part can be written as With direction-dependent masses, e.g.
V (r, ξ) = −g2CF Z d3p (2π)3 (eip·r − 1)1 2 ⇣ D∗L
R + D∗L A + D∗L F
⌘
Re[V (r, ξ)] = −g2CF Z d3p (2π)3 eip·r p2 + m2
α + m2 γ
(p2 + m2
α + m2 γ)(p2 + m2 β) − m4 δ
26
be parameterized as a Debye-screened potential with a direction-dependent Debye mass
imaginary part coming from the Landau damping
exists in the isotropic case
Laine et al hep-ph/0611300
the free energy (F) and also obtained internal energy (U) from this.
Dumitru, Guo, Mocsy, and MS, 0901.1998 MS, 1106.2571; Bazow and MS, 1112.2761 Dumitru, Guo, and MS, 0711.4722 and 0903.4703 Burnier, Laine, Vepsalainen, arXiv:0903.3467 (aniso)
27
Internal Energy
28 cartoon
expect weaker color binding (Debye screening + asymptotic freedom)
the bound states cause them to have shorter lifetimes à larger spectral widths
Υ(1𝑡) Υ(2𝑡) Υ(3𝑡)
29 cartoon
expect weaker color binding (Debye screening + asymptotic freedom)
the bound states cause them to have shorter lifetimes à larger spectral widths
Υ(1𝑡) Υ(2𝑡) Υ(3𝑡)
30 cartoon
expect weaker color binding (Debye screening + asymptotic freedom)
the bound states cause them to have shorter lifetimes à larger spectral widths
Υ(1𝑡) Υ(2𝑡) Υ(3𝑡)
31 cartoon
expect weaker color binding (Debye screening + asymptotic freedom)
the bound states cause them to have shorter lifetimes à larger spectral widths
Υ(1𝑡) Υ(2𝑡) Υ(3𝑡)
32 cartoon
expect weaker color binding (Debye screening + asymptotic freedom)
the bound states cause them to have shorter lifetimes à larger spectral widths
Υ(1𝑡) Υ(2𝑡) Υ(3𝑡)
The CMS (Compact Muon Solenoid) experiment has measured bottomonium spectra for both pp and Pb-Pb collisions. With this we can extract RAA experimentally.
33
pp PbPb
Υ(1𝑡) Υ(2𝑡) Υ(3𝑡)
34
production anisotropies
data explained well by KSU in-medium breakup model
35
For in-medium suppression, given the population of quarkonia states at some t0, we can simply integrate the instantaneous decay/regeneration rate of the state G(t,x,y,h)
36
1 fm/c = 3 x 10-24 seconds
Pb-Pb @ 2.76 TeV T0 = 600 MeV t0 = 0.25 fm/c b = 7 fm
b = 7 fm
Solve the 3d Schrödinger EQ with a complex-valued potential
37
Fold together with the non-EQ spatiotemporal evolution to
Obtain the real and imaginary parts of the binding energies for the U(1s), U(2s), U(3s), cb1, and cb2 as function of energy density and momentum-anisotropy.
Yager-Elorriaga and MS, 0901.1998; Margotta, MS, et al, 1101.4651
Number of nucleon-nucleon collisions per nucleus collision Number produced in a proton-proton collision Number produced in a nucleus-nucleus collision
|| < < < πη/ =
χ χ Υ() Υ() Υ()
39
suppression
(QGP thermometer is continuous!)
U(1s) cb(1p) cb(2p) U(2s) U(3s)
40
|| < < < Υ() Υ() Υ() Υ() πη/ = πη/ = πη/ =
41
|| <
Υ() Υ() Υ() Υ() πη/ = πη/ = πη/ =
< < Υ() Υ() Υ() Υ() Υ() πη/ = πη/ = πη/ =
from CMS and ALICE
data but not perfect
rapidity range 2.5 < y < 4
Y(2s) suppression
|| < < < Υ() Υ()
Υ() Υ() πη/ = πη/ = πη/ =
42
Cold Nuclear Matter Effects
43
2.5 < y < 4.0 0 < pT < 15 GeV
Υ() πη/ = πη/ = πη/ =
|y| < 2.4 0 < pT < 30 GeV Υ(1S) Υ(2S)
Υ() Υ() πη/ = πη/ = πη/ =
covers forward rapidity (2.5 < |y| < 4)
44
4d flow tomography Quark Matter 2019
Cern Courier, Fall 2019
Bhadhuri, Alqahtani, Borghini, Jaiswal, and MS, 2007.03939
45
46
3s, 3p) by solving 3D Schrödinger equation numerically for each trajectory.
constructed from PDG branching
sections.
47
HQQD - Υ(1s) HQQD - Υ(2s) HQQD - Υ(3s)
ALICE - Y(1s) ATLAS - Y(1s) CMS - Y(1s) ATLAS - Y(2s) CMS - Y(2s)
100 200 300 400 0.0 0.2 0.4 0.6 0.8 1.0 Npart RAA
5.02 TeV Pb-Pb ALICE: pT < 15 GeV and 2.5 < y < 4 ATLAS: pT < 15 GeV and |y| < 1.5 CMS: pT < 30 GeV and |y| < 2.4 HQQD - Υ(1s) HQQD - Υ(2s) HQQD - Υ(3s) ALICE - Y(1s) ATLAS - Y(1s) CMS - Y(1s) ATLAS - Y(2s) CMS - Y(2s)
5 10 15 20 25 30 0.0 0.2 0.4 0.6 0.8 1.0 pT RAA
5.02 TeV Pb-Pb CMS: 0-100% and |y| < 2.4 ALICE: 0-100% and 2.5 < y < 4 ATLAS: 0-80% and |y| < 1.5
48
HQQD - Υ(1s) HQQD - Υ(2s) HQQD - Υ(3s)
10 20 30 40 50 60 70 80 90 100 0.00 0.02 0.04 0.06 0.08 Centrality (%) v2[Υ]
5.02 TeV Pb-Pb y=0 pT < 50 GeV
49
CMS 2020 ALICE 2019 HQQD
5 10 15
0.00 0.05 0.10 pT [GeV] v2[Υ(1s)]
5.02 TeV Pb-Pb 5-60% Centrality CMS: |y| < 2.4 ALICE: 2.5 < y < 4
50
CMS 2020 ALICE 2019 HQQD
5 10 15
0.00 0.05 0.10 pT [GeV] v2[Υ(1s)]
5.02 TeV Pb-Pb 5-60% Centrality CMS: |y| < 2.4 ALICE: 2.5 < y < 4
51
52
explain the experimental observations.
suppression seen at LHC à QGP!
today was “regeneration”. This occurs when the density of heavy (anti-)quarks becomes large, making it probable for a pair to recombine in the QGP. At very high temperatures/beam energies this effect is important for charm quarks, but still not so important for bottom quarks.
real-time in-medium Schrödinger equation evolution (student A. Islam) à in-medium quantum regeneration
53