Abstract

Various novel transport phenomena in chiral systems result from the interplay of quantum anomalies with magnetic field and vorticity in high-energy heavy-ion collisions and could survive the expansion of the fireball and be detected in experiments. Among them are the chiral magnetic effect, the chiral vortical effect, and the chiral magnetic wave, the experimental searches for which have aroused extensive interest. The goal of this review is to describe the current status of experimental studies at Relativistic Heavy-Ion Collider at BNL and the Large Hadron Collider at CERN and to outline the future work in experiment needed to eliminate the existing uncertainties in the interpretation of the data.

1. Introduction

High-energy heavy-ion collisions can produce a hot, dense, and deconfined nuclear medium, which is dubbed the quark-gluon plasma (QGP). The thermodynamic states of a QGP can be specified by the axial chemical potential , besides the temperature and the vector chemical potential . characterizes the imbalance of right-handed and left-handed fermions in a system, and a chiral system bears a nonzero . Chiral domains may be created locally in heavy-ion collisions through various mechanisms on an event-by-event basis (e.g., topological fluctuations in the gluonic sector, glasma flux tubes, or fluctuations in the quark sector) [16]. In a noncentral collision, a strong magnetic field ( T) can be produced (mostly by energetic spectator protons) [2, 3] and will induce an electric current along in chiral domains, , which is called the chiral magnetic effect (CME) [1, 2]. On average, is perpendicular to the so-called reaction plane () that contains the impact parameter and the beam momenta, as depicted in Figure 1. Hence the CME will manifest a charge transport across the reaction plane.

In the presence of the CME and other modes of collective motions, we can Fourier decompose the azimuthal distribution of particles of given transverse momentum () and pseudorapidity ():where is the azimuthal angle of a particle and . Here the subscript (+ or −) denotes the charge sign of the particle. Conventionally is called “directed flow” and “elliptic flow” [7]. The parameter (with ) quantifies the electric charge separation due to the CME.

An anomalous transport effect can also occur when a chiral system undergoes a global rotation. The fluid rotation can be quantified by a vorticity , where is the flow velocity field. For a given vorticity , the chiral vortical effect (CVE) induces a vector current [8]. Whereas the CME is driven by , the CVE is driven by in a chiral medium. Here the subscript “” means “vector,” and can be, for example, “” (baryon) or “” (electron). In heavy-ion collisions, is typically larger than by an order of magnitude, making it easier to search for the CVE via the baryonic-charge separation than the electric charge separation. Hence the subscript in (1) represents baryon or antibaryon in the CVE search.

Another complementary transport phenomenon to the CME has been found and named the chiral separation effect (CSE) [9, 10], in which chiral charges are separated along the magnetic field direction in the presence of a finite vector-charge density: . In a chirally symmetric phase, the CME and CSE form a collective excitation, the chiral magnetic wave (CMW), a long wavelength hydrodynamic mode of chiral charge densities [11, 12]. The CMW is a signature of the chiral symmetry restoration and manifests itself in a finite electric quadrupole moment of the collision system, where the “poles” (“equator”) of the produced fireball acquire additional positive (negative) charge [11]. This effect, if present, will be reflected in the measurements of charge-dependent elliptic flow.

There are other chiral magnetic/vortical effects such as the chiral electric separation effect (CESE) [13, 14] and the chiral vortical wave (CVW) [15]; see [16] for a recent review on these effects. This article reviews the experimental results in the past decade to search for the chiral magnetic/vortical effects in high-energy heavy-ion collisions: evidence for the initial magnetic field and vorticity in Section 2, the observation of the electric (baryonic) charge separation in Section 3 (Section 4), and the manifestation of the electric quadrupole moment in Section 5. An outlook for future development is discussed in Section 6.

2. Driving Force

We may intuitively regard the magnetic field (vorticity) as the driving force of the CME (CVE), while the chirality imbalance is the initial condition, and the electric (baryonic) charge separation is the manifestation. A rough estimate of the initial magnetic field gives , where , is the impact parameter, and is the Lorentz factor. Therefore a typical Au+Au collision at = 200 GeV produces . Many computations have attempted to quantify the electromagnetic field on an event-by-event basis (see, e.g., [1720]), in terms of the spatial distribution, the orientation fluctuation as well as the dependence on colliding nuclei, centrality, and beam energy.

