Did we create the QGP at RHIC and the LHC? Scott Pratt Michigan State University. SP, W.P. McCormack and Claudia Ratti, arxiv 1409.

Transcription

1 Did we create the QGP at RHIC and the LHC? Scott Pratt Michigan State University SP, W.P. McCormack and Claudia Ratti, arxiv (214)

2 QGP Chemistry Basics 52 ~massless degrees of freedom Strongly interacting Conserved: up-ness, down-ness, strangeness, color Not conserved: quark number Lattice measures charge fluctuations: χ ab Q a Q b /V Parton gas: χ QGP ab = (n a + n a )δ ab Hadron gas: χ HAD ab = n α q α,a q α,b α a,b = uds α=π +,π -,π,k +...

3 Lattice Charge Fluctuations scaled by entropy Courtesy of Claudia Ratti.12 uu ss us Parton gas: χ ab /s.8.4 Had. Gas S.B. limit χ QGP ab = (n a + n a )δ ab Hadron gas: χ HAD ab = n α q α,a q α,b α α=π +,π -,π,k +... off-diagonal elements (V.Koch, PRL 25) T (MeV)

4 Two waves of quark production hadronization up or down quarks isentropic expansion thermalization strange quarks τ (fm/c)

5 Problems with Comparing Experiment to Lattice 1.Lattice = Grand Canonical (Particle Bath) Experiment = Canonical (net charge = ) 2. Charge created at hadronization 3. One measures hadrons -- not uds 4. One measures momenta, not positions

6 Consider charge correlations gab(δη) is a 3x3 matrix g ud (Δη) Q u (η)q d (η + Δη) [ ] n d (η + Δη) n d (η + δη) = n u (η) n u (η) d Δη g (Δη) = ab charge conservation

7 1. Before hadronization d Δη g ab(δη) = g ud = g us = g ds = σ QGP (QGP) g ab (Δη) = χ ab From lattice! δ (Δη) exp( Δη2 2 / 2σ (QGP ) 2 (2πσ (QGP) ) 1/2 only extra parameter )

8 2. Just after hadronization guu(δη) g ab (Δη) = χ (HAD) ab δ (Δη) ( (QGP) ) exp( Δη2 2 / 2σ ( HAD ) 2 (2πσ (HAD) ) 1/2 χ (HAD) ab χ ab χ ab (QGP) exp( Δη2 2 / 2σ (QGP ) 2 (2πσ (QGP) ) 1/2 ) ) σ QGP σ HAD g ab (Δη) can't change suddenly except at Δη=

9 2nd Bump: Positive or Negative? Crude Expectations guu(δη) g ab (Δη) = χ (HAD) ab δ (Δη) ( (QGP) ) exp( Δη2 2 / 2σ ( HAD ) 2 (2πσ (HAD) ) 1/2 χ (HAD) ab χ ab χ ab (QGP) exp( Δη2 2 / 2σ (QGP ) 2 (2πσ (QGP) ) 1/2 ) ) σ QGP σ HAD Sign depends on χ (HAD) (QGP) ab χ ab

10 Positive or negative (electric charge) χ (HAD) electric = n charged, χ (QGP) electric = 4 9 n u n d n s Narrower bump is same sign, but much taller than wider bump π + π - BF dominated by narrow peak

11 Before and after (STRANGENESS) χ (HAD) strange = n K + n Λ + 4n Ξ + 9n Ω, χ (QGP) strange = n s Narrower bump is very small K + K - BF dominated by broad peak

12 Before and after (BARYON No.) χ (HAD) baryon = n p + n n + n Λ + n Ξ + n Ω, χ (QGP) baryon = 1 9 (n + n + n ) u d s Smaller Bump is smaller with opposite sign p-pbar BF dominated by broad peak, but has dip

13 STAR results verify this qualitatively STAR Preliminary corrected for acceptance/efficiency Noisy in tails

14 Back to being more exact Take into account: hadrons have multiple charges, hyperon decays χ has off-diagonal elements

15 2. Just after hadronization (uds basis) ' g ab (Δη) = χ ab (QGP) e Δη2 2 /2σ (QGP ) 2 2πσ (QGP) + (χ (HAD) ab χ ab (QGP) ) e Δη2 2 /2σ ( HAD ) 2 2πσ (HAD) χ ab n α q α,a q α,b α hadrons χ (QGP) ab 2n a δ ab

16 3. But, we measure G αβ not g ab!!! α,β=π,p,k...a,b=u,d,s G αβ (Δη) [n α n α ][n β n β ] e.g., G pk = [n p n p ][n K n K + ] Generalized Balance Function (aside from factor of <n β >)

17 Analogous problem... Given δρa and n α, find δn α Solution: assign chemical potential δ n α = n α ( e µ aq /T α,a 1) µ n α q αa q b αb T = δρ a = δ n α q α,a α αb b χ (had) µ b ab T µ a T = b δ n α = n α (χ 1 ) ab δρ b b q αa (χ 1 ) ab δρ b

