Quark Gluon Plasma. Rajiv V. Gavai T. I. F. R., Mumbai. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V.

Transcription

1 Quark Gluon Plasma Rajiv V. Gavai T. I. F. R., Mumbai Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 1

2 Quark Gluon Plasma Rajiv V. Gavai T. I. F. R., Mumbai Introduction QGP in Bulk J/ψ suppression Jet Quenching QCD Phase Diagram Summary Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 1

3 Introduction Measurement Fit O meas O fit /σ meas α (5) had (m Z ) ± m Z [GeV] ± Γ Z [GeV] ± σ 0 had [nb] ± R l ± A 0,l fb ± A l (P τ ) ± R b ± R c ± A 0,b fb ± A 0,c fb ± A b ± A c ± A l (SLD) ± sin 2 θ lept eff (Q fb ) ± m W [GeV] ± Γ W [GeV] ± m t [GeV] ± Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 2

4 Measurement Fit O meas O fit /σ meas α (5) had (m Z ) ± m Z [GeV] ± Γ Z [GeV] ± σ 0 had [nb] ± R l ± A 0,l fb ± A l (P τ ) ± R b ± R c ± A 0,b fb ± A 0,c fb ± A b ± A c ± A l (SLD) ± sin 2 θ lept eff (Q fb ) ± m W [GeV] ± Γ W [GeV] ± m t [GeV] ± Introduction Precision tests from LEP Standard Model Very Successful! Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 2

5 Measurement Fit O meas O fit /σ meas α (5) had (m Z ) ± m Z [GeV] ± Γ Z [GeV] ± σ 0 had [nb] ± R l ± A 0,l fb ± A l (P τ ) ± R b ± R c ± A 0,b fb ± A 0,c fb ± A b ± A c ± A l (SLD) ± sin 2 θ lept eff (Q fb ) ± m W [GeV] ± Γ W [GeV] ± m t [GeV] ± Introduction Precision tests from LEP Standard Model Very Successful! All tests based on perturbation theory Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 2

6 Measurement Fit O meas O fit /σ meas α (5) had (m Z ) ± m Z [GeV] ± Γ Z [GeV] ± σ 0 had [nb] ± R l ± A 0,l fb ± A l (P τ ) ± R b ± R c ± A 0,b fb ± A 0,c fb ± A b ± A c ± A l (SLD) ± sin 2 θ lept eff (Q fb ) ± m W [GeV] ± Γ W [GeV] ± m t [GeV] ± Introduction Precision tests from LEP Standard Model Very Successful! All tests based on perturbation theory Need to understand non-perturbative QCD, e.g., to explain baryonic mass in our Universe, i.e., us. Also, total cross sections, matrix elements, structure functions... Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 2

7 Measurement Fit O meas O fit /σ meas α (5) had (m Z ) ± m Z [GeV] ± Γ Z [GeV] ± σ 0 had [nb] ± R l ± A 0,l fb ± A l (P τ ) ± R b ± R c ± A 0,b fb ± A 0,c fb ± A b ± A c ± A l (SLD) ± sin 2 θ lept eff (Q fb ) ± m W [GeV] ± Γ W [GeV] ± m t [GeV] ± Introduction Precision tests from LEP Standard Model Very Successful! All tests based on perturbation theory Need to understand non-perturbative QCD, e.g., to explain baryonic mass in our Universe, i.e., us. Also, total cross sections, matrix elements, structure functions... Lattice QCD only well-understood, viable tool for this: Iconic Results - Wilczek. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 2

8 Measurement Fit O meas O fit /σ meas α (5) had (m Z ) ± m Z [GeV] ± Γ Z [GeV] ± σ 0 had [nb] ± R l ± A 0,l fb ± A l (P τ ) ± R b ± R c ± A 0,b fb ± A 0,c fb ± A b ± A c ± A l (SLD) ± sin 2 θ lept eff (Q fb ) ± m W [GeV] ± Γ W [GeV] ± m t [GeV] ± Introduction Precision tests from LEP Standard Model Very Successful! All tests based on perturbation theory Need to understand non-perturbative QCD, e.g., to explain baryonic mass in our Universe, i.e., us. Also, total cross sections, matrix elements, structure functions... Lattice QCD only well-understood, viable tool for this: Iconic Results - Wilczek It predicts a transition to Quark-Gluon Plasma (QGP), and QGP properties. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 2