A major uncertainty in theoretical calculations of the magnetic field is its duration in the QCD fluid created in the heavy-ion collision (see, e.g., [2124]). The time dependence of after the impact of the two nuclei crucially depends on whether/when/how a conducting medium may form and the lifetime of the magnetic field may be elongated. The electric conductivity [25] and the time evolution of the quark densities [26] can be studied via directed flow of charged hadrons in asymmetric collisions, such as Cu+Au. Figure 2 illustrates the transverse plane for a Cu+Au collision with = 6 fm [27]. The difference in the number of protons creates a strong electric field in the initial stage of the collision, pointing along the arrow in the figure. The lifetime of the electric field might be very short (e.g., from [25, 26]), but the electric charges from quarks and antiquarks that are present in the early stage of the collision would experience the Coulomb force and lift the degeneracy in between positively and negatively charged particles [18, 19, 25]:where denotes the azimuthal angle of the electric field and the coefficient characterizes the strength of the dipole deformation induced by the electric field and is proportional to the electric conductivity of the medium. Here represents the rapidity-even component of directed flow that is dominant in asymmetric collisions, whereas in symmetric collisions conventionally denotes the rapidity-odd component.

Figure 3 shows recent STAR measurements of charge-dependent and the difference as functions of in 10–40% Cu+Au and Au+Au collisions [27]. For < 2 GeV/, seems to increase with . The results from Au+Au collisions have much smaller values, roughly by a factor of 10, than those in Cu+Au. Note that in Au+Au collisions is similarly small [28]. in Au+Au is consistent with zero. Calculations for charged pions from the parton-hadron-string-dynamics (PHSD) model [26], which is a dynamical transport approach in the partonic and hadronic phases, are compared with the data. The PHSD model calculates two cases: charge-dependent with and without the initial electric field (EF). In the case with the EF switched on, the model assumes that all electric charges are affected by the EF and this results in a large separation of between positive and negative particles as shown in Figure 3(a). In Figure 3(b), the calculations of with and without the EF are shown together, but note that the EF-on calculation points are scaled by 0.1. After scaling by 0.1, the model describes rather well the dependence of the measured data for < 2 GeV/. This qualitative evidence for the strong initial electric field in asymmetric collisions provides an indirect evidence for the strong initial magnetic field in heavy-ion collisions that could leave an imprint on the final-stage particles.

The vorticity is induced by the global rotation of the QGP in heavy-ion collisions. In a noncentral collision, the majority of the global angular momentum, , is carried away by spectator nucleons. However, a considerable fraction (about 10–20%) of could remain in the QGP and be approximately conserved in time [29, 30]. This implies a relatively long duration of the vortical effects. On average, the angular momentum is pointing in the out-of-plane direction, so the CME and CVE are very much alike in terms of their experimental observables. Attempts to compute local vorticity and its space-time distribution have also been made [2934].

Experimentally global polarization of hyperons such as provides a measure for both the plasma vorticity and the magnetic field. Whereas the vortical effects will generate a positive polarization for both and , the coupling of the hadronic magnetic dipole moment to the magnetic field will generate a positive contribution for and a negative one for . Therefore, a splitting between and polarization will be a direct evidence for the strong initial magnetic field. Recently, preliminary STAR measurements [35] have reported the first observation of global and polarization in heavy-ion collisions. At < 100 GeV, the signal is on the order of a few percent and displays a weak beam energy dependence. The average polarization over and evidences the plasma vorticity, while the splitting observation requires much higher statistics to be delivered in the second beam energy scan (BES-II) program at RHIC [36] to signify the magnetic field.

The future search for evidence for the initial magnetic field (vorticity) can be conducted via photon (vector meson) polarization measurements [37]. The initial magnetic helicity () of the collision system can be quite large and bears opposite signs in the upper and lower hemispheres. Owing to the chiral anomaly, the helicity can be transferred back and forth between the magnetic flux and fermions as the collision system evolves, so that the magnetic helicity could last long enough to yield photons with opposite circular polarizations in the hemispheres above and below the reaction plane [3841]. A similar asymmetry in photon polarization can also result from the initial global quark polarization [42], which could effectively lead to a polarization of photons [39]. This local imbalance of photon circular polarization could be observed in experiments, for example, by studying the polarization preference with respect to the reaction plane for photons that convert into pairs [37]. Similarly, vector mesons that decay into two daughters can also have their polarization preferences measured with the scheme outlined in [37], and the helicity separation in this case originates from vorticity [4244].

3. Chiral Magnetic Effect

In high-energy heavy-ion collisions, the occurrence of the CME requires the simultaneous existence of ultra-strong magnetic fields, chiral symmetry restoration, and topological charge changing transitions. Therefore, the confirmation of the CME will further our understanding of all these fundamental physics mechanisms. The experimental searches for the CME have been carried out extensively in the past decade at RHIC and the LHC. This section will introduce the observables pertinent to the electric charge separation induced by the CME, present the experimental results, and discuss the background contributions due to the coupling of elliptic flow and other physics mechanisms.

