Abstract

Low-mass dilepton spectra as measured in high-energy heavy-ion collisions are a unique tool to obtain spectroscopic information about the strongly interacting medium produced in these reactions. Specifically, in-medium modifications of the vector spectral function, which is well known in the vacuum, can be deduced from the thermal radiation off the expanding QCD fireball. This, in particular, allows to investigate the fate of the resonance in the dense medium and possibly infer from it signatures of the (partial) restoration of chiral symmetry, which is spontaneously broken in the QCD vacuum. After briefly reviewing calculations of thermal dilepton emission rates from hot QCD matter, utilizing effective hadronic theory, lattice QCD, or resummed perturbative QCD, we focus on applications to dilepton spectra at heavy-ion collider experiments at RHIC and LHC. This includes invariant-mass spectra at full RHIC energy with transverse-momentum dependencies and azimuthal asymmetries, as well as a systematic investigation of the excitation function down to fixed-target energies, thus making contact to previous precision measurements at the SPS. Furthermore, predictions for the energy frontier at the LHC are presented in both dielectron and dimuon channels.

1. Introduction

The exploration of matter at extremes of temperature () and baryon density () is at the forefront of research in contemporary nuclear physics, with intimate connections to high-energy, condensed-matter, and even atomic physics [1]. Theoretical efforts over the last few decades are suggesting an extraordinary richness of the phase diagram of strongly interacting matter, which should ultimately emerge from the underlying theory of quantum chromodynamics (QCD) as part of the standard model. However, several basic questions, both qualitative and quantitative, such as the possible existence of first order transitions and their location as function of baryon-chemical potential () and temperature, remain open to date [2]. A close interplay of experiment and theory is needed to create a robust knowledge about the QCD phase structure. On one hand, naturally occurring matter at temperatures close to or beyond the expected pseudo-critical one,  MeV [3, 4], may last have existed ~14 billion years ago, during the first tens of microseconds of the Universe. On the other hand, at small temperatures, matter with baryon densities close to or beyond the critical one for the transition into quark matter may prevail in the interior of compact stars today, but its verification and exploration from observational data are challenging [5]. It is quite fascinating that tiny man-made samples of hot QCD matter can nowadays be created and studied in the laboratory using ultrarelativistic heavy-ion collisions (URHICs). Significant progress has been made in understanding the properties of this medium through analyses of experiments conducted at the CERN's Super-Proton Synchrotron (SPS), BNL's Relativistic Heavy-Ion Collider (RHIC), and CERN's Large Hadron Collider (LHC) (see, e.g., the recent Quark Matter conference proceedings [6, 7]). For example, systematic investigations of the produced hadron spectra have revealed a hydrodynamic behavior of the bulk matter in the region of low transverse momenta ( GeV) and a strong absorption of hadrons with high transverse momentum (). Even hadrons containing a heavy quark (charm or bottom) exhibit substantial energy loss and collectivity due to their coupling to the expanding fireball. While the total charm and bottom yields are essentially conserved, the production of heavy quark-antiquark bound states (charmonia and bottomonia) is largely suppressed. The relation of the above hadronic observables to spectral properties of the medium is, however, rather indirect. Low-mass dileptons, on the other hand, are radiated from the interior of the medium throughout the fireball's lifetime, as their mean-free path is much larger than the size of the fireball. Thus, their invariant-mass spectra directly measure the in-medium vector spectral function, albeit in a superposition of the varying temperature in the fireball's expansion.

The dilepton program at the SPS has produced remarkable results. The CERES/NA45 dielectron data in Pb-Au collisions, and particularly the NA60 dimuon spectra in In-In collisions, have shown that the -meson undergoes a strong broadening, even complete melting, of its resonance structure, with quantitative sensitivity to its spectral shape, see [810] for recent reviews. The QCD medium at SPS energies is characterized by a significant net-baryon content with chemical potentials of  MeV at chemical freezeout,  MeV [11], and further increasing as the system cools down [12]. Baryons have been identified as a dominant contributor to the medium modifications of the ’s spectral function [10]. The question arises how these develop when moving toward the net baryon-free regime in the QCD phase diagram, . Theoretical expectations based on the hadronic many-body approach [13] suggest comparable medium effects in this regime, since the relevant quantity is the sum of baryon and antibaryon densities, and this turns out to be similar at SPS and RHIC/LHC [12], at least close to . Since at collider energies, the total baryon density at RHIC and LHC in the subsequent hadronic evolution of the fireball will remain similar. We also note that the  MeV regime is amenable to numerical lattice QCD calculations, both for the equation of state of the medium evolution, and in particular for the microscopic dilepton production rate, at least in the QGP phase for now [14, 15]. Furthermore, since the phase transition at  MeV presumably is a continuous crossover [16], a realistic dilepton rate should vary smoothly when changing the temperature through . Thus, after the successful fixed-target dilepton program at the CERN-SPS, the efforts and attention are now shifting to collider energies around experiments at RHIC and LHC.

In the present paper we will focus on the theory and phenomenology of dilepton production at collider energies (for a recent overview including an assessment of SPS data, see, e.g., [17]). The presented material is partly of review nature, but also contains thus far unpublished results, for example, updates in the use of nonperturbative QGP dilepton rates and equation of state, and detailed predictions for invariant-mass and transverse-momentum spectra for ongoing and upcoming experiments at RHIC and LHC, including an excitation function of the beam energy scan program at RHIC.