9 Quest for Quark-Gluon Plasma : Heavy Ion Collisions at SPS, RHIC and LHC. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 3

10 Quest for Quark-Gluon Plasma : Heavy Ion Collisions at SPS, RHIC and LHC. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 3

11 Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 4

12 Another fundamental aspect Critical Point in T -µ B plane; based on symmetries and models. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 5

13 Another fundamental aspect Critical Point in T -µ B plane; based on symmetries and models. Expected QCD Phase Diagram Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 5

14 Another fundamental aspect Critical Point in T -µ B plane; based on symmetries and models. Expected QCD Phase Diagram... but could, however, be... Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 5

15 Another fundamental aspect Critical Point in T -µ B plane; based on symmetries and models. Expected QCD Phase Diagram... but could, however, be... Τ Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 5

16 Transition temperature, Critical energy density, Order of Phase Transition? What are the properties of QGP (EoS, Excitations, Screening..)? Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 6

17 Transition temperature, Critical energy density, Order of Phase Transition? What are the properties of QGP (EoS, Excitations, Screening..)? Need N s N t for thermodynamic limit and large N t for continuum limit. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 6

18 Transition temperature, Critical energy density, Order of Phase Transition? What are the properties of QGP (EoS, Excitations, Screening..)? Need N s N t for thermodynamic limit and large N t for continuum limit. Completely parameter-free : Λ QCD and quark masses from hadron spectrum. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 6

19 The Transition Temperature T c 170 MeV (Bielefeld & CP-PACS, Tsukuba 2001) and Equation of State (EOS), have been predicted by lattice QCD. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 7

20 The Transition Temperature T c 170 MeV (Bielefeld & CP-PACS, Tsukuba 2001) and Equation of State (EOS), have been predicted by lattice QCD. OLD results on order parameters & EoS from Bielefeld 2001 (Coarse N t =4 lattices). Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 7

21 The Transition Temperature T c 170 MeV (Bielefeld & CP-PACS, Tsukuba 2001) and Equation of State (EOS), have been predicted by lattice QCD. OLD results on order parameters & EoS from Bielefeld 2001 (Coarse N t =4 lattices). Other quantities, notably strangeness enhancement in Heavy Ion Physics, the Wróblewski Parameter λ s (RVG & Sourendu Gupta PR D 2002) have also been predicted by lattice QCD. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 7

22 Thrust of new results now on continuum limit, lighter quarks, 2+1 flavours transition temperature T c = 192(7)(4)MeV (Bielefeld-RBC hep-lat/ ). = A gap between freeze-out and transition temperatures? Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 8

23 Thrust of new results now on continuum limit, lighter quarks, 2+1 flavours transition temperature T c = 192(7)(4)MeV (Bielefeld-RBC hep-lat/ ). = A gap between freeze-out and transition temperatures? more complex observables, speed of sound, transport coefficients, Fluctuations, J/ψ-dissolution/persistence, dileptons... etc. and T -µ phase diagram. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 8

24 Thrust of new results now on continuum limit, lighter quarks, 2+1 flavours transition temperature T c = 192(7)(4)MeV (Bielefeld-RBC hep-lat/ ). = A gap between freeze-out and transition temperatures? more complex observables, speed of sound, transport coefficients, Fluctuations, J/ψ-dissolution/persistence, dileptons... etc. and T -µ phase diagram. An interesting theoretical issue Conformal Invariance and AdS/CFT predictions. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 8

25 Thrust of new results now on continuum limit, lighter quarks, 2+1 flavours transition temperature T c = 192(7)(4)MeV (Bielefeld-RBC hep-lat/ ). = A gap between freeze-out and transition temperatures? more complex observables, speed of sound, transport coefficients, Fluctuations, J/ψ-dissolution/persistence, dileptons... etc. and T -µ phase diagram. An interesting theoretical issue Conformal Invariance and AdS/CFT predictions. Lot of activity in Model Building to explain Lattice QCD results: Quasi-particle models, Hadron Resonance Gas, Quarkonia from Lattice Q Q potential, sqgp and coloured states... Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 8

