![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-16314377-1d30-402f-aeb6-60ff9a8d1391.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-c338a33d-8348-4f14-823b-267689f55db1.jpg)

Nuclear Physics A

Nuclear Physics A 00 (2014) 1-4

## Properties of Thermal Photons at RHIC and LHC

## R. Rapp a , H. van Hees b , M. He c

a Cyclotron Institute and Department of Physics and Astronomy, Texas A &amp; MUniversity, College Station, TX 77843, USA b FIAS, Ruth-Moufang Str. 1 and Institute for Theoretical Physics, Max-von-Laue-Str. 1, D-60438 Frankfurt, Germany c Department of Applied Physics, Nanjing University of Science and Technology, Nanjing 210094, China

## Abstract

We study the emission characteristics of thermal photons at RHIC and LHC as a GLYPH&lt;11&gt; ected by both the space-time evolution of the bulk medium and the thermal emission rates. For the former we compare the results of two evolution models (expanding fireball and hydrodynamics). For the latter, we detail the influence of hadronic emission components and study a speculative scenario by upscaling the default QGP and hadronic rates around the pseudo-critical region.

Keywords: QCD Matter, Ultrarelativistic Heavy-Ion Collisions, Thermal Photon Emission

## 1. Introduction

The recent observations of a large enhancement of direct photons in heavy-ion collisions at RHIC [1, 2] and LHC [3] are indicative for thermal radiation o GLYPH&lt;11&gt; an interacting fireball of QCD matter formed in these reactions. It furthermore turns out that these direct photons carry a remarkable elliptic flow ( v 2 ) [4, 5], suggestive for prevalent emission from later stages in the fireball evolution, when most of the v 2 of the bulk has already built up. Typically, this takes about 5 fm / c by which time the system (in semi-central Au-Au collisions at RHIC) has cooled down to near the (pseudo-) critical temperature, T pc ' 170 MeV, cf. Fig. 1 left. Assuming an average transverse-flow velocity of ¯ GLYPH&lt;12&gt; ' 0.4 c around this region yields a schematic estimate of the inverse-slope parameter of such a source of T e GLYPH&lt;11&gt; ' T pc p (1 + h vT i ) = (1 GLYPH&lt;0&gt; h vT i ) ' 255 MeV, which is within the range of values measured by PHENIX [1, 2], see Fig. 1 [6]. Such a source is compatible with what has been deduced from other observables, e.g., for the quark distribution functions used in phenomenological coalescence models to reproduce the constituent-quark scaling of light and multistrange hadron spectra and v 2 , or for thermal dileptons, where large emission contributions from around T pc [7, 8] explain the apparent dissolution of the GLYPH&lt;26&gt; resonance from SPS to RHIC energies. Especially the latter is closely related to photons which correspond to the M ! 0 limit of dileptons. This link has been explored in Ref. [6]; using a schematic fireball model for the bulk evolution with state-of-the-art photon rates, the direct-photon v 2 turns out to be not far from the PHENIX data. In particular the hadronic rates are larger than in other calculations in the literature. Most notably, the M ! 0 limit of the GLYPH&lt;26&gt; spectral function [9] includes a rather extensive set of baryoninduced processes, known to be of critical importance for the low-mass dilepton enhancement [7]. In the following, we scrutinize these findings by testing various ingredients of the thermal photon calculations, both with respect to the space-time evolution (utilizing hydrodynamics with initial flow and sequential chemical and kinetic freezeout) and the photon rates (switching o GLYPH&lt;11&gt; the GLYPH&lt;26&gt; spectral function contribution, and upscaling current QGP and hadronic rates around T pc ) [10].