3.1. Charge Separation Observable

It is tempting to directly measure the event average of . However, from event to event, the signs of the values are equally likely, and the signs of finite and will flip accordingly, which will result in . Figure 4 presents the STAR measurements of with the 1st harmonic event plane reconstructed from spectator neutrons [45]. These results indicate nonsignificant charge dependence in all centrality intervals, and the difference between positive and negative charges is typically less than 10−4.

One solution to search for the CME is via charge separation fluctuations perpendicular to the reaction plane, for example, with a three-point correlator [55], , where the averaging is done over all particles in an event and over all events. In practice, the reaction plane is approximated with the “event plane” () reconstructed with measured particles, and then the measurement is corrected for the finite event plane resolution. The expansion of the γ correlator,reveals the difference between the in-plane and out-of-plane projections of the correlations. The first term () in the expansion provides a baseline unrelated to the magnetic field. The background contribution () is suppressed to a level close to the magnitude of [55].

The γ correlator was first measured by the STAR Collaboration with the 2nd harmonic event plane for Au+Au (shown with crosses in Figure 5) and Cu+Cu (not shown here) collisions at 62.4 and 200 GeV with data from the 2004/2005 RHIC runs [46, 47]. All the results have been found to be in qualitative expectation with the CME: the opposite-charge () is above the same-charge () correlations, displaying a “right” ordering. The opposite-charge correlations in Cu+Cu collisions are stronger than those in Au+Au, possibly because of the suppression of the correlations among oppositely moving particles in a larger system. STAR also presented and dependences of the signal. The signal has a width of about one unit of rapidity, suggesting that short-range correlations are not the dominant contributions. The signal is found to increase with the pair average , and it was later shown [56] that the radial flow expansion can explain this feature.

Similar γ results for 200 GeV Au+Au and 2.76 TeV Pb+Pb were observed by the PHENIX Collaboration [57] and the ALICE Collaboration [49], respectively. PHENIX also employed a multiparticle charge-sensitive correlator, [58], and their preliminary results showed a concave distribution, also evidencing the charge separation effect. Recently this method was also applied to the STAR data, which yielded qualitatively similar results [59]. However, there are a few issues with this method that keep us from further discussing its results. First, this correlator employs charge-sign-randomized particles to provide a baseline for the charge separation, but the randomization will introduce a bias by destroying local charge conservation. Second, there is no correction for the event plane resolution in the measurement, making it impossible to directly compare the measurements between different experiments. Third, the extraction of the charge separation signal relies on the comparison with a simplified Monte Carlo simulation, which has to well describe many experimental phenomena, including , spectra, event plane resolution, transverse momentum conservation, and local charge conservation. Among other things, we are not aware of any simple Monte Carlo simulation that properly handles momentum conservation. In short, the reliability of this method depends on how realistic the simple simulation can be.

The conventional-physics backgrounds were investigated by employing heavy-ion event generators MEVSIM [50], UrQMD [60], and HIJING [61] (with and without an elliptic flow afterburner implemented). None of these generators could achieve reasonable agreement with the data. However, these generators do not well describe main features of heavy-ion collisions, so their failure to describe effects as subtle as charge separation provides only limited support to the CME picture.

The charge separation signal was cross-checked with data from the 2007 RHIC run (shown in Figure 5) [45]. The γ correlations from these data were measured with respect to both the 1st harmonic event plane (of spectators at large rapidity) and the 2nd harmonic event planes at mid-rapidity. The agreement of results using different types of event planes ensures that the signal is not caused by three-particle background correlations and is truly due to reaction-plane-related correlations.

Another test was carried out by replacing one of the two charged particles in γ with a neutral particle, for example, , and the results show no separation between and [62]. Thus the charge separation observed in the γ correlation between two charged particles is indeed due to the electric charge. Short-range correlations are visible in the differential measurements of at low relative momenta, which are related to quantum interference (“HBT”) and final-state-interactions (Coulomb dominated) [45], and they belong to late-stage effects. To suppress short-range correlations, the conditions of > 0.15 GeV/ and were applied to the three-point correlator, shown with the gray bars in Figure 6. When pairs with low relative momenta are excluded, the positive opposite-charge correlations in peripheral collisions are significantly reduced, but the difference between same- and opposite-charge correlations remains largely unchanged.

