# Thermal Dileptons from a Nonperturbative Quark-Gluon Phase

## I Introduction and Summary

Extensive lattice simulations have shown that QCD exhibits a phase transition at a temperature of the order of the infrared scale $`\mathrm{\Lambda }`$. This transition is believed to be a confinement-deconfinement transition for heavy quarks, and a chiral transition for light quarks. The high temperature phase is characterized by a small coupling constant $`g(T)`$. Unfortunately, perturbation theory in $`g(T)`$ is plagued with infrared singularities since the magnetic sector is non-perturbative at all temperatures.

High temperature lattice  and instanton simulations  seem to concur on the fact that in the vicinity of the critical temperature, $`T_c`$, the chromoelectric and chromomagnetic condensates are still substantial. A number of theoretical arguments, both perturbative  and nonperturbative , suggests that in the Euclidean formulation, the high temperature phase of QCD exhibits strong $`A_4A_4`$ correlations. It has even been argued that the global color symmetry might be spontaneously broken with the appearance of an $`A_4`$ condensate . Also, lattice simulations show that chiral symmetry is restored at the critical temperature, so that the quark condensate can be set to zero above $`T_c`$.

In the present analysis we will incorporate these three important observations that are pertinent to QCD near and above $`T_c`$; the quark condensate vanishes, the chromoelectric and chromomagnetic correlations remain non-zero, and strong correlations develop in the static part of the $`A_4`$ gauge field that can be described by a local $`A_4^2`$ condensate. In the spirit of the QCD sum rule method, we will not address the complicated dynamical questions of how the gluon condensates come about, but simply assume that they exist. We will also restrict our discussion to only those with dimension 2, $`\frac{\alpha _s}{\pi }A_4^2`$, and 4, $`\frac{\alpha _s}{\pi }E^2`$ and $`\frac{\alpha _s}{\pi }B^2`$, following an earlier suggestion by two of us . The value of the condensates can only be estimated using non-perturbative calculations such as lattice QCD, and at high temperature, the scale is expected to be set by the mass gap in the magnetic sector which is $`g^2T`$.

That chromoelectric and chromomagnetic condensates could be present in the high temperature phase is not very surprising, and their size can be roughly estimated using the bulk energy and momentum measurements at finite temperature on the lattice. The analysis in  gives the values : $`\frac{\alpha _s}{\pi }E^2\frac{\alpha _s}{\pi }B^2(200\mathrm{MeV})^4`$

The presence of an $`A_4^2`$ condensate is more controversial , so we will give a summary of the arguments presented in  for why such a condensate should be considered. The basic observation is that the Wilson line (or Polyakov loop) operator $`W=\mathrm{TrPexp}(i_0^\beta 𝑑x_4A_4(x_4))`$, is a gauge invariant operator that is a good order parameter for the pure glue theory. Although this operator is non-local, it is possible to relate it to a string of local operators by choosing the static gauge $`_4A_4=0`$. For example, in an SU(2) theory we have $`W=1+\mathrm{cos}(\frac{\beta }{2}A_4)`$, and one could now try to use factorization to extract $`A_4^2`$. Here we will simply equate the condensate to the magnetic mass, i.e. we will typically use, $`\frac{\alpha _s}{\pi }A_4^2/T^2\frac{\alpha _s}{\pi }(g^2)^216\pi \alpha _s^3`$, which is about $`0.4`$ for $`\alpha _s0.2`$ This estimate is also supported by a lattice calculation in pure Yang-Mills . Indeed, an analysis of the gluon-propagator in four-dimensions in Landau-gauge gives $`a^2(gA_4)^21/3`$ on a $`4\times 8^3`$ lattice, which is equivalent to $`\frac{\alpha _s}{\pi }A_4^2/T^20.4`$. A direct reconstruction of this condensate from the dimensionally reduced theory has proven to be more subtle .

