Deconfinement at high temperatures and moderately high baryon densities Péter Petreczky What is the limiting temperature on hadronic matter? What is the nature of the deconfined matter? In this talk: Chiral transition in QCD at T>0 Color screening deconfinement Equation of state Fluctuations of conserved charges How does the transition and EoS changes with baryon density LQCD INT Program INT-16-2b, Phases of Dense Matter, Seattle, July 14, 2016
The temperature dependence of chiral condensate Renormalized chiral condensate introduced by Budapest-Wuppertal collaboration With choice : Bazavov et al (HotQCD), Phys. Rev. D85 (2012) 054503; Bazavov, PP, PRD 87(2013)094505, Borsányi et al, JHEP 1009 (2010) 073 Calculations with chiral (Domain Wall) fermions: disc = V T h( ) 2 i h i 2 Bhattacharya et al (HotQCD), PRL 113 (2014)082001 0.025 0.02 0.015 0.01 0.005 l R HISQ, m l =m s /20 HISQ, m l =m s /27 stout, m l =m s /27 χ MS l,disc. /T 2 40 30 20 48 3 12, HISQ,m l =0.05 m s,m π =161MeV 16 3 8,m π =200MeV 24 3 8,m π =200MeV 32 3 8,m π =200MeV 32 3 8,m π =140MeV 64 3 8,m π =140MeV 0 10-0.005 120 130 140 150 160 170 180 190 200 210 O(4) scaling analysis and continuum limit: 0 120 130 140 150 160 170 180 190 200 210 220 230 T (MeV) Transition is a crossover (no volume dep.) T c = (155 ± 8 ± 1)MeV
Deconfinement and color screening Onset of color screening is described by Polyakov loop (order parameter in SU(N) gauge theory) exp( F Q Q(r, T )/T )= 1 9 htrl(r)trl (0)i F Q Q(r!1,T)=2F Q (T ) 2+1 flavor QCD, continuum extrapolated (work in progress with Bazavov, Weber ) 0.9 L ren (T) T QCD c T PG c free energy of static quark anti-quark pair shows Debye screening SU(N) gauge theory QCD! at high temperatures Similar results with stout action Borsányi et al, JHEP04(2015) 138 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 stout, cont. HISQ, cont. SU(3) SU(2) 150 200 250 300 350 400 450 500
The entropy of static quark 600 500 400 300 200 Bazavov, PP, PRD 87 (2013) 094505 Hadron gas (of static-light) vs. LQCD F Q (T) [MeV] HISQ stout 5 4.5 4 3.5 3 2.5 2 1.5 Bazavov et al, PRD 93 (2016) 114502 S Q S Q = T (N ) @F Q @T local fit global 1/N 4 fit global 1/N 2 fit HRG N =8 100 120 140 160 180 200 1 120 140 160 180 200 220 240 260 280 300 At low T the entropy S Q increases reflecting the increase of states the heavy quark can be coupled to; at high temperature the static quark only sees the medium within a Debye radius, as T increases the Debye radius decreases and S Q also decreases The onset of screening corresponds to peak is S Q and its position coincides with T c
Equation of state in the continuum limit Equation of state has been calculated in the continuum limit up to T=400 MeV using two different discretizations and the results agree well 16 12 8 4 Bazavov et al (HotQCD), PRD 90 (2014) 094503 T c HRG 3p/T 4 /T 4 3s/4T 3 non-int. limit 0 130 170 210 250 290 330 370 Hadron resonance gas (HRG): Interacting gas of hadrons = non-interacting gas of hadrons and hadron resonances ç virial expansion Prakash,Venugopalan, NPA546 (1992) 718 HRG agrees with the lattice for T< 155 MeV 0.34 0.32 0.3 0.28 0.26 0.24 0.22 0.2 0.18 0.16 0.14 0.12 Borsányi et al, PLB 730 (2014) 99 c 2 s T c = (154 ± 9)MeV c ' 300MeV/fm 3 low ' 180MeV/fm 3 HISQ/tree stout HRG 150 200 250 300 350 400 nucl ' 150MeV/fm 3 high ' 500MeV/fm 3 proton ' 450MeV/fm 3
QCD thermodynamics at non-zero chemical potential Taylor expansion : p(t,µ B,µ Q,µ S,µ C ) T 4 = X ijkl 1 i!j!k!l! BQSC ijkl µb i µb j µs k µc T T T T l hadronic p(t,µ u,µ d,µ s,µ c ) T 4 = X ijkl 1 i!j!k!l! udsc ijkl µu i µd j µs k µc T T T T l quark abcd ijkl = T i+j+k+l @ i @µ i b @ j @µ j b @ k @µ i c @ l @µ l d ln Z(T,V,µ a,µ b,µ c,µ d ) µa =µ b =µ c =µ d =0 Taylor expansion coefficients give the fluctuations and correlations of conserved charges, e.g. information about carriers of the conserved charges ( hadrons or quarks ) probes of deconfinement
Deconfinement : fluctuations of conserved charges 1 0.8 B = 1 VT 3 hb2 i hbi 2 Q = 1 VT 3 hq2 i hqi 2 S = 1 VT 3 hs2 i hsi 2 i / i SB baryon number electric charge strangeness Ideal gas of massless quarks : 0.6 0.4 0.2 0 i=b Q S filled : HISQ, N =6, 8 open : stout continuum 150 200 250 300 350 SB S =1 conserved charges carried by light quarks Bazavov et al (HotQCD) PRD86 (2012) 034509 Borsányi et al, JHEP 1201 (2012) 138 conserved charges are carried by massive hadrons
Quark number fluctuations at high T At high temperatures quark number fluctuations can be described by weak coupling approach due to asymptotic freedom of QCD 0.9 0.85 u 4 / ideal 4 quark number fluctuations 0-0.002 ud 11 quark number correlations 0.8 0.75 0.7 0.65 0.6 EQCD 4 u (cont) 3-loop HTL 1.05 0.95 0.55 300 350 400 450 500 550 600 650 700 1 0.9 u 2 / ideal 2 0.85 300 500 700 900 EQCD Lattice results converge as the continuum limit is approached Good agreement between lattice and the weak coupling approach for 2 nd and 4 th order quark number fluctuations as well as for correlations -0.004-0.006-0.008-0.01-0.012 Bazavov et al, PRD88 (2013) 094021, Ding et at, PRD92 (2015) 074043 EQCD N =6 8 10 12 cont -0.014 250 300 350 400 450 500 550 600 650 700
140 180 220 260 300 340 Deconfinement of strangeness Partial pressure of strange hadrons in uncorrelated hadron gas: Strange hadrons are heavy è treat them as Boltzmann gas should vanish! 0.30 non-int. quarks v 1 and v 2 do vanish within errors at low T v 1 and v 2 rapidly increase above the transition region, eventually reaching non-interacting quark gas values Bazavov et al, PRL 111 (2013) 082301 0.25 0.20 0.15 0.10 0.05 0.00 uncorr. hadrons χ 2 B -χ4 B v 1 v 2
Deconfinement of strangeness (cont d) Using the six Taylor expansion coefficients related to strangeness it is possible to construct combinations that give up to terms 1.0 0.5 0.0 0.10 0.05 B 3 (0,0) B 3 (0,1/6) B 3 (1/3,1/6) HRG P S =3,B -0.5-1.0-1.5-2.0 M(0,0) M(0,-3) M(-6,-3) HRG P S =1,M 140 180 220 260 300 340 0.00-0.05 140 180 220 260 300 340 Hadron resonance gas descriptions breaks down for all strangeness sectors above T c Bazavov et al, PRL 111 (2013) 082301
Deconfinement of charm XY C nml = T m+n+l @n+m+l p(t,µ X,µ Y,µ C )/T 4 @µ n X @µm Y @µl C m c T only C =1 sector contributes 3.0 2.5 2.0 In the hadronic phase all BC-correlations are the same! BC BC 13 / 22 BC BC 11 / 13 N : 8 6 non-int. quarks 0.5 0.4 0.3 0.2 0.4 Bazavov et al, PLB 737 (2014) 210 Charm baryon to meson pressure BC C BC 13 /( 4-13 ) BQC QC BQC 0.5 112 /( 13-112 ) QM-HRG-3 QM-HRG PDG-HRG non-int. quarks 1.5 0.3 1.0 un-corr. hadrons 140 160 180 200 220 240 260 280 0.7 0.5 BSC SC BSC - 112 /( 13-112 ) N : 8 6 Hadronic description breaks down just above T 0.3 c open charn deconfines above T c 140 150 160 170 180 190 200 210 The charm baryon spectrum is not well known (only few states in PDG), HRG works only if the missing states are included
0.8 0.6 Quasi-particle model for charm degrees of freedom Charm dof are good quasi-particles at all T because M c >>T and Boltzmann approximation holds p C (T,µ B,µ c )=p C q (T ) cosh(ˆµ C +ˆµ B /3) + p C B(T ) cosh(ˆµ C +ˆµ B )+p C M (T ) cosh(ˆµ C ) C 2, BC 13, BC 22 ) p C q (T ),p C M (T ),p C B(T ) Partial meson and baryon pressures described by HRG at T c then drop gradually, charm quark only dominant dof at T>200 MeV or ε > 6 GeV/fm 3 1.0 p q C /p C p B C /p C ˆµ X = µ X /T and dominate the charm pressure Partial pressures drop because hadronic excitations become broad at high temperatures (bound state peaks merge with the continuum) 0.4 0.2 p M C /p C Mukherjee, PP, Sharma, PRD 93 (2016) 014502 See Jakovác, PRD88 (2013), 065012 Biró, Jakovác, PRD(2014)065012 0.0 150 170 190 210 230 250 270 290 310 330 0.3 Vice versa for quarks
Equation of state at non-zero baryon density At the highest energy heavy ion collisions create a system with almost zero net baryon density. by decreasing the the center of mass energy of the collisions E CM = p s the net baryon number can be increased and the EoS calculated in LQCD using Taylor expansion Taylor expansion up to 4 th order for net zero strangeness n S =0 and r = n Q /n B = Z/A =0.4 Moderate effects due to non-zero baryon density up to µ B /T =2$ p s 20GeV and Taylor expansion works well BNL-Bielefeld-CCNU
Freeze-out and transition at non-zero baryon density T c (µ B ) T c (0) =1 apple µb T c 2 +... apple =0.020(4) Cea et al, PRD 93 (2016) 014507 apple =0.0135(20) Bonati et al, PRD 92 (2015) 054503 apple =0.0149(21) Bellweid et al, PLB 751 (2015) 559 apple =0.0066(7) Kaczmarek et al, PRD 83 (2011) 014504 0.07 n S = 0, n Q = 0.4n B O(µ B ) 0.06 0.05 0.04 N τ = 8 O(µ B 3 ) O(µ B 5 ) n B [fm -3 ] 0.03 0.02 s 1/2 NN [GeV] BNL-Bielefeld-CCNU 15 20 25 30 35 40 45 50 55
Summary The deconfinement transition temperature coincides with the chiral transition temperature T c = 154(9) MeV Equation of state are known in the continuum limit up to T=400 MeV at zero baryon density and the energy density at the transition is small 300 MeV/fm 3 Hadron resonance gas (HRG) can describe various thermodynamic quantities at low temperatures Deconfinement transition can be studied in terms of fluctuations and correlations of conserved charges; deconfinement means an abrupt breakdown of hadronic description that occurs at T c and appearance of quark like excitations Hadronic excitations containing charm quark can exist above T c For T > (300-400) MeV weak coupling expansion may work for certain quantities (e.g. quark number susceptibilities) The region of moderately large baryon densities relevant for RHIC can be accessed in LQCD using Taylor expansion. The effects of the baryon densities are not large and consistent with HRG (no hint for critical point?); Much more work is needed to obtain reliable results, e.g. T c (µ B ), energy density at T c.
How Equation of state changed since 2002? 16 14 12 10 8 6 4 /T 4 ideal gas m q /T=0.4, N =4 (2000) HISQ (2014) 2 0 100 150 200 250 300 350 400 450 Much smoother transition to QGP The energy density keeps increasing up to 450 MeV instead of flattening
Entropy of static quark at high temperature S Q (T) 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 LO NLO, µ=(1-4) T NNLO lattice -0.1 1000 1500 2000 2500 3000 3500 4000 4500 5000 M. Berwein, N. Brambilla, P. Petreczky and A. Vairo, arxiv:1512.08443 [hep-ph]
Does the quasi-particle model makes sense? 4 non-trivial constraints on the model provided by : BC 31, SC 31, BSC 121, BSC 211 Diquark pressure is zero! 1.0 0.5 c 1 /p C c 2 /p C c 3 /p C c 4 /p C 0.0 0.3 p B C,S=2 /p C 0.0-0.5-1.0 150 170 190 210 230 250 270 290 310 330 0.2 p C,S=1 B /p C p C,S=1 M /p C 0.1 0.0 150 170 190 210 230 250 270 290 310 330