Different azimuthal regions of charge separation are weighted differently in the correlator; that is, oppositely charged pairs emitted azimuthally at 90° from the event plane (maximally out-of-plane) are weighted more heavily than those emitted only a few degrees from the event plane (minimally out-of-plane). A test was carried out to modify the γ correlator such that all azimuthal regions of charge separation receive an identical weight. This may be done by first rewriting (3) aswhere and stand for the absolute magnitude () and the sign (±1) of the sine or cosine function, respectively. IN represents the cosine part of (3) (in-plane), and OUT represents the sine part (out-of-plane). A modulated sign correlation (msc) is obtained by reducing the γ correlator [45]:The modulated sign correlations are compared with the three-point correlator for Au+Au collisions at 200 GeV in Figure 6. The msc turns out to well reproduce the same trend as the three-point correlator although their magnitudes differ slightly. STAR also invented another approach called the charge multiplicity asymmetry correlation (CMAC), whose methodology is similar to the msc, and yielded very similar results [63].

The beam energy dependence of the correlation could further our understanding of the origin of the observed charge separation. The charge separation effect depends strongly on the formation of the QGP and chiral symmetry restoration [1], so the signal will diminish if a QGP is short-lived or even absent at low collision energies. The fact that the lifetime of the magnetic field is longer at smaller collision energies could lead to an almost threshold effect: with decreasing collision energy, the signal might slowly increase with an abrupt drop thereafter. Unfortunately, the exact energy dependence of the CME is not calculated yet.

Figure 7 presents and correlators as functions of centrality for Au+Au collisions at = 7.7–200 GeV measured by STAR [48] and for Pb+Pb collisions at 2.76 TeV by ALICE [49]. In most cases, is still above , manifesting extra charge separation fluctuations perpendicular to the reaction plane. With decreased beam energy, both and tend to rise up starting from peripheral collisions. This feature seems to be charge-independent and can be explained by momentum conservation and elliptic flow [45]. Momentum conservation forces all produced particles, regardless of charge, to separate from each other, whereas elliptic flow works in the opposite sense. In peripheral collisions, momentum conservation dominates because the multiplicity () is small. The lower beam energy, the smaller , and the higher and . In more central collisions where the multiplicity is large enough, this type of charge-independent background can be described by [45, 64]. MEVSIM is a Monte Carlo event generator developed for STAR simulations [50]. In Figure 7, we also show the model calculations of MEVSIM with the implementation of and momentum conservation, which qualitatively describe the beam energy dependence of the charge-independent background. The difference between and seems to vanish at low collision energies, but the interpretation involves an ambiguity to be discussed in Section 3.2.

3.2. Flow Backgrounds

By construction, the correlator contains the background terms and , and their difference was originally studied for the “flowing cluster” case [55]:where is the cluster emission azimuthal angle and and are the azimuthal angles of two decay products. The flowing cluster can be generalized to a larger portion of or even the whole event, through the mechanisms of transverse momentum conservation (TMC) [65, 66] and/or local charge conservation (LCC) [67]. We introduce a useful tool to study the background, the two-particle correlator, . Ideally should be proportional to , but in reality it is dominated by backgrounds. For example, the TMC effect leads to the following pertinent correlation terms in δ and γ [66]:where and and represent the - and -weighted moments of , respectively. The subscript “” denotes an average of all produced particles in the full phase space; the actual measurements will be only in a fraction of the full space, denoted by “Ω.” The background contribution due to the LCC effect has a similar characteristic structure as the above [65, 67].

It is convenient to express γ and with a two-component framework [48, 66]:where and are the CME and background contributions, respectively. The background-subtracted correlator, , can be obtained from the ensemble averages of several observables:The major uncertainty in the above expression, the coefficient , depends on particle charge combination, particle transverse momentum, centrality, and collision energy.

Figure 8 shows the values estimated for Au+Au collisions at 200 GeV [51], with the measurements by the PHOBOS collaboration [68, 69], and with the calculations from the AMPT model [7072]. Here only the TMC effect has been taken in account, and is typically within for . The values attained this way will vary slightly if a smaller acceptance is involved. In reality, should be averaged over various mechanisms such as TMC, LCC, and resonance decays. The AMPT model gives a more comprehensive estimate in Figure 8 via , where the numerator is solely due to flow backgrounds. For the centrality range of 10–50%, where the statistical uncertainties are small, the values thus obtained are close to those estimated with the information.

Figure 9 shows () as a function of beam energy for 30–60% Au+Au (Pb+Pb) collisions [48, 49]. is depicted with the solid line. In both cases of , demonstrates a weak energy dependence above 19.6 GeV and tends to diminish from 19.6 to 7.7 GeV, though the statistical errors are large for 7.7 GeV. This may be explained by the probable domination of hadronic interactions over partonic ones at low energies. A more definitive conclusion may be reached with a more accurate estimation of and with higher statistics at lower energies in the proposed phase II of the RHIC Beam Energy Scan program, as illustrated by the shaded band in Figure 9. Recently the CMS collaboration carried out the γ correlations in p+Pb and Pb+Pb collisions at 5.02 TeV and observed similar magnitudes of versus multiplicity between p+Pb and Pb+Pb collisions [73]. This seems to indicate no CME signals in Pb+Pb collisions at this collision energy, which could be along the trend shown in Figure 9: a simplified extrapolation to 5.02 TeV points to zero if is around 1.5. A further measurement of the correlator at 5.02 TeV would help verify/falsify this speculation.

