Hindawi Publishing Corporation Advances in High Energy Physics Volume 2013, Article ID 176578, 41 pages http://dx.doi.org/10.1155/2013/176578

Review Article

Electromagnetic Radiations from Heavy Ion Collision

Payal Mohanty,1,2 Sabyasachi Ghosh,2,3 and Sukanya Mitra2

1 HENPP Division, Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Kolkata 700064, India

2 Theoretical Physics Division, Variable Energy Cyclotron Centre, 1/AF Bidhan Nagar, Kolkata 700064, India

3 Ctr. for Astroparticle Physics and Space Bose Institute, 5/EN Bidhan Nagar, Kolkata 700091, India

Correspondence should be addressed to Payal Mohanty; payal.mohanty@saha.ac.in Received 15 April 2013; Accepted 5 July 2013 Academic Editor: Jan E. Alam

Copyright © 2013 Payal Mohanty et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In this review, we have discussed the different sources of photons and dileptons produced in heavy ion collision (HIC). The transverse momentum (pT) spectra of photons for different collision energies are analyzed with a view of extracting the thermal properties of the system formed in HIC. We showed the effect of viscosity on pT spectra of produced thermal photons. The dilepton productions from hot hadrons are considered including the spectral change of light vector mesons in the thermal bath. We have analyzed the pT and invariant mass (M) spectra of dileptons for different collision energies too. As the individual spectra are constrained by certain unambiguous hydrodynamical inputs, so we evaluated the ratio of photon to dilepton spectra, Rem, to overcome those quantities. We argue that the variation of the radial velocity extracted from Rem with M is indicative of a phase transition from the initially produced partons to hadrons. In the calculations of interferometry involving dilepton pairs, it is argued that the nonmonotonic variation of HBT radii with invariant mass of the lepton pairs signals the formation of quark gluon plasma in HIC. Elliptic flow (v2) of dilepton is also studied at ^sNN = 2.76 TeV for 30-40% centrality using the (2 + 1 )d hydrodynamical model.

1. Introduction

The main objective of relativistic heavy ion collisions is to study the transient phase, that is, quark gluon plasma (QGP) which is believed to permeate the early universe a few microseconds after the Big Bang. Collision between nuclei at ultrarelativistic energies produces charged particles either in hadronic or in partonic state depending on the collision energy. Interaction among these charged particles produces electromagnetic (EM) radiation [1-9]. However, hadrons being strongly interacting objects give snapshot of evolution only from the freezeout surface. So they have hardly any information about the interior of the plasma. Whereas EM radiations, for example, the thermal photons and dileptons, are expected to provide an accurate information about the initial condition and the history of evolution of the plasma. This is possible since photons and dileptons interact only through the EM interaction. The EM interaction strength is small compared to that of strong interaction (a < as) and thus dominates the dynamics of nuclear collision processes. Therefore, its mean free path (A = 1/n<r) is larger than

the size of the system. Because of their negligible finalstate interactions with the hadronic environment, once produced it brings the electromagnetic particles about to escape unscathed carrying the clean information of all stages of the collision. The EM radiations produce all stages of collision process which contribute to the measured photon spectra; in principle, a careful analysis may be useful to uncover the whole space-time history of nuclear collision. Hence EM radiations—real and the virtual photons (dilepton)—are considered as efficient probes to study dynamical evolution of the matter formed in relativistic heavy ion collision. However, as they are emitted continuously, they sense in fact the entire space-time history of the reaction. This expectation has led to an intense and concerted efforts toward the identification of various sources of such radiations. While initially these signals were treated as thermometer of the dense medium created, but later on recent calculations suggest it might serve as chronometer [10] and flow-meter [11-16] ofHIC.

The review is organized as follows. In Section 2, we start with possible sources of photons and dileptons that were produced in HIC. We have discussed the formalism of static

emission rate of photons and dileptons in Section 3. To get total yield, we need concept of hydrodynamics. So, we briefly outlined relativistic hydrodynamics in Section 4 which takes care of the evolution. In Section 5 we have presented the thermal emission rate of photons from QGP (Section 5.1) and hadronic matter (Section 5.2) which is used to produce the results, and total invariant yield of direct photon for different collision energy is shown in Section 5.4. The effect of viscosity on the transverse momentum (pT) spectra of photon is discussed in Section 5.5. Similarly, the details of the emission rate of dileptons from QGP and hot hadrons are given in Sections 6.1 and 6.2, respectively. Using these rates, the results of pT and invariant mass (M) spectra of dileptons are presented in Section 6.4. In Section 7 the radial flow is extracted by simultaneous use of pT spectra of photons and dileptons and ratio of the spectra, and Rem is conferred. The correlation function for dilepton has been calculated and HBT radii are extracted as function of M in Section 8. We have also evaluated the dilepton v2 in Section 9 taking into account the medium effect on the spectral function of the vector mesons. Finally, we have summarized the work in Section 10.

2. Various Sources of EM Radiations

As argued previously that EM radiations emerge out copiously from all stages of collision, so, in order to proceed, it is useful to identify various sources of photons and dileptons produced in the HIC. So the "inclusive" photon spectrum coming from such collision in usual sense can be defined as the unbiased photon spectrum observed in pp, pA, or AA collision. This spectrum is built up from a cocktail of various components.

Depending on their origin, there are two different types of sources which are "direct photons" and "photons from decay of hadrons." The term "direct photons" is meant for those photons and dileptons which produce directly from collision between the particles. One can subdivide this broad category of "direct photons" into "prompt photons," "preequilibrium photons" and "thermal photons," depending on their origin. On the other hand, the decay photons do not come directly from the collision, rather from the decay of hadrons.

2.1. Transverse Momentum (pT) Dependence of EM Radiations. The EM spectra provided by the experimentalist are mingled with various sources of photons and dileptons and it is difficult to distinguish different sources experimentally. However, real interest lies in the thermal photons and dileptons since it is expected to render information about the initial condition and the history of evolution of the plasma while it cools and hadronizes. Thus, theoretical models are used with great advantage to identify these sources of photons and their relative importance and characteristics in the spectrum.

Depending on the process through which photons/dilep-tons produce, they are categorized as follows.

(1) Prompt: the EM radiations produced by hard scattering of the partons inside the nucleons of incoming

nuclei in the initial stage of collision, before the thermalization sets in, are known as prompt photons and dileptons (Drell Yan). This contribution may be evaluated by using pQCD.

(2) Preequilibrium: the preequilibrium photons and dilep-tons are produced in the preequilibrium stage, that is, before the thermalization sets in the system. In such scenario the contribution from preequilibrium stage will be very small and hence neglected.

(3) Thermal: EM radiations which are emitted from the thermalized systems of quarks and gluons or hadronic gas.

(4) Decay: after the freezeout of the fireball, photons and dileptons are also produced from the decays of long-lived (compared to strong interaction time scale) hadrons and known as "photons from decay."

Out of different sources, the thermal photons and dilep-tons are privileged as they carry information about the formation of QGP. As indicated in Figure 1, the hard photons dominate the high pT part of the invariant momentum spectra, and decay photon populates the low pT part and rest over thermal contribution shines in the intermediate domain of the pT spectra ~1-3GeV. And the calculations based on theory infer that the photons and lepton pairs from hadronic matter dominate the spectrum at lower pT (~1-

2 GeV) whereas photons and dileptons form QGP dominate in the intermediate pT range, that is, pT ~ 2-3 GeV (depending on the models) [17]. This small window may help in learning the properties of QGP. Thus one has to subtract out the nonthermal sources to understand the properties of the QGP. However, it is not possible experimentallyto distinguish between different sources. Thus, theoretical models and calculations can be used to great advantage to identify different sources of direct photons and their relative importance and characteristics in the spectrum. The hard photons and dilep-tons are well understood in the framework of pQCD, and decay contributions can be filtered out experimentally using different subtraction methods, like invariant mass analysis, mixed event analysis, internal conversion method, and so forth.

The invariant momentum distribution of photons and dileptons produce from a thermal source depends on the temperature (T) of the source through the thermal phase space distributions of the participants of the reaction that produces the photons and dileptons [18]. As a result the pT spectra of thermal photons and dileptons reflect the temperature of the source through the phase space Hence ideally the photons with intermediate pT values (~2-

3 GeV, depending on the value of initial temperature) reflect the properties of QGP (realized when T > Tc, Tc is the transition temperature). Therefore, one should look into the pT spectra for these values of pT for the detection of QGP. However, for an expanding system the situation is far more complex. The thermal phase space factor changes by flow; for example, the transverse kick received by low pT photons due to flow originating from the low temperature hadronic phase (realized when T < Tc) populates the high pT part of the

Figure 1: Schematic diagram of different sources of photon production in heavy ion collision as function of transverse momentum (Pt)-

; I I I I | I I I I | I I I - n°, ^ Dalitz-decays I | I I I I | I I I I ;

' v /\P,w

Vf : \ \

\ \ 1 f ^ L '■\.DD \ J/V 1

_ V_ \ — \ \ \ V 1 Drell-Yan

r Low- : Intermediate! >10fm >1 fm r , , , , : 1 , , , , 1 , , , High-mass region n <0.1 fm 1 .... 1 ....

Mass (GeV/c2)

Figure 2: Expected different sources of dilepton production in HIC as function of invariant mass [5].

spectra [19]. As a consequence the intermediate or the high pT part of the spectra contains contributions from both QGP and hadrons. Thus, it is not an easy task to disentangle the photons coming from pure partonic phase. Thus photons appear to be a more restrictive probe since they are characterized only by their momentum whereas the dileptons have two kinematic variable, pT and invariant mass (M) to play with. A soft photon (low pT) in one frame of reference can be hard (high pT) in another frame, whereas the pT integrated invariant mass distribution of dileptons is independent of any frame. In addition to it the pT spectra are affected by the flow; however, the pT integrated M spectra remain unaltered by the flow in the system. Also in the M spectra of dileptons, again in M spectra, the dileptons from QGP dominates over its hadronic counterpart above the 0 peak. All these suggests that a judicious choice of pT and M windows will be very useful to characterize the QGP and hadronic phase separately.

2.2. Invariant Mass (M) Dependence of EM Radiations. Being massive, dileptons make situation different from photons. They have two kinematic variables—pT and M. As argued before, the pT spectra are affected by the flow; however, the _pT-integrated M spectra remain unaltered by the flow in the system. It should be mentioned here that for M below p peak and above <p peak dileptons from QGP dominates over its hadronic counterpart (assuming the contributions from hadronic cocktails are subtracted out) if the medium effect of spectral function of the low mass vector mesons are not taken into account. However, the spectral function of low mass vector mesons (mainly p) may shift toward lower invariant mass region due to nonzero temperature and density effects. As a consequence the contributions from the decays of p mesons to lepton pairs could populate the low M window and may dominate over the contributions from the QGP phase [5, 8, 20]. All these suggest that the invariant mass distribution of dilepton can be used as a clock for HIC, and a

judicious choice of pT and M windows will be very useful to characterize the flow in QGP and hadronic phase.

The measured dilepton spectra can be divided into several phases. Depending on the invariant mass of the emitted dileptons, it can be classified into three distinct regimes (discussed below [5]), and a schematic diagram of dilepton mass distribution is shown in Figure 2.

(i) High mass region (HMR):

(M>Mr/w (=3.1 GeV), pT -3-5 GeV). (1)

The HMR region corresponds to early preequilibrium phase (r < Tj), where the lepton pairs are produced with large invariant mass (M > 3 GeV) and the dominant contributions are from the hard scattering between the partons, like Drell Yan annihilation [21, 22]. The final abundance of the heavy quarkonia (//¥, Y) and their contribution to the spectrum is suppressed due to the Debye screening and as a result the bound states are dissolved.

(ii) Intermediate mass region (IMR):

(M0 <M< MJ/W, pT -1-3 GeV) . (2)

Thermalization is achieved in the system after a time scale (> Tj). In this domain, the dileptons from the QGP are produced via quark-antiquark annihilation. In this regime, due to higher temperature the continuum radiation from QGP dominates the dilepton mass spectrum and thus this region is important for the detection of QGP. The decays of "open charm" mesons, that is, pairwise produced DD mesons [23] followed by semileptonic decays, contribute largely in this domain of M. Although an enhanced charm production is interesting in itself—probably related to the very early collision states—it may easily mask the thermal plasma

signal. To somewhat lesser extent, this also holds true for the lower-mass tail of Drell-Yan production [21, 22]. If the heavy quark does not get thermalized, then their contribution maybe estimated from pp collision data with the inclusion of nuclear effects like shadowing, and so forth, and they do not contribute to the flow also [24].

(iii) Low mass region (LMR):

(M < M^ (_ 1.02 GeV) ,pT <1 GeV) .

With subsequent expansion and cooling, the QGP converts into a hot hadron gas at the transition temperature, Tc. At later stages, the dileptons are preferentially radiated from hot hadron gas from the decay of (light) vector meson, such as the p, to, and 0. The low M domain of the lepton pairs is dominated by the decays of p. Medium modification of p will change the yield in this domain of M. The change of p spectral function is connected with the chiral symmetry in the bath; therefore, the measurement of low M lepton pairs has great importance to study the chiral symmetry restoration [25, 26] at high temperature and density. Thus the invariant mass of the lepton pair directly reflects the mass distribution of the light vector mesons. This explains the distinguished role that vector mesons in conjunction with their in-medium modifications play for dilepton measurements in HIC.

So far, we have discussed the different sources of photons and dileptons. As QGP is expected to form in the HIC experiments, the basic intention of the present study is to study the properties of QGP. Therefore, we have emphasized more on the study of thermal photons and dileptons in this review, as they may provide information to understand the formation and unique properties of the novel matter.

The emission of thermal photons and dileptons coming from HIC consists of two important segments:

dN _ J dR

d4p _ J d4p

(1) firstly, static emission rate (dR/d4p) which takes care of the basic interactions in respective phases (QGP or hadronic phase),

(2) Secondly, the space-time integration over four volume (d4x) which takes care of the evolution of the thermal matter created in HIC. As the EM radiations produced from each space-time point of the evolving matter, we need the concept of relativistic hydrodynamics (described in Section 4) for understanding the evolution.

3. Formulation of Thermal Emission Rate of EM Radiations

The importance of the electromagnetic probes for the study of thermodynamic state of the evolving matter was first proposed by Feinberg in 1976 [27]. Feinberg showed that the emission rates can be related to the electromagnetic current-current correlation function in a thermalized system.

3.1. Dilepton Emission Rate from Thermal Medium. Let us consider an initial state |7) which goes to a final state |F) producing a lepton pair l+T with momenta p1 and p2, respectively. The dilepton multiplicity thermally averaged over initial states is given by [4, 28]

N_H\{F>l+l-

i\ Lintd4

e'M d3p1 d3p2

(2n)32E1 (2n)32E2 '

where Lint = efl(x)y.ifi(x)A^(x) + e] (x)A^(x) in which Vi(x) is the lepton field operator and J^(x) is the electromagnetic current and Z = Tr[e-^H]. Following [1, 4, 8] this expression can be put in the form

dN d4qd4:

:L(M2)fBE (q0)g>"(q0,q), (6)

where the factor L(M2) = (1 + 2m2l/M2)(1 - Amj/M2)1'2 is of the order of unity for electrons, M(= ^q2) being the invariant mass of the pair, and the electromagnetic (e.m.) current correlator is defined by

W,v (%, q) = \ d4xe'qx ([j™ (x), C (0)]> . (7)

Here J?(x) is the electromagnetic current and {) indicates ensemble average. The rate given by (6) is to leading order in electromagnetic interactions but exact to all orders in the strong coupling encoded in the current correlator . The q in the denominator indicates the exchange of a single virtual photon and the Bose distribution implies the thermal weight of the source. We can also express the dilepton rate in terms of a photon spectral function . Using the relation [4],

AnaW^ = 2n [q2g^a - q^qa) A°f (-q) {q2gpv - qpqv)

in (6), we have dN

d4qd4x (■■■) q4

where A^{q0,q) _ Im ^m^.

f be (qo) im n™ (qo)'

3.2. Photon Emission Rate from Thermal Medium. The photon emission rate is calculated in the similar way to that of dilepton rate. The photon emission rate differs from the dilep-ton rate in the following way: the factor e2L /q4 appearing in the dilepton rate (in (9)) which is nothing but the product of electromagnetic vertex y* ^ l+T, the leptonic current involving Dirac spinors, and the square of the photon propagator should be replaced by the factor £ e^ev (= ). And the phase space factor d3p1/(2n) E1d3p2/(2n) E2 should be replaced by d3 p/(2n) E. Then the photon emission rate becomes

qo ^ = r /be (qo) Im n

The above emission rate is correct up to order e2 in electromagnetic interaction but exact, in principle, to all order in strong interaction. However, for all practical purposes, one is able to evaluate up to a finite order of loop expansion. Now it is clear from the above results that to evaluate photon and dilepton emission rate from a thermal system we need to evaluate the imaginary part of the photon self-energy. The Cutkosky rules at finite temperature or the thermal cutting rules [29-32] give a systematic procedure to calculate the imaginary part of a Feynman diagram. The Cutkosky rule expresses the imaginary part of the n-loop amplitude in terms of physical amplitude of lower order (n - 1 loop or lower). This is shown schematically in Figure 3. When the imaginary part of the self-energy is calculated up to and including L order loops where L satisfies x + y < L + 1, then one obtains the photon emission rate for the reaction x particles ^ y particles +y, and the above formalism becomes equivalent to the relativistic kinetic theory formalism [2, 3].

3.3. Emission Rate Using Relativistic Kinetic Theory Formalism. According to relativistic kinetic theory formulation, the production of «-type particles from the reaction of type 1(p1) + 2(p2) ^ 3(p3)+4(p) is given as follows:

d3p1 d3p2 d3p3 d3p (2n)32E1 (2n)32E2 (2n)32E3 (2n)32E

x(2n)4S(4) (Pl +P2 -P3 -p)\M\2 X fj 2 (l±f3)(l±fi),

where N is the overall degeneracy for the reaction under consideration, \.M\2 is the square of the invariant amplitude for the process under consideration, pt, Et, and fi(Et) are the three momentum, energy, and thermal distribution functions (Fermi-Dirac or Bose-Einstein) of the incoming and outgoing particles "«."

The transverse momentum (pT) distribution of photons from a reaction of the type: 1 + 2 ^ 3 + y taking place in a thermal bath at a temperature, T is given by [2, 3]:

,dR = N J d3p! d3p2 <£p1 Jd3p 2(2n)8 J 2EX h!2 ( ± J3)

2E2 2E3

xS(4) (Pl +P2 -p3 -P)\.M\2.

Using the Mandelstam variables (s,t,u) we can write the differential photon production rate as [33]

d3p \6(2n)7E .

rœ rt,

I ,ds\

J(m1+m 2 ) Jtn

, t, u)

\dEl \

f(Ei )f(E2)[l + f(E3)]^ ^aE'2 + 2bE2 + c

yy /yy yy \

wwiQww*/ = Im + <wwv<^)ww\/ + ••• ) =

Figure 3: Optical Theorem in Quantum Field Theory.

/ 2 2\2 a = -(s + t - m2 - m3) ,

