# Lattice studies of magnetic phenomena in heavy-ion collisions

###### Abstract

We review some experimental consequences of the presence of superstrong magnetic fields of order of the nuclear scale in noncentral heavy-ion collisions. We present lattice estimates for the strength of the Chiral Magnetic Effect (CME) for different quark flavours and argue that the dependence of the anisotropy of the distribution of emitted hadrons on their flavor content might be used as another experimental evidence of the CME. Another possible effect of superstrong magnetic field might be the observed abnormal enhancement of dilepton yield. We show that the presence of the magnetic field leads to a specific anisotropy of the dilepton emission rate.

## 1 Introduction

It has been realized recently that in heavy-ion collision experiments hadronic matter is affected not only by extremely high temperatures and densities, but also by superstrong magnetic fields with field strength being comparable to hadron masses squared. Such superstrong fields are created due to the relative motion of heavy ions themselves, since they carry large charge Kharzeev:08:1 ().

Obviously, the magnetic field is perpendicular to the collision plane, which can be reconstructed in experiment from the angular distribution of produced hadrons Abelev:09:1 (). There is no direct experimental way to measure the absolute value of the field strength, but it can be estimated in some microscopic transport model, such as the Ultrarelativistic Quantum Molecular Dynamics model (UrQMD) Skokov:09:1 ().

Probably the most notable effect which arises due to magnetic fields in heavy-ion collisions is the so-called Chiral Magnetic Effect. The essence of the effect is the generation of electric current along the direction of the external magnetic field in the background of topologically nontrivial gauge field configurations Kharzeev:08:1 (); Kharzeev:08:2 (). Such generation is not prohibited by -invariance, since topological charge density is a pseudoscalar field and thus nonzero topological charge explicitly breaks parity. However, since QCD is parity-invariant, the net current or the net electric charge should vanish when averaged over multiple collision events. Nevertheless, the nontrivial effect can still be detected if one considers dispersions of electric current or electric charge Buividovich:09:7 (); Shevchenko:10:1 (). Experimentally, the Chiral Magnetic Effect manifests itself as the dynamical enhancement of fluctuations of the numbers of charged hadrons emitted above and below the reaction plane Selyuzhenkov:06:1 (); Abelev:09:1 (); Voloshin:04:1 (); Voloshin:08:1 ().

Another effect, which is closely related to the CME, but has different experimental signatures, is the anisotropic electric conductivity of hadronic matter in the strong magnetic field, which was discovered in lattice simulations Buividovich:10:1 (). Since the conductivity of the hadronic matter is directly related to the lepton emission rate Teryaev:95:1 (); Gupta:04:1 (), such anisotropic conductivity should result in specific anisotropy of the dilepton emission rate w.r.t. the reaction plane. This anisotropy should grow with the centrality of the collision and with the charges of the colliding ions. This effect might also contribute to the observed abnormal dilepton yield in heavy-ion collisions Adare:09:1 (). Some theoretical considerations Chernodub:10:1 () as well as preliminary lattice data Buividovich:11:2 () suggest that at very strong magnetic fields with strength , where is the mass of the charged -meson, the anisotropic conductivity might even turn into the anisotropic superconductivity in the direction of the field. Unfortunately, such extremely large field strengths are hardly reachable with present-day heavy-ion colliders.

In this paper we give some estimates of the expected experimental signatures of superstrong magnetic fields, basing on the lattice data. In Section 2 we consider the Chiral Magnetic Effect and argue that its strength should decrease with increasing quark mass, which can be used to discriminate between the CME and other possible effects which might result in preferential emission of charged hadrons in the direction perpendicular to the reaction plane. In Section 3 we consider the dilepton emission rate, and estimate the contribution of the induced conductivity to the total dilepton yield and dilepton angular distribution in heavy-ion collisions.

## 2 Chiral Magnetic Effect

Chiral Magnetic Effect is usually characterized by the following experimental observables, suggested first in Ref. Voloshin:04:1 ():

(1) |

where denotes hadrons with positive or negative charges, respectively, is the number of events used for data analysis, and are the total multiplicities of positively/negatively charged particles produced in each event, and , are the angles w.r.t. the reaction plane at which the hadrons with labels and are emitted. Summation in (1) goes over all produced hadrons. In practice, only sufficiently energetic particles are considered.