This paper is organized as follows. In Section 2, we briefly review the calculation of the thermal dilepton emission rates from hadronic matter and the quark-gluon plasma (QGP). We elaborate on how recent lattice-QCD results at vanishing three-momentum () may be extended to finite to enable their application to URHICs. In Section 3, we discuss in some detail the calculations of dilepton spectra suitable for comparison with experiment; this involves a brief discussion of the medium evolution in URHICs (including an update of the equation of state) in Section 3.1 and of nonthermal sources (primordial production and final-state decays) in Section 3.2. It will be followed by analyses of mass and momentum spectra, as well as elliptic flow at full RHIC energy in Section 3.3, and of an excitation function as obtained from the RHIC beam energy scan in Section 3.4; predictions for dielectron and dimuon spectra at current (2.76 ATeV) and future (5.5 ATeV) LHC energies are presented in Section 3.5. We end with a summary and outlook in Section 4.

2. Thermal Dilepton Rates in QCD Matter

The basic quantity for connecting calculations of the electromagnetic (EM) spectral function in QCD matter to measurements of dileptons in heavy-ion collisions is their thermal emission rate; per unit phase space, it can be written as where is a lepton phase-space factor (=1 for vanishing lepton mass), denotes the thermal Bose distribution, and is the energy of the lepton pair (or virtual photon) in terms of its invariant mass and 3-momentum. As mentioned above, this observable is unique in its direct access to an in-medium spectral function of the formed system, namely, in the vector (or EM) channel, . It is defined via the correlation function of the EM current, , as transported by the electric-charge carriers in the system. In quark basis, the EM current is given by the charge-weighted sum over flavor: while in hadronic basis, it is in good approximation given by the vector-meson fields: known as vector-dominance model (VDM). Since the significance of thermal dilepton radiation is limited to masses below the mass,  GeV, we will focus on the light- and strange-quark sector in this article.

In the vacuum, the EM spectral function is well known from the annihilation cross section into hadrons, usually quoted relative to the annihilation into dimuons as the ratio (cf. Figure 1). It illustrates that the nonperturbative hadronic description in terms of VDM works well in the low-mass region (LMR),  GeV, while the perturbative partonic description appears to apply for  GeV. Thus, in URHICs, dilepton spectra in the LMR are ideally suited to study the properties of vector mesons in the medium. A central question is if and how these medium modifications can signal (the approach to) deconfinement and the restoration of the dynamical breaking of chiral symmetry (DBCS). After all, confinement and DBCS govern the properties of hadrons in vacuum. At masses  GeV, the perturbative nature of the EM spectral function suggests that in-medium modifications are suppressed, coming in as corrections in powers of and . In this case, invariant-mass spectra of thermal radiation become an excellent measure for the prevalent temperatures of the produced system, free from blue shifts due to the medium expansion which strongly affect spectra.

2.1. Hadronic Matter

Over the last two decades, broad efforts have been undertaken to evaluate the medium modifications of the -meson. The latter dominates in the EM spectral function over the by about a factor of 10 (the appears to be rather protected from hadronic medium effects, presumably due to the OZI rule, at least for its coupling to baryons). Recent overviews of these efforts can be found, for example, in [10, 18, 19]. Most approaches utilize effective hadronic (chiral) Lagrangians and apply them in diagrammatic many-body theory to compute thermal (or density) loop corrections. The generic outcome is that of a substantial broadening of the ’s spectral shape, with little mass shift (in a heat bath, chiral symmetry protects the from mass shifts at order [20]). The magnitude of the ’s in-medium width (and/or its precise spectral shape) varies in different calculations, but the discrepancies can be mostly traced back to the differing contributions accounted for in the Lagrangian (e.g., the set of baryon and/or meson resonance excitations, or medium effects in the ’s pion cloud). Similar findings arise when utilizing empirically extracted on-shell -meson scattering amplitudes off hadrons in linear-density approximation [21]. Since these calculations are restricted to resonances above the nominal (or ) threshold, quantitative differences to many-body (field-theoretic) approaches may arise; in particular, the latter account for subthreshold excitations, for example, , which induce additional broadening and associated enhancement of the low-mass part in the spectral function (also causing marked deviations from a Breit-Wigner shape). Appreciable mass shifts are typically found in mean-field approximations (due to large in-medium scalar fields) or in calculations where the bare parameters of the underlying Lagrangian are allowed to be temperature dependent [22].

An example for dilepton rates following from a spectral function calculated in hot and dense hadronic matter at SPS energies is shown in Figure 2(a). The EM spectral function follows from the -meson using VDM, (3), although corrections to VDM are necessary for quantitative descriptions of the EM couplings in the baryon sector [23, 24]. When extrapolated to temperatures around , the resonance peak has essentially vanished leading to a structureless emission rate with a large enhancement in the mass region below the free mass. The decomposition of the rate into in-medium self-energy contributions illustrates the important role of the pion cloud modifications and of multiple low-energy excitations below the free mass, for example, resonance-hole , that is, for off-shell -mesons. The hadronic medium effects are slightly reduced at collider energies (Figure 2(b)), where a faint resonance structure appears to survive at around (it is significantly more suppressed at  MeV). A recent calculation in a similar framework, combing thermal field theory with effective hadron Lagrangians [25] and including both finite-temperature and -density contributions to the self-energy through baryon and meson resonances, shows fair agreement with the results shown in Figure 2(a).