In this letter, we present a calculation that suggests that the gluon condensates strongly affect the emission rates of low mass dileptons from the plasma. Our calculation will be very naive. We will simply assume that the quark-antiquark annihilation processes that give rise to the dileptons occur in the presence of low momentum fluctuating background gluon fields. By averaging over the background fields, we get leading corrections to the emission rates proportional to the various condensates. Note that in doing so we take the concept of a fluctuating background field at face value and perform the calculations directly in the time-like region assuming that the temperature is high enough for perturbation theory to be useful. Thus our method is very different from the QCD sum rules, where the background field method is just a convenient tool to calculate coefficient functions in the operator product expansion which is then used in the deep Euclidean region.<sup>#1</sup><sup>#1</sup>#1 A work similar in spirit to ours is the calculation by Patkos and Sakai  of the nonperturbative magnetic contribution to the $`e^+e^{}`$ annihilation rate. That calculation was inspired by the so called Copenhagen vacuum based on fluctuating magnetic fluxtubes.  However, since the emission rates can be expressed in terms of current-current correlation functions we can use the standard QCD sum rule *techniques* to calculate the condensate contributions. How this is done in detail, especially when it comes to gauge choice, is explained in . To calculate the contribution from the $`A_4^2`$ condensate we will use the following observation: A constant potential $`A_4`$ is the same as an imaginary (color) chemical potential, $`\mu ^a`$, and thus one can obtain the term $`A_4^2`$ by using a grand canonical ensemble for the quarks with a chemical potential $`\mu ^a=iA_4^a`$, *averaged* over the arbitrary directions, $`a`$, in color space, and pick up the term $`\mu ^2`$. (We have also checked that this procedure gives the same result as the standard fixed point gauge methods used in .)

At this point we want to stress that although our approach is inspired by the above theoretical arguments, it is basically phenomenological, and should primarily be judged by how successfully it can accommodate data. The theoretical uncertainties are not limited to the assumptions we made about the non-perturbative corrections to correlations functions in a high temperature quark gluon plasma. We also have to make several assumptions about the physics of the heavy ion collision, and we will use an oversimplified model, namely that of an expanding fireball containing only the plasma phase. The main motivation for using such a simple model (the details are given in section III below) is to demonstrate that the gluon condensate effects are important. Hopefully more realistic models e.g. including regions with mixed or hadronic phases, would also exhibit this sensitivity to the condensates. Also in the transition regime, the pure and mixed quark condensates may still be relevant.

Having said all this, we will demonstrate that a naive fire ball model combined with a simple parametrization of the production amplitudes based on the above ideas, and with the condensate parameters at the values suggested by the scale arguments, is quite successful in explaining the dilepton data from the CERES and HELIOS relativistic heavy ion experiments. The rest of this paper, where this claim is substantiated, is organized as follows: In section II, we calculate the effects of the leading order condensates on the dilepton emission rates from a thermalized quark-gluon phase. In section III, we use a simple fire-ball model to estimate the dilepton yield, and compare to the CERES and HELIOS experiments. In section IV, we show that our emission rates are also consistent with the recent analysis of the CERES results using different $`p_t`$ cuts.

## II Emission Rates

In a translationally invariant, but non-perturbative, quark-gluon phase, the rate $`𝐑`$ of dileptons produced in a unit four volume is directly related to the electromagnetic current-current correlation function . For a pair of leptons with mass $`m_l`$ and momenta $`p_{1,2}`$, the rate per unit invariant momentum $`q=p_1+p_2`$ is ,

$`{\displaystyle \frac{d𝐑}{d^4q}}`$ $`=`$ $`{\displaystyle \frac{\alpha ^2}{6\pi ^3q^2}}\left(1+{\displaystyle \frac{2m_l^2}{q^2}}\right)\left(1{\displaystyle \frac{4m_l^2}{q^2}}\right)^{1/2}𝐖(q),`$ (1)   
where $`\alpha =e^2/4\pi `$ is the fine structure constant, and

$`𝐖(q)={\displaystyle d^4xe^{iqx}\mathrm{Tr}\left(𝐏e^{iA_4^a𝐐^a}e^{(𝐇𝐅)/T}𝐉^\mu (x)𝐉_\mu (0)\right)}.`$ (2)   
Here $`e𝐉_\mu `$ is the hadronic part of the electromagnetic current, $`𝐇`$ the hadronic Hamiltonian, $`𝐅`$ the free energy, $`T`$ the temperature, and $`𝐐^a`$ a global color generator. (As mentioned above the direction, $`a`$, is arbitrary but fixed.) $`𝐖(q)`$ is related to the retarded current-current correlation function by,

$`𝐖(q)={\displaystyle \frac{2}{e^{q^0/T}1}}\mathrm{Im}𝐖^R(q)`$ (3)   
where

$`𝐖^R(q)=i{\displaystyle d^4xe^{iqx}\mathrm{Tr}\left(𝐏e^{iA_4^3𝐐^3}e^{(𝐇𝐅)/T}𝐑^{}\left(𝐉^\mu (x)𝐉_\mu (0)\right)\right)},`$ (4)   
and $`𝐑^{}`$ is the retarded product. The leading contribution to the retarded function is condensate independent, and is just the “Born” $`q\overline{q}`$ annihilation term, which for massless quarks becomes,