26 QGP in Bulk : EoS, Speed of Sound... Recent results for EoS : N t =6, Smaller quark masses. Bernard et al., MILC hep-lat/ ; Aoki et al., hep-lat/ Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 9

27 Small differences; ɛ(t c ) 6Tc 4 still. Too small volumes = Thermodynamic Limit? Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 10

28 Small differences; ɛ(t c ) 6Tc 4 still. Too small volumes = Thermodynamic Limit? C s Crucial for elliptic flow, hydrodynamical studies... C v Event-by-event temperature/p T fluctuations. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 10

29 Small differences; ɛ(t c ) 6Tc 4 still. Too small volumes = Thermodynamic Limit? C s Crucial for elliptic flow, hydrodynamical studies... C v Event-by-event temperature/p T fluctuations. Can be obtained from ln Z by taking appropriate derivatives which relate it to the temperature derivative of anomaly measure /ɛ. (RVG, S. Gupta and S. Mukherjee, PR D71 (2005) ) New method to obtain these differentially without getting negative pressure. Introduced an improved operator than used in earlier Bielefeld studies. (RVG, S. Gupta and S. Mukherjee, hep-lat/ ) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 10

30 Small differences; ɛ(t c ) 6Tc 4 still. Too small volumes = Thermodynamic Limit? C s Crucial for elliptic flow, hydrodynamical studies... C v Event-by-event temperature/p T fluctuations. Can be obtained from ln Z by taking appropriate derivatives which relate it to the temperature derivative of anomaly measure /ɛ. (RVG, S. Gupta and S. Mukherjee, PR D71 (2005) ) New method to obtain these differentially without getting negative pressure. Introduced an improved operator than used in earlier Bielefeld studies. (RVG, S. Gupta and S. Mukherjee, hep-lat/ ) Using lattices with 8, 10, and 12 temporal sites ( and 38 4 lattices) and with statistics of million iterations, ɛ, P, s, C 2 s and C v obtained in continuum. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 10

31 Ideal Gas Ideal Gas C v /T 3 15 C s ( f ) ( e ) T/T c T/T c (RVG, S. Gupta and S. Mukherjee, hep-lat/ ) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 11

32 Ideal Gas Ideal Gas C v /T 3 15 C s ( f ) ( e ) T/T c T/T c (RVG, S. Gupta and S. Mukherjee, hep-lat/ ) C v 4ɛ for 2T c but No Ideal Gas limit. Specific heat fluctuations in p T? Cs 2 closer to Ideal Gas limit; Any structure near T c?? Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 11

33 Entropy agrees with strong coupling SYM prediction (Gubser, Klebanov & Tseytlin, NPB 98, 202) for T = 2 3T c but fails at lower T, as do various weak coupling schemes : s s 0 = f(g 2 N c ), where f(x) = ζ(3)x 3/2 + and s 0 = 2 3 π2 Nc 2 T 3. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 12

34 Entropy agrees with strong coupling SYM prediction (Gubser, Klebanov & Tseytlin, NPB 98, 202) for T = 2 3T c but fails at lower T, as do various weak coupling schemes : s s 0 = f(g 2 N c ), where f(x) = ζ(3)x 3/2 + and s 0 = 2 3 π2 Nc 2 T T c 2T c 1.5T c 1.25T c ( b ) AdS/CFT (s/s 0 )-(3/4) g 2 N c Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 12

35 Entropy agrees with strong coupling SYM prediction (Gubser, Klebanov & Tseytlin, NPB 98, 202) for T = 2 3T c but fails at lower T, as do various weak coupling schemes : s s 0 = f(g 2 N c ), where f(x) = ζ(3)x 3/2 + and s 0 = 2 3 π2 Nc 2 T 3. 1 ( b ) 3T c 2T c 1.5T c AdS/CFT 1.25T c Ideal gas (s/s 0 )-(3/4) P/T Conformal g 2 N c ε/t 4 Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 12

36 Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 13

37 Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 14

38 QGP - (Almost) Perfect Liquid v 2 (p T ) b 6.8 fm (16-24% Central) STAR Data Γ s /τ o = Γ s /τ o = 0.1 Γ s /τ o = p T (GeV) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 15