2.2. Quark-Gluon Plasma

In a perturbative QGP (pQGP), the leading-order (LO) mechanism of dilepton production is EM quark-antiquark annihilation as following from a free quark current in (2). The corresponding EM spectral function is essentially given by the “naive quark model” curve in Figure 1, extended all the way down to vanishing mass, where (an additional phase-space factor occurs for finite current quark masses) and . Finite-temperature corrections are induced by a quantum-statistical Pauli-blocking factor (written for ) which produces a nontrivial 3-momentum dependence [26]; for , it simplifies to , where is the thermal Fermi distribution. The pertinent 3-momentum integrated dilepton rate is structureless (cf. long-dashed curve in Figure 2(b)). It's finite value at implies that no real photons can be produced from this mechanism.

A consistent implementation of corrections in a thermal QGP at vanishing quark chemical potential has been achieved by resumming the hard-thermal-loop (HTL) action [27]. Quarks and gluons acquire thermal masses , but bremsstrahlung-type contributions lead to a marked enhancement of the rate over the LO pQCD results (cf. the dash-dotted line in Figure 2(b)).

Recent progress in calculating dilepton rates nonperturbatively using thermal lattice QCD (lQCD) has been reported in [14, 15, 28]. The basic quantity computed in these simulations is the Euclidean-time correlation function which is related to the spectral function, , via Results for obtained in quenched QCD for at vanishing (in which case ) are shown by the data points in Figure 3(a), normalized to the free (noninteracting) pQGP limit. At small , corresponding to large energies in the spectral function, this ratio tends to one as expected for the perturbative limit. For larger , a significant enhancement develops which is associated with a corresponding enhancement in the low-energy (or low-mass) regime of the spectral function (and thus dilepton rate). This enhancement may be quantified by making an ansatz for the spectral function in terms of a low-energy Breit-Wigner part plus a perturbative continuum [14], (note that ). The strength () and width () of the Breit-Wigner, as well as a perturbative correction (), are then fit to the Euclidean correlator. The large- enhancement in the correlator generates an appreciable low-energy enhancement in the spectral function (cf. Figure 3(b)). The zero-energy limit of the spectral function defines a transport coefficient, the electric conductivity, . Similar to the viscosity or heavy-quark diffusion coefficient, a small value for , implied by a large value for , indicates a strong coupling of the medium; for example, in pQCD, [29]. The results for the dilepton rate (or spectral function) at a smaller temperature of are found to be similar to the ones at [28], suggesting a weak temperature dependence in this regime. Note, however, that the phase transition in quenched QCD is of first order; that is, a stronger variation is expected when going across . Recent results for two-flavor QCD [15] also indicate rather structureless spectral functions similar to the quenched results. Ultimately, at sufficiently small temperatures, the lattice computations should recover a -meson resonance peak; it will be interesting to see at which temperatures this occurs.

For practical applications, a finite 3-momentum dependence of the lQCD dilepton rate is needed, which is currently not available from the simulations. We here propose a “minimal” construction which is based on a matching to the 3-momentum dependence obtained from the LO pQCD photon rate [30]. The latter reads The idea is now to adopt the transverse part of the EM spectral function as given by (7) for the 3-momentum dependence of the spectral function in (6) by replacing the Breit-Wigner part with it; that is, Here, we have introduced a factor into , which serves two purposes: (i) with , it rather accurately accounts for the enhancement of the complete LO photon rate calculation [31] over the rate in (7); (ii) it better reproduces the low-energy regime of the lQCD spectral function; for example, for , the electric conductivity following from (8) is , not far from the lQCD estimate with the fit ansatz (6), (also compatible with [32]; the systematic uncertainty in the lattice result, due to variations in the ansatz, is significantly larger). The resulting spectral function (upper dashed line in Figure 3(b)) somewhat overestimates the lQCD result at high energies, where the latter coincides with the annihilation term. This can be improved by an additional form factor, , resulting in the lower dashed line in Figure 3(b) (using ).

Finally, care has to be taken to include a finite longitudinal part which develops in the timelike regime. Here, we employ a dependence that follows, for example, from standard constructions of gauge-invariant -wave -baryon interactions, yielding [34]. Thus, we finally have

The lQCD results for the isovector spectral function are compared to hadronic calculations in Figure 3(b). Close to the phase transition temperature, the “melting” of the in-medium spectral function suggests a smooth transition from its prominent resonance peak in vacuum to the rather structureless shape extracted from lQCD, signaling a transition from hadronic to partonic degrees of freedom. It would clearly be of interest to extract the conductivity from the hadronic calculations, which currently is not well resolved from the , limit of the spectral function. The mutual approach of the nonperturbative hadronic and lQCD spectral functions is also exhibited in the 3-momentum integrated dilepton rate shown in Figure 2(b), especially when compared to the different shapes of the LO pQCD and vacuum hadronic rates. Arguably, the in-medium hadronic rate still shows an indication of a broad resonance. A smooth matching of the rates from above and below might therefore require some additional medium effects in the hot and dense hadronic medium and/or the emergence of resonance correlations in the correlator in the QGP. Unless otherwise noted, the thermal emission rates used in the calculations of dilepton spectra discussed below will be based on the in-medium hadronic rates of [35] and the lQCD-inspired QGP rates [14], extended to finite 3-momentum as constructed above (with and form factor).