Collisions of Uranium nuclei have been utilized at RHIC to study the dependence of multiplicity production, flow, and the CME on the initial overlap geometry [74, 75]. It was initially thought that the prolate shape of Uranium nuclei would facilitate the selection of nearly fully overlapping events with large elliptic flow values, but with small magnetic fields. However, the square of the magnetic field may be not particularly small, owing to fluctuations. Measurements of very central collisions also demonstrated that the multiplicity does not depend on the configuration of the collisions as strongly as anticipated in the two-component multiplicity model, which significantly reduces the feasibility to independently manipulate the flow and the magnetic field [76].

Figure 10 shows measurements of for the 0–10% centrality range in 193 GeV U+U (a) and 200 GeV Au+Au collisions (b) [5254]. In both U+U and Au+Au collisions, the signal increases almost linearly with , suggesting that the charge separation observable is dominated by a -related background. However, the charge separation goes to zero while is still sizable in central Au+Au and central U+U collisions. Model calculations show that the quantity as a function of eccentricity exhibits the same trend [53]: although remains large owing to fluctuations, goes to zero as and become decorrelated in very central collisions. In other words, while fluctuations in central collisions keep the participant eccentricity positive-definite, the decorrelation of and drives to zero. Therefore, the data are in better agreement with a CME picture than a flow background interpretation. A phenomenological study for extrapolating both signal and background from Au+Au to U+U collisions was done in [77].

Another measurement that exploits the 4th harmonic event plane could clarify the origin of the charge-dependent correlations and the role of the background [78]. The correlations, , should include the effect of the flow-related background, but not any contribution from the CME. The correlations due to the background in this case are expected to be smaller in magnitude since the 4th harmonic flow is usually weaker than the elliptic flow. The preliminary results of such measurements by ALICE are presented in Figure 11 [79] with separate same-charge and opposite-charge correlations relative to the 4th harmonic event plane as functions of centrality (a) and the difference between same-charge and opposite-charge correlations with respect to the 2nd and 4th harmonic event planes (b). The correlation difference relative to the 4th harmonic event plane is very small and suggestive of small background contributions. Detailed simulations have to be performed to draw more definite conclusion from this measurement.

4. Chiral Vortical Effect

The chiral vortical effect (CVE) is in analog to the CME, and its experimental manifestation is the baryonic-charge separation, instead of the electric charge separation, perpendicular to the reaction plane. Therefore, the three-point γ correlator is still applicable, only now between two (anti)baryons. However, if both particles are (anti)protons that carry also electric charges, there will be an ambiguity due to the possible CME signal. The study of the γ correlator with an electrically neutral baryon, such as Λ, will provide more conclusive evidence for the baryonic-charge separation effect.

Although (anti)Λ’s are electrically neutral, we still need to answer whether strange quarks behave the same way as up/down quarks in the chiral dynamics during the collision. If the answer is no, then (anti)Λ’s may still act like electrically charged particles in the γ correlation. Figure 12 shows the γ correlation of () and () as functions of centrality in Au+Au collisions at 200 GeV [62]. In this analysis, (anti)protons are excluded from the charged hadrons in the correlator to avoid any possible CVE contribution. Tentatively assigning a positive (negative) charge to Λ (), we observe that the “same-charge” and “opposite-charge” correlations are consistent with each other; hence there is no charge-dependent effect. As mentioned in Section 3.1, the null correlations [62] tell that the difference between same-charge and opposite-charge correlations as shown in Figure 5 is truly due to the electric charge, so the null charge separation effect in manifests no electric charges of (anti)Λ’s in the γ correlation. Thus strange quarks seem to behave the same way as up/down quarks in the chiral dynamics. Furthermore, the correlation provides a baseline for the correlation, and any possible signal in the latter should not come from the CME contribution.