$`\mathrm{Im}𝐖^R`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}\left(N_c{\displaystyle \underset{q=u,d,s}{}}e_q^2\right)q^2\left[1+{\displaystyle \frac{2T}{|\stackrel{}{q}|}}\mathrm{ln}\left({\displaystyle \frac{n_+}{n_{}}}\right)\right],`$ (5)   
where $`N_c`$ is the number of colors, $`e_q`$ the charge of the quarks, and $`n_\pm =1/(e^{(q_0\pm |\stackrel{}{q}|)/2T}+1)`$. At $`T=0`$ this is just the $`R_3`$ ratio well known from $`e^+e^{}`$ annihilation

$`R_3(q^2)={\displaystyle \underset{i}{}}{\displaystyle \frac{\sigma _i(3\pi )}{\sigma (e^+e^{}\mu ^+\mu ^{})}}={\displaystyle \frac{4\pi }{q^2}}\mathrm{Im}𝐖_0^R|_{T=0}=N_c{\displaystyle \underset{q}{}}e_q^2.`$ (6)   
Using only the three light flavors, $`R_3`$ is $`2`$.

As mentioned above, the contribution of the $`A_4A_4`$ condensate to the retarded function $`𝐖^R`$ can easily be calculated by picking up the term quadratic in the imaginary colored chemical potential $`A_4^a`$. The result is

$`\mathrm{Im}𝐖_{11}^R`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}\left(N_c{\displaystyle \underset{q}{}}e_q^2\right)q^2\left[{\displaystyle \frac{\alpha _s}{\pi }}A_4^2\left({\displaystyle \frac{4\pi ^2}{T|\stackrel{}{q}|}}\right)\left(n_+(1n_+)n_{}(1n_{})\right)\right].`$ (7)   
According to the estimates in the previous section, $`\beta \frac{\alpha _s}{\pi }A_4^2/T^2`$ is of the order of $`0.4`$. Note that (7) has the same sign as (5) and will thus increase the dilepton yield. In the next section, we will see that this effect can be rather dramatic in the invariant mass region studied in the CERES experiments.

The contribution of the chromoelectric and chromomagnetic condensates to the retarded function can be calculated using the methods in , and the result for $`q^2>0`$ is,

$`\mathrm{Im}𝐖_{12}^R`$ $`=`$ $`{\displaystyle \frac{1}{4\pi }}\left(N_c{\displaystyle \underset{q}{}}e_q^2\right)\left[\left({\displaystyle \frac{1}{6}}{\displaystyle \frac{\alpha _s}{\pi }}E^2+{\displaystyle \frac{1}{3}}{\displaystyle \frac{\alpha _s}{\pi }}B^2\right)\left({\displaystyle \frac{4\pi ^2}{T|\stackrel{}{q}|}}\right)\left(n_+(1n_+)n_{}(1n_{})\right)\right].`$ (8)   
The temperature dependence of (8) is the same as that of (7) and can be understood from the condition of detailed balance. Note that the emission rate is enhanced by a chromomagnetic field and depleted by a chromoelectric field as expected<sup>#2</sup><sup>#2</sup>#2 Roughly speaking, the $`e^+e^{}`$ pairs are repelled from each other in the CM frame in the presence of a chromoelectric field, and deflected in the same direction in the presence of a chromomagnetic field., and also that all condensate contributions vanish at $`T=0`$.

## III Results

In Fig. 1 we show the dilepton emission rates at temperature $`T=150`$ MeV, using the parameters, $`\delta \frac{\alpha _s}{\pi }E^2/(200\mathrm{MeV})^4=\frac{\alpha _s}{\pi }B^2/(200\mathrm{MeV})^4`$ . The contribution from the non-perturbative plasma, i.e. perturbative plus condensate contributions (dotted and dash-dotted lines) should be compared to the purely perturbative part (thick solid line) and, as a reference, to that of a gas of mesons (thin solid line)  . The condensate contribution to the rates are sizable in the 500 MeV region, and the effect of the $`A_4^2`$ condensate, $`\beta `$, dominates that of the chromoelectric and chromomagnetic ones, $`\delta `$, for masses above 250 MeV. The dependence on both condensates $`\beta `$ and $`\delta `$ is negligible at energies higher than 2 GeV, and above 2.5 GeV the perturbative plasma gives the same result as that of a hadron gas, as expected from duality. In the Dalitz region ($`<200`$MeV) the effects of both the electric and magnetic condensates are appreciable and dominates the perturbative part by orders of magnitude. This observation maybe important for photon emission as well.