3. Dilepton Spectra at RHIC and LHC

The calculation of dilepton mass and transverse-momentum () spectra, suitable for comparison to data in heavy-ion collisions, requires an integration of the thermal rates of hadronic matter and QGP over a realistic space-time evolution of the AA reaction: In addition to the thermal yield, nonthermal sources have to be considered, for example, primordial Drell-Yan annihilation and electromagnetic final-state decays of long-lived hadrons. We will briefly discuss space-time evolutions in Section 3.1 and nonthermal sources in Section 3.2, before proceeding to a more detailed discussion of thermal spectra and comparisons to data, as available, in Sections 3.3, 3.4, and 3.5 for full RHIC energy, the beam-energy scan, and LHC, respectively.

3.1. Medium Expansion

The natural framework to carry out the space-time integral over the dilepton rate in URHICs is relativistic hydrodynamics. The application of this approach to AA collisions at RHIC and LHC works well to describe bulk hadron observables (e.g., spectra and elliptic flow) up to momenta of  GeV, which typically comprises more than 90% of the total yields. Some uncertainties remain, for example, as to the precise initial conditions at thermalization, viscous corrections, or the treatment of the late stages where the medium becomes dilute and the hadrons decouple (see, e.g., [38] for a recent review). Another key ingredient is the equation of state (EoS) of the medium, , which drives its collective expansion. Figure 4(a) illustrates the effects of updating a previously employed bag-model EoS (a quasiparticle QGP connected to a hadron resonance gas via a first-order phase transition) [39] by a recent parametrization of a nonperturbative QGP EoS from lQCD data [40, 41] (continuously matched to a hadron-resonance gas at  MeV) [42]; within a 2+1-D ideal hydro calculation, the most notable change is a significant increase of the temperature (at fixed entropy density) in the regime just above the transition temperature (up to ca. 30 MeV at the formerly onset of the first-order transition). Together with the fact that the hadronic portion of the formerly mixed phase is now entirely associated with the QGP, this will lead to an increase (decrease) of the QGP (hadronic) contribution to EM radiation relative to the first-order scenario. In addition, the harder lattice EoS induces a stronger expansion leading to a slightly faster cooling and thus reduction in the lifetime by about 5%. This effect becomes more pronounced when modifying the initial conditions of the hydrodynamic evolution, for example, by introducing a more compact spatial profile (creating larger gradients) and/or initial transverse flow (associated with interactions prior to the thermalization time, ) [42] (cf. the solid line in Figure 4(a)). The resulting more violent expansion plays an important role in understanding the HBT radii of the system [43]. The relevance for EM radiation pertains to reducing the fireball lifetime by up to ~20%.

More simplistic parametrizations of the space-time evolution of AA collisions have been attempted with longitudinally and transversely expanding fireballs. With an appropriate choice of the transverse acceleration (in all applications below it is taken as at the surface), an approximate reproduction of the basic features (timescales and radial flow) of hydrodynamic evolutions can be achieved, see Figure 4(b). Most of the results shown in the remainder of this article are based on such simplified fireball parametrizations, utilizing the EoS of [42] where a parametrization of lQCD results is matched with a hadron resonance gas at  MeV and subsequent chemical freezeout at  MeV (see also [44]). We note that the use of this EoS, together with the lQCD-based QGP emission rates, constitutes an update of our earlier calculations [45] where a quasiparticle bag-model EoS was employed in connection with HTL rates in the QGP. We have checked that the previous level of agreement with the acceptance-corrected NA60 spectra is maintained, which is essentially due to the duality of the QGP and hadronic rates around (a more detailed account in the context of the SPS dilepton data will be given elsewhere [46]). For our discussion of collider energies below, the initialization (or thermalization times) are chosen at at full RHIC energy (increasing smoothly to at  GeV) and in the LHC regime. This results in initial temperatures of 225 MeV and 330 MeV in minimum-bias (MB) Au-Au collisions at 20 and 200 GeV, respectively, increasing to ~380 MeV in central Au-Au(200 GeV) and ~560(620) MeV in central Pb-Pb at 2.76(5.5) ATeV. These values differ slightly from previous calculations with a quasiparticle EoS; they are also sensitive to the initial spatial profile (cf. Figure 4(a)). However, for our main objective of calculating low-mass dilepton spectra, the initial temperature has little impact.

3.2. Nonthermal Sources

In addition to thermal radiation from the locally equilibrated medium, dilepton emission in URHICs can arise from interactions prior to thermalization (e.g., Drell-Yan annihilation) and from EM decays of long-lived hadrons after the fireball has decoupled (e.g., Dalitz decays or ). Furthermore, paralleling the structure in hadronic spectra, a nonthermal component from hard production will feed into dilepton spectra, for example, via bremsstrahlung from hard partons traversing the medium [47] or decays of both short- and long-lived hadrons which have not thermalized with the bulk (e.g., “hard” -mesons or long-lived EM final-state decays). Hadronic final-state decays (including the double semileptonic decay of two heavy-flavor hadrons originating from a or pair produced together in the same hard process) are commonly referred to as the “cocktail,” which is routinely evaluated by the experimental collaborations using the vacuum properties of each hadron with spectra based on measured spectra or appropriately extrapolated using thermal blast-wave models. In URHICs, the notion of the cocktail becomes problematic for short-lived resonances whose lifetime is comparable to the duration of the freezeout process of the fireball (e.g., for , , etc.). In their case, a better approximation is presumably to run the fireball an additional ~ to treat their final-decay contribution as thermal radiation including medium effects. However, care has to be taken in evaluating their dilepton -spectra, as the latter are slightly different for thermal radiation and final-state decays (cf. [45] for a discussion and implementation of this point). For light hadrons at low , the cocktail scales with the total number of charged particles, , at given collision energy and centrality, while for hard processes, a collision-number scaling is in order (and compatible with experiment where measured, modulo the effects of “jet quenching”). The notion of “excess dileptons” is defined as any additional radiation observed over the cocktail, sometimes quantified as an “enhancement factor” in a certain invariant-mass range. The excess radiation is then most naturally associated with thermal radiation, given the usual limitation where hard processes take over, that is,  GeV.