39 QGP - (Almost) Perfect Liquid v 2 (p T ) b 6.8 fm (16-24% Central) Elliptic flow consistent with Ideal Hydrodynamics bound on Shear viscosity. (D. Teaney, nucl-th/ ; PRC 2003) STAR Data Γ s /τ o = Γ s /τ o = Γ s /τ o = p T (GeV) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 15

40 QGP - (Almost) Perfect Liquid v 2 (p T ) b 6.8 fm (16-24% Central) STAR Data Γ s /τ o = 0 Γ s /τ o = 0.1 Elliptic flow consistent with Ideal Hydrodynamics bound on Shear viscosity. (D. Teaney, nucl-th/ ; PRC 2003) 4 3 Γ s = η st, (1) where η is Shear Viscosity and s is entropy density; τ = t 2 z 2 is the time scale of expansion Γ s /τ o = p T (GeV) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 15

41 QGP - (Almost) Perfect Liquid v 2 (p T ) b 6.8 fm (16-24% Central) STAR Data Γ s /τ o = 0 Γ s /τ o = 0.1 Elliptic flow consistent with Ideal Hydrodynamics bound on Shear viscosity. (D. Teaney, nucl-th/ ; PRC 2003) 4 3 Γ s = η st, (1) where η is Shear Viscosity and s is entropy density; τ = t 2 z 2 is the time scale of expansion Γ s /τ o = 0.2 Perturbation theory Large η/s Small η/s Strongly Coupled Liquid p T (GeV) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 15

42 η s Perturbative Theory KSS bound Nakamura and Sakai, PRL 94 (2005). T Tc Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 16

43 η s Kubo s Linear Response Theory : Transport Coefficients in terms of equilibrium correlation functions. 1.0 Perturbative Theory 0.5 KSS bound Nakamura and Sakai, PRL 94 (2005). T Tc Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 16

44 η s Perturbative Theory KSS bound Kubo s Linear Response Theory : Transport Coefficients in terms of equilibrium correlation functions. Obtain Energy-Momentum Correlation functions on Lattice (at discrete Matsubara frequencies) T Tc Continue them to get Retarded ones Shear, Bulk Viscosities. Nakamura and Sakai, PRL 94 (2005). Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 16

45 η s Perturbative Theory KSS bound Kubo s Linear Response Theory : Transport Coefficients in terms of equilibrium correlation functions. Obtain Energy-Momentum Correlation functions on Lattice (at discrete Matsubara frequencies) Nakamura and Sakai, PRL 94 (2005). T Tc Continue them to get Retarded ones Shear, Bulk Viscosities. Larger lattices and inclusion of dynamical quarks in future. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 16

46 Wróblewski Parameter Measure of Strangeness produced. Quark number susceptibilities, χ ij ln Z/ µ i µ j, can provide a handle; QNS also useful theoretical check on models. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 17

47 Wróblewski Parameter Measure of Strangeness produced. Quark number susceptibilities, χ ij ln Z/ µ i µ j, can provide a handle; QNS also useful theoretical check on models. Fluctuation-Dissipation Theorem Production of Strange quark-antiquark pair imaginary part of generalized strange quark susceptibility. Kramers - Krönig relation can be used to relate it to the real part of the susceptibility, which we obtain from lattice QCD simulations. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 17

48 Wróblewski Parameter Measure of Strangeness produced. Quark number susceptibilities, χ ij ln Z/ µ i µ j, can provide a handle; QNS also useful theoretical check on models. Fluctuation-Dissipation Theorem Production of Strange quark-antiquark pair imaginary part of generalized strange quark susceptibility. Kramers - Krönig relation can be used to relate it to the real part of the susceptibility, which we obtain from lattice QCD simulations. Finally, make a relaxation time approximation (ωτ 1) ratio of real parts is the same as the ratio of imaginary parts. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 17

49 We use m/t c = 0.03 for u, d and m/t c = 1 for s quark; At each T, ratio of χ s and χ ud λ s (T ). Extrapolate it to T c. (RVG & Sourendu Gupta, PRD 2002, PRD 2003 and PRD 2006) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 18