Figure 13 shows γ correlations of () and () as functions of centrality in Au+Au collisions at 200 GeV [62]. The same-baryonic-charge correlation is above the opposite-baryonic-charge correlation from mid-central to peripheral collisions. This baryonic-charge separation with respect to the event plane is consistent with the CVE picture. More investigations into the background contribution are needed. For example, in analog to the local charge conservation, the local baryonic-charge conservation could play a similar role as LCC when coupled to . The magnitudes of the correlations are much larger than those of the correlations. This is partially because of baryons is higher than that of mesons, and the correlation strength increases with the average of the particle pair in the correlation. Future differential measurements versus the pair average and further correlations between identified particles may provide a better insight into the correlation strength between the CME- and CVE-related correlations.

5. Chiral Magnetic Wave

The CMW is a signature of the chiral symmetry restoration in the QGP and consists of actually two chiral seamless modes traveling at the same speed [11]: the right-handed (left-handed) wave transports the right-handed (left-handed) density and current in the direction parallel (antiparallel) to the direction. A more general analysis [84] studied various possible collective modes based on a non-neutral-background QGP (i.e., with nonzero and/or ) in external electric and/or magnetic fields and found a new type of collective motion, the chiral electric wave (CEW), arising from CESE and propagating in parallel/antiparallel to the field. In symmetric collisions there should be no net electric field on average, but asymmetric collisions like Cu+Au could provide a test ground for the CEW measurements.

5.1. Electric Quadrupole Observable

The CMW will induce a finite electric quadrupole moment of the collision system, with additional positive (negative) charge at the “poles” (“equator”) of the produced fireball [11]. This electric quadrupole, if boosted by radial flow, will lead to charge-dependent elliptic flow. Let us take pions as an example: on top of the baseline , the CMW will lead to [11] where , , and are the quadrupole moment, the net charge density, and the charge asymmetry of the collision event, respectively. is always positive, so the -integrated of () should be above (below) the baseline owing to the CMW. However, the baseline may be different for and because of several other physics mechanisms [85, 86]. Therefore, it is less ambiguous to study the CMW via the dependence of pion than via the -integrated .

Figure 14 [80] takes 30–40% 200 GeV Au+Au, for example, and shows pion as a function of in panel (a).   increases with , whereas decreases with a similar magnitude of the slope. Note that was integrated over a narrow low- range (0.15 < < 0.5 GeV/) to focus on the soft physics of the CMW. This selection also ensures that the is independent of and is the same for and , so that the splitting is not a trivial effect. This splitting was also confirmed by ALICE results for Pb+Pb collisions at 2.76 TeV [82]. The difference between and is fitted with a straight line in panel (b). The slope parameter , or presumably / from (10), is positive and qualitatively consistent with the expectation of the CMW picture. The fit function is nonzero at , indicating the -integrated for and are different, which was also observed in [87].

STAR followed the same procedure as shown in Figure 14 to retrieve the slope parameter as a function of centrality for Au+Au collisions at 200, 62.4, 39, 27, 19.6, 11.5, and 7.7 GeV, as shown in Figure 15 [80]. A similar rise-and-fall trend is observed in the centrality dependence of the slope parameter for all the beam energies except 11.5 and 7.7 GeV, where the slopes are consistent with zero with large statistical uncertainties. It was argued [85] that at lower beam energies the -integrated difference between particles and antiparticles can be explained by the effect of quark transport from the projectile nucleons to mid-rapidity, assuming that the coalescence mechanism still holds and that of transported quarks is larger than that of produced ones. The same model, however, suggested a negative slope for as a function of [88], which is contradicted by the data.

As was done in the CME case, contributions from conventional physics were checked via the same analysis of the Monte Carlo events from UrQMD [60]. For Au+Au collisions at 200 GeV, the slopes extracted from UrQMD events are consistent with zero for the 10–70% centrality range, where the signal from the real data is prominent. Similarly, the AMPT event generator [70, 71] also yields slopes consistent with zero (not shown here). On the other hand, the simplified CMW calculations [81] demonstrate a centrality dependence of the slope parameter similar to the data. Recently a more realistic implementation of the CMW [89] confirmed that the CMW contribution to is sizable, and the centrality dependence of is qualitatively similar to the data. A quantitative comparison between data and theory requires further work on both sides. For example, the STAR experiment only measured in a pseudorapidity slice () of an event, instead of that of the whole collision system. for these two cases are expected to be proportional to each other, but the determination of the ratio will be model dependent.

Figure 16 shows a comparison in the slope parameter between STAR results for 200 GeV Au+Au [80] and ALICE results for 2.76 TeV Pb+Pb [82]. Overall, the slopes are surprisingly similar when considering the different collision energies and multiplicities, as well as the different kinematic acceptance: the STAR data estimated for charged pions with 0.15 < < 0.5 GeV/ and , while the ALICE data are for unidentified hadrons with 0.2 < < 5 GeV/ and . Note that, in central and peripheral collisions, the STAR slopes at 200 GeV are consistent with zero, while the ALICE slopes at 2.76 TeV still remain positive-definite. This difference leaves room for potential backgrounds in the measurements.