3.3. RHIC-200

We start our discussion of low-mass dilepton spectra at full RHIC energy where most of the current experimental information at collider energies is available, from both PHENIX [48] and STAR [49] measurements.

3.3.1. Invariant-Mass Spectra

Figure 5 shows the comparison of thermal fireball calculations with low-mass spectra from STAR [49]. As compared to earlier calculations with a bag-model EoS [13], the use of lQCD-EoS and emission rates for the QGP enhances the pertinent yield significantly. It is now comparable to the in-medium hadronic contribution for masses below  GeV and takes over in the intermediate-mass region ( GeV). The hadronic part of the thermal yield remains prevalent in a wide range around the free mass, with a broad resonance structure and appreciable contributions from annihilation for  GeV. Upon adding the thermal yield to the final-state decay cocktail by the STAR collaboration (without decay), the MB data are well described. For the central data, a slight overestimate around  GeV and around the peak is found. A similar description [51] of the STAR data arises in a viscous hydrodynamic description of the medium using the spectral function from on-shell scattering amplitudes [21] (see also [52]) and in the parton-hadron string dynamics transport approach using a schematic temperature- and density-dependent broadening in a Breit-Wigner approximation of the spectral function [53]. More studies are needed to discern the sensitivity of the data to the in-medium spectral shape, as the latter significantly varies in the different approaches. For the PHENIX data (not shown), the enhancement as recently reported in [54] for noncentral collisions (carried out with the hadron-blind detector (HBD) upgrade) agrees with earlier measurements [48], is consistent with the STAR data, and thus should also agree with theory. For the most central Au-Au data, however, a large enhancement was reported in [48], which is well above theoretical calculations with broad spectral functions [13, 53, 55, 56], even in the MB sample. More “exotic” explanations of this effect, which did not figure at the SPS, for example, a Bose-condensed like glasma in the pre-equilibrium stages [57], have been put forward to explain the “PHENIX puzzle.” However, it is essential to first resolve the discrepancy on the experimental side, which is anticipated with the HBD measurement for central collisions.

To quantify the centrality dependence of the thermal radiation (or excess) yield, one commonly introduces an exponent as , which describes how the excess (or thermal) yield in a given mass range scales relative to the charged-particle multiplicity. For full RHIC energy, the theoretical calculation gives (with a ca. 10% error), similar to what had been found for integrated thermal photon yields [58].

3.3.2. Transverse-Momentum Dependencies

When corrected for acceptance, invariant-mass spectra are unaffected by any blue-shift of the expanding medium, which renders them a pristine probe for in-medium spectral modifications. However, the different collective flow associated with different sources may be helpful in discriminating them by investigating their spectra, see, for example, [26, 5964]. As is well known from the observed final-state hadron spectra, particles of larger mass experience a larger blue-shift than lighter particles due to collective motion with the expanding medium. Schematically, this can be represented by an effective slope parameter which for sufficiently large masses takes an approximate form of , where and denote the local temperature and average expansion velocity of the emitting source cell. Dileptons are well suited to systematically scan the mass dependence of by studying spectra for different mass bins (provided the data have sufficient statistics). At the SPS, this has been done by the NA60 collaboration [65], who found a gradual increase in the slope from the dimuon threshold to the mass characteristic for a source of hadronic origin (a.k.a in-medium mesons), a maximum around the mass (late decays), followed by a decrease and leveling off in the intermediate-mass region (IMR,  GeV) indicative for early emission at temperatures  MeV (where at the SPS the collective flow is still small).

Figure 6 shows the spectra for thermal radiation from hadronic matter and QGP in MB Au-Au(200 AGeV) in two typical mass regions where either of the two sources dominates. In the low-mass region (LMR), both sources have a surprisingly similar slope ( MeV), reiterating that the emission is from mostly around where the slope of both sources is comparable (also recall from Figure 5 that in the mass window  GeV the QGP emission is largest at the lower mass end, while the hadronic one is more weighted toward the higher end). For definiteness, assuming  MeV and  GeV, one finds that , which is right in the expected range [42]. On the other hand, in the IMR, where the QGP dominates, the hadronic slope has significantly increased to ca. 360 MeV due to the larger mass in the collective-flow term. On the other hand, the slope of the QGP emission has only slightly increased over the LMR, indicating that the increase in mass in the flow term is essentially offset by an earlier emission temperature, as expected for higher mass (for hadronic emission, the temperature is obviously limited by ). Consequently, at RHIC, the effective slope of the total thermal radiation in the IMR exceeds the one in the LMR, contrary to what has been observed at SPS. Together with blue-shift free temperature measurements from slopes in invariant-mass spectra, this provides a powerful tool for disentangling collective and thermal properties through EM radiation from the medium.