b = E1 (s + t - m"2 - m^) (m>2 - t)

+ E[(s + t - m2 - m2 )(s - m - m ) - 2mx (m\ - i)],

2i 2 \2 2\( 2 2\2 2 2~l c = -E1 (m2 -t) -E \(s - m1 - m2) - 4m1m2]

2 i 2 2\2 2(2 \ 2 + m2(s + t-m2 - m3) + m1(m2 -t) - 2E1E

x [2m2 (s + t - m"2 - m^) - (m\ -t)(s-m^ - m2)]

- (s + t - m^ - m"2) (m2 - t)(s -m2 - m,

s + t - m

Em^ m^ -t

s + t - m2 - m2

v = _^ +

min m2 -t+ 4E '

b Vb2 - ac

= -- +

In a similar way the dilepton emission rate for a reaction aa ^ l+T can be obtained as

^ = N J. d4p N J

3pa d3pâ d3pi d3p2

d4p " J 2Ea(2n)3 2E^(2n)3 2E1(2n)3 2E2(2n)3

x(2n)4S(4 (pa + p, -P1 - P2)S(4 (15)

X{p-pa - Pa) VKa^l-f (Pa) f (Pa) >

where f(pa) is the appropriate occupation probability for bosons or fermions.

4. Relativistic Hydrodynamics

To evaluate the photon and dilepton production from HIC we need to convolute the static rate over space-time integration. Thus, we need to know hydrodynamics which takes care of the evolution of the matter. In this section, we briefly discuss the relativistic hydrodynamics for an ideal as well as viscous medium formed in HIC. Ideally, one cannot describe heavy ion experimental data from the first principle, that is, quantum chromodynamics (QCD) due to its complexity which mainly arises from nonlinearity of interactions of gluons, strong coupling, dynamical many body system, and color confinement. One promising strategy to connect the first

principle with phenomena is to introduce hydrodynamics as a phenomenological theory. Relativistic hydrodynamics [18,3443] plays an important role for an expanding system where pressure, temperature, and so forth vary with space and time. It is assumed that, due to intense rescatterings among the produced secondaries, the system reaches a state of local thermal equilibrium and then the evolution of the system is described by relativistic fluid dynamics. To describe the space-time evolution of such expanding system during the collision, the prescription of relativistic hydrodynamics is essential which assumes the system to be in local thermodynamic equilibrium, which means that pressure and temperature are not constant but rather are the function of space and time. This prescription is valid in the regime where the mean-free path in this "thermalised" system A is much smaller than the characteristic dimensions of the system (L); that is, L> A.

4.1. Basic Equations of Ideal Hydrodynamics. The space-time evolution of the pressure, energy density, particle densities, and the local fluid velocities is controlled by energy momentum conservation equations from hydrodynamics. The basic equations of relativistic hydrodynamics which result from applying constraints of energy-momentum conservations relevant for heavy ion collision at relativistic energies are expressed in

du_ 0,

where T^ is the energy-momentum tensor of fluid element, and in its local rest frame it is given by

Jf = diag [e (x), -P (x), -P (x), -P (*)]. (17)

Local rest frame is the frame in which the velocity of the fluid element is zero. In such a frame the becomes diagonal since the energy flux of the fluid Ti0 and the momentum density T°i turns to be zero. In absence of any dissipative processes the T component becomes the energy density e and T,} = PS,} since T,}dSj is the ith component of force acting on the surface element which according to Pascal's law is isotropic and perpendicular to the surface. P is the pressure of the fluid element in the local rest frame. Isotropy implies that the energy flux Ti0 and the momentum density T0j vanish in the rest frame of fluid. In addition, it implies that the pressure tensor is proportional to the identity matrix, that is, T^ = PS¡j, where P is the thermodynamic pressure.

By doing a proper Lorentz transformation, the energy-momentum tensor T^ in a moving frame, where the fluid moves with an arbitrary four-velocity, u^ = y(1, v) where

Y = (1 - v2)

2\-H2 ■

is given in T^ _ [e (x) + P (x)] uV - P (x) g1

where = diag(1,-1,-1,-1) is the Mankowski metric tensor and u!" is thefluid4-velocityreferredtoas"collectivity" of the system which can be defined as u^ = y(1, V) with y = 1/^1 - v2 and u^u^ = 1, where V is the velocity of fluid element. In the above equation, the e and P are the energy

density and pressure, respectively, in the fluid rest frame, and both are functions of space time coordinate x^.

Apart from the energy-momentum conservation, a fluid may contain several conserved charges, such as total electric charge, and net baryon number. The conserved charges obey the following continuity equation given in (19):

^ _ 0.

N^ = nBu^ is the conserved net baryonic current and nB is baryon number density. For the present work the net baryon number is assumed to be negligible small, so (16) is the only relevant equation to deal with. In addition to it, the total entropy of an inviscid fluid is conserved throughout (S = constant). If we define the entropy current: s^ = su^, then the conservation of entropy results in d s^ = 0 [44].

4.2. Basic Equations ofViscous Hydrodynamics. In the above discussion we considered an idealized situation of a perfect fluid with no internal friction or energy dissipation. But in practice most of the times we have to deal with a system of imperfect fluid in which the density, pressure, and fluid velocity changes over a distance of the order of mean-free path. Such presence of a space-time gradient of those thermodynamic quantities results in modifying the energy momentum tensor T^ and the conserved current N^ to the first-order gradient of these quantities:

_(e + P)u^uv -Pgr + ,

One thing should be mentioned here that for a relativistic fluid it is necessary to specify whether u^ is the velocity of energy transport or velocity of particle transport. In the approach of Landau and Lifshitz, u^ is taken to be the velocity of energy transport and so T'0 vanishes in a comoving frame. In the approach of Eckart, u^ is taken to be the velocity of particle transport and so N* _ nu* _ 0 in a comoving frame. The second approach is adopted here to obtain the following assumptions. The modification in the energy momentum tensor ATand conserved current AN^ is such that in a comoving frame:

u^AT^ _ 0, AN1' _ 0.

With these assumptions we need to construct AT^ to quantify the dissipative processes within the system. This has to be done in such a way that the rate of entropy production per unit volume is positive, which is again required from second law of thermodynamics. To accomplish this task some guidelines are to be followed.

(1) The thermodynamic quantities e, P, and u^ vary slightly over the mean-free path of the particles within the fluid; that is, the system is only very slightly away from equilibrium. So the dissipative term in energy momentum tensor AT^ must be a linear combination of space-time derivatives of e, P, u^, and so forth.

(2) Only the space time derivative of T and u^ can occur in ATV because if derivative of e, P, or n appeared in AT^V, then d^S^ would contain pressure or density gradient, with velocity or temperature gradient and these products are not always positive for all fluid configurations.

The entropy production rate comes out to be

dS* dx*

Lu' -Ldx

AT'0 + L (djUj) AT'1. (22)

From the condition that d s^ > 0 for all fluid configuration we obtain

/ du' = n I -

at11 = nI— + d—-2-^- US'1 ) + (V • US'1, '\dx¡ dx' 3 1

'0 dT du' AT = + C^—

Adx' S dt

AT00 = 0,

where q is the coefficient of shear viscosity, Z is coefficient of bulk viscosity, and \ is thermal conductivity.

Generalizing this expression comes out to be [45]

AT^ = q (vV + VV - ^A^y) + ZA^vy. (24)

Here we have ignored the terms related to thermal conductivity since we are not showing any effect of that on any observables. We have defined = A*vdv, where A*v =

- u*uv is the projection operator.

For the present study, the evaluation of matter from QGP (initial) to the hadronic system (final) via an intermediate quark-hadron transition is studied by applying relativistic hydrodynamics.

4.3. Space-Time Evolution. Hydrodynamics is a general framework to describe the space-time evolution of locally thermalized matter for a given equation of state (EoS). The basic ingredients required to solve the ideal hydrodynamic equations are EoS and initial conditions. As the system expands from its initial state, the mean-free path between particles within the system increases. At certain stage, the mean-free path becomes comparable to the system size, and then the hydrodynamic description breaks down and the phase space distribution of the particle gets fixed by the temperature of the system at this stage. This stage of evolution is called freezeout state and the corresponding temperature of the system is called thermal freezeout temperature (Tj). The hydrodynamic evolution stops at the freezeout point.

4.3.1. Initial Condition. The initial conditions are crucial to the description of space-time evolution. Initial conditions in hydrodynamics may be constrained in the following ways to reproduce the measured final multiplicity. We assume that the system reaches equilibration at a time (called initial thermalization time) after the collision. The T¡ can be

related to the measured hadronic multiplicity (dN/dy) by the following relation [46]:

1 1 dN

45C (3) 4aeff nR2. dy

where RA - 1.1Np/3t is the radius of the system, £(3) is the

Riemann zeta function, and aeff = n2geff/90, geff(= 2x8 + (7/8)x2x2xNc xNp) is the degeneracy of quarks and gluons in QGP, Nc = number of colors, Np = number of flavors. The factor "7/8" originates from the difference between the Bose-Einstein and the Fermi-Dirac statistics. Tt depends on the centrality through the multiplicity, dN/dy. The value of dN/dy for various beam energies and centralities can be obtained directly from experiment or calculated using the following relation [47]:

dN (Npart)

-dy=(1-x)n??^r-

(Ncoll) ,

where npp is the multiplicity per unit rapidity measured in pp collisions: npp = 2.5 - 0.25 ln(s) + 0.023ln2(s), the fraction x of npp is due to "hard" processes, with the remaining fraction (1 - x) being "soft" processes. The multiplicity in nuclear collision has then two components: "soft," which is proportional to number of participants, Npart and "hard," which is proportional to number of binary collision, NcoU.

After the initial thermalization time, Tj, the system can be treated hydrodynamically. The initial conditions to solve the hydrodynamic equations are given through the energy density and velocity profile:

e(T',r) =

1 + exp ((r-RA) IS)' v (T',r) = 0,

where e0 is the initial energy density which is related to initial ( T;), Ra is the nuclear radius, and S is the diffusion parameter taken as 0.5 fm.

4.3.2. Equation of State (EoS). The set of hydrodynamic equations are not closed by itself; the number of unknown variable exceeds the number of equations by one. Thus a functional relation between any two variables is required so that the system become deterministic. The most natural course is to look for such relation between the pressure P and the energy density e. Under the assumption of local thermal equilibrium, this functional relation between P, e, and nB is the EoS:

P = P(e,nB)

which expresses the pressure as function of energy density, e, and baryon density, nB. This can be obtained by exploiting numerical lattice QCD simulation [48].

Different EoSs (corresponding to QGP vis-a-vis that of hadronic matter) will govern the hydrodynamic flow quite differently. It is thus imperative to understand in what

(d) (e)

Figure 4: Partonic processes for production of photons.

respects the two EoSs differ and how they affect the evolution in space and time. The role of the EoS in governing the hydrodynamic flow lies in the fact that the velocity of sound, c2 = (dP/de) sets an intrinsic scale in hydrodynamic evolution. One can thus write simple parametric form of the EoS: P = c2se, for baryon-free system which is relevant for the present study.

4.3.3. Freeze-Out Criteria. The expansion persists as long as the fluid particles interact. At sufficiently longer A when it is comparable to system size the particles decouple to behave as free particles which is called "freeze-out" stage. This freeze-out scenario is characterized by a system temperature Tj which is of the order of pion mass and defines a space-time surface T(x, t) = Tf which serves as the boundary of the hydrodynamical flow [49].

5. Emission of Thermal Photons from Heavy Ion Collision

The thermal photons emerge just after the system thermalizes (t > t¡) from both QGP due to partonic interactions and hot hadrons due to interactions among the hadrons. Now with the formalism discussed in Section 3, the production of thermal photons from QGP and hot hadronic gas is given in Sections 5.1 and 5.2, respectively. And using the hydrodynamic equations, we have convoluted these static rates by space-time integration (discussed in Section 4) and obtained the total invariant yield of photon for different collision energies. The space-time integration is constrained to the hydrodynamical inputs which has been discussed elaborately in this section.

5.1. Photons Emission from Quark Gluon Plasma. The contribution from QGP to the spectrum of thermal photons due to annihilation (qq ^ gy) and Compton (q(q)g ^ q(q)y) processes has been calculated in [17, 50] using hard thermal loop (HTL) approximation [51, 52]. The rate of hard photon emission is then obtained as [17]

^ T2e-E/T ln

f 2n2 '

2.912E)

where as is the strong coupling constant. Later, it was shown that photons from the processes [53]: gq ^ gqy, qq ^ qqy, qqq ^ qy, and gqq ^ gy contribute in the same order 0(aas) as Compton and annihilation processes (shown in Figure 4). The complete calculation of emission rate from QGP to order as has been performed by resuming ladder diagrams in the effective theory [54, 55]. In the present work this rate has been used. The temperature dependence of the strong coupling, as, has been taken from [56].

5.2. Photons Emission from Hot Hadronic Gas. For the photon spectra from hadronic phase we consider an exhaustive set of hadronic reactions and the radiative decay of higher resonance states [33, 57, 58].

To evaluate the photon emission rate from a hadronic gas we model the system as consisting of n, p, to, and The relevant vertices for the reactions nn ^ py and np ^ ny and the decay p ^ nny are obtained from the following Lagrangian [57] (see Figure 5):

L = -g^ • (n x d^n) - eJ"A„ + fa x pv)3,

л+ У л+ p°

------VWWWWVí ------j-

0 - ! л;

Л ' P Л 1 У ------1--------Кллллллллл,

+ л+ У

Sp0 Л+^ — 0 ^ ~ ~

r p Л

5.3. Total Invariant Momentum Spectra of Thermal Photons. In this section we evaluate photon spectrum from a dynamically evolving system. The evolution of the system is governed by relativistic hydrodynamic. The photon production from an expanding system can be calculated by convoluting the static thermal emission rate with the expansion dynamics, which can be expressed as follows:

^Л- p0

dNy d2pTdy

d2pTdy

where the d4x is the four volume. The energy, E*, appearing

Figure 5: Photon-producing reactions and decays in hadronic gas. in (33) should be replaced by u^p^ f°r a system expanding with space-time-dependent four-velocity u^. Under the assumption of cylindrical symmetry and longitudinal boost

where = d^Av - dvA is the Maxwell field tensor and J^ is the hadronic part of the electromagnetic current given by

j" = (Pr

(n x (&

with B^ = d^pv - dvpp - gpnn(p^ X pv).

For the sake of completeness we have also considered the photon production due to the reactions nq ^ ny, nn ^ qy, and the decay to ^ ny using the following interaction:

Зшрп v-JX в ~""P л "

+ ш д УЛ —А?Р ■

0рлл

The last term in the above Lagrangian is written down on the basis ofvector meson dominance (VMD) [60,61]. To evaluate the photon spectra, we have taken the relevant amplitudes for the abovementioned interactions from [33, 57]. The effects of hadronic form factors [62] have also been incorporated in the present calculation. The reactions involving strange mesons: nK* ^ Ky, nK ^ K*y, pK ^ Ky, and KK* ^ ny [62, 63] have also been incorporated in the present work. Contributions from other decays, such as K* (892) ^ Ky, $ ^ ny,bl(1235) ^ ny,a2(1320) ^ ny,andK1(1270) ^ ny, have been found to be small [63] for pT >1 GeV.

With all photon-producing hadronic reaction, the static thermal emission rate ofphotons for hadronic phase has been evaluated [17, 33, 54, 55, 57, 62]. The reaction involving p mesons has dominant contribution. The rate at low photon energy is dominated by reaction with p in final state, because these reactions are endothermic with most of the available energy going into rho mass. At high photon energy reactions with the p in initial state are dominant because these reactions are exothermic; most of the rho mass is available for the production of high energy photons. Similar remarks can be made concerning reactions involving q mesons, but as the value of gpm is smaller, thus so are the rates. All the isospin combinations for the above processes have properly been implemented.

invariance, и" can be written as

и = yT (т, r) vr (т,г),

= yT (MT cosh ц, их, иу, MT sinh ц) = yT (MT cosh ц, vr cos ф, vr sin ф, MT sinh ц),

where vr(T, r) is the radial velocity, yr(T, r) = (1 - vr(T, r)) ' and, therefore, for the present calculations,

U"pu = Yr (MT cosh (У-П)- vrpT cos ф) ■

For massless photon the factor и^р^ can be obtained by replacing MT in (35) by pT. For the system produced in QGP phase reverts to hot hadronic gas at a temperature T ~ Tc. Thermal equilibrium may be maintained in the hadronic phase until the mean-free path remains comparable to the system size. The term "(dR/d2 pTdy)¡ = [(•••)fbe\" is the static rate of photon production where i stands for quark matter (QM), mixed phase (M) (in a lst-order phase transition scenario), and hadronic matter (HM), respectively. The pT dependence of the photon and dilepton spectra originating from an expanding system is predominantly determined by the thermal factor fBE. The total momentum distribution can be obtained by summing the contribution from QM and HM, where the distribution for both the phases can be obtained by choosing the phase space appropriately.

The d4x integration has been performed by using relativistic hydrodynamics with longitudinal boost invariance [41] and cylindrical symmetry [64] along with the inputs (given in Table 1) as the initial conditions for SPS and RHIC energies.

To estimate dN/dy for RHIC, we have taken dnpp/dy =

2.43 and x = 0.1 at

' = 200 GeV. It should be mentioned

here that the values of dN/dy (through N"part and Ncon in (26)) and hence the Tt (through dN/dy in (25)) depend on the centrality of the collisions. For SPS, dN/dy is taken from experimental data [65]. We use the EoS obtained from the lattice QCD calculations by the MILC collaboration [66]. We consider kinetic freeze-out temperature, Tj = 140 MeV for all the hadrons. The ratios of various hadrons measured experimentally at different -^sNN indicate that the system formed in heavy ion collisions chemically decouple at Tch

Table 1: The values of various parameters—thermalization time {ri), initial temperature (T), and hadronic multiplicity dN/dy (the value of dN/dy for various beam energies and centralities is calculated from (26))—used in the present calculations.

4s nn Centrality dN/dy Г (fm) T (MeV)

17.3 GeV 0-06% 700 1.0 200

0-20% 496 0.6 227

200 GeV 20-40% 226 0.6 203

min. bias 184 0.6 200

2.76 TeV 0-40% 1212 0.1 553

which is higher than Tj which can be determined by the transverse spectra of hadrons [67, 68]. Therefore, the system remains out of chemical equilibrium from Tch to Tj. The deviation of the system from the chemical equilibrium is taken into account by introducing chemical potential for each hadronic species. The chemical nonequilibration affects the yields through the phase space factors of the hadrons which in turn affects the productions of the EM probes. The value of the chemical potential has been taken into account following [69].

5.4. Results and Discussion on pT Distributions of Photons. For comparison with direct photon spectra as extracted from HIC two further ingredients are required. With all the ingredients we have reproduced the pT spectra of direct photon for both SPS and RHIC energies. The prompt photons are normally estimated by using perturbative QCD. However, to minimize the theoretical model dependence here, we use the available experimental data from p-p collisions to estimate the hard photon and normalized it to A-A data with TAA(b) for different centrality; that is, the photon production from A-A collision and p-p collision are related to the following relation:

Ncoii (b) da

d2pTdy

d2pTdy

= Taa (b)

d2pTdy'

where NcoU(b) is taken for the corresponding experiments and the the typical app (app 41 mb for RHIC and 30 mb for SPS).

5.4.1. Photon Spectrum for WA98 Collaboration. The WA98 photon spectra from Pb+Pb collisions are measured at Vnn = 17.3 GeV. However, no data at this collision energy is available for pp interactions. Therefore, prompt photons for p+p collision at Vsnn = 19.4 GeV have been used [70] to estimate the hard contributions for nuclear collisions at Vnn = 17.3 GeV. Appropriate scaling [65] has been used to obtain the results at Vsnn = 17.3 GeV. For the Pb+Pb collisions the result has been appropriately scaled by the number of collisions at this energy (this is shown in Figure 6 as prompt photons). The high pT part of the WA98 data is reproduced by the prompt contributions reasonably well. At low pT the hard contributions underestimate the data indicating the presence of a thermal source. The thermal photons with initial temperature = 200 MeV along with the

— 1 n-2 ^ 10

10-4 г

10 5 r

10-6 r

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 pT (GeV)

• Pb+Pb (WA98)

- Thermal + pQCD

- Prompt photons

Figure 6: Transverse momentum spectra of photon at SPS energy for Pb+Pb collision at midrapidity [12].

prompt contributions explain the WA98 data well (Figure 6), with the inclusion of nonzero chemical potentials for all hadronic species considered [25, 26, 69, 71, 72]. In some of the previous works [73-78] the effect of chemical freezeout is ignored. As a result either a higher value of or a substantial reduction of hadronic masses in the medium was required [73]. In the present work, the data has been reproduced without any such effects.

5.4.2. Photon Spectrum for PHENIX Collaboration. In Figure 7, transverse momentum spectra of photons at RHIC energy for Au-Au collision for three different centralities (0-20%, 20-40%, and min. bias.) at midrapidityi shown, where the red tangles are the direct photon data measured by PHENIX collaboration [79] from Au-Au collision at ■^sNN = 200 GeV, blue-dashed line is the contribution of the prompt photons and the black solid line is thermal + prompt photons. For the prompt photon contribution at ■^sNN = 200 GeV, we have used the available experimental data from pp collision and normalized it to Au-Au data with TAA(b) for different centrality [80] (using (36)). At low pT the prompt photons underestimate the data indicating the presence of a possible thermal source. The thermal photons along with the prompt contributions explain the data [79] from Au-Au collisions at -^sNN = 200 GeV reasonably well. The reproduction of data is satisfactory (Figure 7) for all the centralities with the initial temperature shown in Table 1 [81].

5.4.3. Photon Spectrum for ALICE Collaboration. The direct photon spectra from Pb+Pb collisions are measured at ^sNN = 2.76 TeV for 0-40% centrality by ALICE collaboration. However, no data at this collision energy is available

¡1 1 1 1 1 1 1 1 1 1 1 x 100 =

''^„, 20-40% x 10-2 j

_ ~„_Mm. bias x 10-4

! RHIC ........ .................."

pr (GeV)

A PHENIX data

...... Prompt photon

- Prompt + thermal

Figure 7: Transverse momentum spectra ofphotons at RHIC energy for Au-Au collision for different centralities at midrapidity [12].

- Thermal + prompt photon

---Prompt photon from pp data

...... Thermal photon

Figure 8: Transverse momentum spectra of photons at LHC energy for Pb-Pb collision for different 0-40% centrality at midrapidity [59].

for pp interactions. Therefore, prompt photons from p+p collision at —sNN = 7 TeVhavebeenusedtoestimatethe hard contributions for nuclear collisions at -sNN = 2.76 TeV by using the scaling (with -sNN) procedure used in [65]. For the Pb+Pb collisions the result has been scaled up by the number of collisions at this energy (this is shown in Figure 8 as prompt photons). The high pT part of the data is reproduced by the prompt contributions reasonably well. At low pT the hard contributions underestimate the data indicating the presence of a possible thermal source.

The thermal photons with initial temperature -553 MeV along with the prompt contributions explain the data well (Figure 8), with the inclusion of nonzero chemical potentials for all hadronic species considered [69] (see also [71, 72]).

It is well known that transverse momentum spectra of photons act as a thermometer of the interior of the plasma. The inverse slope of the thermal distribution is a measure of the average (over evolution) effective (containing flow) temperature of the system. We have extracted the average effective temperature (- (pT)) from the thermal distributions of photons at different collision energies—that is, for SPS, RHIC, and LHC energies. Figure 9 shows the variation of (pT) with multiplicity for different collision energies. To minimize the centrality dependence of the results the dNch/dq is normalized by Npart. The results clearly indicate a significant rise in the average pT ((pT)) while going from SPS to RHIC to LHC. The values of (pT) for different collision energies are given in Table 1. Since photons are emitted from each space time point of the system, therefore, the measured slope of the pT spectra represents the average effective temperature of the system.

The quantity, pff(= 1/NpartdNch/d^), is proportional to the entropy density. Therefore, /(pr)3 a the average

ALICE @ 2.76 TeV

sr 3 "в 3

PHENIX @ 200 Ge"^

WA98 @ 17.3 GeV

» WA98 I ALICE

■•V*..... "

pr (GeV) _l_ï_ï_ï_L_

0.22 0.24 0.26 0.28 0.3 (pr) (GeV)

Figure 9: The variation of (pT) with the increase in multiplicity for different collision energies.

effective statistical degeneracy, a quantity which changes drastically if the colour degrees of freedoms deconfined; that is, if a phase transition takes place in the system. We find that the entropy density (s ~ geffT ) at LHC increases by almost 96% compared to RHIC, and there is an enhancement of 46% at RHIC compared to SPS. However, part of this increase is due to the increase in the temperature and part is due to increase in degeneracy. To estimate the increase in the degeneracy we normalize the quantity pff by (pT) . Therefore, we estimate pffKpr)3 from the analysis of the experimental data and found that there is a 15% increase in this quantity from SPS to RHIC and 35% increase from RHIC to LHC.

5.5. Total Invariant Momentum Spectra of Thermal Photons in Viscous Medium. Effects of viscosity on the transverse momentum distribution of photons were earlier considered in [82, 83] and recently the interest in this field is renewed [84-86]. The measured photon spectra (dN/d2 pTdy) are the yield obtained after performing the space time integration over the entire evolution history—from the initial state to the freezeout point using (33). Beyond a certain threshold in collision energy the system is expected to be formed in QGP phase which will inevitably make a transition to the hadronic matter later. The measured spectra contain contributions from both QGP and hadronic phases. Therefore, it becomes imperative to estimate the photon emission with viscous effects from QGP as well as hadrons and identify a kinematic window where photons from QGP dominate. While in some of the earlier works [84-86] contributions from hadrons were ignored, in others [82, 83] the effects of dissipation on the phase space factors were omitted. In the present work we study the effects of viscosity on the thermal photon spectra originating from QGP and hadronic matter and argue that photons can be used as a very useful tool to estimate j/s and hence characterize the matter.

Equation (12) can be simplified to the following form (see the appendix) [87]:

d2pTdy 16(2n)i

J plTdplTdp2Td^ldyldy2

"fl (Pl)f2 (p2)(l±fs (p3))*\M\2

X \PiT sin {$1 -$2) + pT sin 02|-21=0O •

The effects of viscosity on the photon spectra resulting from HIC enter through two main factors: (i) the modification of the phase space factor due to the deviation of the system from equilibrium and (ii) the space time evolution of the matter governed by dissipative hydrodynamics. One more important issue deserves to be mentioned here. Normally, the initial temperature (T;) and the thermalization time (r;) are constrained by the measured hadron multiplicity (dN/dy). This approach is valid for a system where there is no viscous loss and the time reversal symmetry is valid. However, for a viscous system the entropy at the freezeout point (which is proportional to the multiplicity) contains the initially produced entropy as well as the entropy produced during the space time evolution due to nonzero shear and bulkviscosity. Therefore, the amount of entropy generated during the evolution has to be subtracted from the total entropy at the freezeout point, and the remaining part which is produced initially should be used to estimate the initial temperature. Therefore, for a given dN/dy (which is associated with the freezeout point) and ri the magnitude of Tt will be lower in case of viscous dynamics compared to ideal flow.

5.5.1. Viscous Correction to the Distribution Function. We assume that the system is slightly away from equilibrium which relaxes back to equilibrium through dissipative processes. Here we briefly recall the main considerations leading

to the commonly used form for the first viscous correction, Sf, to the phase space factor, f, defined as follows [88]:

f, (p)=/¡a (i+m

= f,o (l + pp [c(Vauß)+AAaß W-u]

where fi0 is the equilibrium distribution function of "ith" particle, (yaup) = VaUp + Vf}ua - (2/3)AapVyuy, Aa<3 =

Bap - uaup, = (gap - uaup)dp, up being the four-velocity of the fluid. The coefficients C and A can be determined in the following way. Substituting f in the expression for stress-energy tensor Tpv we get

T"=\éJhp'p"fo (1+Sf

(2n) = T**v + AT^,

where T^v = (e + P)u^uv - g^P is the energy momentum tensor for ideal fluid. From general considerations [44] the dissipative part can be written as

AT= n (W) + ÇA*™W - u.

Equating the part containing Sf from (38) with (40), C and A can be expressed in terms of the coefficients of shear (j) and bulk (0 viscosity, respectively, in terms of which the phase space distribution for the system can be written as

f==fo (l + ^papß (VaUß)-5TT3papß A aß W-u).

For a boost invariant expansion in (1 + 1) dimension this can be simplified to get

f = fo [l+Sf-Sf]

*fn ==3Ïi-T(p2r-2p'z2),

*=5Tk&+p*

where p'z = mT sinh(y-q) is the z-component of the momentum in the fluid comoving frame. The phase space distribution with viscous correction (42) thus enters the production rate of photon through (37).

5.5.2. Viscous Correction to the Expansion Dynamics. As mentioned before the pT distribution of thermal photons is obtained by integrating the emission rate over the evolution history of the expanding fluid. Relativistic viscous hydrodynamics can be used as a tool for the space-time dynamics of the fluid.

For a (1 + 1) dimensional boost invariant expansion [41] the evolution equation, d^T^ = 0, can be written as [89]

de e + P _ (4/3) q+ C dr r r2

where P is the pressure and e is the energy density. We assume that the baryonic chemical potential is small in the central rapidity region for RHIC/LHC collision energies. Therefore, the equation corresponding to the net baryon number conservation need not be considered in these situations.

We assume that the system achieves thermal equilibrium at a time ri after the collision at an initial temperature Ti. With this initial condition and equation of state (EoS) P = e/3 the solution of (44) can be written as [82]

T==nT. )U3 + AQ

where Aq = ((4/3)^ + Cqo), 1qo = 1q/T3 = 4aQ(q/s)Q, and

Cqo = Cq/T3 = 4aQ(C/s)Q.

Equation (45) dictates the cooling of the QGP phase from its initial state to the transition temperature, Tc, at a time, rq, when the QGP phase ends.

In a first-order phase transition scenario, the pure QGP phase is followed by a coexistence phase of QGP and hadrons. The energy density, shear, and bulk viscosities in the mixed phase can be written in terms of the corresponding quantities of the quark and hadronic phases at temperature Tc as follows [82]:

eM (r) = fQeQ (Tc) + (l-fQ (r)) eH (Tc) >

1m (r) = fQVQ (Tc) + (l-fQ (r)) Vh (Tc) , (46)

Cm (r) = fQCQ (Tc) + (l-fQ (t))Ch (Tc),

where fQ(r)(fH(r)) indicates the fraction of the quark (hadronic) matter in the mixed phase at a proper time r. We have eQ(Tc) = 3aQTc4 + B, e^T) = 3asTc4, aQ = gQn2/ 90, aH = gHn2/90, B is the bag constant, gQ (gH) denote statistical degeneracy for the QGP (hadronic) phase. In the mixed phase the temperature remains constant but the energy density varies with time as the conversion of QGP to hadrons continues. This time variation is executed through fQ(r). Substituting (46) in (44) and solving for fQ(r) we get [82]

„-b/r fT'

dr' + rQe(b/rQ-b/r)> r

where a = 4eH/(3Ae), b = [4(iq - 1h)/3 + 2(Cq - Ch)]/&, c = ((4/3)qH + 2CH)/Ae, and Ae = eQ - eH. Equation (47) indicates how the fraction of QGP in the coexistence phase evolves with time.

The variation of T with r in the hadronic phase can be obtained by solving (44) with the boundary condition T = Tc and r = rH, where rH is the (proper) time at which the mixed

Table 2: The values of various q/s and parameters—initial temperature (Tj), starting time of mixed phase (rq), and ending time of mixed phase (rh)—used in the present calculations.

q/s Tt (MeV) rq (fm) rh (fm)

0 328 3.95 8.8

1/4n 315 3.93 8.78

2/4n 302 3.92 8.76

phase ends; that is, when the conversion of QGP to hadronic matter is completed,

T = Tc(TH

Similar to QGP, P = e/3 has been used for hadronic phase. For a vanishing bulk viscosity (C = 0) the cooling of the QGP is dictated by

1/3 2 + —

r ) 3ri\s Jq

Similarly the time variation of temperature in the hadronic phase is given by

T = Tc

In a realistic scenario the value of q/s may be different for QGP [90-94] and hadronic phases [95-98]. However, in the present work we take the same value of q/s both for QGP and hadronic matter as shown in Table 2.

5.5.3. Results and Discussion on Viscous Effect on pT Distributions of Photons. In case of an ideal fluid, the conservation of entropy implies that the rapidity density dN/dy is a constant of motion for the isoentropic expansion [41]. In such circumstances, the experimentally observed (final) multiplicity, dN/dy, may be related to a combination of the initial temperature Ti and the initial time ri as Tj3ri. Assuming an appropriate value of I (taken to be ~0.6 fm/c in the present case), one can estimate Ti.

For dissipative systems, such an estimate is obviously inapplicable. Generation of entropy during the evolution invalidates the role of dN/dy as a constant of motion. Moreover, the irreversibility arising out of dissipative effects implies that estimation of the initial temperature from the final rapidity density is no longer a trivial task. We can, nevertheless, relate the experimental dN/dy to the freezeout temperature, Tf, and the freezeout time, if, by the relation:

nR2A4asT3frf

where RA is the radius of the colliding nuclei (we consider AA collision for simplicity) and k is a constant ~3.6 for massless bosons.

To estimate the initial temperature for the dissipative fluid we follow the following algorithm. We treat Ti as a parameter;

for each Tt, we let the system evolve forward in time under the condition of dissipative fluid dynamics (44) till a given freezeout temperature Tj is reached. Thus rj is determined. We then compute dN/dy at this instant of time from (51) and compare it with the experimental dN/dy. The value of Tt for which the calculated dN/dy matches the experimental number is taken to be the value of the initial temperature. Once Tt is determined, the evolution of the system from the initial to the freezeout stage is determined by (45), (47), and (48).

In Figure 10 we display the variation of temperature with proper time. It is clear from the results shown in the inset (Figure 10) that initial temperature for system which evolves with nonzero viscous effects is lower compared to the ideal case for a fixed dN/dy. Because of a nonviscous isentropic evolution scenario the multiplicity (measured at the freezeout point) is fixed by the initial entropy. However, for a viscous evolution scenario the generation of entropy due to dissipative effects contributes to the multiplicity. Therefore, for a given multiplicity (which is proportional to the entropy) at the freezeout point one requires lower initial entropy; hence, initial temperature will be lower. It is also seen (Figure 10) that the cooling of the system is slower for viscous dynamics because of the extra heat generated during the evolution.

In this section we present the shift in the pT distribution of the photons due to viscous effects. The integrand in (33) is a Lorentz scalar; consequently the Lorenz transformation of the integrand from the laboratory to the comoving frame of the fluid can be effected by just transforming the argument; that is, the energy of the photon (E = pT cosh(y)) in the laboratory frame should be replaced by u^p^ in the comoving frame of the fluid, where p^ is the four momenta of the photon.

The results presented here are obtained with vanishing bulk viscosity. The effects of viscosity enter into the photon spectra through the phase space factor as well as through the space time evolution. We would like to examine these two effects separately. For convenience we define two scenarios:

(i) the effects of viscosity on the phase space factor are included (Sf, = 0) in (42), but the viscous effects on the evolution are neglected (j = 0 )in (44),

(ii) the effects of j = 0 are taken into account in the phase space factors as well as in the evolution dynamics.

The space time-integrated photon yield originating from the QGP in scenario (i) is displayed in Figure 11. Note that the value of the initial temperatures for the results displayed in Figure 11 is the same (for all j/s) because the viscous effects on the evolution are ignored in scenario (i). The viscous effects on the pT distribution of the photons are distinctly visible. The higher values of j/s make the spectra flatter through the pT dependence of the correction, Sf ,'.

Next we assess the effects of viscosity on photon spectra for scenario (ii). In Figure 12 we depict the photon spectra for various values of j/s. In this scenario the value of Tt is lower for higher j/s for reasons described above. As a result the enhancement in the photon production due to change

0.34....................... :

0.32 ^ \ : ii

0.28 ^ \

O 0.24 ^ \

0.12 r r

0 2 4 6 8 10 12 14 t (fm)

- fj/s = 0

---tj/s = l/4n

...... tj/s = 2/4n

Figure 10: Variation of temperature with proper time for different phases for various values of the shear viscosities. Inset shows the effect of viscosity on the cooling of the QGP phase (in an amplified scale) for different values of y/s.

- q/s = 0

--- q/s = l/4n

...... q/s = 2/4n

Figure 11: Transverse momentum distribution of thermal photons from QGP for various values of j/s in the scenario (i).

in phase space factor, Sf^, is partially compensated by the reduction in Tt for nonzero y, which is clearly seen in the results displayed in Figures 11 and 12.

In Figures 13 and 14 we exhibit results for the hadronic phase for scenarios (i) and (ii), respectively. The effects of dissipation on the pT distribution of photons from hadronic phase are qualitatively similar to the QGP phase; that is,

pT (GeV)

- ^/s = 0

---^/s = 1/4re

...... ^/s = 2/4«

Figure 12: Transverse momentum distribution of thermal photons from QGP for various values of q/s in the scenario (ii).

- ^/s = 0

---^/s = 1/4rc

...... ^/s = 2/4«

Figure 13: Transverse momentum distribution of photons from thermal hadrons for various values of q/s in the scenario (i).

the effect is more prominent in scenario (i) than in (ii). It is also clearly seen that the effects of viscosity though the effect is stronger in the QGP phase than in the hadronic phase. It is expected that the observed shift in the photon spectra due to viscous effects may be detected in future high precision experiments.

Finally in Figures 15 and 16 we plot the pT spectra of photons for the entire life time of the thermal system; that is, the photon yield is obtained by summing up contributions from QGP, mixed and hadronic phases for different values of

- ^/s = 0

---^/s = 1/4«

...... ^/s = 2/4«

Figure 14: Transverse momentum distribution of photons from thermal hadrons for various values of q/s in the scenario (ii).

- ^/s = 0

---^/s = 1/4«

...... ^/s = 2/4«

Figure 15: Transverse momentum distribution of thermal photons from the entire evolution history of the system for various values of q/s in the scenario (i).

q/s for scenario (i) and (ii), respectively. The effect ofviscosity for the scenario (i) is stronger than (ii).

6. Emission of Thermal Dileptons from Heavy Ion Collision

Unlike real photon, dilepton is massive. Thus dilepton has two kinematic variables, invariant mass (M) and transverse momentum (pT). Again, the pT spectra are affected due to

10 10 10 10 10 ■а ю-6 10 10

convenient to break up the quark current Jp into parts with definite isospin:

2.5 3 3.5 pT (GeV)

- ^/s = 0

---^/s = 1/4«

...... ^/s = 2/4«

Figure 16: Transverse momentum distribution of thermal photons from the entire evolution history of the system for various values of q/s in the scenario (ii).

flow, whereas the _pT-integrated M spectra remain unaltered by flow. By tuning these two parameters, different stages of expanding fireball can be understood. Dileptons having large M and high pT are emitted early from the hot zone of the system. On the other hand, those having lower M and pT produced at later stage of the fireball when the temperature is low. Because of an additional variable, the invariant pair mass M, dileptons have the advantage over real photons [99].

The production of thermal dileptons from QGP (Section 6.1) and hot hadronic gas (Section 6.2) is described below.

6.1. Dileptons Emission from QGP. In the QGP, where quarks and gluons are the relevant degrees of freedom, the Wpv can be directly evaluated by writing the electromagnetic current in terms of quarks of flavor f, that is, Jpm = ^fCfffY^yf. Confining to the leading order contribution we obtain

The rate in this case corresponds to dilepton production due to process qq ^ y* ^ l+l~. The static thermal emission rate of dilepton from QM is given by (qq ^ y* ^ l+T)

[100,101] (also [102,103]),

ff = ef

1 + 2L in ("+

%=2 (иУии - dbd)+6 +dbd)+•••

where V and S denote iso-vector and iso-scalar currents and the dots denote currents comprising of quarks with strangeness and heavier flavors. These currents couple to individual hadrons as well as multiparticle states with the same quantum numbers and are usually labeled by the lightest meson in the corresponding channel [104]. We thus identify the isovector and isoscalar currents with the p and w mesons, respectively. Defining the correlator of these currents W^'^ analogously as in (7), we can write

¡ЛУ

= Wp +WL

The correlator of vector-isovector currents W?v has in fact been measured [105, 106] in vacuum along with the axial-vector correlator by studying r decays into even and odd number of pions. The former is found to be dominated at lowerenergiesbythe prominentpeakofthe p meson followed by a continuum at high energies. The axial correlator, on the other hand, is characterized by the broad hump of the a1 . The distinctly different shape in the two spectral densities is an experimental signature of the fact that chiral symmetry of QCD is dynamically broken by the ground state [107]. It is expected that this symmetry may be restored at high temperature and/or density and will be signaled by a complete overlap of the vector and axial-vector correlators [17].

In the medium, both the pole and the continuum structure of the correlation function gets modified [8,108]. We will first evaluate the modification of the pole part due to the self-energy of vector mesons in the following. Using vector meson dominance the isovector and scalar currents are written in terms of dynamical field operators for the mesons allowing us to express the correlation function in terms of the exact (full) propagators or the interacting spectral functions of the vector mesons in the medium. To reach that goal we have to specify the coupling of the currents to the corresponding vector fields. For this purpose we write, in the narrow width approximation [104],

(0\j;m (0)\R)=FRmRei,

where e f is the charge of the quark and n± = l/(e(P°±lPl/2T)+1). e+e through the relation

where R denotes the resonance in a particular channel and ep is the corresponding polarization vector. The coupling constants Fr are obtained from the partial decay widths into

6.2. Dileptons Emission from Hot Hadronic Gas. To obtain the rate of dilepton production from hadronic interactions it is

г ' e

yielding FR = 0.156 GeV, 0.046 GeV, and 0.079 GeV for p, to, and <p, respectively. Equation (56) suggests the operator relations:

Jp. (x) = FpmpVP (x)

J" (x) = 3FamaV^ (x), and so forth,

where V."^ ( x) denotes the field operator for the p(u>) meson. So using the above relations connecting currents to fields (so-called field-current identity), the current commutator becomes

W = I FRmR f d4xeiqx ([vRR (x),vR])

D— J

R=Prr>...

= 1 FRmRARv (%>$)

= 2e i%) Z FRmR Im {clo> q) >

where A^v are the spectral functions of corresponding vector _g

meson resonances (R) and D is the diagonal element of the thermal propagator matrix. The form of the diagonal element of the exact thermal propagator matrix for the spin 1 particle is given by

uvy'1J q2mR '

—R , , —R

Duv (q) = -

Qr /q2

q2 -mR - nR (q) q2 - mR - q2n" (q)

ittR ,

The imaginary part is then put in (59) and then in (6) to arrive at the dilepton emission rate:

d4qd4x n3q2

L(q2)fBE (qo

x [Fpm2pAp (qo, q) + FlmlAa (qo, q} + -"]

where, for example, Ap(= -guv Im D /3) is given by

A0 = —

2^ Im nt

q2 -m2. -I Re < ) + (l Im ^

i —R

1 Im n,_

(q2-m.-q2l Re nR) + q4(l Im n

the sum running overall meson loops nH and baryon loops _g

NB. Here ntl is the diagonal element of vector meson (R = p and to) self-energy at finite temperature and density which

V(q) f

\ v(q)

H(q - k)

v(q) v(q)

B(k + q)

B(k - q)

Figure 17: One-loop Feynman diagrams for p or m self-energy involving mesons (first figure) and baryons (second and third figures). V stands for the p or m in the external line. In the internal lines of meson loops, H = n, m, a1, and h1 for V = p whereas H = p for V = m. For the baryonic loops, N and B indicate, respectively, nucleon and baryonic internal lines.

is also a matrix in the real-time formalism. We have taken H = n, m, a1,h1 [109] and B = N(940), N'(1520), N'(1650), N'(1700), N'(1720), A(1230), and A'(1620) [110] for p meson whereas for m meson, H = p (with nn folding), and B = N(940), N'(1440), N'(1520), N'(1535), N'(1650), and N'(1720) [111] are taken. These self-energy graphs are

diagrammatically represented in Figure 17.

The general expression of ntl for meson loop nH (representing the first diagram of Figure 17) is given by [109,111]

(2n)3 4w„wH

(l + r?)Lti + nHL3tl qo + im (qo)

+ -n7Llti + nHL4ti qo + ^h + im (qo) n7IL2t,l -nHL%

% +Mn -mh + im

+ -n"L2ti -(1 + nH)L{l

% + Mn + mh + im {<io)

where n's are Bose-Einstein distribution functions for the internal meson lines and m's are their on-shell energies. In the above expression L't=j1'""4 denote the values of Lt l (k0) for k0 = Mn, -mu, q0 - mh, and q0 + mh, respectively. The corresponding expression for the baryon loop NB (second diagram of Figure 17) is given by [110,111]

(l-n^L), -nBL3t, qo + im (qo)

nNjl -nBj4 n+ L t,l n+L t,l

qo + ^b + im (qo) + -nNL2ti + nBL3ti qo + ^n + im (qo)

nNL2t, + (-1+nB)L

qo + toN + toB + iye (qo

e e(M -60

«X JS 20

0.4 0.5 0.6 M (GeV)

0.1 0.2 0.3

0.4 0.5 0.6 M (GeV)

0.7 0.8 0.9

B(^ = 0 MeV) B (^ = 150 MeV)

----- Total (^ = 0 MeV)

- Total (^ = 150 MeV)

-20 -15 -10 -5 0

10 5 0 -5 -10

1 ' 1» ' 1 A /» -

.... /i

_____ .......

—" /'•.

1 — r i^' . 1 ______

0.1 0.2 0.3 0.4

0.5 0.6 M (GeV)

0.7 0.8

----- M

■-■- B(№ = 0) ---Total (^B = 0)

0.5 0.6 M (GeV)

...... B (^B = 150 MeV)

- Total (^B = 150 MeV)

Figure 18: The imaginary (upper) and real (lower) part of self-energy function of p (a) and w (b) at different chemical potential (p or pB). The contribution coming from meson, baryon loops, and their sum is shown.

where n+(-_)'s are Fermi-Dirac distribution functions for the internal baryon (antibaryon) lines. Here, L'f}"""4 denote the values of Lti(k0) for k0 = wN, -wN, q0 - wB, q0 + wB, respectively. The expression for the third diagram of Figure 17 can be obtained by changing the sign of the external momentum q in (65).

The numerical results for the p and w meson self-energy are, respectively, shown in Figures 18(a) and 18(b). The individual contribution from the meson and baryon loops is also shown for two values of the baryon chemical potential. For both p and w mesons, the small positive contribution from the baryon loops to the real part is partly compensated by the negative contributions from the meson loops which can be clearly seen in the lower panels of Figure 18.

We now use these self-energy functions in the expression for the exact propagator (60) to obtain an explicit results of in-medium spectral functions for p and w meson. In view of the fact that the p and w peaks are close to each other, it is worthwhile to compare their relative spectral strengths below their nominal masses. This is shown in Figure 19 for two values of the chemical potential. The characteristic 2n and 3n thresholds for the p and w in the vacuum case are also visible. At fixed temperature and density, the w contribution is lower than p but of comparable magnitude below their nominal masses. However, the fact that the w is suppressed by a factor ~10 (- FpF2) compared to the p in the dilepton emission rate makes a quantitative study of the w difficulty. In the above expressions the meson (H) and baryon resonances (£) have been treated in the narrow width approximation. These have then been folded with the width of the resonances as shown in [110].

----W (^B = 0) — p (^B = 0)

- w (^B = 400 MeV) p (^B = 400 MeV)

Figure 19: The w spectral function seen in comparison with the p.

Thus, the dilepton emission rate in the present scenario actually boils down to the evaluation of the self-energy graphs of p and w as a function of q0, q, temperature (T), and net baryon density (pB). Using those functions in (63) we can get a numerical estimation of dilepton static rates. With all the ingredients discussed previously, we have calculated the static emission rate of dilepton from QGP and hadronic matter. The emission rate from both the phases is plotted in Figure 20 for a given temperature of 175 MeV and baryonic chemical potential of 30 MeV. We observe significant enhancement in the dilepton yield in the mass region below the p pole

10 ^ 10

i-5 10-6 -

« ■a

10-' r

10-8 r

10-9 r

0 0.2 0.4 0.6 0.8 1 1.2 1.4 M (GeV)

---QGP

- Had (medium)

...... Had (vacuum)

Figure 20: The dilepton emission rate from different sources at T -175 MeV and uK = 30 MeV.

compared to vacuum. This rate has been used in the analysis of the dimuon spectra obtained from In-In collisions at 17.3 GeV at CERN SPS [112, 113] (discussed in Section 6.4.1). The calculations show a reasonable agreement with the invariant mass spectra for different pT ranges as well as the MT spectra for different M bins.

As indicated earlier, coupling of the hadronic current to multiparticle states gives rise to a continuum structure in the current correlation function . Following [104], we take a parameterized form for this contribution and augment the dilepton emission rate with

dN d4qd4.

= ^2)fBE (%) I A

cont V >

V=p,œ

n ) 1 + exp (w0 - q0) /S

with w0 = 1.3,1.1 GeV for p, w, and S = 0.2 for both p and w. The continuum contribution for the w contains an additional factor of 1/9.

6.3. Invariant Mass and Momentum Spectra ofDileptons. The total invariant transverse momentum distribution of thermal dileptons (l+T) is obtained as follows:

d2pTdy

i=Q,M,H

d2pTdydM2 ^

MdMd4x. (68)

In a similar manner, the invariant transverse mass distribution of thermal dileptons (l+l~) can be obtained by integrating static emission rate over certain pT window and by convoluting that by four volume and expressed as follows:

2 Md Mdy

i=Q,M,H

d2pTdydM2

pTdpTd x. (69)

The limits for integration over pT and M can be fixed judiciously to detect contributions either from quark matter or hadronic matter. Experimental measurements [79,114,115] are available for different M window.

6.4. Results and Discussion on pT and M Distributions of Dileptons. Thus far we have discussed the dilepton emission rate for a given temperature. In HIC the dilepton yield is obtained by convoluting the static emission rate over space and time.

6.4.1. Dileptons at SPS Energy. With all these ingredients the MT and M spectra of dileptons measured by NA60 collaboration at SPS energy are reproduced as follows.

Finally, we have obtained the dimuon yield (dN/dM) in In-In collisions at SPS at a center of mass energy of 17.3 GeV. The initial energy density is taken as 4.5 GeV/fm3 corresponding to a thermalisation time r' = 0.7 fm, the QGP to hadronic matter transition temperature Tc = 175 MeV, and the freezeout temperature Tj = 120 MeV (fixed from the slope of the hadronic spectra measured by the NA60 Collaboration) has been taken to compare the data measured by NA60 collaboration. In Figure 21 we have shown the invariant mass spectra for different transverse momentum (pT) windows calculated for 17.3 GeV energies.

The theoretical curves agree quite well with the experimental data [114, 115] for all the pT ranges. The strong enhancement in the low M domain is clearly due to the large broadening of the p in the thermal medium which comes entirely from the Landau cut in the self-energy diagrams. In the last panel, the blue-dashed line curve is the result of a previous calculation [113] where the self-energy due to baryons has been evaluated following the approach of [116]. In the present work [112] we have included an exhaustive set of baryon loops using the real-time thermal field theoretical approach where we have employed the full relativistic baryon propagators in which baryons and anti-baryons appear on an equal footing. This [112] is seen to be in better agreement with the experimental data [114, 115] than [113] in the range 0.35 <M< 0.65 GeV.

Apart from the M spectra, we have also evaluated the transverse mass spectra of dimuon pairs at SPS energy measured by NA60 collaboration [114,115].

The results are compared with the data obtained by NA60 collaborations [114, 115, 117, 118] at SPS energy (Figure 22). Theoretical results contain contributions from the thermal decays of light vector mesons (p, m, and $) and also from the decays of vector mesons at the freezeout [10,49] of the system has also been considered. The nonmonotonic variation of the effective slope parameter extracted from the MT spectra of the lepton pair with (M) evaluated within the ambit of the present model [113] reproduces the NA60 [114, 115] results reasonably well.

6.4.2. Dileptons at RHIC Energy. For Au+Au collisions at ■^sNN = 200 GeV, we have evaluated the dilepton spectra for different invariant mass bins with the initial condition (min bias) shown in Table 1 and lattice QCD equation of

0.6 0.8 M (GeV)

In-In semicentral = 0.6 < pT (GeV) < 0.8

0.2 0.4 0.6 0.8 1

M (GeV)

1 - 10

In-In semicentral 1 < pT (GeV) < 1.2

0.6 0.8 M (GeV)

1.2 1.4

2 •a

0.6 0.8 1 M (GeV)

2" 10-

1 •a

: 1 1 1 • i | i | i | i | i .

I For all pT

e' = 4.5, Tj = 0.7

i , i p(pole + cont) .........

0.2 0.4

0.6 0.8 M (GeV)

1.2 1.4

• NA60 dimuon data

- This work

- Eletsky et al.

• NA60 dimuon data

- This work

- Eletsky et al.

Figure 21: Dilepton invariant mass spectra for different pT bins compared with the NA60 data.

^ 10 ■H 10

^ 10 U 10

Jr 106

MT - (M) (GeV)

• 0.2 < M (GeV) < 0.4 ■ 0.6 < M (GeV) < 0.9 ▲ 1 < M (GeV) < 1.4 ♦ 0.4 < M (GeV) < 0.6

Figure 22: Transverse mass spectra of dimuons in In+In collisions at SPS energy. Solid lines denote the theoretical results [12].

x 10-4

2? 10-8 •a

0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 pt (GeV)

• 0.1 < M (GeV) < 0.2 • 0.5 < M (GeV) < 0.75 ▲ 0.2 < M (GeV) < 0.3 ▼ 0.81 < M (GeV) < 0.99 ■ 0.3 < M (GeV) < 0.5

Figure 23: Transverse momentum spectra of dileptons for different invariant mass windows for minimum bias Au-Au collisions at RHIC energy [12].

state. The results are displayed in Figure 23. The slopes of the experimental data on pT distribution of lepton pairs for different invariant mass windows measured by the PHENIX collaboration [119, 120] could be reproduced well with the same initial condition that reproduces photon spectra [79]. In fact, the reproduction of data for the mass bins 0.5 < M(GeV) < 0.75 and 0.81 < M(GeV) < 0.99 does not need any normalization factors (Figure 23). For lower mass windows slopes are reproduced well but fail to reproduce the absolute normalization. Therefore, it should be clarified here that the theoretical results shown in Figure 22 for lower mass windows (to be precise for 0.1 < M(GeV) < 0.2, 0.2 < M(GeV) < 0.3 and 0.3 < M(GeV) < 0.5) contain arbitrary normalization constant.

Assuming 10% hard (i.e., x = 0.10) and 90% soft collisions for initial entropy production the value of dNpjj/dy turns out to be about 2.43 at /s = 200 GeV. For RHIC energy, we take T' = 320 MeV with initial time r' = 0.2 fm/c which acts as inputs to the hydrodynamic evolution.

For studying thermal dileptons at the RHIC energy (as well as the LHC energy) we have included the vacuum spectral function of $ meson because its mass appears at a boundary between quark and hadronic sources of dileptons.

We begin by plotting the space-time integrated invariant mass spectra of dileptons. In Figure 24 we plot the yield of lepton pairs from the hadronic matter (HM), evaluated with and without the modified p spectral function for RHIC energy. The enhancement in the region 0.1 < M < 0.7 GeV is purely a medium effect and is a contribution from the Landau cut of the meson and baryon loops. In contrast, the vacuum spectral function naturally starts from the 2mn threshold coming from the unity in the unitary cut contribution. The (small) kink at 0.42 GeV in this curve is due to the 3mn

10 E—i-'-1-'-1-'-1-'-1-'-r

- 10 3

0.2 0.4 0.6 0.8 1 1.2 1.4 M (GeV)

HM vacuum HM medium

Figure 24: Invariant mass distribution of dileptons from hadronic matter (HM) for modified and unmodified p meson for RHIC energy.

threshold for w production. Theenhancement in theyield due to medium effects is ~20 for M around 400 MeV.

In Figure 25, we have shown the dependence of the yield from the two phases on the EoS. Dilepton radiation from hadronic phase outshines the emission from quark matter for M up to 0 mass. Since the internal loops of p self-energy contain aln and w7 interactions, we ignore the four pion annihilation process [121] to avoid double counting. The contributions from quark matter phase dominate over its hadronic counter part for both the EoS for M beyond 0

- 10-3

0.2 0.4 0.6 0.8 1 1.2 1.4 M (GeV)

- Total LQCD ...... QM HRG

----Total HRG ---Freeze out LQCD

---QM LQCD ----- Freeze out HRG

Figure 25: Freeze-out, QGP, and total contribution for EoS (a) HRG (dash-double dotted, dotted, and dash-dotted lines) and EoS (b) LQCD (bold dashed, dashed, and solid lines) at RHIC energy.

0.2 0.4 0.6 0.8 1 1.2 1.4 M (GeV)

- LQCD total ---QM LQCD

----HRG total ...... QM HRG

Figure 26: QGP and total contribution for EoS (a) HRG (dotted and dash-dotted lines) and EoS (b) LQCD (dashed and solid lines) at LHC energy.

peak. This fact may be used to extract various properties, that is, average flow, temperature, and so forth of quark matter and hadronic matter by selecting M windows judiciously. The dilepton yield from hadronic matter is observed to be larger when the HRG EoS is employed in comparison with LQCD. This can be understood in terms of the velocity of sound c^(= dP/de evaluated at constant entropy) which controls the rate of expansion. For EoS of the type (a) c2s ~ 1/3 in the QGP phase which is larger than the value of the corresponding quantity for EoS of the type (b). Therefore, the rate of expansion in the scenario (b) is comparatively slower, allowing the QGP to emit lepton pairs for a longer time resulting in greater yield for LQCD EoS. In contrast, for the EoS (a), the lower value of cs for the hadronic phase results in a slower cooling and hence a larger yield. Also shown for comparison is the yield from the decays of p mesons at the freezeout for the two types of EoS used. The yield from this source is much smaller and we will not consider it any further.

6.5. Dileptons at LHC Energy. At LHC the measured values of dN^/dy for ^sNN = 900 GeV, 2.36 TeV, and 7 TeV are 3.02, 3.77, and 6.01, respectively [122]. The value dNf^/dy at nn = 5-25 TeV is obtained by interpolating the above experimental data mentioned above. Assuming x = 0.2 in (26) we obtain dN/dy = 2607 in Pb+Pb collision for 0-10% centrality. For i = 0.1 fm/c we get = 756 MeV.

The invariant mass spectra of lepton pairs are displayed for LHC initial conditions in Figure 26. Although the results are qualitatively similar to RHIC, quantitatively the yield at LHC is larger by an order of magnitude, primarily because of the large initial temperature. This enhancement is also seen in the transverse mass distributions of the lepton pairs at LHC.

All the results presented above for photon and dilepton production are reproduced using hydrodynamical model

with cylindrical symmetry [64] and boost invariance along the longitudinal direction [41]. The two approximations, such as cylindrical symmetry and boost invariance, remain intact at very high pT and in central collision. So our results will not differ much even if (3+1) hydrodynamics is used at very high pT and in central collision. In [123], the transverse momentum spectra of photons and v2 of photons are calculated at RHIC energy using (3 + 1) hydrodynamical model. Recently, in [124, 125], a realistic (3+1) hydrodynamical model is established which can be used further to calculate the dilepton production, higher harmonics of flow of dileptons, interferometry with dileptons, and many more.

7. Radial Flow of Thermal Photons and Dileptons

The average magnitude of radial flow can be extracted from the transverse mass spectra mT(= ^pT + m2) spectra of the hadrons only at freezeout surface. However, hadrons being strongly interacting objects can bring the information of the state of the system when it is too dilute to support collectivity; that is, the parameters of collectivity extracted from the hadronic spectra are limited to the evolution stage where the collectivity ceases to exist. These collective parameters have hardly any information about the interior of the matter. On the other hand, electromagnetic (EM) probes; that is, photons and dileptons are produced and emitted [1-8] from each space time point. Therefore, estimating radial flow from the EM probes will shed light on the time evolution of the collectivity in the system.

The calculations of EM probes from thermal sources depend on the parameters such as Tt, rt, Tch, Tf, are EoS, which are not known uniquely. These abovementioned uncertainties have been used in the evaluation of individual single spectra of photon as well as dilepton. In order to

0 0.5 1 1.5 2 2.5 3 3.5 pt (GeV)

- 0.2 < M (GeV) <0.4 ---0.6 < M (GeV) < 0.9

...... 0.4 < M (GeV) < 0.6 --- 1 < M (GeV) < 1.4

Figure 27: Variation thermal photon to dilepton ratio, Rem, with pT f (b) (see text).

...... 0.1 < M (GeV) <0.2 ---0.5 < M (GeV) < 0.75

-o- 0.2 < M (GeV) < 0.3 - 0.81 < M (GeV) < 0.99

----0.3 <M (GeV) < 0.5

different invariant mass windows at SPS energy (a) and RHIC energy

overcome the uncertainties and minimize the dependence of thermal sources on these parameters, the importance of the ratio of the transverse momentum spectra of photon to dilepton (Rem) has been emphasized in the present study (see [12,16,126,127]), where the uncertainities are canceled out partially. It may be mentioned here that in the limit of M —> 0 the lepton pairs (virtual photons) emerge as real photons. Therefore, the evaluation of the ratio of the pT spectra of photons to dileptons for various invariant mass bins along with a judicious choice of the pT and M windows will be very useful to extract the properties of QGP as well as those of hadronic phase. This will be demonstrated in the present work by analyzing WA98 and PHENIX photons (results are shown in Section 5.4) and NA60 and PHENIX dilepton (results are shown in Section 6.4) spectra.

The pT spectra of photon and dilepton can be parametrized as follows:

d2pTdyJ y

= Ai( — i exp [-ciPt] ; ci = —>

( dN \

\d2pTdy)t

1 \B2 M

exp [-C2Mt] ; C2 =

where Tefli = Tav^(1 + vr)/(1 - vr) is the blue-shifted effective temperature for massless photons and Tefl, = Tav + Mvl is the effective temperature for massive dileptons. Tav is theaverage temperatureand vr is the average radial flow of the system. The Teff can be obtained by parameterizing the pT spectra of photons and dileptons (see Sections 5.4 and 6.4), respectively, with the expressed form of (70). The ratio, Rem,

for different M windows (Figure 27) can be parametrized as follows:

A(M )B exp [-c(Mt -Pt)]; c=± (71)

V PT / T eff

with different values of Teff for different invariant mass windows. The argument of the exponential in (71) can be written as [126]

MT - pT MT pT

Tav + MV22 TJ(1+ Vr)/(1 - Vr

As mentioned before some of the uncertainties prevailing in the individual spectra may be removed by taking the ratio, Rem, of the pT distribution of thermal photon to dileptons. In the absence of experimental data for both photon and dilepton from the same colliding system for SPS energies, we have calculated the ratio Rem for Pb+Pb system, where the initial condition and the EoS are constrained by the measured WA98 photon spectra. The results are displayed in Figure 27.

Also we evaluate the ratio of the thermal photon to dilepton spectra constrained by the experimental data from Au+Au collisions measured by PHENIX collaboration. The results for the thermal ratio, Rem, displayed in Figure 27(b) are constrained by the experimental data on the singlephoton and -dilepton spectra. The behavior of Rem with pT for different invariant mass windows which is extracted from the available data is similar to the theoretical results obtained

0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19

Tav (GeV)

---<M> = 0.75 GeV

- <M>= 1.2 GeV

0.14 0.16

Tv (GeV)

- <M> = 0.625 GeV

---<M> = 0.9 GeV

Figure 28: The variation of radial flow velocity with average temperature of the system for (M) = 0.75 GeV and 1.2 GeV at SPS energy (a) and for (M) = 0.625 GeV and 0.9 GeV at RHIC energy (b).

in [16,126]. It is observed that the ratio decreases sharply and reaches a plateau beyond pT > 1.5 GeV.

This behavior of Rem as a function of pT can be understood as follows: (i) for pT » M, MT - pT and consequently Rem ~ A giving rise to a plateau at large pT. The height of the plateau is sensitive to the initial temperature of the system [16,126] and (ii) for pT < M, Rem - exp(-pT/Teff)/pT indicating a decrease of the ratio with pT (at low pT) as observed in Figure 27.

7.1. Variation of Radial Flow with Average Temperature. For a given pT and M, (72) can be written as vr = f(Tav). The Teff's are obtained from the ratio of individual spectra of photon and dilepton (by parametrising the Rem in Figure 27 using (71)). Thus we obtained the variation of radial flow with average temperature (vr(Tav)) for SPS (a) and RHIC (b) that has been depicted in Figure 28. It is clear from Figure 28 that the magnitude of the flow is larger in case of RHIC than SPS because of the higher initial pressure. Because of the larger initial pressure and QGP life time the radial velocity for QGP at RHIC is larger compared to SPS.

The Teff obtained from the parametrization of ratio at SPS energy is 263 MeV and 243 MeV for M = 0.75 and 1.2 GeV, respectively. The average flow velocity vr versus Tav has been displayed for M = 0.75 GeV and 1.2 GeV in Figure 28(a). The hadronic matter (QGP) dominates the M-0.75(1.2) GeV region. Therefore, these two mass windows are selected to extract the flow parameters for the respective phases. The vr increases with decreasing Tav (increase in time) and reaches its maximum when the temperature of the system is minimum; that is, when the system attains Tf, the freezeout temperature. Therefore, the variation of vr with Tav may be

treated as to show how the flow develops in the system. The vr is larger in the hadronic phase because the velocity of sound in this phase is smaller, which makes the expansion slower, as a consequence system lives longer—allowing the flow to fully develop. On the other hand, vr is smaller in the QGP phase because it has smaller life time where the flow is only partially developed. In Figure 28(b) the variation of average transverse velocity with average temperature for RHIC initial conditions is depicted.

7.2. Variation of Radial Flow with Invariant Mass. Obtaining Teffi and Teff2 from the individual spectra and eliminating Tav one gets the variation of vr with M. Figure 29(a) shows the variation of vr with M for SPS conditions. The radial flow velocity increases with invariant mass M up to M = Mp and then drops. How can we understand this behavior? From the invariant mass spectra, it is well known that the low M (below p mass) and high M (above 0 peak) pairs originate from a partonic source [16]. The collectivity (or flow) does not develop fully in the QGP because of the small life time of this phase, which means that the radial velocity in QGP will be smaller for both low and high M. Whereas the lepton pairs with mass around p-peak mainly originate from a hadronic source (at a late stage of the evolution of system) are largely affected by the flow resulting in higher values of flow velocity. In summary, the value of vr for M below and above the p-peak is small but around the p peak is large— with the resulting behavior displayed in Figure 29. Similar nonmonotonic behavior is observed in case of elliptic flow of photon as a function of pT [128]. The variation of vr with M in RHIC (Figure 29(b)) is similar to SPS though the values

0.1.............

0 0.3 0.6 0.9 1.2

M (GeV)

0 0.2 0.4 0.6 0.8 1

M (GeV)

(a) (b)

Figure 29: The variation of radial flow with invariant mass pairs for SPS (a) and RHIC (b) energies.

of vr at RHIC are larger than those of SPS as expected due to higher initial pressure.

It is shown that simultaneous measurements of photon and dilepton spectra in HIC will enable us to quantify the evolution of the average radial flow velocity for the system, and the nature of the variation of radial flow with invariant mass indicates the formation of partonic phase at SPS and RHIC energy. The stronger radial flow at RHIC compared to SPS is due to higher initial energy densities and a longer lifetime of the reaction zone.

8. Dilepton Interferometry

The two-particle intensity interferometry, commonly known as Hanbury Brown Twiss (HBT) interferometry [129], is considered as one of the efficient methods to extract the information of space-time structure of the fireball formed in HIC. The utility of the intensity interferometry with dileptons [14,15] for extracting fireball properties is as follows. As EM radiation produces from each stages of HIC, it retains information of the fireball at each space-time point. In contrast to correlation studies with hadrons which give information of the system when the system has frozen out, two-particle intensity interferometry using lepton pairs [14,15], or photon [130, 131], can provide the information on the history of evolution of hot matter efficiently because EM probes do not rescatter after its production. As argued previously, photons appear to be more restrictive compared to dilepton, and we have attempted to do the correlation calculations with dileptons. So that with judicious choice of pT and M windows we can get the spatial as well as temporal information of QGP and the hadronic phases separately by making use of correlations between two dilepton pairs.

8.1. Formalism and Equation of Bose-Einstein Correlation Function (BECF). As interferometry of the dilepton pairs actually reflects correlations between two virtual photons, the analysis then concentrates on computing the Bose-Einstein correlation function (BECF) for two identical particles defined as

C2 (p1,p2) =

P2 (pi,p2) P (pl)Pl (p2)'

where pt is the three momenta of the particle i, and Pi(pi) and P2(pi,p2) represent the one- and two-particle inclusive lepton pair spectra, respectively, and is expressed as follows:

Pi (p)= Jd4xw(x,K),

P2 (pi,p2) = Pi (pi) Pi (p2)

+ J d4xid4x2w (xit K)

x m (x2, K) cos (Axpqp) ,

where K = (pi + p2)/2, q^ = p- p2^ = Ax^ = x1p - x2fi, x^ and p¡p are four coordinates for position and momentum variables, respectively, and m(x, K) is the source function related to the thermal emission rate of lepton pairs per unit four volume, expressed as as follows:

w (x, K) =

dM2d2KTdy'

Table 3: Values of the various parameters used in the relativistic hydrodynamical calculations.

Input RHIC LHC

dN/dy 1100 2376

T 290 MeV 640 MeV

ri 0.6 fm 0.1 fm

Tc 175 MeV 175 MeV

Tch 170 MeV 170 MeV

f 120 MeV 120 MeV

EoS 2 + 1 Lattice QCD 2 + 1 Lattice QCD

With further simplification, the C2 can be redefined as

C2 (pl,p2) = 1 +

A [J d4xw (x, K) cos (Aa) ] 3 p1 (p1)p1 (h)

[J d4xw (x, K) sin (Aa) ]

+ Pi (hi)Pi (?i)

where Aa = ai - a2, ai = TMiT cosh(y; -q)- rpir cos(d -f), MiT = ^p2j, + M2 is the transverse mass, yi is the rapidity, and f i's are the angles made by piT with the x-axis.

The inclusion of the spin of the virtual photon will reduce the value of C2 - 1 by 1/3. The correlation functions can be evaluated for different average mass windows, (M)(= Ml+l-) = (Mi + M2)l2. The leading order process through which lepton pairs are produced in QGP is qq ^ l+l~ [100, 101]. For the low M dilepton production from the hadronic phase the decays of the light vector mesons p, w, and <p have been considered including the continuum [1, 5, 7, 8, 104]. Since the continuum part of the vector meson spectral functions is included in the current work, the processes like four pions annihilations [121] are excluded to avoid double counting.

For the space time the evolution of the system relativistic hydrodynamical model with cylindrical symmetry [64] and boost invariance along the longitudinal direction [41] has been used. The values of the parameters required for spacetime evolution are displayed in Table 3. With all these ingredients we evaluate the correlation function C2 for 0-5% Au+Au collisions centrality for RHIC at ^sNN = 200 GeV

[132] and Pb+Pb collisions at for LHC at ^Nn = 2.76 TeV

[133] for different invariant mass windows as a function of qside and qout which are related to transverse momenta of individual pair [134, 135]. By choosing appropriate phase space for the QGP and hadron gas and performing the space time integration using the initial condition tabulated in Table 3, the C2 for different phase has been evaluated. We have evaluated the C2 for (M) = 0.3, 0.5, 0.7, 1.2, 1.6, and 2.5 GeV. In Figure 30 the results for only three values of (M) corresponding to low and high mass which are expected to be dominated by radiations from QGP ((M) ~ 1.6 GeV) and hadronic phase ((M)~0.77 GeV), respectively, are displayed.

In Figure 30, we plot the C2 as a function of qside and qout for RHIC initial conditions as tabulated in Table 3. A clear difference of dilepton pair mass dependence of the BEC studied as a function of qside is observed for the contributions from different M domains. The differences are however small when BEC is studied as a function of qout.

8.2. Source Dimension. The source dimensions can be obtained by parameterizing the calculated correlation function of the dilepton pairs with the empirical (Gaussian) form:

C2 (q,K)=1 + A exp (-R2 (K) q, (77)

where i stands for side, out, and long. Thus Rside, Rout, and Rlong appearing in (77) are commonly referred to as HBT radii, which is measure of Gaussian widths of source size and can be expressed as follows:

4e (K)=(f), (K) = ((x- vrt)2),

(K) = {(z- vzt)2).

The qout, qside, and qlong can be expressed in terms of individual particle momenta as [136]

qside =

q-r- qout

2pirp2r cos2 (fi -f2)

qout =

f (kir, k2r)

q-r ■ K- (pir - P2r) K

f (kir, k2r)

qlong = Piz - p2z = Pi- sinh 7i - p2r sinh y2>

where f(kir, k2r) = ^plr + p\r + 2pirp2r cos(fi - f2) and A (in this review, A = 1/3) represents the degree of chaotic of the source. The deviation of A from 1/3 will indicate the presence of nonthermal sources. A representative fit to the correlation functions is shown in Figure 30 (solid lines). While the radius (Rside) corresponding to qside is closely related to the transverse size of the system and considerably affected by the collectivity, the radius ( Rout) corresponding to qout measures both the transverse size and duration ofparticle emission [136-139]. The extracted Rside and Rout for different (M) are shown in Figures 31(a) and 31(b), respectively.

8.2.1. Variation ofRside and Rout with (M)

Rside. The variation of Rside for QGP, hadronic and QGP + hadronic phase, is obtained from the respective C2 in that phase with an appropriate selection of phase space in space-time integration. Figure 31(a) shows nonmonotonic dependence of Rside on M, starting from a value close to QGP value (indicated by the dashed line); it drops with increase in M finally again approaching the QGP value for (M) > m^. It can be shown that R^ - 1/(1 + £collectivel£thermal) [134,135].

0.05 0.1 0.15 qside (GeV/c)

o <M> = 0.3 Gev/c2 A <M> = 0.77 Gev/c2 □ <M> = 1.6 Gev/c2

0.05 0.1

qout (GeV/c)

o <M> = 0.3 Gev/c2 A <M> = 0.77 Gev/c2 □ <M> = 1.6 Gev/c2

Figure 30: Correlation function for dilepton pairs as a function of <jside ((a), for pJT = p2T = 2 GeV and y2 = 0) and qout ((b), for yj = y2 = 0 and pJT = 2 GeV) for three values of (M) [14,15]. The solid lines show the parameterization of C2 using (77).

<M> (GeV)

(fm) 4.2

1 1 1 1 1 1 1 1

it-^ Q /

- ° O /

..........

O No p and o>

...... QGP

...... Hadronic

<M> (GeV)

Figure 31: (a) Rside is evaluated with pJT = p2T = 2 GeV and y2 =0 as a function of (M) for RHIC energy and (b) Rout is evaluated with fj =V2 =0 and pJT = 2 as a function of (M) for RHIC [14,15].

In the absence of radial flow, ,Rside is independent of ^side. With the radial expansion of the system a rarefaction wave moves toward the center of the cylindrical geometry; as a consequence the radial size of the emission zone decreases with time. Therefore, the size of the emission zone is larger at

early times and smaller at late time. The high (M) regions are dominated by the early partonic phase where the collective flow has not been developed fully; that is, the ratio of collective-to-thermal energy is small, and hence the source has larger ^side. In contrast, the lepton pairs with M ~ mp

— LHC ---RHIC

(M> (GeV)

LHC --- RHIC

(M> (GeV)

Figure 32: Rside (a) is evaluated with pJT = p2T = 1 GeV and y2 = 0 and Rout (b) is evaluated with pJT = 1 GeV and yj = y2 = 0 as a function of (M) for RHIC (dashed line) and LHC (solid line) energies [14,15].

are emitted from the late hadronic phase where the size of the emission zone is smaller due to larger collective flow giving rise to a smaller .Rside. The ratio of collective to thermal energy for such cases is quite large, which is reflected as a dip in the variation of ^side with (M) around the £>-mass region (Figure 31(a)). Thus the variation of ^side with M can be used as an efficient tool to measure the collectivity in various phases of matter. The dip in Rside at (M) ~ mp is due to the contribution dominantly from the hadronic phase. We observe that, by keeping the p and w contributions and setting radial velocity, vr = 0,the dipin ^side vanishes, confirming the fact that the dip is caused by the radial flow of the hadronic matter. Therefore, the value of ^side at (M) ~mp maybe used to estimate the average vr in the hadronic phase.

Rout. The Rout probes both the transverse dimension and the duration of emission, and unlike ^side it does not remain constant even in the absence of radial flow. As a result its variation with M is complicated. The values Rout for different phases are obtained in a similar fashion as followed for obtaining the ^side values for the different phases. The large M regions are populated by lepton pairs from early partonic phase where the effect of flow is small and the duration of emission is also small—resulting in smaller values of Rout. For lepton pair from M ~ mp the flow is large which could have resulted in a dip as in ^side in this M region. However, Rout probes the duration of emission too which is large for hadronic phase because the expansion is slower in this phase for the EoS used in the present work The velocity of sound which controls the rate of expansion and hence the

duration of the phase has larger value in hadronic phase than in the partonic phase. Thus resulting in the larger Rout in the hadronic phase than that in partonic phase, the larger duration compensates the reduction of Rout due to flow in the hadronic phase resulting is a bump in Rout in this region of M (Figure 31(b)). Again the duration of particle emission from both the phases obviously is larger than that from the individual phases.

Both ^side and Rout approach QGP values for (M) ~ 2.5 GeV implying dominant contributions from partonic phase.

8.3. Comparison ofHBT Radii with Different Collision Energies. Now we study the sensitivity of the HBT radii on the different collision energy. The ^side and Rout extracted from the C2's evaluated for 0-5% centrality in Au+Au collisions for RHIC at = 200 GeV [132] and Pb+Pb collisions

for LHC at Vsnn = 2.76 TeV [133] for different invariant mass windows as a function of qside and qout are shown in Figure 32. The change of ^side with (M) for RHIC and LHC is qualitatively similar but quantitatively different. The smaller values of ^side for LHC are due to the larger radial expansion which can be understood from the fact that the quantity -Ecollective/-Ethermal is larger at LHC than RHIC. So, the dip in the ^side variation at LHC is below than that at RHIC confirming a larger flow at LHC than RHIC. As the Rout probes both the transverse size and the duration of emission, from the previous discussion in Section 8.2.1, the larger duration compensates the reduction of Rout due to that flow resulting is a bump in Rout for M ~ mp.

--- RHIC

Figure 33: Variation of (pT) as function of M (see (80)).

Though the duration of particle emission is more at LHC compared to RHIC (shown in Figure 36), the larger flow (corresponds to smaller size) at LHC [126] than that of RHIC compensates other factor (like duration of emission) which has an enhancing effect on Rout. So the value Rout at LHC is smaller than that of RHIC.

8.4. Radial Flow from HBT Radii. According to the discussion given in the Section 8.2, Rside is independent of ^side in the absence of radial flow. Rside is related to radial flow as follows:

Rside (M) =

{pT (M))'

(Pt (M))=Tav + Mv2r. (80)

The values of Rside(M) are obtained for different (M) windows (shown in Figure 31). The higher mass, that is, (M) = 2.5 GeV, corresponds to the initial stage of collision where the flow is not developed fully. So assume for (M) = 2.5, vr = 0, and T = Ti the value of K = Ti x Rside | (m)=25. Once we know the value of K, we can calculate the (pT(M))(= K/Rside(M)) The variation of (pT) with (M) has been displayed in Figure 33.

The high (M) regions are dominated by the early partonic phase where the collective flow has not been developed fully and hence show smaller (pT). In contrast, due to larger collective flow for the lepton pairs with M~mp, emitted from the late hadronic phase, (pT) is larger. The larger value of (pT) around the ^-mass region is due to the contribution of large flow in the hadronic phase. Thus the variation of Rside with M (Figure 31) can be used as an efficient tool to measure the collectivity in various phases of matter.

8.5. Sensitivity of HBT Radii on pir. In this section, the sensitivity of the HBT radii for different values of the individual transverse momentum of the pairs is described. In Figure 34, the variation of Rside and Rout with (M) is shown for pir = 1 and 2 GeV. The lepton pairs coming from higher pr and high mass region enable us to quantify the size of hotter zone. As mentioned before, the pr contains the effect of flow as well as thermal motion. Hence the larger Rside at M - Mp for pr = 2 GeV is associated with longer flow and hence smaller source size. The observed bump in Rout (Figure 34(b)) is resulted from the fact that it contains both the size of the system as well as the duration of dilepton emission as discussed earlier.

8.6. Duration of Particle Emission from HBT Radii. The HBT radii, Rout and Rside, provide the information of average source size. However, in the ratio, Rout|Rside, some of the uncertainties associated with the space time evolution get canceled out. The quantity, Rout|Rside gives the duration of particle emission [134,135,140,141] for various domains of M. The difference between R2ide(K) and Rout(K) at nonzero in K is then only due to the explicit K dependence in (78), that is, the term vr (t2). This implies that the explicit K dependence dominates if the emission duration is sufficiently large or if the position-momentum correlations in the source are sufficiently weak:

Rdiff =R2out (K)-Rlde (K) = Vr (t2).

In this case, the difference between these two HBT radius parameters gives direct access to the average emission duration (t2) of the source and allows to partially disentangle the spatial and temporal information contained in (78).

Figure 35 shows the Rout|Rside and the difference

^R2out - R2de as a function of (M) for Au+Au collisions at ■^sNN = 200 GeV. Both show a nonmonotonic dependence on (M) .Thesmaller values ofboth thequantities, particularly at high mass region, reflect the contributions from the early partonic phase of the system. The peak around ^-meson mass reflects dominance of the contribution from hadronic phase as discussed before. Figure 36 shows a comparative study of the above two quantities (the ratio and the difference of Rout and Rside) for RHIC and LHC energies. They reflect a larger life time of thermal system for LHC than RHIC.

9. Elliptic Flow of Thermal Dileptons

It has been argued that the anisotropic momentum distribution of the hadrons can bring the information on the interaction of the dense phase of the system [142] despite the fact that the hadrons are emitted from the freezeout surfaces when the system is too dilute to support collectivity. Therefore, a suitable dynamical model is required to extrapolate the final hadronic spectra backward in time to get the information about the early dense phase. Such an extrapolation is not required for lepton pairs because they are emitted from the entire space-time volume of the system. Therefore, the v2 of lepton pairs provides information of the hot and dense

1 1.5 2

CM) (GeV)

Pit = p2T = 1 GeV Pit = p2T = 2 GeV

1 1.5 2

CM) (GeV)

- p1T = 2 GeV

---piT = 1 GeV

Figure 34: (a) Shows Rside as a function of (M) which is evaluated with pJT = p2T = 1 and 2 GeV and y2 = 0 and similarly (b) shows Rout as a function of for yj = y2 = 0 and pJT = 1 and 2 GeV (M). The pJT = 1 and 2 GeV results are shown as dashed line and solid line, respectively.

phase directly. The v2 of dileptons can also be used to to test the validity and efficiency of the extrapolation required for hadronic v2. The v2 of real photons and dileptons [128,143146] has been evaluated for RHIC energies and shown that it can be used as effective probes to extract the properties of the partonic plasma. The sensitivity of the v2 of lepton pairs on EoS has been elaborated in [145] for RHIC collision conditions. The lepton pairs are produced from each space time point of the system and hence the study of v2 of lepton pairs will shed light on the time evolution of collectivity in the system [12,147]. The radial flow alters the shape of the pT spectra of dileptons; it kicks the low pT pairs to the higher pT domain, making the spectra flatter. Therefore, the presence of large radial flow may diminish the magnitude of v2 at low pT [37, 38], and this effect will be larger when the radial flow

is large, that is, in the hadronic phase which corresponds to lepton pairs with M~mp.

9.1. Formalism ofElliptic Flow ofDilepton. The elliptic flow of dilepton, v2, can be defined as

V2 (pT>M)

= (cos 2<f)

L=q,h J cos (2cf) (dNy* /d2pTdM2dy\y=0)

1=00 \ (dNr*/d2pTdM2dy\y=0)d$ '

---RHIC

(M> (GeV)

--- RHIC

(M> (GeV)

Figure 36: The ratio Rout/Rside (a) and the difference - R2ide (b) as a function of (M) for RHIC (dashed line) and LHC (solid line) energies are shown.

where the £ stands for summation over quark matter (QM) and hadronic matter (HM) phases. The quantity dN/d2 prdM2 dy\y =0 appearing in (82) can be obtained from the dilepton production per unit four volume, dN/d4pd4x in a thermalized medium by integrating over the space-time evolution of the system. x for lepton pairs

for QGP and hadrons are discussed in Sections 6.1 and 6.2, respectively (see [13, 109] for details).

To evaluate v2 from (82) one needs to integrate the fixed temperature production rate given by (9) over the space time evolution of the system—from the initial QGP phase to the final hadronic freezeout state through a phase transition in the intermediate stage. The space-time evolution is done over the 4-volume, which is defined as d x(= rdrdxdydq) are expressed in terms of x^ = (r,x,y,q). We assume that the matter is formed in QGP phase with zero net baryon density in Pb+Pb collision at ^sNN = 2.76 TeV. The energy of the lepton pair (p0) should be replaced by its value in the comoving frame of the expanding system which is given by [p ■ u = yr(Mr cosh(y - q) - vxpr cos0 -vypr sin <f)],p^ = (M- cosh y,pr cos ^,pr sin$,Mr sinhy), and U = yr(coshq, vx, vy, sinhq). The EoS required to close the hydrodynamic equations is constructed by complementing Wuppertal-Budapest lattice simulation [48] with a hadron resonance gas comprising all the hadronic resonances up to mass of 2.5 GeV [148, 149]. The necessary initial conditions to solve the hydrodynamic equations are Ti = 456 MeV, the value of the temperature corresponding to the maximum of the initial energy profile for 30-40%

centrality at -^sNN = 2.76 TeV, with ri = 0.6 fm/c, the thermalization time. The transition temperature, Tc,for quark hadron conversion is taken as 175 MeV. The system is assumed to get out of chemical equilibrium at T = Tch = 170 MeV

[69]. The kinetic freezeout temperature TF = 130 MeVisfixed

from the pr spectra of the produced hadrons at the same collision energy of Pb+Pb system. The EoS and the values of the parameters mentioned above are constrained by the pr spectra (for 0-5% centrality) and elliptic flow (for 1050% centrality) of charged hadrons [148] measured byALICE collaboration [150,151].

9.2. Results and Discussion. In Figure 37 we depict the constant temperature contours corresponding to Tc = 175 MeV and Tj = 130 MeV in the r-x plane (at zero abscissa) indicating the boundaries for the QM and HM phases, respectively.

The life time of the QM phase -6 fm/c and the duration of the HM are -6-12 fm/c. Throughout this work by early and late will approximately mean the duration of the QM and HM, respectively.

With all the ingredients mentioned above we evaluate the pr integrated M distribution of lepton pairs originating from QM and HM (with and without medium effects on the spectral functions of p and w). The results are displayed in Figure 38 for the initial conditions and centrality mentioned above. We observe that for M > M^ the QM contributions dominate. For Mp < M < M^ the HM shines brighter than QM. For M < Mp, the HM (solid line) over shines

O.B GeV

O.l75 GeV

0 2 4 6 8

x (fm)

Figure 37: Constant temperature contours denoting space-time boundaries of the QGP and hadronic phase.

Tav (fm)

---(M) = O3 GeV (no medium)

- (M) = O3 GeV

-•- (M) = O.77 GeV

(M) = l.O2 GeV (M) = 2.5 GeV

Figure 39: Fractional contribution of lepton pairs for various invariant mass windows as a function of average proper time (see text for details) [13].

M (GeV)

---Quark matter

- Hadronic matter (medium)

...... Hadronic matter (no medium)

Figure 38: Invariant mass distribution of lepton pairs from quark matter and hadronic matter [13].

the QM due to the enhanced contributions primarily from the medium-induced broadening of p spectral function. However, the contributions from QM and HM become comparable in this region of M if the medium effects on p spectral function are ignored (dotted line). Therefore, the results depicted in Figure 38 indicate that a suitable choice of M window will enable us to unravel the contributions from a particular phase (QM or HM).

To further quantify these issues we evaluate the following quantity:

J' (dN/d4xd2pTdM2dy) dxdyd^TdTd2pTdM2 J (dN/d4xd2pTdM2dy) dxdyd^rdrd2pTdM2

where the M integration in both the numerator and denominator is performed for selective M windows from M1 to M2 with mean M defined as (M) = (M1 + M2)/2. The prime in J in the numerator indicates that the t integration in the numerator is done from t1 = Tt to t2 = Tt + At with progressive increment of At, while in the denominator the integration is done over the entire lifetime of the system. In Figure 39, F is plotted against rav(= (t1 + t2)/2). The results substantiate the fact that pairs with high (M) ~ 2.5 GeV originate from QM (rav < 6 fm/c, QGP phase) and pairs with (M) ~ 0.77 GeV mostly emanate from the HM phase (rav > 6 fm/c). The change in the properties of p due to its interaction with thermal hadrons in the bath is also visible through F evaluated for (M) ~ 0.3 GeV with and without medium effects. This clearly indicates that the (M) distribution of lepton pairs can be exploited to extract collectivity of different phases of the evolving matter.

Figure 40(a) shows the differential elliptic flow, v2(pT), of dileptons arising from various (M) domains in quark matter. Similarly Figure 40(b) shows the differential elliptic flow, v2(pT), of dileptons arising from various (M) domains from hadronic matter. The individual v2 for QM and HM is obtained by doing an integration over specific invariant masses (M) window as well as space time integration over the regime where Tc < T(t, x, y) < Tt and Tj < T(t, x, y) < Tc, respectively. The v2 is small at low pT and gradually increases

- (M> = 140 MeV

-o- (M> = 300 MeV

...... (M> = 500 MeV

-*- (M> = 770 MeV

----(M> = 940 MeV

-•- (M> = 1020 MeV

---(M> = 1800 MeV

-o- (M> = 2500 MeV

i i i i i i i i i ................ Pb+Pb @ 2.76 TeV I

J 30-40%

Hadron I

L .........

0 1 2 pr (GeV)

- (M> = 140 MeV -.- (M> = 1020 MeV

-*- (M> = 770 MeV ---(M> = 1800 MeV

----(M> = 940 MeV -0- (M> = 2500 MeV

(a) (b)

Figure 40: Elliptic flow of quark matter (a) and hadronic matter (b) as a function of pr for various mass windows.

i-1-1-1-1-1-1-1-1-1-1-1-1-r

(M> = 0.3 GeV - (M> = 0.49 GeV - (M> = 0.77 GeV

l-1-1-1-1-1-1-1-1-1-1-1-1-r

■■*■■ (M> = 0.94 GeV (M> = 1.8 GeV

■-■- (M> = 1.02GeV -o- (M> = 2.5GeV

(a) (b)

Figure 41: Total elliptic flow as function of pr for various mass windows.

and attains large value around pr - 2-3 GeV/c. Also there is clear mass ordering that has been observed for v2(pr) for QM; that is, v2 decreases with increase in M. This is because dileptons come from high M region, M > M^, come mostly from hot partonic phase where the fluid velocity is not strong to support the collectivity but the spatial eccentricity

of the source is large. On the other hand dileptons that come from low M region, M below 0 peak dominantly come from late hadronic matter where the collectivity is strong and the spatial asymmetry dissolve into momentum asymmetry.

Figure 41 shows the differential elliptic flow, v2(pr) of dileptons arising from various (M) domains. We observe

that, for {M) = 2.5 GeV, v2 is small for the entire pT range because these pairs arise from the early epoch (see Figure 39) when the flow is not developed entirely. However, the v2 is large for {M) = 0.77 GeV as these pairs originate predominantly from the late hadronic phase when the flow is fully developed.

It is also interesting to note that the medium-induced enhancement of p spectral function provides a visible modification in v2 for dileptons below p peak. Figure 42 shows the comparison between v2(pT) of dilepton at {M) = 300 MeV with and without medium effects.

In Figure 43 we depict the variation of Rq with pT for {M) = 0.3 GeV (line with solid circle) 0.77 GeV (solid line) and 2.5 GeV (line with open circle). The quantity Rq and RH is defined as

where v2QM and v2HM are the elliptic flow of QM and HM phases, respectively. The results clearly illustrate that v2 of lepton pairs in the large (M)(= 2.5 GeV) domain (open circle in Figure 43) originates from QM for the entire pT range considered here. The value of Rq is large in this domain because of the large (negligibly small) contributions from QM (HM) phase. /qm is large here. It is also clear that the contribution from QM phase to the elliptic flow for (M)(= 0.77 GeV) is very small (solid line in Figure 43). The value of RH for (M) = 0.77 GeV is large (not shown in the figure). The v2 at the (late) hadronic phase (either at p or <p peak) is larger than the (early) QGP phase (say at (M) = 2.5 GeV) for the entire pT range considered here. Therefore, the pT-integrated values of v2 should also retain this character of v2 at the corresponding values of (M). It is also important to note that the differential elliptic flow, v2(pT), obtained here at LHC is larger than the values obtained at RHIC [143-145] for all the invariant mass windows. The value of RH for (M) = 0.77 GeV is large (not shown in the figure).

The v2 at the HM phase (either at p or <p peak) is larger than its value in the QGP phase (say at (M) = 2.5 GeV) for the entire _pT-range considered here. Therefore, the pT integrated values of v2 should also retain this character at the corresponding values of (M), which is clearly observed in Figure 44 which displays the variation of v2((M)) with (M). The v2(<x ep) of QM is small because of the small pressure gradient in the QGP phase. The v2 resulting from hadronic phase has a peak around p pole indicating the full development of the flow in the HM phase. For (M) > m^ the v2 obtained from the combined phases approachs the value corresponding to the v2 for QGP. Therefore, measurement of v2 for large (M) will bring information of the QGP phase at the earliest time of the evolution. It is important to note that the _pT-integrated v2((M)) of lepton pairs with (M) ~mn, mK is close to the hadronic vj and v^ (symbol * in Figure 44) if the thermal effects on p properties are included. Exclusion of

Pt (GeV)

w/o medium w medium

Figure 42: The figure displays the effect of the broadening of p spectral function on the elliptic flow for (M) = 300 MeV.

12 Pt (GeV)

---CM) = 2.5 GeV

-•- (M) = 0.3 GeV - (M) = 0.77 GeV

Figure 43: It shows the variation of Rq (see text) with pT for (M) 0.3 GeV, 0.77 GeV and 2.5 GeV

medium effects gives lower v2 for lepton pairs compared to hadrons. The fact that the v2 of the (penetrating) lepton pairs is similar in magnitude to the v2 of hadrons for ((M) ~ mn, mK, mproton, etc.) ascertains that the anisotropic momentum distribution of hadrons carries the information of the HM phase with duration ~6-12fm/c [13]. We also observe that the variation of v2((M)) with (M) has a structure similar to dN/dM versus M. As indicated by (82) we can write

0.36 0.32 0.28 0.24 0.2 0.16 0.12 0.08 0.04 0

1 1 1 1 1 1 1 1 1 1 1 1 1 : Pb+Pb @ 2.76 TeV ; 30-40% 1 1 1 1 1 1 1 i i | i i i

- S'C"" m _ H \\ ----- _____0 -

VP? - n J*' T - i-Vi V.

- 0...o-............. ....... ......

-o- HM -o - QM

1.2 1.6 CM) (GeV)

Total (no medium)

Figure 44: (Color online) Variation of dilepton elliptic flow as function of (M) for QM, HM (with and without medium effects), and for the entire evolution. The symbol * indicates the value of v2 for hadrons, for example, n, kaon, proton, and (.

V2((M)) - Z;=qm,hm v2 x fi, where f is the fraction of QM or HM from various space-time regions. The structure of dN/dM is reflected in v2((M)) through f We find that the magnitude of v2((M)) at LHC is larger than its value at RHIC [143-145].

10. Summary

In this review, we have attempted to summarize the photon and dilepton spectra measured at different beam energies by different experimental collaborations which have been analyzed to understand the evaluation of collectivity in the system. The initial conditions of the evolving matter required to calculate the photon and dilepton spectra have been constrained to reproduce the measured multiplicity in these collisions. The EoS, the other crucial input to the calculations, has been taken from lattice QCD calculations. The deviation of the hadronic phase from chemical equilibrium is taken into account by introducing nonzero chemical potential for each hadronic species. For dilepton emission from hot hadrons, the spectral function of p, w at finite temperature and ( in vacuum is considered. We have also observed the dissipative effects on the photon spectra by taking into account the viscous contribution to phase space factors of all the participating partons/hadrons in the photon production rate as well as to the expansion dynamics with finite shear viscosity.

It is shown that simultaneous measurements of photon and dilepton spectra in heavy ion collisions will enable us to quantify the evolution of the average radial flow velocity for the system, and the nature of the variation of radial flow with invariant mass will indicate the formation of partonic phase. All this study suggests with judicious choice of pT and

M window will infer information of partonic and hadronic phase separately.

In this work, we present a new proposal for carrying out an experimental measurement of dilepton interferometry both for RHIC and LHC. We establish through a hydrody-namical model-based space-time evolution the promise of such a dilepton interferometry analysis that will be useful to understand the properties of the partonic phase. We have evaluated the correlation function, C2, for two dilepton pairs for various invariant mass domains and extracted the HBT radii, that is Rside and Rout as a function of M. These HBT radii show a nonmonotonic dependence on the invariant mass, reflecting the evolution of collective flow in the system which can be considered as a signal of the QGP formation in heavy ion collisions. The M dependence of the Rout/Rside

and ^R2out - R2ide which can be experimentally measured could be used to characterize the source properties at various instances of the evolution.

We have evaluated the v2 of dileptons originating from the Pb+Pb collisions at ^s^ = 2.76 TeV for 30-40% centrality. Our study shows that v2(M) provides useful information on the collective motion of the evolving QCD matter formed in high energy heavy-ion collisions. The present work indicates that experimental observation of the reduction of v2(M) with increasing M beyond ( mass would reflect the presence of small momentum space anisotropy through small collective motion in the partonic phase. We observe that v2 ((M)) of the penetrating probe (lepton pairs) for (M) = mn and mK is similar to the hadronic vj and vf when the medium-induced change in the p spectral function is included in evaluating the dilepton spectra. Since the medium effects are large during the dense phase of the system, therefore, this validates the statement that the hadronic v2 carries the information of the early dense phase of the collisions. Our study also establishes the fact that the invariant mass dependence of dilepton v2 can in principle act as a clock for the space time evolution of the system formed in HIC.

Some comments on effect of magnetic field on photon and dilepton production are in order here. The magnetic contribution is significant when photons propagate perpendicular to magnetic field or in reaction plane in noncentral collision [152]. Recently, Tuchin has implemented the magnetic field effect on the photon as well as dilepton production which is summarized here. Recently in [153, 154], it has been argued that photons, in particular thermal and direct photons from earliest times of QGP, will be affected by the magnetic field. In dilepton sector, the magnetic contribution is important at lower electron energies (below -0.5 GeV at midrapidity at RHIC). In fact, it becomes dominant source of dileptons at lower electron energies. For LHC energies, the effect of magnetic field on dilepton should be taken seriously because 40% of energy density of the QGP resides in the strong magnetic field (-1018 Gauss [152]), which can be generated in heavy ion collision at LHC. Whereas in RHIC energy this fraction is within 3% so can be safely ignored. Some more realistic calculation and hydrodynamic model with magnetic field are required to conclude anything in future.

Appendix Phase Space

In this appendix we derive (37) from (12). The photon production rate from the process, 1 + 2 ^ 3 + y is given by

„dR 1 N

d3p 2 (2n)8

jf Îf if ^ ^ M

x[1±f3 (£3)] \M\2S X (Pl +P2 -P3 -P).

Performing the d3p3 integration using the delta function and using d3p/E = pTdpTdydfa we get

„dR 1 N

' d3p 16 (2n)8

j p1Tdp1Tdy1d<t>1p2Tdp2Tdy2dÏ2

x±fi (E.)f2 (E2)[l±f3 (E3)] (A.2) h

x \M\2S(El+E2 -E3 -E),

where fa and fa2 are the angles made by the transverse momenta of first and second particles with the transverse momentum of the emitted photon. The momentum conservation along the z-direction: p3z = p.z + p2z - pz can be written in terms of rapidity as

m3T sinh y3 = m1T sinh y. + m2T sinh y2 - pT sinh y.

Now the energy, E3, can be written as

E3 = m3T cosh y3 = vm3T + m^sinh Substituting (A.3) in (A.4) we get

r, \2 2 11/2

E3 = [ (m1T sinh y1 + m2T sinh y2 - pT sinh y) + m3T \ .

Considering the energy conservation (E3 = E1 + E2 - E) and writing the energies in terms of rapidity (Et = miT cosh yt) we get

E3 = m1T cosh y1 + m2T cosh y2 - pT cosh y. (A.6) Equating (A.5) and (A.6) we have m3T = [m2^ + + p? + 2m1Tm2T cosh (y1 - y2)

-2m1TpT cosh (yl - y) - 2m2TpT cosh (y2 -

However, we also have

m3T = {pit + md

= [(pit + P2T -pt)2 +ml]112 = [P21T + pit + pt + 2pitp2t cos (fai2)

-2PtPit cos (fa.) - 2pTp2T cos ($2) + m23]112,

cos ($12) = cos (faj) cos (fa) + sin (fa) sin (fa). (A.9) Equating (A.7) with (A.8) leads to the expression: [(p1T cos fa1 - pi) cos fa + p1T sin fa1 sin fa2]

1— [(ml + m2 - m23) + 2m1Tm2T cosh (y1 - y2)

- 2mlTpT cosh (y. - y) - 2m2TpT cosh (y2 - y) +2pTp.T cos fa. ].

(A.10)

Solving (A.10) for fa one gets

i0 t- -1

fa2 = tan

p1T sin fa. p1T cos fa. - pT

(A.11)

R = ^p2t + pI - 2p1TpT cos fa1, (A.12)

H = (ml + m2 - m2) + 2m1Tm2T cosh (y1 - y2)

- 2m1TpT cosh (y1 - y) - 2m2TpT cosh (y2 - y) + 2pTp1T cos fa1.

(A.13)

Now we express the argument of the delta function in (A.2) as function of fa2 as

— Ej + E

= m1T cosh y1 + m2T cosh y2 - pT cosh y

+ (m1T sinh y. + m2T sinh y2 - pT sinh y)

2]1I2 (A.14)

and performing the integration in (A.2) we get [13

„ dR 1 N (œ (œ ,

dp = Ï6{2nf Jo ^^ Jo dp2T

\p1T sin ($1 -$2) + pT sin $2 \f2

with the constraint \H/2Rp2T\ < 1 originating from \ cos($)\ < 1.

References

Îœ rœ r 2n

dyi \ dy2 \ d$i -œ J-œ Jo (A.15)

xfi (Ei)f2 (E2)[l±fs (E3)] [15

[16 [17

Acknowledgments [18

The authors thank Sourav Sarkar and Jan E. Alam for their [19 major collaborative contribution in these works presented here and also for comments and suggestions to complete the paper. Payal Mohanty and Sabyasachi Ghosh thank Victor Roy, Santosh K. Das, Jajati K. Nayak, Bedangadas Mohanty, [20 and Asis K. Chaudhuri for their assistance in the work presented in the review.

[1] L. D. McLerran and T. Toimela, "Photon and dilepton emission [22 from the quark-gluon plasma: some general considerations," Physical Review D, vol. 31, pp. 545-563,1985.

[2] C. Gale and J. Kapusta, "Dilepton radiation from high tempera- [23 ture nuclear matter," Physical Review C, vol. 35, no. 6, pp. 21072116, 1987.

[3] C. Gale and J. I. Kapusta, "Vector dominance model at finite temperature," Nuclear Physics B, vol. 357, no. 1, pp. 65-89,1991. [24

[4] H. A. Weldon, "Reformulation of finite-temperature dilepton production," Physical Review D, vol. 42, pp. 2384-2387,1990.

[5] R. Rapp and J. Wambach, "Chiral symmetry restoration and [25 dileptons in relativistic heavy-ion collisions," Advances in Nuclear Physics, vol. 25, pp. 1-205, 2000.

[6] "Elementary particles, nuclei and atoms," Relativistic Heavy Ion [~6 Physics, vol. 23, pp. 134-175, 2010. [

[7] J.-E. Alam, S. Raha, and B. Sinha, "Electromagnetic probes of quark gluon plasma," Physics Reports, vol. 273, no. 5-6, pp. 243362, 1996.

[8] J.-E. Alam, S. Sarkar, P. Roy, T. Hatsuda, and B. Sinha, "Thermal photons and Lepton pairs from Quark Gluon plasma and hot hadronic matter," Annals of Physics, vol. 286, no. 2, pp. 159-248, 2000.

[9] R. Chatterjee, L. Bhattacharya, and D. K. Shrivastava, "Electro- [29 magnetic probes," Lecture Notes in Physics, vol. 785, pp. 219-264,

2010. [30

[10] H. van Hees and R. Rapp, "Dilepton radiation at the CERN super-proton synchrotron," Nuclear Physics A, vol. 806, no. 14, pp. 339-387, 2008. [31

[11] R. Rapp, "Theory of soft electromagnetic emission in heavy-ion collisions," Acta Physica Polonica B, vol. 42, no. 12, p. 2823, 2011.

[12] P. Mohanty, J. K. Nayak, J.-E. Alam, and S. K. Das, "Radial flow [32 from electromagnetic probes and signal of quark gluon plasma," Physical Review C, vol. 82, no. 3, Article ID 034901, 2010.

P. Mohanty, V. Roy, S. Ghosh, S. K. Das, B. Mohanty et al., "Elliptic flow of thermal dileptons as a probe of QCD matter," Physical Review C, vol. 85, Article ID 031903, 2012. P. Mohanty, J.-E. Alam, and B. Mohanty, "Evolution of collectivity as a signal of quark gluon plasma formation in heavy ion collisions," Physical Review C, vol. 84, no. 2, Article ID 024903, 2011.

P. Mohanty, J.-E. Alam, and B. Mohanty, "Characterizing quark gluon plasma by dilepton interferometry," Nuclear Physics A, vol. 862, pp. 301-303, 2011.

J. K. Nayak, J.-E. Alam, S. Sarkar, and B. Sinha, "Thermalphoton to dilepton ratio in high energy nuclear collisions," Physical Review C, vol. 78, no. 3, Article ID 034903, 2008. J. I. Kapusta and E. V. Shuryak, "Weinberg-type sum rules at zero and finite temperature," Physical Review D, vol. 49, pp. 4694-4704, 1994.

C. Y. Wong, Introduction of High Energy Heavy Ion Collisions, World Scientific, Singapore, 1994.

J.-E. Alam, D. K. Srivastava, B. Sinha, and D. N. Basu, "Transverse flow effects on high-energy photons emitted by expanding quark-gluon plasma," Physical Review D, vol. 48, pp. 1117-1131, 1993.

J. Ruppert, C. Gale, T. Renk, P. Lichard, and J. I. Kapusta, "Low mass dimuons produced in relativistic nuclear collisions," Physical Review Letters, vol. 100, Article ID 162301, 2008. S. D. Drell and T.-M. Yan, "Massive lepton-pair production in hadron-hadron collisions at high energies," Physical Review Letters, vol. 25, no. 5, pp. 316-320,1970. S. D. Drell and T. M. Yan, "Massive Lepton-pair production in Hadron-Hadron collisions at high energies," Physical Review Letters, vol. 25, p. 902, 1970.

S. Gavin, P. L. McGaughey, P. V. Ruuskanen, and R. Vogt, "Lepton production from charm decay in nuclear collisions at ^s = 200 GeVand5.5TeVpernucleon," Physical Review C,vol. 54, no. 5, pp. 2606-2623,1996.

N. Armesto, N. Borghini, S. Jeon, U. A. Wiedemann et al., "Heavy-ion collisions at the LHC-Last call for predictions," Journal of Physics G, vol. 35, Article ID 054001, 2008. R. Rapp, "Duality and chiral restoration from low-mass dilep-tons at the CERN-SpS," Nuclear Physics A, vol. 661, no. 1-4, Article ID 990734, pp. 33-44,1999.

R. Rapp and J. Wambach, "Low-mass dileptons at the CERN-SpS: evidence for chiral restoration?" The European Physical Journal A, vol. 6, pp. 415-420,1999.

E. L. Feinberg, "Direct production of photons and dileptons in thermodynamical models of multiple hadron production," Nuovo Cimento A, vol. 34, pp. 391-412, 1976.

S. Sarkar, "Photons from ultra-relativistic heavy ion collisions," http://arxiv.org/abs/nucl-th/0007011.

A. Das, Finite Temperature Field Theory, World Scientific, Singapore, 1997.

R. Kobes and G. Semenoff, "Discontinuities of green functions in field theory at finite temperature and density," Nuclear Physics

B, vol. 260, no. 3-4, pp. 714-746,1985.

R. Kobes and G. Semenoff, "Discontinuities of green functions in field theory at finite temperature and density (II)," Nuclear Physics B, vol. 272, no. 2, pp. 329-364,1986.

F. Gelis, "Cutting rules in the real-time formalisms at finite temperature," Nuclear Physics B, vol. 508, no. 1-2, pp. 483-505, 1997.

[33] S. Sarkar, J.-E. Alam, P. Roy, A. K. Dutt-Mazumder, B. Dutta-Roy, and B. Sinha, "Photons from hadronic matter at finite temperature," Nuclear Physics A, vol. 634, no. 1-2, pp. 206-230, 1998.

[34] T. Hirano and Y. Nara, "Hydro dynamic afterburner for the color glass condensate and the parton energy loss," Nuclear Physics A, vol. 743, no. 4, pp. 305-328, 2004.

[35] D. Teaney, J. Lauret, and E. V. Shuryak, "Flow at the SPS and RHIC as a Quark-Gluon plasma signature," Physical Review Letters, vol. 86, pp. 4783-4786, 2001.

[36] P. F. Kolb, P. Huovinen, U. W. Heinz, and H. Heiselberg, "Elliptic flow at SPS and RHIC: from kinetic transport to hydrodynamics," Physics Letters B, vol. 500, no. 3-4, pp. 232240, 2001.

[37] P. Huovinen, P. F. Kolb, U. Heinz, P. V. Ruuskanen, and S. A. Voloshin, "Radial and elliptic flow at RHIC: further predictions," Physics Letters, Section B, vol. 503, no. 1-2, pp. 5864, 2001.

[38] T. Hirano and M. Gyulassy, "Perfect fluidity of the quark-gluon plasma core as seen through its dissipative hadronic corona," Nuclear Physics A, vol. 769, no. 1-4, pp. 71-94, 2006.

[39] C. Nokana and S. A. Bass, "Space-time evolution of bulk QCD matter," Physical Review C, vol. 75, Article ID 014902, 2007.

[40] L. D. Landau, "On the multiple production of particles in high energy collisions," Izvestiya Akademii Nauk Azerbaidzhanskoi SSR. Seriya Fiziko-Tekhnicheskikh, vol. 17, p. 51,1953.

[41] J. D. Bjorken, "Highly relativistic nucleus-nucleus collisions: the central rapidity region," Physical Review D, vol. 27, no. 1, pp. 140151, 1983.

[42] T. Hirano, N. van der Kolk, and A. Bilandzic, "Hydrodynamics and flow," Lecture Notes in Physics, vol. 785, pp. 139-178, 2010.

[43] J.-Y. Ollitrault, "Relativistic hydrodynamics for heavy-ion collisions," European Journal of Physics, vol. 29, no. 2, pp. 275-302, 2008.

[44] S. Weinberg, "Entropy generation and the survival of protogalaxies in an expanding universe," The Astronomical Journal, vol. 168, pp. 175-194,1971.

[45] P. Danielewich and m. Gyulassy, "Dissipative phenomena in quark-gluon plasmas," Physical Review D, vol. 31, pp. 53-62, 1985.

[46] R. C. Hwa and K. Kajantie, "Diagnosing quark matter by measuring the total entropy and the photon or dilepton emission rates," Physical Review D, vol. 32, pp. 1109-1118,1985.

[47] D. Kharzeev and M. Nardi, "Hadron production in nuclear collisions at RHIC and high-density QCD," Physics Letters B, vol. 507, no. 1-4, pp. 121-128, 2001.

[48] S. Borsanyi, G. Endrodi, Z. Fodor, A. Jakovac et al., "The QCD equation of state with dynamical quarks," Journal ofHigh Energy Physics, vol. 1011, p. 077, 2010.

[49] F. Cooper and G. Frye, "Single-particle distribution in the hydrodynamic and statistical thermodynamic models of mul-tiparticle production," Physical Review D, vol. 10, no. 1, pp. 186189, 1974.

[50] R. Bair, H. Nakkagawa, A. Niegawa, and K. Redlich, "Production rate of hard thermal photons and screening of quark mass singularity," Zeitschrift fur Physik C Particles and Fields, vol. 53, no. 3, pp. 433-438, 1992.

[51] E. Braaten and R. D. Pisarski, "Soft amplitudes in hot gauge theories: a general analysis," Nuclear Physics B, vol. 337, no. 3, pp. 569-634,1990.

[52] E. Braaten and R. D. Pisarski, "Deducing hard thermal loops from Ward identities," Nuclear Physics B, vol. 339, no. 2, pp. 310324,1990.

[53] P. Aurenche, F. Gelis, R. Kobes, and H. Zaraket, "Bremsstrahlung and photon production in thermal QCD," Physical Review D, vol. 58, Article ID 085003,1998.

[54] P. Arnold, G. D. Moore, and L. G. Yaffe, "Photon emission from ultrarelativistic plasmas," Journal of High Energy Physics, vol. 0111, Article ID 057, 2001.

[55] P. Arnold, G. D. Moore, and L. G. Yaffe, "Photon emission from Quark-Gluon plasma: complete leading order results," Journal ofHigh Energy Physics, vol. 0112, Article ID 009, 2001.

[56] O. Kaczmarek and F. Zantow, "Static quark-antiquark interactions in zero and finite temperature QCD: I. Heavy quark free energies, running coupling, and quarkonium binding," Physical Review D, vol. 71, Article ID 114510, 2005.

[57] P. Roy, S. Sarkar, J.-E. Alam, and B. Sinha, "Electromagnetic radiation from hot and dense hadronic matter," Nuclear Physics A, vol. 653, no. 3, pp. 277-300,1999.

[58] J.-E. Alam, P. Roy, and S. Sarkar, "Comment on "Hadronic production of thermal photons"," Physical Review C, vol. 71, Article ID 059802, 2005.

[59] S. Mitra, P. Mohanty, S. Ghosh, S. Sarkar, and J.-E. Alam, "Characterizing quark gluon plasma by thermal photons and lepton pairs," http://arxiv.org/abs/1303.0675.

[60] J. J. Sakurai, Currents and Mesons, The University of Chicago Press, Chicago, Ill, USA, 1969.

[61] R. K. Bhaduri, Models of Nucleon—From Quarks to Soliton, Addison-Wesley, 1988.

[62] S. Turbide, R. Rapp, and C. Gale, "Hadronic production of thermal photons," Physical Review C, vol. 69, no. 1, Article ID 014903, 2004.

[63] K. L. Haglin, "Photon production from hot hadronic matter," Journal of Physics G, vol. 30, p. L27, 2004.

[64] H. Von Gersdorff, L. D. McLerran, M. Kataja, and P. V. Ruuskanen, "Studies of the hydrodynamic evolution of matter produced in fluctuations in pp collisions and in ultrarelativistic nuclear collisions," Physical Review D,vol. 34, pp. 794-810,1986.

[65] M. M. Aggarwal, A. Agnihotri, Z. Ahammed, WA98 Collaboration et al., "Observation of direct photons in central 158A GeV 208P208b+Pb collisions," Physical Review Letters, vol. 85, pp. 3595-3599, 2000.

[66] C. Bernard, T. Burch, C. DeTar et al., "QCD equation of state with 2+1 flavors of improved staggered quarks," Physical Review D, vol. 75, Article ID 094505, 2007.

[67] I. Arsene, I. G. Bearden, D. Beavis, BRAHMS Collaboration et al., "Quark-gluon plasma and color glass condensate at RHIC? The perspective from the BRAHMSexperiment," Nuclear Physics A, vol. 757, no. 1-2, pp. 1-27, 2005.

[68] B. B. Back, M. D. Baker, M. Ballintijn, PHOBOS Collaboration et al., "The PHOBOS perspective on discoveries at RHIC," Nuclear Physics A, vol. 757, no. 1-2, pp. 28-101, 2005.

[69] T. Hirano and K. Tsuda, "Collective flow and two-pion correlations from a relativistic hydrodynamic model with early chemical freeze-out," Physical Review C, vol. 66, Article ID 054905, 2002.

[70] D. L. Adams, N. Achurin, N. I. Belikov, E704 Collaboration et al., "Measurement of single spin asymmetry for direct photon production in pp collisions at 200 GeV/c," Physics Letters B, vol. 345, no. 4, pp. 569-575,1995.

[98 [99

[71] T. Renk, "A comprehensive description of multiple observables [89 in heavy-ion collisions at SPS," Journal of Physics G, vol. 30, Article ID 064905, p. 1495, 2004.

[72] T. Renk, "Photonic measurements of the longitudinal expansion [90 dynamics in relativistic heavy-ion collisions," Physical Review C,

vol. 71, Article ID 064905, 2005.

[73] J.-E. Alam, S. Sarkar, T. Hatsuda, T. K. Nayak, and B. Sinha, [91 "Photons from Pb-Pb collisions at ultrarelativistic energies," Physical Review C, vol. 63, Article ID 021901, 2001.

[74] D. K. Srivastava and B. Sinha, "Radiation of single photons from [92 Pb+Pb collisions at relativistic energies and the quark-hadron phase transition," Physical Review C, vol. 64, Article ID 034902,

2001. [93

[75] P. Huovinen, P. V. Ruuskanen, and S. S. Räsanen, "Photon emission in heavy ion collisions at the CERN SPS," Physics [94 Letters, Section B, vol. 535, no. 1-4, pp. 109-116, 2002.

[76] K. Gallmeister, B. Kämpfer, and O. P. Pavelenko, "Unique large [95 thermal source of real and virtual photons in the reactions Pb(158A GeV) + Pb, Au," Physical Review C, vol. 62, Article ID 057901,2000. [96

[77] F. D. Steffen and M. H. Thoma, "Hard thermal photon production in relativistic heavy ion collisions," Physics Letters B, vol. [97 510, no. 1-4, pp. 98-106, 2001.

[78] R. Chatterjee, D. K. Srivastava, and S. Jeon, "Single photons from relativistic collisions of lead nuclei at energies available at the CERN Super Proton Synchrotron (SPS): a reanalysis," Physical Review C, vol. 79, Article ID 034906, 2009.

[79] A. Adare, S. Afanasiev, C. Aidala, PHENIX collaboration et al., "Enhanced production of direct photons in Au+Au collisions

at = 200 GeV and implications for the initial tempera- [100 ture," Physical Review Letters, vol. 104, Article ID 132301, 2010.

[80] S. S. Adler, S. Afanasiev, C. Aidala, PHENIX collaboration et al., "Suppressed n0 Production at Large Transverse Momentum in Central Au+Au Collisions at ^sNN = 200 GeV"' Physical Review [101 Letters, vol. 91, Article ID 072301, 2003.

[81] F. M. Liu, T. Hirano, K. Werner, and Y. Zhu, "Centrality-dependent direct photon pt spectra in Au+Au collisions at the

BNL Relativistic Heavy Ion Collider (RHIC) energy = [102

200 GeV" Physical Review C, vol. 79, Article ID 014905, 2009.

[82] S. Sarkar, P. Roy, J.-E. Alam, S. Raha, and B. Sinha, "Dissipative effects in photon diagnostics of quark—gluon plasma," Journal of Physics G, vol. 23, p. 469,1997.

[83] A. K. Chaudhuri, "Effect of viscosity on one dimensional hydrodynamic flow and direct photons from 200 AGeV S+Au collisions at CERN SPS," Physica Scripta, vol. 61, p. 311, 2000.

[84] A. K. Chaudhuri and B. Sinha, "Direct photon production from viscous quark-gluon plasma," Physical Review C, vol. 83, Article ID 034905, 2010.

[85] J. R. Bhatt, H. Mishra, and V Sreekanth, "Thermal photons in [106 QGP and non-ideal effects," Journal of High Energy Physics, vol.

1011, p. 106, 2010.

[86] K. Dusling, "Photons as a viscometer of heavy ion collisions," Nuclear Physics A, vol. 839, no. 1-4, pp. 70-77, 2010. [107

[87] S. Mitra, P. Mohanty, S. Sarkar, and J.-E. Alam, "Thermal radiation from an evolving viscous Quark Gluon plasma," [108 International Journal of Modern Physics E, vol. 22, Article ID 1350004, 2013.

[88] D. Teaney, "Effect of shear viscosity on spectra, elliptic flow, and [109 Hanbury Brown-Twiss radii," Physical Review C, vol. 68, Article

ID 034913, 2003.

[104 [105

P. Danielewicz and M. Gyulassy, "Dissipative phenomena in quark-gluon plasmas," Physical Review D, vol. 31, no. 1, pp. 5362, 1985.

H. B. Meyer, "Calculation of the shear viscosity in SU(3) gluodynamics," Physical Review D, vol. 76, Article ID 101701, 2007.

P. Arnold, G. D. Moore, and L. G. Yaffe, "Transport coefficients in high temperature gauge theories: (II) Beyond leading log," Journal ofHigh Energy Physics, vol. 0305, Article ID 051, 2003. J.-W. Chen, H. Dong, K. Ohnishi, and Q. Wang, "Shearviscosity of a gluon plasma in perturbative QCD," Physics Letters, Section

B, vol. 685, no. 4-5, pp. 277-282, 2010.

S. K. Das and J. Alam, "Transport coefficients of gluonic fluid,"

Physical Review D, vol. 83, Article ID 114011, 2011.

Z. Xu and C. Greiner, "Shear viscosity in a Gluon gas," Physical

Review Letters, vol. 100, Article ID 172301, 2008.

K. Itakura, O. Morimatsu, and H. Otomo, "Shear viscosity of a

hadronic gas mixture," Physical Review D, vol. 77, no. 1, Article

ID 014014, 2008.

A. Dobado and I. J. Llanes-Estrada, "Viscosity of meson matter," Physical Review D, vol. 69, Article ID 116004, 2004. P. Chakraborty and J. I. Kapusta, "Quasiparticle theory of shear and bulk viscosities of hadronic matter," Physical Review C, vol. 83, no. 1, Article ID 014906, 2011.

A. S. Khvorostukhin, V. D. Toneev, and D. N. Voskresensky, "Viscosity coefficients for hadron and quark-gluon phases," Nuclear Physics A, vol. 845, no. 1-4, pp. 106-146, 2010.

I. Tserruya, "Summary of low-mass dilepton and direct photon results," Nuclear Physics A, vol. 590, no. 1-2, pp. 127-138,1995. J. Cleymans, J. Fingberg, and K. Redlich, "Transverse-momentum distribution of dileptons in different scenarios for the QCD phase transition," Physical Review D, vol. 35, pp. 2153-2165, 1987.

K. Dusling and I. Zahed, "Transverse momentum spectra of dileptons measured by the NA60 Collaboration in In+In collisions at 158 GeV/nucleon," Physical Review C, vol. 80, no. 1, Article ID 014902, 2009.

T. Altherr and P. V. Ruuskanen, "Low-mass dileptons at high momenta in ultra-relativistic heavy-ion collisions," Nuclear Physics B, vol. 380, no. 3, pp. 377-390,1992. M. H. Thoma and C. T. Traxler, "Production of energetic dileptons with small invariant masses from the quark-gluon plasma," Physical Review D, vol. 56, no. 1, pp. 198-202,1997. E. V. Shuryak, "Correlation functions in the QCD vacuum," Reviews of Modern Physics, vol. 65, pp. 1-46,1993. R. Barate, D. Buskulic, D. Decamp, ALEPH Collaboration et al., "Measurement of the axial-vector t spectral functions and determination of a3(M2) from hadronic t decays," The European Physical Journal C, vol. 4, no. 3, pp. 409-431,1998. K. Ackerstaff, G. Alexander, J. Allison et al., "Measurement of the strong coupling constant and the vector a3 and axial-vector spectral functions in hadronic tau decays," The European Physical Journal C, vol. 7, pp. 571-593,1999. S. Sarkar, "Hadrons at finite temperature," Nuclear Physics A, vol. 862-863, pp. 13-19, 2011.

S. Mallik and S. Sarkar, "Vector and axial-vector mesons at finite temperature," The European Physical Journal C, vol. 25, pp. 445452, 2002.

S. Ghosh, S. Mallik, and S. Sarkar, "Analytic structure of p meson propagator at finite temperature," The European Physical Journal

C, vol. 70, pp. 251-262, 2010.

40 [110

[114 [115 [116

[118 [119 [120

[121 [122

[124 [125 [126 [127

S. Ghosh and S. Sarkar, " p self-energy at finite temperature and [128 density in the real-time formalism," Nuclear Physics A, vol. 870871, no. 1, pp. 94-111, 2011.

S. Ghosh and S. Sarkar, "Analysis of self-energy at finite temperature and density in the real-time formalism," The European [129 Physical Journal A, vol. 49, article 97, 2013. S. Sarkar and S. Ghosh, "In-medium vector mesons and low mass lepton pairs from heavy ion collisions," Journal ofPhysics, [130 vol. 374, Article ID 012010, 2012.

J. K. Nayak, J.-E. Alam, T. Hirano, S. Sarkar, and B. Sinha, "Muon pairs from In + In collision at energies available at the CERN [131 Super Proton Synchrotron," Physical Review C, vol. 85, Article ID 064906, 2012.

R. Arnaldi, R. Averbeck, K. Banicz et al., "First measurement [132 of the p spectral function in high-energy nuclear collisions," Physical Review Letters, vol. 96, Article ID 162302, 2006. R. Arnaldi, R. Averbeck, K. Banicz et al., "Evidence for radial flow of thermal dileptons in high-energy nuclear collisions," Physical Review Letters, vol. 100, Article ID 022302, 2008.

V. L. Eletsky, M. Belkacem, P. J. Ellis, and J. I. Kapusta, "Properties of p and œ mesons at finite temperature and density as inferred from experiment," Physical Review C, vol. 64, Article ID 035202, 2001. H. Van Hees and R. Rapp, "Comprehensive interpretation of thermal dileptons measured at the CERN super proton [136 synchrotron," Physical Review Letters, vol. 97, no. 10, Article ID 102301, 2006.

R. Rapp, G. Chanfray, and J. Wambach, "Rho Meson propagation and dilepton enhancement in hot hadronic matter," Nuclear Physics A, vol. 617, pp. 472-495,1997.

A. Toia and PHENIX Collaboration, "e+e- pairs: a clock and a thermometer of heavy ion collisions," Journal ofPhysics G, vol. 35, Article ID 104037, 2008.

A. Adare, S. Afanasiev, C. Aidala, PHENIX Collaboration et al., "Detailed measurement of the e+e- pair continuum in p+p and Au+Au collisions at ^sNN = 200 GeV and implications for direct photon production," Physical Review C, vol. 81, Article ID 034911, 2010.

P. Lichard and J. Jura, "Electron-positron annihilation into four charged pions and the a. pn Lagrangian," Physical Review D, vol. [1 76, Article ID 094030, 2007. [

K. Aamodt, N. Abel, U. Abeysekara, ALICE collaboration et al., "Charged-particle multiplicity measurement in proton-proton collisions at ^s = 0.9 and 2.36 TeV with ALICE at LHC," The European Physical Journal C, vol. 68, pp. 89-108, 2010. F. -M. Liu, T. Hirano, K. Werner, and Y. Zhu, "Centrality- [143 dependent direct photon pt spectra in Au+Au collisions at the BNL Relativistic Heavy Ion Collider (RHIC) energy ^sNN = 200 GeV" Physical Review C, vol. 79, Article ID 014905, 2009. [144

B. Schenke, S. Jeon, and C. Gale, "(3+1)D hydrodynamic simulation of relativistic heavy-ion collisions," Physical Review

C. vol. 82, Article ID 014903, 2010.

B. Schenke, S. Jeon, and C. Gale, "Higher flow harmonics from (3+1)D event-by-event viscous hydrodynamics," Physical Review C, vol. 85, Article ID 024901, 2012. J. K. Nayak and J. Alam, "Measuring radial flow ofpartonic and hadronic phases in relativistic heavy-ion collisions," Physical Review C, vol. 80, Article ID 064906, 2009.

B. Sinha, "Universal signals of a quark-gluon plasma," Physics [147 Letters B, vol. 128, no. 1-2, pp. 91-94,1983.

[134 [135

[137 [138 [139 [140

[145 [146

R. Chatterjee, E. S. Frodermann, U. Heinz, and D. K. Srivastava, "Elliptic flow of thermal photons in relativistic nuclear collisions," Physical Review Letters, vol. 96, no. 20, Article ID 202302, 2006.

R. Hanbury Brown and R. Q. Twiss, "A test of a new type of stellar interferometer on sirius," Nature, vol. 178, pp. 1046-1048, 1956.

D. Peressounko, "Hanbury Brown-Twiss interferometry of direct photons in heavy ion collisions," Physical Review C, vol. 67, Article ID 014905, 2003.

Y. M. Sinyukov, "Hanbury Brown/Twiss correlations for expanding hadron and quark-gluon matter," Nuclear Physics A, vol. 498, pp. 151-159, 1989.

A. Adare, S. Afanasiev, C. Aidala et al., "Energy loss and flow of heavy Quarks in Au+Au collisions at ^sNN = 200 GeV," Physical Review Letters, vol. 98, Article ID 172301, 2007.

B. Abelev, B. Abelev, A. A. Quintana, The ALICE Collaboration et al., "Charged-particle multiplicity density at Midrapidity in central Pb-Pb collisions at ^sNN = 2.76 TeV,"' Physical Review Letters, vol. 105, Article ID 252301, 2010.

S. Pratt, "Pion interferometry for exploding sources," Physical Review Letters, vol. 53, no. 13, pp. 1219-1221, 1984. S. Pratt, "Pion interferometry of quark-gluon plasma," Physical Review D, vol. 33, pp. 1314-1327,1986.

U. A. Widemann and U. Heinz, "Particle interferometry for relativistic heavy-ion collisions," Physics Reports, vol. 319, no. 45, pp. 145-230,1999.

U. Heinz and B. V. Jacak, "Two-particle correlations in relativistic heavy-ion collisions," Annual Review of Nuclear and Particle Science, vol. 49, no. 1, pp. 529-579,1999.

T. Csörgö and B. Lörstad, "Bose-Einstein correlations for three-dimensionally expanding, cylindrically symmetric, finite systems," Physical Review C, vol. 54, pp. 1390-1403,1996. D. H. Rischke and M. Gyulassy, "The time-delay signature of quark-gluon plasma formation in relativistic nuclear collisions," Nuclear Physics A, vol. 608, no. 4, pp. 479-512,1996. S. Chapman, P. Scotto, and U. Heinz, "New cross term in the two-particle Hanbury-Brown-Twiss correlation function in ultrarelativistic heavy-ion collisions," Physical Review Letters, vol. 74, no. 22, pp. 4400-4403, 1995.

M. Herrmann and G. F. Bertsch, "Source dimensions in ultrarelativistic heavy-ion collisions," Physical Review C, vol. 51, pp. 328-338, 1995.

P. Huovinen and P. V. Ruuskanen, "Hydrodynamic models for heavy ion collisions," Annual Review of Nuclear and Particle Science, vol. 56, pp. 163-206, 2006.

R. Chatterjee, D. K. Srivastava, U. W. Heinz, and C. Gale, "Elliptic flow of thermal dileptons in relativistic nuclear collisions," Physical Review C, vol. 75, Article ID 054909, 2007. R. Chatterjee, D. K. Srivastava, U. W. Heinz, and C. Gale, "Elliptic flow of thermal dileptons in relativistic nuclear collisions," Physical Review C, vol. 75, Article ID 054909, 2007. J. Deng, Q. Wang, N. Xu, and P. Zhuang, "Dilepton flow and deconfinement phase transition in heavy ion collisions," Physics Letters B, vol. 701, pp. 581-586, 2011.

H. Van Hees, C. Gale, and R. Rapp, "Thermal photons and collective flow at energies available at the BNL relativistic heavy-ion collider," Physical Review C, vol. 84, no. 5, Article ID054906, 2011.

T. Renk and J. Ruppert, "Dimuon transverse momentum spectra as a tool to characterize the emission region in heavy-ion

collisions," Physical Review C, vol. 77, no. 2, Article ID 024907, 2008.

[148] V. Roy and A. K. Choudhuri, "Charged particle's elliptic flow in (2+12+1)D viscous hydrodynamics at LHC (^s = 2.76 TeV) energy in Pb + Pb collision," Physics Letters B, vol. 703, no. 3, pp. 313-317, 2011.

[149] B. Mohanty and J. Alam, "Velocity of sound in relativistic heavy-ion collisions," Physical Review C, Article ID 064903, 2003.

[150] K. Aamodt, A. Abrahantes Quintana, D. Adamova, ALICE Collaboration et al., "Centrality dependence of the charged-particle multiplicity density at midrapidity in Pb-Pb collisions at ^sNN = 2.76TeV' Physical Review Letters, vol. 106, Article ID 032301, 2011.

[151] K. Aamodt, A. Abrahantes Quintana, D. Adamova', ALICE Collaboration et al., "Suppression of charged particle production at large transverse momentum in central Pb-Pb collisions at ^s^ = 2.76 TeV"' Physics Letters B, vol. 696, p. 30, 2011.

[152] K. Itakura, "Strong-field physics as a new probe of early-time dynamics in heavy-ion collisions," Journal of Physics, vol. 422, Article ID 012029, 2013.

[153] K. Tuchin, "Photon decay in a strong magnetic field in heavy-ioncollisions," Physical Review C, vol. 83, Article ID 017901,2011.

[154] K. Tuchin, "Electromagnetic radiation by quark-gluon plasma in a magnetic field," Physical Review C, vol. 87, Article ID 024912, 2013.

Copyright of Advances in High Energy Physics is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.