The main signatures of the CME is the growth of with impact parameter, as well as the negativity of and and the positivity of Selyuzhenkov:06:1 (); Abelev:09:1 (); Voloshin:04:1 (); Voloshin:08:1 (). However, it has been pointed out recently that these results can also be explained by other effects, such as the influence of nuclear medium on jet formation Petersen:11:1 () and other in-medium effects Ma:11:1 (). It is therefore important to think about more refined experimental tests of the CME. Our main message in this paper is that the dependence of the charge fluctuations on quark mass can be used to discriminate between the CME and other possible phenomena which contribute to the observed asymmetry of charge fluctuations. Indeed, the CME emerges due to fluctuations of quark chirality Kharzeev:08:1 (), which are suppressed when the quark mass is increased.

The dependence of the observables (1) on the quark mass can be studied if one sums separately over charged mesons with different quark content, e.g. and (charged pions), and (charged kaons) or , , , (-mesons). The observables should then decrease with the meson mass. The dependence of on the centrality of the collision should also become weaker. Here we give a rough estimate of this effect basing on the results of lattice simulations.

Let us first note that the observables can be expressed in terms of the differences of multiplicities of charged hadrons emitted above and below the reaction plane Kharzeev:08:1 ():

(2) |

where , are the differences of the multiplicities of hadrons with positive or negative charges above or below the reaction plane, respectively. The factor depends on the hydrodynamical evolution of hadronic matter, and is usually close to unity. For multiplicities one can also neglect the initial charge of heavy ions with a good precision, and assume that , where is the mean multiplicity of the same-charge hadrons per event.

Lattice results can be compared to experimental data by considering the quantity , which can be expressed solely in terms of the difference of net charges of hadrons emitted above and below the reaction plane:

(3) |

In turn, the dispersion of the charge difference can be related to the vacuum expectation values of the squared current densities Buividovich:09:7 (). The contribution of each quark flavor to the total electromagnetic current is . We do not consider here the third-generation quarks, which are extremely rarely produced in heavy-ion collisions.

The simplest model which allows to express in terms of is the model of spherical fireball, which emits positively and negatively charged hadrons from its surface with intensity proportional to . This leads to the following relation Buividovich:09:7 ():

(4) |

where and are, respectively, the currents along the magnetic field and perpendicular to it, is some characteristic collision time, is the fireball radius and is some typical correlation length for electric charge density in the fireball. In our estimates, we take (this is a typical decay time for the magnetic field in heavy-ion collisions Skokov:09:1 ()), , which are reasonable parameters for, say, gold-gold collisions at nuclon. We also assume that the fireball is a sphere with radius within the overlapping region between the two heavy ions of radius which collide at impact parameter . The net multiplicity and the impact parameter as the functions of collision centrality can be found in Table 1 in Kharzeev:08:1 (). For simplicity, we assume that are approximately constant on the surface of the fireball. Note also, that in order to exclude the effects related to the dependence of multiplicities of strange and charmed mesons on the collision centrality, which might be different from that of light mesons, we normalize the charge of emitted hadrons by the square of the total multiplicity of all hadrons.

Several technical remarks are in order. First, we assume that the matter within the fireball is in the state of thermal equilibrium, and thus the expectation values can be calculated from gauge theory in Euclidean space. We also assume that the magnetic field is uniform and nearly time-independent. Of course, these are rather rough approximations, but we are aiming here at qualitative rather than quantitative estimates. We have calculated the currents in lattice gauge theory with background magnetic field both in the confinement phase, neglecting the contribution from the virtual quark loops (quenched approximation). A comparison with gauge theory suggests that this is a reasonable approximation Kalaydzhyan:10:1 (). A more detailed study has shown also that in the deconfinement phase the dispersions of local current densities are practically independent of the magnetic field Buividovich:09:7 (), thus we do not consider here this case. The expectation values contain also the ultraviolet divergent part, which we have removed by subtracting the corresponding expectation values at zero temperature and zero magnetic field.