50 We use m/t c = 0.03 for u, d and m/t c = 1 for s quark; At each T, ratio of χ s and χ ud λ s (T ). Extrapolate it to T c. (RVG & Sourendu Gupta, PRD 2002, PRD 2003 and PRD 2006) AGS Si-Au Quenched QCD (T ) c RHIC Au-Au SpS S-S SpS S-Ag SpS Pb-Pb AGS Au-Au λs Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 18

51 We use m/t c = 0.03 for u, d and m/t c = 1 for s quark; At each T, ratio of χ s and χ ud λ s (T ). Extrapolate it to T c. (RVG & Sourendu Gupta, PRD 2002, PRD 2003 and PRD 2006) Quenched QCD (T ) c 1 RHIC Au-Au 0.8 SpS S-S SpS S-Ag 0.6 SpS Pb-Pb AGS Au-Au λ s 0.4 Nt = 4 Nf = 2 QCD AGS Si-Au λs T/Tc Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 18

52 Baryon-Strangeness Correlation Correlation between quantum numbers K and L can be studied through the ratio C (KL)/L = KL K L. L 2 L 2 These are robust : theoretically & experimentally. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 19

53 Baryon-Strangeness Correlation Correlation between quantum numbers K and L can be studied through the ratio C (KL)/L = KL K L. L 2 L 2 These are robust : theoretically & experimentally. Baryon Number(Charge) Strangeness correlation : C (BS)/S (C (QS)/S ) (Koch, Majumdar and Randurp, PRL 95 (2005); RVG & Sourendu Gupta, PR D 2006; S. Mukherjee, hep-lat/ ); u-d Correlation. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 19

54 Baryon-Strangeness Correlation Correlation between quantum numbers K and L can be studied through the ratio C (KL)/L = KL K L. L 2 L 2 These are robust : theoretically & experimentally. Baryon Number(Charge) Strangeness correlation : C (BS)/S (C (QS)/S ) (Koch, Majumdar and Randurp, PRL 95 (2005); RVG & Sourendu Gupta, PR D 2006; S. Mukherjee, hep-lat/ ); u-d Correlation X=Q 1 XS C X=B T/Tc Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 19

55 Baryon-Strangeness Correlation Correlation between quantum numbers K and L can be studied through the ratio C (KL)/L = KL K L. L 2 L 2 These are robust : theoretically & experimentally. Baryon Number(Charge) Strangeness correlation : C (BS)/S (C (QS)/S ) (Koch, Majumdar and Randurp, PRL 95 (2005); RVG & Sourendu Gupta, PR D 2006; S. Mukherjee, hep-lat/ ); u-d Correlation. XS C X=Q X=B C (ud)/u T/Tc T/Tc Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 19

56 Anomalous J/ψ Suppression : CERN NA50 results Matsui-Satz idea J/ψ suppression as a signal of QGP. Deconfinement Screening of coloured quarks, which cannot bind. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 20

57 Anomalous J/ψ Suppression : CERN NA50 results Matsui-Satz idea J/ψ suppression as a signal of QGP. Deconfinement Screening of coloured quarks, which cannot bind. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 20

58 Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 21

59 J/ψ Suppression Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 22

60 J/ψ Suppression or Not? Original Matsui-Satz idea Based on Quarkonium potential model calculations and an Ansatz for temperature dependence dissolution of J/ψ and χ c by 1.1T c. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 22

61 J/ψ Suppression or Not? Original Matsui-Satz idea Based on Quarkonium potential model calculations and an Ansatz for temperature dependence dissolution of J/ψ and χ c by 1.1T c. Impressive NA50 results from CERN and now PHENIX at RHIC. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 22

62 J/ψ Suppression or Not? Original Matsui-Satz idea Based on Quarkonium potential model calculations and an Ansatz for temperature dependence dissolution of J/ψ and χ c by 1.1T c. Impressive NA50 results from CERN and now PHENIX at RHIC. A critical assessment of the original theoretical argument: Made feasible by the recognition of MEM technique as a tool to extract spectral functions from the temporal correlators computed on the Euclidean lattice. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 22