The measurement of requires a correction to the observed because of the finite detector tracking efficiency. A novel correlator [90] that is less dependent on efficiency was proposed:Here and are the azimuthal angles of particles 1 and 2, and is the charge (±1) of particle 3. The single brackets represent the average over particles and events, and the double bracket denotes the cumulant. In the absence of charge-dependent correlations, this correlator should be equal to zero. When the charge of the third particle is averaged over all particles in a specified kinematic acceptance of an event, the mean is equal to the charge asymmetry; that is, .

The three-particle correlator (as in (11)) for the 2nd harmonic was measured by ALICE and is presented in Figure 17 as a function of centrality in Pb+Pb collisions at 2.76 TeV [82]. The correlation strength substantially increases in more peripheral collisions, which can be caused by a combination of several factors. For example, the magnetic field strength increases with the impact parameter and causes the stronger correlations due to the CMW. In addition, the LCC effect could play a role [90], together with other unknown sources of correlation, which is a similar case as the slope parameter discussed above.

5.2. Possible Backgrounds

It was argued [91] that LCC at freeze-out, when convoluted with the characteristic shape of and , may provide a qualitative explanation for the finite slope observed from data. A realistic estimate of the contribution of this mechanism turns out to be smaller than the measurement by an order of magnitude [80]. Reference [91] also proposes a test with the measurement, and the corresponding slope parameters for were reported by STAR to be consistent with zero [92], which further suggests the smallness of this effect. ALICE measured the three-particle correlator multiplied by for the 3rd harmonic and the 4th harmonic [82]. In both cases, the centrality dependence of the charge dependence is flat, in contrast to the 2nd harmonic that has a significant centrality dependence. This may indicate a different nature of the correlation or reflect a weaker centrality dependence of compared with that of . Future measurements of these higher harmonics with better precision will shed light on the true origin of this correlator.

Recently a hydrodynamic study [93] claimed that simple viscous transport of charges, combined with certain specific initial conditions, will lead to a sizable splitting of charged pions. To have their results of pion splitting resemble data, the authors had to invoke a crucial relation between isospin chemical potential and the electric charge asymmetry, which needs to be verified. Furthermore, certain predictions of this model (e.g., splitting for kaons) appear to be contradicted by experimental data [94]. The validity of such an idea needs to be thoroughly vetted by realistic viscous hydrodynamic simulations. Nevertheless, this study raises an important point: the normal (viscous hydrodynamical) transport of charges should be quantitatively understood to make a firm case for the observation of anomalous charge transport via the CMW.

6. Future Measurements

The confirmation of the experimental observation of several chiral anomalous effects will bring forth an exciting program to directly study the nonperturbative sector of QCD. Future experimental measurements should aim at more detailed study of the observed signals as well as understanding the background effects. Previous sections have covered a few such topics: initial magnetic field and vorticity, correlations with identified particles, higher-harmonic correlations, BES-II, and U+U collisions. In the following, we will focus on the event shape engineering (ESE) and isobaric collisions.

6.1. Event Shape Engineering

Flow-related backgrounds could be potentially removed via ESE [75, 95], with which spherical events or subevents are selected, so that the particles of interest therein carry zero . A previous attempt was made with the charge separation observable of CMAC (roughly equivalent to γ), as a function of event-by-event “observed ” [63]. However, there are several issues in this approach that prevent a clear interpretation of the result. Reference [51] studied the flow vector of the subevent of interest, :and found that is a good handle on event shape.

Figure 18 shows the event plane resolution () for the subevents and (a) and the true elliptic flow and the corrected as functions of (b), from AMPT simulations of 20–60% Au+Au collisions at = 200 GeV [51]. Each AMPT event has been divided into three subevents according to pseudorapidity, : subevent contains particles of interest with , and subevent () serves as a subevent plane using particles with (). Flow fluctuation causes a positive correlation in flow between subevents in the same event, and as a result, for subevent () increases with for subevent . Figure 18(b) displays a discrepancy between and the corrected , owing to the difference between the reaction plane and the participant plane [96], in terms of nonflow and flow fluctuation. What matters more is the fact that both values decrease with and drop to (), which demonstrates ’s capability of selecting spherical subevents in the second harmonic.