The masses of the valence quarks took the values (for and quarks), (for -quark) and (for -quark). The lowest value of the quark mass is dictated by the numerical stability of our algorithm. Thus our calculations of at should be considered as only the lower bound for at realistic masses of and quarks.

On Fig. 2 we compare the quantity calculated for the experimental data by the STAR collaboration Voloshin:08:1 () with the estimate (4) based on the results of lattice simulations at different quark masses and different temperatures. In order to match the value of the magnetic field strength in experiment and in simulations, we take the rough estimate from eq. (A.12) in Kharzeev:08:1 ():

(5) |

where is the impact parameter. This rough fit also agrees by the order of magnitude with the results of more sophisticated calculations of the magnetic field within the UrQMD model Skokov:09:1 ().

One can see that the best agreement with the STAR data is obtained at the smallest quark mass. The combination quickly decreases as the quark mass is increased - approximately by a factor of as the quark mass changes from to , and by a factor of almost as the quark mass further increases to . Since all observables are typically of the same order, one can expect that each such observable will also decrease with the quark mass. This dependence of asymmetry of angular distributions of mesons of different flavors on their mass can be used to discriminate between the CME and other phenomena which might cause such asymmetry.

The result may be compared with the perturbative analog of CME Teryaev:11:1 () resulting from straightforward generalization of Heisenberg-Euler Lagrangian depending on quark mass as

(6) |

The correspondence of this perturbative and Abelian effect to CME is manifested by the substitution

(7) |

As a result, the perturbative mechanism may become essential when quark mass is exceeding the inverse correlation length of topological charge density. One may expect that the transition point from non-perturbative to perturbative mass dependence is not too far from the strange quark mass, like it happen for vacuum quark condensates and strangeness polarization in nucleons (see Teryaev:11:1 () and Ref. therein) so that the non perturbative lattice results are applicable for experimentally important case of strangeness separation.

## 3 Induced conductivity and abnormal dilepton yield

Another phenomenon which might be caused by superstrong magnetic fields acting on the hadronic matter is the induced anisotropic conductivity along the magnetic field Buividovich:10:1 (). While the Chiral Magnetic Effect is related to the local fluctuations of current density, the induced conductivity reflects the fact that these fluctuations also have long-range correlation in time. Indeed, by virtue of the Green-Kubo relations the conductivity is related to the zero-frequency limit of the spectral function which corresponds to the correlator in Minkowski space Kadanoff:63:1 ():

(8) |

This spectral function can also be extracted from the results of lattice simulations using the so-called Maximal Entropy Method Aarts:07:1 (); Asakawa:01:1 ().

The spectral function determines also the dilepton emission rate from either cold or hot hadronic matter McLerran:85:1 (); Teryaev:95:1 ():

(9) |

where and are the momenta of the leptons, ,

is the dilepton tensor ( is the Minkowski metric), is the lepton mass. Thus the low-momentum limit of is related, on the one hand, to the emission rate of soft dileptons, and, on the other hand, to the conductivity of hadronic matter.

The enhancement of the conductivity due to the magnetic field should thus lead to the enhancement of the dilepton emission rate. This might provide a viable explanation of the abnormal soft dilepton yield observed in heavy-ion collisions Adare:09:1 (). Moreover, the anisotropy of the conductivity should lead to specific correlations between the momenta of the dileptons and the direction of the magnetic field (in other words, with the orientation of the reaction plane).

In a similar way, one can also consider the correlators of charged vector currents , where are the flavour indices. The correlators of charged vector currents have been studied recently in Buividovich:11:2 () in order to test the conjectured -meson condensation in superstrong magnetic fields Chernodub:10:1 (). It was found that as the magnetic field is switched on, this correlator decays slower, which corresponds to the the increase of the corresponding spectral function at lower frequencies. Since charged currents are coupled to leptons and neutrinos, this effect might be manifested in the enhanced lepton and neutrino yield. Excitations over the -meson condensate might show up then in the usual decay channel .

Lattice simulations show that the electric conductivity is nonzero only in the direction of the magnetic field and depends linearly on Buividovich:10:1 (); Buividovich:10:4 (). For sufficiently small momenta and (and thus for small ) by virtue of the Green-Kubo theorem one has Kadanoff:63:1 (); Aarts:07:1 (). Let us also neglect the lepton masses and go to the rest frame of the dilepton pair, where , and the spatial components of the dilepton tensor are . The dilepton emission rate is therefore proportional to