63 J/ψ Suppression or Not? Original Matsui-Satz idea Based on Quarkonium potential model calculations and an Ansatz for temperature dependence dissolution of J/ψ and χ c by 1.1T c. Impressive NA50 results from CERN and now PHENIX at RHIC. A critical assessment of the original theoretical argument: Made feasible by the recognition of MEM technique as a tool to extract spectral functions from the temporal correlators computed on the Euclidean lattice. Caution : nonzero temperature obtained by making temporal lattices shorter : to Lattices used. (S. Datta et al., Phys. Rev. D 69, (2004).) Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 22

64 0.16 η c J/ψ T c 1.1T c 1.5T c 0.12 χ c0 χ c1 ρ(ω) 0.08 ρ(ω) T c 1.5T c 2.25T c 3T c ω[gev] ω[gev] χ c seems to indeed dissolve by 1.1T c, however, J/ψ and η c persist up to 2.25 T c and are gone at 3T c ; Similar results by Asakawa-Hatsuda and Matsufuru. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 23

65 0.16 η c J/ψ T c 1.1T c 1.5T c 0.12 χ c0 χ c1 ρ(ω) 0.08 ρ(ω) T c 1.5T c 2.25T c 3T c ω[gev] ω[gev] χ c seems to indeed dissolve by 1.1T c, however, J/ψ and η c persist up to 2.25 T c and are gone at 3T c ; Similar results by Asakawa-Hatsuda and Matsufuru. Since about % observed J/ψ come through χ and ψ decays, expect changes of suppression patterns as a function of T or s. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 23

66 0.16 η c J/ψ T c 1.1T c 1.5T c 0.12 χ c0 χ c1 ρ(ω) 0.08 ρ(ω) T c 1.5T c 2.25T c 3T c ω[gev] ω[gev] χ c seems to indeed dissolve by 1.1T c, however, J/ψ and η c persist up to 2.25 T c and are gone at 3T c ; Similar results by Asakawa-Hatsuda and Matsufuru. Since about % observed J/ψ come through χ and ψ decays, expect changes of suppression patterns as a function of T or s. No Significant Effect of inclusion of dynamical fermions? Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 23

67 Jet Quenching Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 24

68 Jet Quenching Rare, Highly Energetic Scatterings produce jets of particles : g + g g + g. Quark-Gluon Plasma, any medium in general, interacts with a jet, causing it to lose energy Jet Quenching. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 24

69 On-Off test possible Compare Collisions of Heavy-Heavy nuclei with Light-Heavy or Light-Light. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 25

70 On-Off test possible Compare Collisions of Heavy-Heavy nuclei with Light-Heavy or Light-Light. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 25

71 Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 26

72 QCD Phase Diagram Expected QCD Phase Diagram and Lattice Approaches to unravel it. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 27

73 QCD Phase Diagram Expected QCD Phase Diagram and Lattice Approaches to unravel it. Lee-Yang Zeroes and Two parameter Re-weighting (Z. Fodor & S. Katz, JHEP 0203 (2002) 014 ). Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 27

74 QCD Phase Diagram Expected QCD Phase Diagram and Lattice Approaches to unravel it. Lee-Yang Zeroes and Two parameter Re-weighting (Z. Fodor & S. Katz, JHEP 0203 (2002) 014 ). Imaginary Chemical Potential (Ph. Forcrand & O. Philipsen, NP B642 (2002) 290; M.-P. Lombardo & M. D Elia PR D67 (2003) ). de Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 27

75 QCD Phase Diagram Expected QCD Phase Diagram and Lattice Approaches to unravel it. Lee-Yang Zeroes and Two parameter Re-weighting (Z. Fodor & S. Katz, JHEP 0203 (2002) 014 ). Imaginary Chemical Potential (Ph. Forcrand & O. Philipsen, NP B642 (2002) 290; M.-P. Lombardo & M. D Elia PR D67 (2003) ). de Taylor Expansion (C. Allton et al., PR D66 (2002) & D68 (2003) ; R.V. Gavai and S. Gupta, PR D68 (2003) ). Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 27

76 Critical Point Estimate B /T µ 2 1 RVG & S. Gupta, PR D /T B µ n n Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 28

77 Critical Point Estimate B /T µ 2 1 RVG & S. Gupta, PR D /T B µ n n Radii of convergence as a function of the order of expansion at T = 0.95T c on N s = 8 (circles) and 24 (boxes). Left panel for ρ n and right one for r n. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 28

