Quarkonium as Probe of Hot QCD Medium · 2013-01-31 · Ágnes Mócsy, XQCD July 18-20 2011, San...
Transcript of Quarkonium as Probe of Hot QCD Medium · 2013-01-31 · Ágnes Mócsy, XQCD July 18-20 2011, San...
Quarkonium as Probe of Hot QCD Medium
Ágnes Mócsy
Quarkonium as Probe of Hot QCD Medium
Ágnes Mócsy
in collaboration with Peter Petreczky (BNL), Chuan Miao (Uni Mainz)
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
QCD Expectations
2
J/ψ suppression proposed
signal of deconfined QGP
- Temperature effects in deconfined medium: screening, Landau damping, ...
- In-medium properties encoded in quarkonium spectral functions
- Theory predicts the J/ψ disappears in the plasma
Dissolution (“melting”) seen as progressive broadening and disappearance of bound-state peaks
T=0 T>0 T>Td
Matsui, Satz PLB 1986
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
What We Know from Lattice
3
Deconfinement Color screening
F1(r,T)[GeV]
T=0
- Rapid rise of the energy density: liberation of new degrees of freedom
- Strong screening of static Q-Qbar free energy - at shorter distances with increasing T
r[fm]Cheng et al (RBC-Bielefeld) PRD 77, 014511 (2008)
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
What We Know from Lattice
4
Deconfinement Color screening
F1(r,T)[GeV]
T=0
- Rapid rise of the energy density: liberation of new degrees of freedom
- Strong screening of static Q-Qbar free energy - at shorter distances with increasing T
rscr < rJ/ψ “melting” of the J/ψ
J/ψ
Y
r[fm]Cheng et al (RBC-Bielefeld) PRD 77, 014511 (2008)
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 5
- J/ψ-suppression pattern observed at SPS and RHIC and LHC
- J/ψ nuclear modification factor: yield in AA collisions relative to yield in pp (where no QGP formation expected) scaled with number of binary NN collisions
Experimental RAA
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Experimental RAA
6
It is difficult to unambiguously interpret - we are still not there
The J/ψ story: Two decades worth of data Modest theory advancement
Lots of ad-hoc phenomenological modeling
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Experimental RAA
7
It is difficult to unambiguously interpret - we are still not there
The J/ψ story: Two decades worth of data Modest theory advancement
Lots of ad-hoc phenomenological modeling
The Υ story: just started !
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
What’s the Physics Behind ?
8
To answer this question we need to know:
- How the properties of J/ψ change in a deconfined medium Determine the spectral function
- Relate an equilibrium spectral function to RAA Through real-time dynamics
- Identify what physics might contribute to RAA for example: Suppression is seen in pA, dA data as well (where no QGP formation expected) Cold nuclear matter effects - could be relevant to AA
But is the J/ψ RAA a signal for deconfinement and screening ?!
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Determine the Spectral Function
9
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Theory Progress
10
- Lattice QCD, potential models, effective field theories (EFT)
- Recently: Spectral functions are calculated Cabrera, Rapp, Mócsy, Petreczky, Alberico, Beraudo,Laine et al Kaczmarek et al, Morita et al, ...
True when Ebin is the smallest scale. - Potential model assumes: most medium effects on quarkonium described by a T-dependent potential.
Internal energy U or free energy F or combination of these
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Theory Progress
11
- Lattice QCD, potential models, effective field theories (EFT)
- Recently: Spectral functions are calculated Cabrera, Rapp, Mócsy, Petreczky, Alberico, Beraudo,Laine et al Kaczmarek et al, Morita et al, ...
True when Ebin is the smallest scale. - Potential model assumes: most medium effects on quarkonium described by a T-dependent potential.
Internal energy U or free energy F or combination of these
Which one of these is the correct potential has become an irrelevant question!
Thermodynamic quantities calculated on the lattice related to potential only in very special cases (ex. LO weak coupling)
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Theory Progress
12
- Recently: Potential model (can be placed on more solid grounds) appears as the tree-level approximation of the EFT and can be systematically improved
Isotrop: Laine et al, Brambilla et al, Blaizot et al, Escobado, Soto, ...Anisotrop: Strickland et al, Philipsen et al,
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Theory Progress
12
- Recently: Potential model (can be placed on more solid grounds) appears as the tree-level approximation of the EFT and can be systematically improved
T-effects: screening, Landau-damping, singlet-octet transition...
Isotrop: Laine et al, Brambilla et al, Blaizot et al, Escobado, Soto, ...Anisotrop: Strickland et al, Philipsen et al,
complex potential
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Theory Progress
12
- Recently: Potential model (can be placed on more solid grounds) appears as the tree-level approximation of the EFT and can be systematically improved
T-effects: screening, Landau-damping, singlet-octet transition...
Isotrop: Laine et al, Brambilla et al, Blaizot et al, Escobado, Soto, ...Anisotrop: Strickland et al, Philipsen et al,
complex potentialbonus: the applicability of potential model approach is apparent
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Theory Progress
12
- Recently: Potential model (can be placed on more solid grounds) appears as the tree-level approximation of the EFT and can be systematically improved
T-effects: screening, Landau-damping, singlet-octet transition...
Isotrop: Laine et al, Brambilla et al, Blaizot et al, Escobado, Soto, ...Anisotrop: Strickland et al, Philipsen et al,
complex potentialbonus: the applicability of potential model approach is apparent
drawback: weak coupling assumption. how to include perturbative effects?
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Complex Potential
13
Mócsy, Petreczky, PRL 99 (07) 211602 Burnier, Laine, Vepsalainen JHEP 0801 (08) 043 Beraudo, arXiv:0812.1130
Maximal value (upper bound) Minimal value (lower bound)
Constrain ReVs(r) by lattice QCD data on the singlet free energy
Take ImVs(r) from HTL resummedpQCD calculations
Thermal contributions come from Re (→ binding energy) and Im (→ width) part
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 14
1.2Tc1.2Tc
T=0 potential
Internal energyTS
- lower limit
- upper limit
Free energy
Mocsy, Petreczky, PRD 77 (2008)
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 14
1.2Tc1.2Tc
T=0 potentialMost confining potential
Internal energyTS
- lower limit
- upper limit
Free energy
Mocsy, Petreczky, PRD 77 (2008)
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Complex Potential
15
Mócsy, Petreczky, PRL 99 (07) 211602 Burnier, Laine, Vepsalainen JHEP 0801 (08) 043 Beraudo, arXiv:0812.1130
Maximal value (upper bound) Minimal value (lower bound)
Constrain ReVs(r) by lattice QCD data on the singlet free energy
Take ImVs(r) from HTL resummedpQCD calculations
Thermal contributions come from Re (→ binding energy) and Im (→ width) part
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Complex Potential
15
Mócsy, Petreczky, PRL 99 (07) 211602 Burnier, Laine, Vepsalainen JHEP 0801 (08) 043 Beraudo, arXiv:0812.1130
Maximal value (upper bound) Minimal value (lower bound)
Constrain ReVs(r) by lattice QCD data on the singlet free energy
Take ImVs(r) from HTL resummedpQCD calculations
Thermal contributions come from Re (→ binding energy) and Im (→ width) part
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Complex Potential
15
Mócsy, Petreczky, PRL 99 (07) 211602 Burnier, Laine, Vepsalainen JHEP 0801 (08) 043 Beraudo, arXiv:0812.1130
Maximal value (upper bound) Minimal value (lower bound)
Constrain ReVs(r) by lattice QCD data on the singlet free energy
Take ImVs(r) from HTL resummedpQCD calculations
Thermal contributions come from Re (→ binding energy) and Im (→ width) part
Study the effect of color screening and of dissipation on the quarkonium spectral functions.
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
16
- 1S peak with reduced binding energy - residual c-cbar correlations persist
Mocsy, Petreczky, PRL 2007, PRD 2008,
residual correlations
Take the upper limit for the real part of the potential allowed by lattice calculations
Im Vs(r) =0
Calculated in full QCD
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
2.6 2.8 3 3.2 3.4 3.6 3.8 4![GeV]
"(!)/!2T=0
245MeV326MeV449MeV
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
16
- 1S peak with reduced binding energy - residual c-cbar correlations persist
Mocsy, Petreczky, PRL 2007, PRD 2008,
residual correlations
Take the upper limit for the real part of the potential allowed by lattice calculations
Im Vs(r) =0
Calculated in full QCD
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
2.6 2.8 3 3.2 3.4 3.6 3.8 4![GeV]
"(!)/!2T=0
245MeV326MeV449MeV
Take the perturbative imaginary part of the potential and the code from Burnier, Laine, Vepsalainen JHEP 0801 (08) 043
Petreczky, Miao, Mocsy, Nucl Phys A 2011
Im Vs(r) ≠0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
2.6 2.8 3 3.2 3.4 3.6 3.8 4![GeV]
"(!)/!2T=0
245MeV326MeV449MeV
- dramatic changes with Im part, peaks strongly broaden
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
16
- 1S peak with reduced binding energy - residual c-cbar correlations persist
Mocsy, Petreczky, PRL 2007, PRD 2008,
residual correlations
Take the upper limit for the real part of the potential allowed by lattice calculations
Im Vs(r) =0
No charmonium state could survive for T > 240 MeV
Calculated in full QCD
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
2.6 2.8 3 3.2 3.4 3.6 3.8 4![GeV]
"(!)/!2T=0
245MeV326MeV449MeV
Take the perturbative imaginary part of the potential and the code from Burnier, Laine, Vepsalainen JHEP 0801 (08) 043
Petreczky, Miao, Mocsy, Nucl Phys A 2011
Im Vs(r) ≠0
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
2.6 2.8 3 3.2 3.4 3.6 3.8 4![GeV]
"(!)/!2T=0
245MeV326MeV449MeV
- dramatic changes with Im part, peaks strongly broaden
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Bottomonium Spectral Function
17
Im Vs(r) ≠0Im Vs(r) =0
2S1S 3S 1S 2S 3S
Calculated in full QCD
Mocsy, Petreczky, PRL 2007, PRD 2008
Take the upper limit for the real part of the potential allowed by lattice calculations
Take the perturbative imaginary part of the potential and the code from Burnier, Laine, Vepsalainen JHEP 0801 (08) 043
Petreczky, Miao, Mocsy, Nucl Phys A 2011
- 1S peak and remnant of 2S state there - Binding energies reduced - Threshold enhancement
- Changing Re part has little effect - Dramatic changes with Im part peaks significantly broaden
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Bottomonium Spectral Function
17
No bottomonium state could survive for T > 450 MeV
Im Vs(r) ≠0Im Vs(r) =0
2S1S 3S 1S 2S 3S
Calculated in full QCD
Mocsy, Petreczky, PRL 2007, PRD 2008
Take the upper limit for the real part of the potential allowed by lattice calculations
Take the perturbative imaginary part of the potential and the code from Burnier, Laine, Vepsalainen JHEP 0801 (08) 043
Petreczky, Miao, Mocsy, Nucl Phys A 2011
- 1S peak and remnant of 2S state there - Binding energies reduced - Threshold enhancement
- Changing Re part has little effect - Dramatic changes with Im part peaks significantly broaden
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Melting” Temperatures
18
T/TC 1/〈r〉
ϒ(1S)
J/ψ(1S)
χc(1P)≤ 1
2
1.2
χb(1P)
- Quantitative estimates of peak disappearance Thermometer of upper limits Tdiss
Note - Ebin would be smaller with other potentials -- lower melting temperatures
- Charmonium sensitive to Re V Bottomonium sensitive to Im V
Note: At finite temperature Ebinding=0 is
an overkill for dissociation !
Could be different mechanism behind melting J/ψ and Y
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Lattice says so ... “
19
- Ugly rumor : “Lattice tells J/ψ survives to 2Tc” There is no evidence for that.
- Neither the unchanged Euclidean correlator ratios, nor the extracted spectral functions suggest charmonium survival.
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Lattice says so ... “
20
Correlator ratios not sensitive to spectral function changes.
€
G τ,T( ) = σ ω,T( )K τ,ω,T( )dω∫
€
Grec τ,T( ) = σ ω,T = 0( )K τ,ω,T( )dω∫
€
G τ, r p ,T( ) = d3∫ xeir p r x jH τ, r x ( ) jH
+ 0,r 0 ( )
- Correlation function of mesonic currents in Euclidean time are calculated
- Correlators ratios studied
- We now understand that changes come from zero modes - Correlator ratios do not change in all channels- Potential model agree with lattice
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
21
MEM
€
G τ,T( ) = σ ω,T( )K τ,ω,T( )dω∫
€
Grec τ,T( ) = σ ω,T = 0( )K τ,ω,T( )dω∫€
G τ,T( ) = σ ω,T( )K τ,ω,T( )dω∫
€
Grec τ,T( ) = σ ω,T = 0( )K τ,ω,T( )dω∫
Extracted from correlation function of mesonic currents in Euclidean time
€
G τ,T( ) = σ ω,T( )K τ,ω,T( )dω∫
€
Grec τ,T( ) = σ ω,T = 0( )K τ,ω,T( )dω∫spectral function extractednot directly calculated
correlator directly calculated
Shortcomings: limited # of data points limited extent in taudefault model dependence large
O (10) data but O (100) degrees of freedom to reconstruct
Extracted in quenched QCD
Umeda et al, EPJ C39S1 (05) 9, Asakawa, Hatsuda,PRL 92 (2004) 01200, Datta et al,PRD 69 (04) 094507,Jakovac et al PRD 2007, ...
maximizes the conditional probability of having the spf given the data and
some prior knowledge
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
22
Jakovac et al PRD 2007
Extracted in quenched QCD
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
22
Jakovac et al PRD 2007
Extracted in quenched QCD
σ(ω)/ω2
ω[GeV
Mocsy, Petreczky, PRD 2008
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
22
Jakovac et al PRD 2007
Extracted in quenched QCD
σ(ω)/ω2
ω[GeV
Mocsy, Petreczky, PRD 2008
Peak (previously regarded as bound state) consistent with
threshold enhancement
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
23
Jakovac et al PRD 2007
Extracted in quenched QCD
Strong default model dependence
prior: from T=0 data at high energiesprior: perturbative continuum spf (no lattice effects)
peak - no peak ?!
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
24
Comparing low resolution confined phase (blue) to low resolution deconfined phase (red) and getting an agreement does not imply the agreement will hold at high resolution
Extracted in quenched QCD
Jakovac et al PRD 2007
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
agreement in low resolution doesnot imply agreement at high resolution
25
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
agreement in low resolution doesnot imply agreement at high resolution
26
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
agreement in low resolution doesnot imply agreement at high resolution
27
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
agreement in low resolution doesnot imply agreement at high resolution
28
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Charmonium Spectral Function
29
Jakovac et al PRD 2007
Extracted in quenched QCD
(extracted) Lattice spectral functions do not suggest Jpsi survival well into the deconfined medium
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
30
Most recent extracted charmonium spectral functionsTalk by H-T Ding (Bielefeld)
at BNL in June 2011
prior: from T=0 data at high energies
Ding, Kaczmarek, Karsch, Satz 2010
“Our analysis suggests that J/ψ is melted already by 1.46 Tc”
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
31
Unlike temporal correlators, spatial correlators can be studied at arbitrary large distances
Charmonium screening in quark gluon plasmaTalk by Swagato Mukherjee (BNL)
at BNL in June 2011
Screening mass C(z) ~ exp[−Mz ]inverse screening length of charmonium in the medium
dynamical lattice QCD simulations
freely propagating c & cbar in the medium
charmonium does not feel T
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
31
Unlike temporal correlators, spatial correlators can be studied at arbitrary large distances
Charmonium screening in quark gluon plasma
In-medium charmonia differs from vacuum for T ~ 1.2Tc
Talk by Swagato Mukherjee (BNL)at BNL in June 2011
Screening mass C(z) ~ exp[−Mz ]inverse screening length of charmonium in the medium
dynamical lattice QCD simulations
freely propagating c & cbar in the medium
charmonium does not feel T
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
32
Compared to lattice: No reduction of data points, allowing a direct comparison of T=0 and T≠0 spectral functions.
Sum rule approach to quarkonium at finite temperature
P. Gubler and M. Oka, Prog. Theor. Phys. 124, 995 (2010). P. Gubler, K. Morita and M. Oka, arXiv:1104.4436 [hep-ph].
Talk by Philipp Gubler (TokyoTech)at BNL in June 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
32
Compared to lattice: No reduction of data points, allowing a direct comparison of T=0 and T≠0 spectral functions.
Sum rule approach to quarkonium at finite temperature
P. Gubler and M. Oka, Prog. Theor. Phys. 124, 995 (2010). P. Gubler, K. Morita and M. Oka, arXiv:1104.4436 [hep-ph].
Talk by Philipp Gubler (TokyoTech)at BNL in June 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
32
Compared to lattice: No reduction of data points, allowing a direct comparison of T=0 and T≠0 spectral functions.
Sum rule approach to quarkonium at finite temperature
P. Gubler and M. Oka, Prog. Theor. Phys. 124, 995 (2010). P. Gubler, K. Morita and M. Oka, arXiv:1104.4436 [hep-ph].
Talk by Philipp Gubler (TokyoTech)at BNL in June 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
32
Compared to lattice: No reduction of data points, allowing a direct comparison of T=0 and T≠0 spectral functions.
Sum rule approach to quarkonium at finite temperature
P. Gubler and M. Oka, Prog. Theor. Phys. 124, 995 (2010). P. Gubler, K. Morita and M. Oka, arXiv:1104.4436 [hep-ph].
Talk by Philipp Gubler (TokyoTech)at BNL in June 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Comparison to Other Works
32
Compared to lattice: No reduction of data points, allowing a direct comparison of T=0 and T≠0 spectral functions.
Sum rule approach to quarkonium at finite temperature
P. Gubler and M. Oka, Prog. Theor. Phys. 124, 995 (2010). P. Gubler, K. Morita and M. Oka, arXiv:1104.4436 [hep-ph].
Both ηc and J/ψ melt between T ~ 1.0 TC and T ~ 1.1 TC
Talk by Philipp Gubler (TokyoTech)at BNL in June 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Melting” Temperatures
33
T/TC 1/〈r〉
ϒ(1S)
J/ψ(1S)
χc(1P)≤ 1
2
1.2
χb(1P)
Tc used in quarkonium studies by everyone MEM, Cabrera, Rapp, Zhao, Mocsy, Petreczky, Young, Shuryak, ...
Tc = 270MeV, 204MeV, 196MeV, 190MeV , ...
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Melting” Temperatures
33
T/TC 1/〈r〉
ϒ(1S)
J/ψ(1S)
χc(1P)≤ 1
2
1.2
χb(1P)
- Be aware of the meaning of Tc ! In pure gauge theory Tc = 270 MeV
In full QCD there is Tchiral ≅157 MeV and Tonset of screening ≅190-200 MeVBazavov, Petreczky, 2010Budapest-Wuppertal 2010, HotQCD 2011
Boyd et al 1996
Tc used in quarkonium studies by everyone MEM, Cabrera, Rapp, Zhao, Mocsy, Petreczky, Young, Shuryak, ...
Tc = 270MeV, 204MeV, 196MeV, 190MeV , ...
Chiral vs Deconfinement
34
defining transition temperature: peak of chiral susceptibility
Everyone agrees
Chiral vs Deconfinement
35
Deconfinement: - liberation of many degrees of freedom, which can be understood as a transition from hadronic degrees of freedom to partonic ones- the onset of screening
Cheng et al (RBC-Bielefeld) PRD 77, 014511 (2008)
pion gas: 3 free quark-gluon gas: ~ 8x2 + 3x2x2x3 = 52
dof ?(160-200) MeV
T > 200MeV Most relevant for
quarkonium physics at finite temperature
Debye screening, Landau damping - are effects from deconfined quarks and gluons
Chiral vs Deconfinement
36
Deconfinement: - liberation of many degrees of freedom, which can be understood as a transition from hadronic degrees of freedom to partonic ones- the onset of screening
Bazavov and P.P., arXiv:1005.1131, arXiv:1009.4914, arXiv:1012.1257, Söldner, arXiv:1012.4484
Partonic screening becomes relevant at about ~1.3Tc, where Tc is chiral transition
So we are back to T > 200MeV region
It is most relevant for quarkonium
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Melting” Temperatures
37
T/TC 1/〈r〉
ϒ(1S)
J/ψ(1S)
χc(1P)≤ 1
2
1.2
χb(1P)
Tc used in quarkonium studies by everyone MEM, Cabrera, Rapp, Zhao, Mocsy, Petreczky, Young, Shuryak, ...
Tc = 270MeV, 204MeV, 196MeV, 190MeV , ...
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Melting” Temperatures
38
New Rules: Let’s talk about absolute temperatures instead in Tc units
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
“Melting” Temperatures
38
T/TC 1/〈r〉
ϒ(1S)
J/ψ(1S)
χc(1P)≤ 1
2
1.2
χb(1P)
T
450 MeV
240 MeV
200 MeV
New Rules: Let’s talk about absolute temperatures instead in Tc units
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
From Theory to Data
39
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Famous Plot
40
How can we relate these?
Quarkonium spectral functionsin equilibrated plasma Experimental data
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Famous Plot
40
How can we relate these?
Quarkonium spectral functionsin equilibrated plasma Experimental dataD
ynam
ical
mod
els
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Connection
41
- Modeling the motion of c-cbar in the evolving fireball - stochastic force from the heat bath
- attractive interaction between c-cbar
- lifetime of the plasma
- Spectral function calculation → no bound states only correlated c-cbar pairs
- Input: charm diffusion constant if small enough that attraction between c and cbar may survive
- What is the probability that c-cbar find themselves in proximity at the hadronization time?
illustration by Alex Doig
- Input: heavy quark potential
Zhao, Rapp, van Hees, Fries, Young, Shuryak, Strickland ...
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Connection
42
- Direct J/ψ suppressed, but ~ 50% of correlated c-cbar recombine- Coalescence gives relative small contribution- Quite good agreement with data for small charm diffusion and Tc=190 MeV
- statistical recombination c and cbar originate fromdifferent hard processes
- suppression + correlated regeneration
- total yield
Note: Tc = Tdeconfinement
Young, Shuryak, PRC 2010The 1st microscopic calculationof non-correlated recombination!
- Comparison to PHENIX data
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Connection
42
- Direct J/ψ suppressed, but ~ 50% of correlated c-cbar recombine- Coalescence gives relative small contribution- Quite good agreement with data for small charm diffusion and Tc=190 MeV
- statistical recombination c and cbar originate fromdifferent hard processes
- suppression + correlated regeneration
- total yield
Note: Tc = Tdeconfinement
Young, Shuryak, PRC 2010The 1st microscopic calculationof non-correlated recombination!
Outlook: Petreczky, Mocsy (BNL/Pratt) and Young, Schenke (McGill/BNL) using MUSIC
- Comparison to PHENIX data
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Connection
43
- Note There are effects not included in this model: initial state effects and absorption in the crossover-hadronic region (CNM effects)- A quantitative comparison with data is difficult
- statistical recombination c and cbar originate fromdifferent hard processes
- suppression + correlated regeneration
- total yield
Young, Shuryak, PRC 2010The 1st microscopic calculationof non-correlated recombination!
- Comparison to PHENIX data
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Isolate Hot Effects from Other Physics
44
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Can we isolate hot effects?
45
c and cbar are produced at
early times ... ... end up in hidden (J/ψ) or
in open (D) charm
... go through the entire evolution ...
illustration by Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 46
Can we isolate hot effects?
Cold Nuclear Matter
Initial: PDF’s modification (shadowing)Final: nuclear absorption
Hot Matter
screeninggluo-dissociation Landau damping threshold enhancement
Coalescence(regeneration)
coalescence of single c quarks in the plasma
Feed-down
J/ψ from decays:ψ’, χc → J/ψ χb, ϒ’, ϒ’’→ ϒ
keywords
Hadronic Absorption
illustration by Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 47
Can we isolate hot effects?new way of organizing
Deconfined Matter
Crossover/Hadronic Region
Initial State
illustration by Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Can we isolate hot effects?
48
“Anomalous suppression” With CNM effects divided out
This makes sense if - all CNM effects are initial state (shadowing), or - absorption in the crossover region (“mixed phase”) is similar to absorption in nuclear matter
M. Leitch, INT Quarkonium Workshop 2008
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Can we isolate hot effects?
49
“Anomalous suppression” With CNM effects divided out
- SPS described well with hadronic - reaches into crossover region - RHIC reaches into deconfined region onset of plasma effects (onset of screening) ~ 3.5 GeV/fm3 ~ 200 MeV
energy-density
M. Leitch, INT Quarkonium Workshop 2008
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Other Controls
50
Great for separating the different contributions At high pT
- CNM effects are less important: larger x- Statistical recombination has little effect
RAA versus pT
At RHIC: no suppression for J/y at high pT (~5 GeV) in 200GeV Cu+Cu and peripheral Au+Au collisions, but suppression at high pT in central Au At LHC suppression persists to higher pT
- A suppression at high pT would indicate suppression of direct J/ψ by the hot medium
- If no suppression then J/ψ forms outside of (or after) the hot plasma. Formation time!
- larger high pT suppression at LHC: can be from smaller x and/or longer lifetime
Zebo Tang (STAR), QM 2011 STAR CuCu: PRC80, 014922(R) PHENIX: PRL98, 232301
Silvestre (CMS), QM 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Other Controls
51
ϒpsilon
- Initial state effects not very relevant (mb >> Q2)- Absorption is small in the crossover/hadronic region- No recombination: number of b and bbar is negligible (at RHIC) - Easy to calculate spectral function, but dynamical modeling harder- Ground state can survive at RHIC and be suppressed at LHC ?!
- Y (theoretically) is a much cleaner signal
Reed (STAR), QM 2011
Silvestre (CMS), QM 2011
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Summary• The imaginary part of the potential plays a prominent role as a quarkonium dissolution mechanism. • Spectral functions determined with complex potential (with effects of color screening and dissipation) • Microscopic mechanism behind the melting of Jpsi and Upsilon could be quite different
• All signs point towards the melting of Jpsi early on in the plasma phase✴ Even the most binding potential allowed by lattice QCD leads 1S charmonium & excited bottomonium melting by T ≈ 240 MeV and of the 1S bottomonium states for T ≈ 450 MeV ✴ Threshold enhancement has phenomenological consequences
• Dynamical bridging necessary to compare to data - heavy quark diffusion, strength of Q-Qbar correlation, lifetime of deconfined medium
• High pT - clear hot matter signal ?! Y much clearer signal
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 53
The Sound of the Little Bangs
illustration: Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 53
The Sound of the Little Bangs
illustration: Alex Doig
The Rise and Fall of the Ridge in Heavy Ion Collisions - Sorensen, Bolliet, Mocsy, Pandit, Puthri, 2011Analyzing the Power Spectrum of the Little Bangs - Mocsy, Sorensen 2011The Sound of the Little Bangs - Mocsy, Sorensen 2010
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 54
Ágnes Mócsy with Alex Doig and Paul Sorensen
www.soundofthelittlebang.com
http://www.youtube.com/watch?v=jF8QO3Cou-Q&feature=player_embeddedYouTube Video (~ 4 min) :
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Media Attention
55
http://www.newscientist.com/article/dn18998-what-does-the-hottest-matter-ever-made-sound-like.html
Ágnes Mócsy, Pratt Institute, HARD PROBES 2010, Eilat, Israel
Analogy with the Early Universe
We determine the power-spectrum in heavy-ion collisions.
illustration: Alex Doig
Ágnes Mócsy, Pratt Institute, HARD PROBES 2010, Eilat, Israel
Analogy with the Early Universepicture credits NASA/WMAP Science Team
We determine the power-spectrum in heavy-ion collisions.
Quantum fluctuations from the start of the universe that show up as temperature-fluctuations (hotspots) in the CMB.
Power spectrum extracted from CMB.Most power is in l~200 corresponding to small distances.
CMB temperature map: fluctuations ~10-5
illustration: Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Analogy with the Early Universe
We determine the power-spectrum in heavy-ion collisions.
Quantum fluctuations from the start of the universe that show up as temperature-fluctuations (hotspots) in the CMB.
What is the analogy for HIC?
illustration: Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Analogy with the Early Universe
59
! ""#
/ $# r
ef (G
eV/c
)2
%"
0-5%
-2-1
01
2
024
-0.0010
0.0010.0020.0030.0040.0050.0060.007
Au+Au at 200 GeV
Adams, [STAR Collaboration]:J.Phys.G34:S679-684,2007
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Correlation Measurements
0.0
0.002
0.004
0.006
! ""#
/ $# r
ef (G
eV/c
)2
%"
0-5%
-2-1
01
2
024
-0.0010
0.0010.0020.0030.0040.0050.0060.007
€
Δρρref
=ρ − ρrefρref
(GeV/c)2
Peak indicates that particles with above average momentum tend to come out together, suggesting they are born out of the same high T lump
pt-pt correlations
illustration: Alex DoigSTAR Collaboration J. Phys. G 32, L37 (2006)
60
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Analogy with the Early Universe
61
! ""#
/ $# r
ef (G
eV/c
)2
%"
0-5%
-2-1
01
2
024
-0.0010
0.0010.0020.0030.0040.0050.0060.007
n0 2 4 6 8 100
0.002
0.004
]2 [(GeV/c)na
Au+Au at 200 GeV
Power Spectrum Adams, [STAR Collaboration]:J.Phys.G34:S679-684,2007
Mocsy, Sorenen, 2010
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Temperature Fluctuations in HIC
x [fm]-50
5
y [fm] -5
0
5
200
300
400
= 0.6 fm/c!Temperature [MeV]
x [fm]-50
5
y [fm] -5
0
5
170180190200210
= 4.6 fm/c!Temperature [MeV]
HIC analog of the CMB map
Au+Au collisions may initially contain hotspots of ~ 1.5 fm and remnants of these persist
Mocsy, Sorenen, 2010
62
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Length Scales and Power Spectrum How much of the initial inhomogeneity is transferred to the final state?
Spherical harmonic expansion of CMB
Fourier expansion of HIC
Higher harmonics probes smaller length-scales.
Efficiency of conversion depends on relation of lmfp to the scale probed at n:
Power spectrum tells the strength transferred into each harmonic n -- ? get lmfp
n=2 n=3 n=4 n=10 n=15
<R> average radial position of participants
sum l
63
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Length Scales and Power Spectrum How much of the initial inhomogeneity is transferred to the final state?
lmfp
Spherical harmonic expansion of CMB
Fourier expansion of HIC
Higher harmonics probes smaller length-scales.
Efficiency of conversion depends on relation of lmfp to the scale probed at n:
Power spectrum tells the strength transferred into each harmonic n -- ? get lmfp
n=2 n=3 n=4 n=10 n=15
<R> average radial position of participants
sum l
63
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Length Scales and Power Spectrum How much of the initial inhomogeneity is transferred to the final state?
lmfp
Spherical harmonic expansion of CMB
Fourier expansion of HIC
Higher harmonics probes smaller length-scales.
Efficiency of conversion depends on relation of lmfp to the scale probed at n:
Power spectrum tells the strength transferred into each harmonic n -- ? get lmfp
n=2 n=3 n=4 n=10 n=15
lmfp
<R> average radial position of participants
sum l
63
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Acoustic Horizon
64
Horizon means long wavelengths cannot be heard at first. Time-ordering of frequencies
Alex Doig 2010Alex Doig
another important lengthscale H = distance density perturbations can travel
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Acoustic Horizon
65
Horizon means long wavelengths cannot be heard at first. Time-ordering of frequencies
Alex Doig 2010Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Acoustic Horizon
66
Horizon means long wavelengths cannot be heard at first. Time-ordering of frequencies
Alex Doig 2010Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
Acoustic Horizon
67
Horizon means long wavelengths cannot be heard at first. Time-ordering of frequencies
Alex Doig 2010Alex Doig
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico
The Sound of the Little Bang
www.soundofthelittlebang.comillustration: Alex Doig
68
Ágnes Mócsy, XQCD July 18-20 2011, San Carlos, Mexico 69
The End