Figure 1. Left panel: time dependence of inverse slope parameter of the photon qt spectrum (solid blue line), rest frame temperature (red dashed line) and average transverse expansion velocity (dash-dotted line) from a fireball evolution, compared to the measured inverse slope of the excess radiation in 0-20% Au-Au( p sNN = 200 GeV), T slope = (239 GLYPH&lt;6&gt; 26) MeV [2] (pink horizontal lines). Right panel: thermal photon production rates from the QGP (solid blue line: full leading order [11], dash-dotted line: Born rate [13]) and hadronic matter (red solid and dashed line: GLYPH&lt;26&gt; spectral function contribution with and without baryons; long-dashed line: tree-level meson gas reactions [9]).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-23c91c2d-84c2-42df-81e8-0ca718a71287.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-baa7ef31-464b-4787-95b5-8c8d13ed6e5e.jpg)

## 2. Thermal Emission Rates

In thermal field theory the equilibrium emission rate of photons,

$$q _ { 0 } \frac { d N _ { \gamma } } { d ^ { 4 } x d ^ { 3 } q } & = - \frac { \alpha _ { \text{EM} } } { \pi ^ { 2 } } \ f ^ { \text{B} } ( q _ { 0 } ; T ) \text{Im} \, \Pi _ { \text{EM} } ^ { T } ( q _ { 0 } = q ; \mu _ { B }, T ),$$

is determined by the in-medium photon selfenergy, GLYPH&lt;5&gt; T EM , due to the coupling to the partonic or hadronic constituents of the thermal bath. Its imaginary part represents the di GLYPH&lt;11&gt; erent cuts of the corresponding Feynman diagrams which characterize the scattering processes for thermal photon production.

For photon emission from the QGP we employ the complete leading-order results of Ref. [11]. For emission from hadronic matter we employ the zero-mass limit of the GLYPH&lt;26&gt; spectral function of Ref. [12] which includes baryon and meson resonances resulting in pertinent Dalitz decays, as well as pion-exchange reactions with baryons generated by the in-medium pion cloud; in addition, pion, kaon, a 1 and ! t -channel exchange reactions in a GLYPH&lt;25&gt;GLYPH&lt;26&gt; KK GLYPH&lt;3&gt; gas (not included in the GLYPH&lt;26&gt; spectral function) have been added through pertinent Born diagrams using the standard kinetic theory expression [9]. These rates are summarized for a temperature of T = 170 MeV and vanishing baryon chemical potential ( GLYPH&lt;22&gt; B = 0) in Fig. 1 right. Contributions from anti- baryons in hadronic matter are significant up to rest-frame momenta / of q 0 GLYPH&lt;20&gt; 1 GeV; for typical hadronic flow velocities at RHIC and LHC, this translates into lab-frame momenta of up to GLYPH&lt;24&gt; 2 GeV. Meson gas t -channel exchange reactions take over for q 0 &gt; 1 GeV. The sum of all hadronic processes is comparable to the perturbative QGP rates.

## 3. Space-Time Evolution

To model the space-time evolution of the thermal medium, over which the rates have to be integrated, we have considered two approaches: (i) an expanding fireball with ellipsoidal transverse area, assumed to be isotropic at each time step and with conserved total entropy [6]; (ii) a (2 + 1)-D ideal hydrodynamic model [15] based on the original code of Ref. [14]. Both models utilize the same equation of state (lattice EoS for the QGP with a near-smooth transition to a hadron resonance gas at T pc = 170 MeV, followed by hadro-chemical freezeout at T ch = 160MeV), and have been adjusted to the same set of light and multistrange hadron spectra and v 2 . Kinetic freezeout for multistrange hadrons has been imposed at T ch , facilitated by a suitable transverse acceleration in the fireball model, and by a non-vanishing initial flow in the hydro model. For both models this implies that the bulk-medium v 2 levels o GLYPH&lt;11&gt; close to T ch , which, in turn, also helps in a better description of the pion and proton v 2 at kinetic freezeout at T kin ' 110 MeV.

Figure 2. Left: temperature emission profile of thermal photon yields from fireball and hydro evolutions at qT = 2 GeV. Right: reduction of the total direct photon v 2 from the fireball when switching o GLYPH&lt;11&gt; the contributions from either the in-medium GLYPH&lt;26&gt; spectral function or strong final-state decays.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-281933ac-fe5c-4834-b3e5-490339b71b4d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-4d688249-475c-4ac0-89e2-f7bb9b2cda94.jpg)