Alternatively, one can investigate the mass spectra in different momentum bins, possibly revealing a -dependence of the spectral shape, as was done for both data in Pb-Au [66] and in In-In [65] at SPS. Calculations for thermal radiation in Au-Au at full RHIC energy are shown in Figure 7 for four bins from  GeV. One indeed recognizes that the resonance structure becomes more pronounced as transverse momentum is increased. In the lowest bin, the minimum structure around  GeV is caused by the experimental acceptance, specifically the single-electron  GeV, which for vanishing suppresses all dilepton yields below  GeV.

3.3.3. Elliptic Flow

Another promising observable to diagnose the collectivity, and thus the origin of the EM emission source, is its elliptic flow [64, 69, 70]. The latter is particularly useful to discriminate early from late(r) thermal emission sources; contrary to the slope parameter, which is subject to an interplay of decreasing temperature and increasing flow, the medium's ellipticity is genuinely small (large) in the early (later) phases. Figure 8(a) shows hydrodynamic calculations of the inclusive thermal dilepton as a function of invariant mass (using the same emission rates and EoS as in the previous figures) [67]. One nicely recognizes a broad maximum structure around the mass, indicative for predominantly later emission in the vicinity of its vacuum mass, a characteristic mass dependence (together with an increasing QGP fraction) below, and a transition to a dominant QGP fraction with reduced above. All these features are essentially paralleling the mass dependence of the slope parameter at SPS, while the latter exhibits a marked increase at RHIC in the IMR due to the increased radial flow in the QGP and early hadronic phase. Rather similar results are obtained in hydrodynamic calculations with in-medium spectral functions from on-shell scattering amplitudes [51]. When using a less pronounced in-medium broadening, the peak structure in tends to become narrower [64, 69, 70]. First measurements of the dilepton- have been presented by STAR [68], see Figure 8(b). The shape of the data is not unlike the theoretical calculations, while it is also consistent with the simulated cocktail contribution. Note that the total is essentially a weighted sum of cocktail and excess radiation. Thus, if the total were to agree with the cocktail, it would imply that the of the excess radiation is as large as that of the cocktail. Clearly, future measurements with improved accuracy will be a rich source of information.

Significant measurements of EM excess radiation have recently been reported in the limit, that is, for direct photons, by both PHENIX [71, 72] and ALICE [73, 74]. A rather large signal has been observed in both experiments [72, 74], suggestive for rather late emission [75] (see also [7679]). In addition, the effective slope parameters of the excess radiation have been extracted,  MeV [71] at RHIC-200 and  MeV at LHC-2760 [73], which are rather soft once blue-shift effects are accounted for. In fact, these slopes are not inconsistent with the trends in the LMR dileptons when going from RHIC (Figure 6) to LHC (Figure 12). This would corroborate their main origin from around .

3.4. RHIC Beam Energy Scan

A central question for studying QCD phase structure is how the spectral properties of excitations behave as a function of chemical potential and temperature. With the EM (or vector) spectral function being the only one directly accessible via dileptons, systematic measurements as a function of beam energy are mandatory. At fixed target energies, this is being addressed by the current and future HADES efforts ( AGeV) [80, 81], by CBM for (Au) up to 35 AGeV [2], and has been measured at SPS energies at  AGeV, as well as in a CERES run at 40 AGeV [82].

At collider energies, a first systematic study of the excitation function of dilepton spectra has been conducted by STAR [68] as part of the beam-energy scan program at RHIC. The low-mass excess radiation develops smoothly when going down from  GeV via 62 GeV to 20 GeV (cf. Figure 9). Closer inspection reveals that the enhancement factor of excess radiation over cocktail in the region below the mass increases as the energy is reduced [68]. An indication of a similar trend was observed when comparing the CERES measurements in Pb-Au at  GeV and 8.8 GeV. Theoretically, this can be understood by the importance of baryons in the generation of medium effects [24], specifically, the low-mass enhancement in the spectral function. These medium effects become stronger as the beam energy is reduced since the hadronic medium close to becomes increasingly baryon rich. At the same time, the cocktail contributions, which are mostly made up by meson decays, decrease. The hadronic in-medium effects are expected to play a key role in the dilepton excess even at collider energies. The comparison with the STAR excitation function supports the interpretation of the excess radiation as originating from a melting resonance in the vicinity of .

A major objective of the beam-energy program is the search for a critical point. One of the main effects associated with a second-order endpoint is the critical slowing down of relaxation rates due to the increase in the correlation length in the system. For the medium expansion in URHICs, this may imply an “anomalous” increase in the lifetime of the interacting fireball. If this is so, dileptons may be an ideal tool to detect this phenomenon, since their total yield (as quantified by their enhancement factor) is directly proportional to the duration of emission. The NA60 data have shown that such a lifetime measurement can be carried out with an uncertainty of about [45]. In the calculations shown in Figure 9, no critical slowing down has been assumed; as a result, the average lifetime in MB Au-Au collisions increases smoothly from ca. 8 to . Thus, if a critical point were to exist and lead to a, say, 30% increase in the lifetime in a reasonably localized range of beam energies, dilepton yields ought to be able to detect this signature. This signal would further benefit from the fact that the prevalent radiation arises from around where the largest effect from the slowing down is expected.

