Did we create the QGP at RHIC and the LHC? Scott Pratt Michigan State University SP, W.P. McCormack and Claudia Ratti, arxiv 149.2164 (214)
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 +...
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) 1 15 2 25 3 35 4 T (MeV)
Two waves of quark production hadronization up or down quarks isentropic expansion thermalization strange quarks 1 2 3 4 5 6 7 8 9 τ (fm/c)
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
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
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 )
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 Δη=
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
Positive or negative (electric charge) χ (HAD) electric = n charged, χ (QGP) electric = 4 9 n u + 1 9 n d + 1 9 n s Narrower bump is same sign, but much taller than wider bump π + π - BF dominated by narrow peak
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
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
STAR results verify this qualitatively STAR Preliminary corrected for acceptance/efficiency Noisy in tails
Back to being more exact Take into account: hadrons have multiple charges, hyperon decays χ has off-diagonal elements
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
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 β >)
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
3. Apply to our problem abcd δ n α ()δ n β (Δη) = n α n β q αa χ (HAD) 1 ' ac g cd (Δη)χ (HAD) 1 db q βb
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 αβ
3. Prefactors... p + + K + p.441,-.66.485,-.162.491,-.146.479,-.178.535,-.242.529,-.258.578,-.338.6,.16 -.44,.96.183,-.61.242,-.94.242,-.94.242,-.94.32,-.128.32,-.128.361,-.161.,-. -.59,.33.74,-.22.97,-.38.99,-.33.95,-.43.122,-.49.12,-.54.144,-.64.2,.5 -.23,.16 +.72,-.27.97,-.38.95,-.43.99,-.33.12,-.54.122,-.49.144,-.64 -.2,-.5 -.25,.11.46,-.21.69,-.29.7,-.28.69,-.31.93,-.36.92,-.38.115,-.45.1,.1 -.23,.8 +.46,-.22.69,-.29.69,-.31.7,-.28.92,-.38.93,-.36.115,-.45 -.1,-.1 -.23,.7 +.9,-.5.15,-.7.15,-.7.15,-.7.21,-.8.21,-.8.27,-.9 -.,-. -.6,.1.119,.318.,-..239,.636 -.239,-.636.119,.318 -.119,-.318 -.,-..239,.636.119,.318 K -.175,.384 -.627,.352 -.63,.417 -.651,.288-1.55,.385-1.79,.321-1.57,.354.24,.64.452,.31 (QGP,HAD) prefactors completely determined by χqgp and final-state hadronic yields (hadron yields from thermal model with B-reduction)
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 )
1 2 3 π + π - G π + π ( y)/(dn ch /dy).6.4.2.6.4 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.2 1 2 3 y
1 2 3 G + K K -( y)/(dn ch /dy).6.4.2.6.4 hadronic contr. QGP contr. total default half strangeness K + K - Little hadronic contribution Can test whether QGP is rich in strangeness.2 1 2 3 y
1 2 3 1 2 3.3.2 QGP contribution hadronic contr. total p-pbar G pp ( y)/(dn ch /dy).1 -.1.3.2.1.5.4 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.3.2.1 1 2 3 y
1 2 3 G pk -( y)/(dn ch /dy).3.2.1 -.1.1 QGP contribution hadronic contr. total default with off-diagonal elements pk - 1 2 3 y QGP contribution negative dips negative too narrow for one source
.5 1 1.5 2.5 1 1.5 2.4 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).5.1.5 K + K - pp -.5 1 1.5 2 y
.5 1 1.5 2.5 1 1.5 2.4 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).2.1.5 π + π - K + K - Blast wave.1.5 pp -.5 1 1.5 2 y
4 Parameter MCMC!.3 < σqgp < 1.5 < σhad/σqgp <1.75< quark/hadron <2" < s/(ud) <1
MCMC Results
π + π - K + K - pp _
MCMC Results σ qgp.9 σ had.25 ud/hadron and s/ud ratios consistent with equilibratedqgp
Conclusions STAR results validate 2-wave picture uds densities of QGP appear close to lattice values
To-Do List: (Experiment) Other charge combinations, e.g. pk " σφ vs ση" B(Qout,Qside,Qlong) binned by kt, reaction plane pp, pa collisions
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?
.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6 1.1 1 1.1 1 1 1.1.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6.4.8 1.2 1.6-5% cent 1 Sample Spectra from Prior and Posterior dn/(2πp t dp t dy) (GeV/c) -2.1 1 1.1 protons kaons 2-3% cent protons kaons 1 1 1.1 pions pions.4.8 1.2 1.6.4.8 1.2 1.6 p t (GeV/c) p t (GeV/c)
Two Calculations 1. J.Novak, K. Novak, S.P., C.Coleman-Smith & R.Wolpert, ArXiv: 133.5769 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 1.25.25.5 η/s η/s.2 5 T dep. of η Rel. Probability.85 1.25 energy norm. 1.2 3 4 σsat (mb) 5.5 W.N./Sat. frac. 1.25.75 Init. Flow 1.25.2.26 η/s.5 2.5 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
Constraining the Eq. of State with RHIC/LHC Data