## 4. Photon Spectra and Flow

A detailed comparison of the photon spectra from the fireball and hydro evolution (using the same default rates, as plotted in Fig. 1) reveals a reasonable agreement (within ca. 20%) of the QGP spectra; both space-time models also exhibit a pronounced maximum of the temperature di GLYPH&lt;11&gt; erential emission around T pc up to momenta of at least qT = 2 GeV, see Fig. 2 left. This is a consequence of the remnant of the latent heat in the transition region, i.e., the large variation in entropy density over a small range in T (similar to low-mass dileptons). However, in the hadronic phase, the hydro emission falls significantly below the one from the fireball, more so at higher momenta (up to ca. 50% at qT = 2GeV in Au-Au at RHIC). Part of the reason is the continuous freeze-out of the hydro medium, while the entire fireball volume only freezes out at the end of the evolution. While the latter is clearly an over-simplification, one should note that hadrons which are frozen out on the hydro hypersurface never re-thermalize, although, in principle, they could be reabsorbed by the freeze-out front. Therefore, we expect the most realistic hadronic emission in between the hydro and fireball results. We believe that another part of the discrepancy in the hadronic photon emission spectra, which becomes more pronounced toward higher qT , is the 'inward' burning of the freeze-out front. This leads to rather low-flow hadronic freeze-out cells especially toward the end of the evolution. On the contrary, the outward moving fireball front implies a monotonously increasing average flow (cf. Fig. 1 left). The 'inward burning' of the hydro was identified as a possible problem in reproducing the measured HBT radii of pions. The initial flow introduced into our hydro evolution helps with this, but possibly not enough.

As a result of the discrepancy in the hadronic emission, the total direct-photon spectra from the hydro evolution (plus primordial photons) in 0-20% Au-Au at RHIC are below the fireball results by a few tens of percent in the region where thermal radiation is prevalent. Likewise, the maximum in the photon v 2 is about 25% smaller for the hydro evolution. Yet our hydro results for the photon v 2 are not far from the lower limit of the error bars of the PHENIX data, especially when allowing for moderate modifications of the primordial contribution, see Fig. 3 right. This constitutes an improvement over existing calculations. We attribute the main di GLYPH&lt;11&gt; erences to the larger hadronic emission rates, the initial flow in the hydro evolution, and the strong-decay feeddown (e.g., GLYPH&lt;1&gt; ! N GLYPH&lt;13&gt; , a 1 ! GLYPH&lt;25&gt;GLYPH&lt;13&gt; , etc.) in our calculation. For example, when switching o GLYPH&lt;11&gt; the (hadronic) photon contributions from the in-medium GLYPH&lt;26&gt; spectral function, the maximum total direct-photon v 2 in the fireball calculation decreases by about 35%, while switching o GLYPH&lt;11&gt; the strong feeddown reduces it by nearly 10%, cf. Fig. 2 right. At the LHC, the stronger QGP dominance renders closer agreement between hydro and fireball results. Both calculations tend to underestimate the preliminary ALICE qT spectra, but agree with the preliminary v 2 data within experimental errors.

Since the default thermal emission rates do not fully account for the observed direct-photon enhancement, we investigated an extreme scenario of upscaling our rates by a factor reaching a maximum of 3 around T pc ('pseudocritical enhancement'), and convoluted these rates over the hydro medium. This basically resolves the discrepancies with both RHIC and LHC spectra and v 2 , see, e.g., the dash-dotted lines in Fig. 3. However, it is presently un-

Figure 3. Direct photon spectra (left) and elliptic flow (right) from the hydrodynamic evolution using the default emission rates with the PHENIX pp reference for primordial emission (solid lines), with a pQCD-motivated primordial emission (dashed line in the right panel), and with 'pseudocritically' enhanced rates (see text for details) with pQCD primordial (dash-dotted lines). PHENIX data are from Refs. [1, 4].

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-a8528cd1-989a-4505-9fb6-3ad84b778f24.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0612v2-figure-d7239326-76a2-4ed5-8026-25352ace8e8b.jpg)