3.5. LHC

The previous section raises the question whether the smooth excitation function of dilepton invariant-mass spectra in the RHIC regime will continue to LHC energies, which increase by another factor of 20. On the other hand, the dilepton spectra, especially their inverse-slope parameters, indicate an appreciable variation from SPS to RHIC, increasing from ca. 220 to 280 MeV in the LMR, and, more pronounced, from ca. 210 to about 320 MeV in the IMR. This is a direct consequence of the stronger (longer) development of collective flow in the QGP phase of the fireball evolution. This trend will continue at the LHC, as we will see below. In the following Section 3.5.1, we will first discuss the dielectron channel at LHC and highlight the excellent experimental capabilities that are anticipated with a planned major upgrade program of the ALICE detector [83]. In addition, ALICE can measure in the dimuon channel, albeit with somewhat more restrictive cuts whose impact will be illustrated in Section 3.5.2.

3.5.1. Dielectrons

The invariant-mass spectra of thermal radiation at LHC energies show a very similar shape and hadronic/QGP composition as at RHIC energy, see Figure 10. This is not surprising given the virtually identical in-medium hadronic and QGP rates along the thermodynamic trajectories at RHIC and LHC (where at chemical freezeout). It implies that the thermal radiation into the LMR is still dominated by temperatures around , with little (if any) sensitivity to the earliest phases. The total yield, on the other hand, increases substantially due to the larger fireball volumes created by the larger multiplicities. More quantitatively, the (-normalized) enhancement around, for example,  GeV, approximately scales as with relative to central Au-Au at full RHIC energy. This is a significantly stronger increase than the centrality dependent enhancement at fixed collision energy, as quoted in Section 3.3.1.