78 Critical Point Estimate B /T µ 2 1 RVG & S. Gupta, PR D /T B µ n n Radii of convergence as a function of the order of expansion at T = 0.95T c on N s = 8 (circles) and 24 (boxes). Left panel for ρ n and right one for r n. Extrapolation in n µ E /T E = 1.1 ± 0.2 at T E = 0.95T c. Finite volume shift consistent with Ising Universality class. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 28

79 Summary Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 29

80 Lattice QCD predicts transition to Quark-Gluon Plasma and several of its properties, T c, EoS, λ s, η... RHIC data exhibits tell-tale signs of QGP : Flow, Jet Quenching, J/ψ... BUT agreement with theory still not very quantitative, e.g., v 2. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 30

81 Lattice QCD predicts transition to Quark-Gluon Plasma and several of its properties, T c, EoS, λ s, η... RHIC data exhibits tell-tale signs of QGP : Flow, Jet Quenching, J/ψ... BUT agreement with theory still not very quantitative, e.g., v 2. Our results on correlations of quantum numbers suggest quark-like excitations in QGP. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 30

82 Lattice QCD predicts transition to Quark-Gluon Plasma and several of its properties, T c, EoS, λ s, η... RHIC data exhibits tell-tale signs of QGP : Flow, Jet Quenching, J/ψ... BUT agreement with theory still not very quantitative, e.g., v 2. Our results on correlations of quantum numbers suggest quark-like excitations in QGP. Phase diagram in T µ B plane has begun to emerge: Our estimate for the critical point is µ B /T 1 2. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 30

83 Lattice QCD predicts transition to Quark-Gluon Plasma and several of its properties, T c, EoS, λ s, η RHIC data exhibits tell-tale signs of QGP : Flow, Jet Quenching, J/ψ... BUT agreement with theory still not very quantitative, e.g., v 2. T/Tc GeV 20 GeV 18 GeV (CERN) 10 GeV Our results on correlations of quantum numbers suggest quark-like excitations in QGP. Phase diagram in T µ B plane has begun to emerge: Our estimate for the critical point is µ B /T Freezeout curve µ B /T Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 30

84 1.0 NL 0.9 HTL χ 3 /Τ T/T c Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 31

85 Quarkonia moving in the Heat Bath Should see more energetic gluons. More Dissociation at the same T as momentum of J/ψ increases? Datta et al. SEWM 2004, PANIC Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 32

86 Quarkonia moving in the Heat Bath Should see more energetic gluons. More Dissociation at the same T as momentum of J/ψ increases? Datta et al. SEWM 2004, PANIC ρ(ω) η c ρ(ω) J/ψ T p[gev] T c 1.1T c p [GeV] ω[gev] T c 1.1T c ω[gev] Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 32

87 Quarkonia moving in the Heat Bath Should see more energetic gluons. More Dissociation at the same T as momentum of J/ψ increases? Datta et al. SEWM 2004, PANIC ρ(ω) η c ρ(ω) J/ψ T p[gev] T c 1.1T c p [GeV] ω[gev] T c 1.1T c ω[gev] Both J/ψ and η c do show this trend. The effect is significant at both 0.75 and 1.1T c. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 32

88 m ρ /T c m π /m ρ m N /m ρ N s m π flavours T E /T c µ E B /T E (5) (2) (2) 2.2 (2) 5.12 (8) (6) (3) 4.5 (2) 5.4 (2) 0.31 (1) 1.8 (2) (2) 1.1 (2) 5.4 (2) 0.31 (1) 1.8 (2) (1) 0.70 (1) Table 1: Summary of critical end point estimates the lattice spacing is a = 1/4T. N s is the spatial size of the lattice and N s m π is the size in units of the pion Compton wavelength, evaluated for T = µ = 0. The ratio m π /m K sets the scale of the strange quark mass. Results are sequentially from Fodor-Katz 04, Fodor-Katz 02, Gavai-Gupta, de Forcrand- Philipsen and Bielefeld-Swansea. Workshop on LHC Physics 2006, T. I. F. R., Mumbai, September 7, 2006 R. V. Gavai Top 33