Figure 19(a) presents the γ correlators multiplied by the number of participating nucleons, , as functions of , for 20–60% AMPT events of Au+Au collisions at 200 GeV [51]. For both the same-charge and the opposite-charge correlators, the true and the corrected are consistent with each other within the statistical uncertainties. This indicates that, compared with , γ is less sensitive to nonflow or flow fluctuation. At larger , the opposite-charge correlators are above the same-charge correlators, suggesting a finite flow-related background. The opposite- and same-charge correlators converge at small . Figure 19(b) shows versus , and again, the two observables seem to coincide. Linear fits to both observables yield small intercepts that are consistent with zero. The finite values in AMPT events are solely due to background contributions, so the disappearance of background is demonstrated when the “correctable” observable () is projected to zero . Reference [51] has designed a promising recipe for future measurements to effectively remove flow backgrounds and restore the ensemble average of the CME signal.

6.2. Isobaric Collisions

To disentangle the possible CME signal and the flow-related backgrounds, one can utilize experimental setups to either vary the backgrounds with the signal fixed or vary the signal with the backgrounds fixed. The former approach was carried out by exploiting the prolate shape of the Uranium nuclei [75]. However, it was found that the total multiplicity of detected hadrons is far less dependent on the number of binary collisions than expected [76], so it is very hard to isolate tip-tip collisions (that generate small ) from body-body collisions (that generate large ). This significantly reduces the lever arm available to manipulate in order to separate flow backgrounds from the CME.

The latter approach (with the -driven backgrounds fixed) can be realized, especially for mid-central or mid-peripheral events, with collisions of isobaric nuclei, such as and [75]. Ru+Ru and Zr+Zr collisions at the same beam energy are almost identical in terms of particle production [83], while the charge difference between Ru and Zr nuclei provides a handle on the initial magnetic field. Our current knowledge of the deformity () of Ru and Zr is incomplete: scattering experiments (case 1) [97, 98] state that Ru is more deformed () than Zr (), while comprehensive model deductions (case 2) [99] tell the opposite that is smaller than . This systematic uncertainty has different impacts on the signal (via the initial magnetic field) and the background (via the initial eccentricity) to be discussed later. As a by-product, measurements in central collisions will discern which information source is more reliable regarding the deformity of the Zr and Ru nuclei.

Figure 20(a) presents the theoretical calculation [83] of the initial magnetic field squared with correction from azimuthal fluctuation of the magnetic field orientation, (with being the pion mass and the azimuthal angle of the magnetic field), for the two collision systems at 200 GeV, using the HIJING model [19, 100]. quantifies the magnetic field’s capability of driving the CME signal in the γ correlator. For the same centrality bin, the Ru+Ru collision produces a significantly stronger magnetic field than Zr+Zr. Figure 20(b) shows that the relative difference in between Ru+Ru and Zr+Zr collisions is approaching 15% (case 1) or 18% (case 2) for peripheral events and reduces to about 13% (both cases) for central events. Figure 20(b) shows the relative difference in the initial eccentricity, , obtained from the Monte Carlo Glauber simulation. is highly consistent with 0 for peripheral events and goes above (below) 0 for the parameter set of case 1 (case 2) in central collisions, because the Ru (Zr) nucleus is more deformed. The relative difference in should closely follow that in eccentricity, so for the centrality range of interest, 20–60%, the -related backgrounds stay almost the same for Ru+Ru and Zr+Zr collisions. Reference [83] further carried out the projection for the γ measurements in Ru+Ru and Zr+Zr at 200 GeV (400 million events for each collision type) and concluded that a significance can be achieved for the relative difference in the observable between the two collision systems, assuming the flow backgrounds take up to two-thirds of the observable. The results strongly suggest that the isobaric collisions can serve as an ideal tool to disentangle the signal of the chiral magnetic effect from the -driven backgrounds. The isobaric collisions may also be used to disentangle the signal of the CMW from background effects.

7. Summary

The physics of anomalous transport is an underexplored territory of QCD as a non-Abelian gauge theory. The interplay of quantum anomalies with a strong magnetic field or vorticity enables various novel transport phenomena in chiral systems. In heavy-ion collisions, these phenomena provide a unique probe to the topological properties of the QGP and manifest the charge dependence of the azimuthal distributions of the produced hadrons. In this article, we reviewed the experimental data from RHIC at BNL and the LHC at CERN and provided an evidence for the predicted CME, CVE, and CMW, with magnitudes similar to the theoretical estimates. Although all of these experimental observables contain known conventional backgrounds, no compelling alternative explanation at present can describe all of the data without invoking the anomalous chiral effects. Nevertheless, with still miles to go in both experiment and theory to substantiate the existing evidence, we outlined a few such programs that hopefully will be accomplished in the near future.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Huan Huang and other members of the UCLA Heavy-Ion Physics Group for discussions. This work is supported by a Grant (no. DE-FG02-88ER40424) from US Department of Energy, Office of Nuclear Physics.