In a heavy-ion collision the electromagnetic emission occurs from various stages of the collision process, and hence at various temperatures. In this part, we will focus on the emission rates in the low and intermediate mass region (up to 1.5 GeV). We will assume that the dileptons emanate from a simple fire ball composed of the nonperturbative plasma phase only. The fire ball expansion is assumed to be homogeneous, and we will use a standard parametrization of volume and temperature,

$`V(t)`$ $`=`$ $`V_0\left(1+{\displaystyle \frac{t}{t_0}}\right)^3`$ (9)   
$`T(t)`$ $`=`$ $`(T_iT_{\mathrm{}})e^{t/\tau }+T_{\mathrm{}}.`$ (10)   
The parameters in (10) regarding the initial and freeze-out temperature will be taken from present transport calculations . We will not try to justify this point except by saying that it allows us to define a space-time volume for the “hot” spot in the heavy-ion collision that is consistent with the one used for an interacting meson-phase . Hence $`t_0=10`$ fm/c, $`T_i=170`$ MeV, $`T_{\mathrm{}}=110`$ MeV, $`\tau =8`$ fm/c, and the value of $`V_0`$ is absorbed into the over-all normalization constant $`N_0V_0=6.76\times 10^7`$ fm<sup>3</sup>. The freeze-out time will be set to $`t_{f.o.}=10`$ fm/c.

Using (10) and the above rate, the final expression for the integrated emission rate per unit rapidity $`\eta `$ and invariant mass $`M`$ is

$`{\displaystyle \frac{dN/d\eta dM}{dN_{ch}/d\eta }}(M)`$ $`=`$ $`N_0M{\displaystyle _0^{t_{f.o.}}}𝑑tV(t){\displaystyle \frac{d^3q}{q_0}A(q_0,q^2)\frac{d𝐑}{d^4q}}.`$ (11)   
The acceptance function $`A`$ enforces the detector cut $`p_t>200MeV,2.1<\eta <2.65`$, and $`\mathrm{\Theta }_{ee}>35`$ mrad for CERES. In Fig. 2, we show the evolved emission rates for CERES. The enhancement is sizable and goes in the direction of the data in the 500 MeV range in CERES. Again, the nonperturbative QCD results are indicated by the dotted line for the condensates we quoted, the dashed-dotted line (without the E- and B-condensates), while the perturbative QCD result is indicated by the solid line and the meson result by the thin solid line.

For HELIOS-3, (11) can be used by integrating over $`\eta `$ with the cut $`m_t4(72\eta ),m_t\sqrt{(2m_\mu )^2+(15/\mathrm{cosh}(\eta ))^2}`$. In Fig. 3, we show the evolved emission rates for HELIOS. For this energy range the effect of chromoelectric and chromomagnetic condensates are small irrespective of the specific value for the thermal condensates as indicated in Figs. 1 and 2. The enhancement is compatible with the HELIOS data, and slightly overshoots both the data and the meson gas in the 1 GeV region.

## IV $`p_t`$ Cuts

Another important check on the nonperturbative emission rates can be obtained by using the recent $`p_t`$ cut analysis by the CERES collaboration to probe the transverse momentum content of the dilepton emission rates in the fire-ball. Our results versus the data are shown in Figs. 4-6, with and without the chromoelectric and chromomagnetic condensates. The $`p_t`$ dependent effects caused by the thermal condensates are most dramatic in the low $`p_t`$ region (Fig. 4), where there is a sizable condensate effect, mainly due to $`\beta `$. The condensate contributions are in the right direction around the low mass dielectron enhancement, and give important additions to the perturbative result.

In the high $`p_t`$ region, Fig. 5, the resulting condensate effects are rather small, and the curve is consistent with the present data. Therefore, the total transverse momentum content (Fig. 6) is quite accurately reproduced by the addition of the gluon condensates.

In summary, we find it encouraging that not only the dilepton emission rates, but also the observed low mass-region $`p_t`$ enhancement is, at least partially, explained by our simple model, and we think that this further strengthens our proposal that nonperturbative condensate effects play a role in the plasma phase of heavy ion collisions.

## Acknowledgements

We thank a number of our colleagues at Stony-Brook for discussions. We also thank Keijo Kajantie and Su Houng Lee for their comments on the manuscript. This work was supported in part by the U.S. Department of Energy under Grant No. DE-FG02-88ER40388, and the Swedish Natural Science Research Council.