clear whether such an enhancement can be substantiated by microscopic rate calculations, and, if confirmed, what ramifications this would have for dilepton observables (where current theory and data agree [8]).

## 5. Conclusions

We have computed thermal photon spectra and v 2 in heavy-ion collisions at RHIC and LHC, using state-of-the-art thermal emission rates from QGP and hadronic matter including the e GLYPH&lt;11&gt; ect of baryons which importantly figure in the dilepton context. We have compared two di GLYPH&lt;11&gt; erent evolution models based on the same EoS and fits to bulk-hadron data, as well as the concept of sequential freezeout. Their main di GLYPH&lt;11&gt; erence appears in the hadronic emission spectra, which are larger in the thermal fireball than in the hydro evolution. Consequently, the fireball does slightly better in describing the experimental excess spectra and v 2 of direct photons at RHIC, even though the hydro results are not far from the data either given their current uncertainties. At LHC, both models tend to underestimate the spectra but are compatible with the v 2 from preliminary ALICE measurements. Important ingredients in our approach are large hadronic rates and a rapid built-up of radial flow (necessary to implement sequential freeze-out). An upscaling of the thermal rates around T pc further improves the agreement, but a microscopic calculation of such an e GLYPH&lt;11&gt; ect remains an open issue.

Acknowledgemts . This work has been supported by the grants U.S.-NSF PHY-1306359, NSFC 11305089, BMBF 05P12RFFTS, and by LOEWE through HIC for FAIR.

## References

- [1] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 104 (2010) 132301.
- [2] A. Adare et al. [PHENIX Collaboration], arXiv:1405.3940 [nucl-ex].
- [3] M. Wilde [ALICE Collaboration], Nucl. Phys. A 904-905 , 573c (2013).
- [4] A. Adare et al. [PHENIX Collaboration], Phys. Rev. Lett. 109 , 122302 (2012).
- [5] D. Lohner et al. [ALICE Collaboration], J. Phys. Conf. Ser. 446 , 012028 (2013).
- [6] H. van Hees, C. Gale and R. Rapp, Phys. Rev. C 84 , 054906 (2011).
- [7] R. Rapp, J. Wambach and H. van Hees, in Relativistic Heavy-Ion Physics , edited by R. Stock and Landolt B¨ ornstein (Springer), New Series I 23A / , 4-1 (2010); [ arXiv:0901.3289[hep-ph] ].
- [8] R. Rapp, Adv. High Energy Phys. 2013 (2013) 148253.
- [9] S. Turbide, R. Rapp and C. Gale, Phys. Rev. C 69 , 014903 (2004).
- [10] H. van Hees, M. He and R. Rapp, submitted to Nucl. Phys. A (2014); arXiv:1404.2846 [nucl-th].
- [11] P.B. Arnold, G.D. Moore and L.G. Ya GLYPH&lt;11&gt; e, JHEP 0112 , 009 (2001).
- [12] R. Rapp and J. Wambach, Eur. Phys. J. A 6 (1999) 415.
- [13] J.I. Kapusta, P. Lichard and D. Seibert, Phys. Rev. D 44 (1991) 2774 [Erratum-ibid. D 47 (1993) 4171].

[14]

P.F. Kolb and U.W. Heinz, In R. C. Hwa (ed.) et al.: Quark gluon plasma, 634, [nucl-th 0305084].

/

- [15] M. He, R.J. Fries and R. Rapp, Phys. Rev. C 85 , 044911 (2012).
R. Rapp H. van Hees M. He
Properties of Thermal Photons at RHIC and LHC
We study the emission characteristics of thermal photons at RHIC and LHC as
affected by both the space-time evolution of the bulk medium and the thermal
emission rates. For the former we compare the results of two evolution models
(expanding fireball and hydrodynamics). For the latter, we detail the influence
of hadronic emission components and study a speculative scenario by upscaling
the default QGP and hadronic rates around the pseudo-critical region.