18 3. Apply to our problem abcd δ n α ()δ n β (Δη) = n α n β q αa χ (HAD) 1 ' ac g cd (Δη)χ (HAD) 1 db q βb

19 3. Putting this together ' (QGP) e Δη2 2 /2σ (QGP ) G αβ (Δη) = w αβ 2 2πσ (QGP) (HAD) e Δη2 2 /2σ ( HAD ) + w αβ 2 2πσ (HAD) abcd w (QGP) αβ = 2 n α q α,a χ 1(HAD) ab χ (QGP) 1(HAD) bc χ cd ab w (HAD) 1(HAD) αβ = 2 n α q α,a χ ab n β prefactors depend only on yields and χab from lattice n β q β,d (QGP) q β,b w αβ

20 3. Prefactors... p + + K + p.441, , , , , , , , , , , , , , , ,-.161., ,.33.74, , , , , , ,-.64.2, , , , , , , , , , ,.11.46, ,-.29.7, , , , ,-.45.1, , , , ,-.31.7, , , , , ,.7 +.9,-.5.15,-.7.15,-.7.15,-.7.21,-.8.21,-.8.27,-.9 -.,-. -.6,.1.119,.318.,-..239, , , , ,-..239, ,.318 K -.175, , , , , , , , ,.31 (QGP,HAD) prefactors completely determined by χqgp and final-state hadronic yields (hadron yields from thermal model with B-reduction)

21 4. Use blast-wave to go from coordinate space η to momentum-space rapidity (Monte Carlo + decays) Use STAR parameters fit to spectra (T and u )

22 1 2 3 π + π - G π + π ( y)/(dn ch /dy) hadronization contr. QGP contribution total σ QGP =1.5 σ QGP =.5 σ QGP =1. Hadronization part narrower Can't well separate components due to thermal smearing acceptance narrows with centrality y

23 1 2 3 G + K K -( y)/(dn ch /dy) hadronic contr. QGP contr. total default half strangeness K + K - Little hadronic contribution Can test whether QGP is rich in strangeness y

24 QGP contribution hadronic contr. total p-pbar G pp ( y)/(dn ch /dy) default half quark density with off-diagonal elements hadron contribution negative tests two-wave nature no narrowing with centrality sensitive to quark density of QGP y

25 1 2 3 G pk -( y)/(dn ch /dy) QGP contribution hadronic contr. total default with off-diagonal elements pk y QGP contribution negative dips negative too narrow for one source

26 STAR preliminary from STAR (balance functions) thesis of Hui Wang (212).2.1 π + π - Single source won't work KK broader than ππ! pp - broader than both -- and has dip!! B( y) K + K - pp y

27 STAR preliminary Compare to Blast Wave 7 parameters: T blast =12 MeV u =.732 σ qgp =1., σ had =, quark/hadron=1, s/(ud)=.95 B( y) π + π - K + K - Blast wave.1.5 pp y

28 4 Parameter MCMC!.3 < σqgp < 1.5 < σhad/σqgp <1.75< quark/hadron <2" < s/(ud) <1

29 MCMC Results

30 π + π - K + K - pp _

31 MCMC Results σ qgp.9 σ had.25 ud/hadron and s/ud ratios consistent with equilibratedqgp

32 Conclusions STAR results validate 2-wave picture uds densities of QGP appear close to lattice values

33 To-Do List: (Experiment) Other charge combinations, e.g. pk " σφ vs ση" B(Qout,Qside,Qlong) binned by kt, reaction plane pp, pa collisions

34 To-Do List: (Theory) ' g ab Continuous creation/annihilation" (x 1, x 2 ) = τ d 4 x'd a (x', x 1 )D b ( x, x 2 ) d(χ / s) ab (x') d τ! non-zero baryon density pp/e + e - collisions need theory How do quarks arise from gluons/string/fluxtubes?

35 % cent 1 Sample Spectra from Prior and Posterior dn/(2πp t dp t dy) (GeV/c) protons kaons 2-3% cent protons kaons pions pions p t (GeV/c) p t (GeV/c)

36 Two Calculations 1. J.Novak, K. Novak, S.P., C.Coleman-Smith & R.Wolpert, ArXiv: energy norm. Rel. Probability RHIC Au+Au Data 6 parameters 5 σsat (mb) σsat (mb) 3 1 W.N./Sat. frac. W.N./Sat. frac. Init. Flow Init. Flow η/s η/s.2 5 T dep. of η Rel. Probability energy norm σsat (mb) 5.5 W.N./Sat. frac Init. Flow η/s T dep. of η 5 2. S.P., E.Sangaline, P.Sorensen & H.Wang, in progress RHIC Au+Au and LHC Pb+Pb Data 14 parameters, include Eq. of State

37 Constraining the Eq. of State with RHIC/LHC Data