Detailed simulation studies of a proposed major upgrade of the ALICE detector have been conducted in the context of a pertinent letter of intent [83]. The final results after subtraction of uncorrelated (combinatorial) background are summarized in Figure 11, based on an excess signal given by the thermal contributions in Figure 10. (The thermal yields provided for the simulations were later found to contain a coding error in the author's implementation of the experimental acceptance; the error turns out to be rather inconsequential for the shape and relative composition of the signal, as a close comparison of Figures 10(b) and 11(b) reveals; the absolute differential yields differ by up to 20–30%.) Figure 10(a) shows that the thermal signal is dominant for the most part of the LMR (from ca. 0.2–1. GeV), while in the IMR it is outshined by correlated heavy-flavor decays. However, the latter can be effectively dealt with using displaced vertex cuts, while the excellent mass resolution, combined with measured and/or inferred knowledge of the Dalitz spectra of (from charged pions), (from charged kaons), and (from direct dilepton decays), facilitates a reliable subtraction of the cocktail. The resulting excess spectra shown in Figure 10(b) are of a quality comparable to the NA60 data. This will allow for quantitative studies of the in-medium EM spectral function in the LMR which are critical for being able to evaluate signatures of chiral restoration (as discussed elsewhere, see, e.g., [17, 85]). In addition, the yield and spectral slope of the dominantly QGP emission in the IMR will open a pristine window on QGP lifetime and temperature (recall that the spectra, which are little affected by the acceptance cuts in the IMR, are unaffected by blue shifts).

Let us turn to the dilepton spectra at full LHC energy, displayed again for two mass bins representing the LMR and IMR in Figure 12. Compared to RHIC, the LHC fireball is characterized by a marked increase in QGP lifetime and associated build-up of transverse flow by the time the system has cooled down to . Consequently, the spectra exhibit an appreciable increase in their inverse-slope parameters, by about 60% in the LMR (for both hadronic and QGP parts) and for the QGP part in the IMR, and up to 80% for the hadronic IMR radiation (recall that in a scenario with chiral mixing, the hadronic radiation for = 1.1–1.5 GeV is expected to increase by about a factor of 2, so that its larger slope compared to the QGP will become more significant for the total).

3.5.2. Dimuons

Low-mass dilepton measurements are also possible with ALICE in the dimuon channel at forward rapidities, , albeit with somewhat more restrictive momentum cuts [86]. The charged-particle multiplicity in this rapidity range is reduced by about 30% compared to midrapidity [87] but, at 2.76 ATeV, is still ca. 30% above central rapidities in central Au-Au at RHIC.

Figure 13 illustrates the expected thermal mass spectra in central Pb-Pb(2.76 ATeV). For “conservative” cuts on the di-/muons ( GeV,  GeV), their yield is substantially suppressed (see Figure 13(a)), by about one order of magnitude, compared to a typical single- cut of  GeV. In addition, the spectral broadening of the in-medium meson is less pronounced, a trend that was also observed in the -sliced NA60 dimuon spectra. Here, it is mostly due to the suppression of medium effects at larger -meson momentum relative to the heat bath, caused by hadronic form factors (analogous to RHIC, recall Figure 7). It is, in fact, mostly the pair cut which is responsible for the suppression, since is significantly larger than . If the former can be lowered to, say, 1 GeV, the thermal yield of accepted pairs increases by about a factor 3 in the IMR and 2 in the LMR (see dashed line in Figure 13(b)). The LMR acceptance is now mainly limited by the single- cut, as the latter suppresses low-mass pairs whose pair momentum is not at least (the same effect is responsible for the rather sharp decrease in acceptance for low-momentum electron pairs below  GeV in Figure 7(a), leading to a dip toward lower mass in the thermal spectra, even though the thermal rate increases approximately exponentially). This could be much improved by lowering the single- cut to, for example, 0.5 GeV, which would increase the low-mass yield by about a factor of 3. At the same time, the spectral broadening of the becomes more pronounced in the accepted yields; that is, the data would be more sensitive to medium effects.

4. Summary and Outlook

In this article, we have given an overview of medium modifications of the electromagnetic spectral function under conditions expected at collider energies (high temperature and small baryon chemical potential) and how these medium effects manifest themselves in experimental dilepton spectra at RHIC and LHC. For the emission rates from the hadronic phase, we have focused on predictions from effective hadronic Lagrangians evaluated with standard many-body (or thermal field-theory) techniques; no in-medium changes of the parameters in the Lagrangian (masses and couplings) have been assumed. Since this approach turned out to describe dilepton data at the SPS well, providing testable predictions for upcoming measurements at RHIC and LHC is in order. As collision energy increases, the QGP occupies an increasing space-time volume of the fireball evolution. To improve the description of the pertinent dilepton radiation, information from lattice-QCD has been implemented on (i) the equation of state around and above and (ii) nonperturbative dilepton emission rates in the QGP. The latter have been “minimally” extended to finite 3-momentum to facilitate their use in calculations of observables. Since these rates are rather similar to previously employed perturbative (HTL) rates, an approximate degeneracy of the in-medium hadronic and the lQCD rates close to still holds. This is welcome in view of the smooth crossover transition as deduced from bulk properties and order parameters at = 0.

The main features of the calculated thermal spectra at RHIC and LHC are as follows. The crossover in the lQCD EoS produces a noticeable increase of the QGP fraction of the yields (compared to a bag-model EoS), while the hadronic portion decreases (its former mixed-phase contribution has been swallowed by the QGP). However, due to the near-degeneracy of the QGP and hadronic emission rates near , both the total yield and its spectral shape change little; the hadronic part remains prevalent in an extended region around the mass at all collision energies. The very fact that an appreciable reshuffling of hadronic and QGP contributions from the transition region occurs indicates that the latter is a dominant source of low-mass dileptons at both RHIC and LHC. This creates a favorable setup for in-depth studies of the chiral restoration transition in a regime of the phase diagram where quantitative support from lQCD computations for order parameters and the EM correlator is available. Current ideas of how to render these connections rigorous have been reported elsewhere. Phenomenologically, it turns out that current RHIC data for LMR dilepton spectra are consistent with the melting scenario (with the caveat of the central Au-Au PHENIX data), including a recent first measurement of an excitation function all the way down to SPS energies. However, the accuracy of the current data does not yet suffice to discriminate in-medium spectral functions which differ considerably in detail. These “details” will have to be ironed out to enable definite tests of chiral restoration through the EM spectral function.

While the low-mass shape of the spectra is expected to be remarkably stable with collision energy, large variations are predicted in the excitation function of other dilepton observables. First, the total yields increase substantially with collision energy. In the LMR, the dependence on charged-particle rapidity density, , is estimated to scale as from RHIC to LHC, significantly stronger than as function of centrality at fixed . This is, of course, a direct consequence of the longer time it takes for the fireball to cool down to thermal freezeout. For the RHIC beam-energy scan program, this opens an exciting possibility to search for a nonmonotonous behavior in the fireball's lifetime due to a critical slowing down of the system's expansion. If the LMR radiation indeed emanates largely from the transition region, a slowed expansion around would take maximal advantage of this, thus rendering an “anomalous dilepton enhancement” a promising signature of a critical point.

Second, the transverse-momentum spectra of thermal dileptons are expected to become much harder with collision energy, reflecting the increase in the collective expansion generated by the QGP prior to the transition region. This “QGP barometer” provides a higher sensitivity than final-state hadron spectra which include the full collectivity of the hadronic evolution. The inverse-slope parameters for spectra in the LMR are expected to increase from ~220 MeV at SPS to ~280 MeV at RHIC-200 and up to ~450 MeV at LHC-5500. Even larger values are reached in the IMR, although the situation is a bit more involved here, since (a) the QGP emission is increasingly emitted from earlier phases and (b) the hadronic emission, while picking up the full effect of additional collectivity at , becomes subleading relative to the QGP. The trend in the LMR seems to line up with the recent slope measurements in photon excess spectra at RHIC and LHC. A similar connection exists for the elliptic flow; pertinent data will be of great interest. Invariant-mass spectra in the IMR remain the most promising observable to measure early QGP temperatures, once the correlated heavy-flavor decays can be either subtracted or reliably evaluated theoretically.

The versatility of dileptons at collider energies comprises a broad range of topics, ranging from chiral restoration to direct-temperature measurements, QGP collectivity, and fireball lifetime. Experimental efforts are well underway to exploit these, while sustained theoretical efforts will be required to provide thorough interpretations.

Acknowledgments

The author gratefully acknowledges the contributions of his collaborators, in particular Charles Gale, Min He, and Hendrik van Hees. He also thanks the STAR and ALICE collaborations for making available their plots shown in this article. This work has been supported by the U.S. National Science Foundation under Grant nos. PHY-0969394 and PHY-1306359 and by the A.-v.-Humboldt Foundation (Germany).