(10) |

where is the angle between the spatial momentum of the outgoing leptons and the magnetic field. Therefore, there should be more soft dileptons emitted perpendicular to the magnetic field than parallel to it. As a result, they are to large extend hidden inside the hadrons in the scattering plane which should lead to the difficulty in their experimental observation. Similar predictions can be also made for the angular distributions of pairs in the case of charged meson decays (see above).

In order to estimate the effect of the magnetic field on the total dilepton yield in heavy-ion collisions, we normalize the conductivity induced by the magnetic field to the conductivity at zero magnetic field and at the temperature close to the deconfinement phase transition, . By virtue of (9), the ratio of these conductivities should be equal to the ratio of dilepton emission rates in the low-momentum region. The relevant lattice data is summarized in Buividovich:10:1 (); Buividovich:10:4 (). In these works it was found that in the deconfinement phase the conductivity is practically independent of the magnetic field. Therefore here we will try to estimate possible contribution of the magnetic field to the dilepton emission rate from hadronic matter in the confinement phase. From the data presented in Buividovich:10:1 (); Buividovich:10:4 (), we estimate the ratio of the induced conductivity to the conductivity of quark-gluon plasma as:

We now use the estimate (5) for the magnetic field strength and take into account that the ratio of conductivities should be equal to the ratio of dilepton emission rates. We thus obtain for the contribution of the magnetic field to the dilepton emission rate:

(11) |

where the factor appears after averaging the expression (10) over the angle .

We conclude that at large impact parameters () the dilepton yield can increase by up to due to the influence of the magnetic field. This factor should be essentially reduced because the magnetic-induced dileptons are hidden inside the scattering plane. The observed abnormal dilepton yield is, however, maximal for central collisions (where it reaches several hundred percent as compared to the hadron resonance model) and decreases as the impact parameter grows Adare:09:1 (). Such behavior might be caused by several factors, such as the change of the temperature within the fireball or the change of the fireball volume with impact parameter. A proper investigation of such factors is out of the scope of this paper, and cannot be undertaken using the methods of lattice gauge theory.

## 4 Conclusions

In this paper we have summarized the main experimental signatures of the effects caused by superstrong magnetic fields in heavy ion collisions, namely, the Chiral Magnetic Effect and the abnormal dilepton yield.

The Chiral Magnetic Effect Kharzeev:08:1 () results in preferential emission of charged hadrons in the direction perpendicular to the reaction plane. The origin of the Chiral Magnetic Effect is the fluctuations of chirality, which are suppressed as the quark mass grows. Thus this asymmetry in angular distributions of charged hadrons, characterized by the coefficients (1), should be strongly suppressed for strange or charmed hadrons.

The abnormal dilepton yield with specific angular dependence (10) is the consequence of electric conductivity of the hadronic matter induced by the magnetic field. In this case more dileptons are emitted in the direction perpendicular to the magnetic field. Let us also note that since according to our lattice data the magnetic field influences the conductivity only in the confinement phase Buividovich:10:1 (), the significant change of dilepton yield in noncentral heavy-ion collisions might be a signature of the confinement-deconfinement phase transition.

More generally, lattice data suggests that the influence of the magnetic field on the properties of hadronic matter is stronger in the confinement phase. Thus heavy-ion collision experiments on colliders with lower beam energy but with larger luminosity (such as FAIR in Darmstadt, Germany or NICA in Dubna, Russia) might be more advantageous for studying magnetic phenomena.

###### Acknowledgements.

The authors are grateful to M. N. Chernodub, A. S. Gorsky, V.I. Shevchenko, M. Stephanov and A. V. Zayakin for interesting and useful discussions. We are also deeply indebted to D. E. Kharzeev for very valuable and enlightening remarks on the present work. The work was supported by the Russian Ministry of Science and Education under contract No. 07.514.12.4028, by the Grant RFBR-11-02-01227-a of the Russian Foundation for Basic Research and by the Heisenberg-Landau program of JINR. P.B. was supported by the Sofja Kowalewskaja award from the Alexander von Humboldt Foundation. Numerical calculations were performed at the ITEP computer systems ”Graphyn” and ”Stakan” (authors are much obliged to A. V. Barylov, A. A. Golubev, V. A. Kolosov, I. E. Korolko, M. M. Sokolov for the help), the MVS 100K at Moscow Joint Supercomputer Center and at Supercomputing Center of the Moscow State University.## Bibliography

- (1) D.E. Kharzeev, L.D. McLerran, H.J. Warringa, Nucl. Phys. A 803, 227 (2008). ArXiv:0711.0950
- (2) B. Abelevet al. [STAR Collaboration], Phys. Rev. Lett. 103, 251601 (2009). ArXiv:0909.1739
- (3) V. Skokov, A. Illarionov, V. Toneev, Int. J. Mod. Phys. A 24, 5925 (2009). ArXiv:0907.1396
- (4) K. Fukushima, D.E. Kharzeev, H.J. Warringa, Phys. Rev. D 78, 074033 (2008). ArXiv:0808.3382
- (5) P.V. Buividovich, M.N. Chernodub, E.V. Luschevskaya, M.I. Polikarpov, Phys. Rev. D 80, 054503 (2009). ArXiv:0907.0494
- (6) V. Orlovsky, V. Shevchenko, Phys. Rev. D 82, 094032 (2010). ArXiv:1008.4977
- (7) I.V. Selyuzhenkov [STAR Collaboration], Rom. Rep. Phys. 58, 049 (2006). ArXiv:nucl-ex/0510069
- (8) S.A. Voloshin, Phys. Rev. C 70, 057901 (2004). ArXiv:hep-ph/0406311
- (9) S. A. Voloshin for the STAR Collaboration. Probe for the strong parity violation effects at RHIC with three-particle correlations. Proceedings of Quark Matter 2008 (2008). ArXiv:0806.0029
- (10) P.V. Buividovich, M.N. Chernodub, D.E. Kharzeev, T.K. Kalaydzhyan, E.V. Luschevskaya, M.I. Polikarpov, Phys. Rev. Lett. 105 (2010). ArXiv:1003.2180
- (11) E.L. Bratkovskaya, O.V. Teryaev, V.D. Toneev, Phys. Lett. B 348, 283 (1995).
- (12) S. Gupta, Phys. Lett. B 597, 57 (2004). ArXiv:hep-lat/0301006
- (13) A.A. et al. (PHENIX Collaboration), Phys. Rev. C 81, 034911 (2010). ArXiv:0912.0244
- (14) M.N. Chernodub, Phys. Rev. D 82, 085011 (2010). ArXiv:1008.1055
- (15) V.V. Braguta, P.V. Buividovich, M.N. Chernodub, M.I. Polikarpov. Electromagnetic conductivity of vacuum induced by strong magnetic field: numerical evidence in lattice gauge theory (2011). ArXiv:1104.3767
- (16) H. Petersen, T. Renk, S.A. Bass, Phys. Rev. C 83, 014916 (2011). ArXiv:1008.3846
- (17) G. Ma, B. Zhang, Phys. Lett. B 700, 39 (2011). ArXiv:1101.1701
- (18) V. Braguta, P. Buividovich, T. Kalaydzhyan, S. Kuznetsov, M. Polikarpov, PoS LAT2010, 190 (2010). ArXiv:1011.3795
- (19) O.V. Teryaev, Nucl.Phys.Proc.Suppl. 86, 219 (2011).
- (20) L.P. Kadanoff, P.C. Martin, Ann. Phys. 24, 419 (1963).
- (21) G. Aarts, C. Allton, J. Foley, S. Hands, S. Kim, Phys. Rev. Lett. 99, 022002 (2007). ArXiv:hep-lat/0703008
- (22) M. Asakawa, T. Hatsuda, Y. Nakahara, Prog. Part. Nucl. Phys. 46, 459 (2001). ArXiv:hep-lat/0011040
- (23) L.D. McLerran, T. Toimela, Phys. Rev. D 31, 545 (1985).
- (24) P.V. Buividovich, M.I. Polikarpov, Phys. Rev. D 83, 094508 (2011). ArXiv:1011.3001