Quantifying the Chiral Magnetic Effect from Anomalous-Viscous Fluid Dynamics

In this contribution we report a recently developed Anomalous-Viscous Fluid Dynamics (AVFD) framework, which simulates the evolution of fermion currents in QGP on top of the bulk expansion from data-validated VISHNU hydrodynamics. With reasonable estimates of initial conditions and magnetic field lifetime, the predicted CME signal is quantitatively consistent with change separation measurements in 200GeV Au-Au collisions at RHIC. We further develop the event-by-event AVFD simulations that allow direct evaluation of two-particle correlations arising from CME signal as well as the non-CME backgrounds. Finally we report predictions from AVFD simulations for the upcoming isobaric (Ru-Ru v.s. Zr-Zr ) collisions that could provide the critical test of the CME in heavy ion collisions.


Introduction -Chiral Magnetic Effect & Anomalous-Viscous Fluid Dynamics
The Chiral Magnetic Effect (CME) [1,2] refers to the generation of an electric current J Q along the magnetic field B applied to a system of chiral fermions with chirality imbalance, i.e.
where σ 5 = C A µ 5 is the chiral magnetic conductivity, with the chiral chemical potential µ 5 that quantifies the imbalance between fermion densities of opposite (right-handed, RH versus left-handed, LH) chirality.
Given the magnificent physics embodied in the Chiral Magnetic Effect, it is of the utmost interest to search for its manifestation in the quark-gluon plasma (QGP) created in relativistic heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC).Dedicated searches for potential CME signals have been ongoing at RHIC and the LHC [3,4,5,6], with encouraging evidences 1 Presenter arXiv:1704.05531v1[nucl-th] 14 Apr 2017 reported through measuring the charge separation signal induced by the CME current (1).The interpretation of these data however suffers from backgrounds arising from the complicated environment in a heavy ion collision (see e.g.[2,7,8]).Currently the most pressing challenge for the search of CME in heavy ion collisions is to clearly separate background contributions from the desired signal.A mandatory and critically needed step, is to develop state-of-the-art modeling tools that can quantify CME contribution in a realistic heavy ion collision environment.
To address this challenge, we've recently developed a simulation framework, the Anomalous-Viscous Fluid Dynamics (AVFD) [9], focusing on describing anomalous chiral transport in heavy ion collisions at high beam energy (such as the top energy RHIC collisions).The bulk evolution in such collisions is well described by boost-invariant 2+1D 2nd-order viscous hydrodynamics (e.g.VISHNU simulations [10]) where net charge densities are small enough and typically neglected without much influence on bulk evolution.However to study the CME, one needs to accurately account for the evolution of fermion currents.Our approach is to solve the following fluid dynamical equations for the chiral fermion currents (RH and LH currents for u and d flavors respectively) as perturbations on top of the bulk fluid evolution: where χ = ±1 labels chirality for RH/LH currents and f = u, d labels light quark flavor with electric charge Q f and color factor N c = 3.The E µ = F µν u ν and B µ = 1 2 µναβ u ν F αβ are external electromagnetic fields in fluid rest frame.The derivative Dµ is covariant derivative and d = u µ Dµ , with projection operator ∆ µν = (g µν − u µ u ν ).Viscous parameters σ and τ r are diffusion constant and relaxation time respectively.In the above equations the fluid four-velocity field u µ , temperature T as well as all other thermodynamic quantities are determined by background bulk flow.Furthermore the (small) fermion densities n χ, f and corresponding chemical potential µ χ, f are related by lattice-computed quark number susceptibilities c f 2 (T ).

Results and Discussion
As the key ingredients for driving the CME current, the magnetic field and the initial axial charge density are the most important inputs for the AVFD simulation.For the magnetic field B = B(τ)ŷ (with ŷ the eventwise out-of-plane direction), we use a plausible parametrization (see e.g.[11]) as B(τ) = B 0 1+(τ/τ B ) 2 .The peak value B 0 (for each centrality) at the collision point has been well quantified with event-by-event simulations and we use the most realistic values from [12].For the lifetime of the B field we use a reasonable estimate of τ B = 0.6 fm/c which is comparable to the onset time of hydrodynamic evolution.For the initial axial charge density arising from gluonic topological charge fluctuations, one could make the following estimate based on the strong chromo-electromagnetic fields in the early-stage glasma similarly to the recent study in [13]: In the above ρ tube 1fm is the transverse extension of glasma flux tube, A overlap is the geometric overlapping area of the two colliding nuclei, and N coll.the binary collision number for a given centrality.Such axial charge density depends most sensitively upon the saturation scale Q s , with a reasonable range of Q 2 s 1 ∼ 1.5GeV 2 for RHIC 200AGeV collisions.After the preceding discussions on the various aspects of the AVFD tool, let us now proceed to utilize this tool for quantifying CME signal to be compared with available data.The measurement of a CME-induced charge separation is however tricky, as this dipole flips its sign from event to event depending the sign of the initial axial charge arising from fluctuations, thus with a vanishing event-averaged mean value.What can be measured is its variance, through azimuthal correlations for same-sign (SS) and opposite-sign (OS) pairs of charged hadrons.The so-called γ S S /OS ≡ cos(φ 1 + φ 2 ) observables measure a difference between the in-plane versus out-of-plane correlations and are indeed sensitive to potential CME contributions.They however suffer from considerable flow-driven background contributions that are not related to CME (see e.g.[2,8]).One plausible approach to separate background and CME signal is based on a two-component scenario [8], which was recently adopted by the STAR Collaboration to suppress backgrounds and extract the flow-independent part (referred to as H S S /OS ) [4].We consider H S S /OS as our "best guess" thus far for potential CME signal to be compared with AVFD computations.Specifically a pure CME-induced charge separation will contribute as (H S S − H OS ) → 2 a ch 1 2 .The AVFD results for various centrality bins are presented in Fig. 1, with the green band spanning the range of key parameter Q 2 s in 1 ∼ 1.5GeV 2 reflecting uncertainty in estimating the initial axial charge.Clearly the CME-induced correlation is very sensitive to the amount of initial axial charge density as controlled by Q 2 s .The comparison with STAR data [4] shows very good agreement for the magnitude and centrality trend for choices with relatively large values of Q 2 s .Quantitative predictions from Anomalous-Viscous Fluid Dynamics simulations for the CMEinduced H-correlations, in comparison with STAR measurements [4].The green bands reflect current theoretical uncertainty in the estimates of initial axial charge generated by gluonic field topological fluctuations.
While the AVFD simulations based on smooth average hydro profile provide a quantitative account of the "pure" CME signal, it is imperative to perform event-by-event simulations that could allow direction evolution of two-particle correlations and thus ultimate comparison with experimental measurements.We've recently implemented such simulations with event-wise fluctuating initial conditions as well as a hadronic stage after hadronization using UrQMD to include the effects of hadron cascade and resonance decay.To demonstrate the CME-driven correlations, we run simulations for 50 − 60% 200GeV Au-Au collisions with three different setups: (a) a "null" case with no magnetic field or chiral imbalance; (b) eB = 5m 2 π and n A /s = 0.1 (corresponding to estimates with Q 2 s 1GeV 2 ); (c) eB = 5m 2 π and n A /s = 0.2 (corresponding to estimates with Q 2 s 1.5GeV 2 ).For each setup, we accumulated ∼ 10 7 events and measured the correlator γ OS −S S versus different event-wise v 2 bins, shown in Fig. 2 (left).An obvious linear dependence γ OS −S S on v 2 is observed for all these three cases.The extracted slope and intercept of such linear dependence are shown in Fig. 2 (right): while the slope is basically independent of initial axial charge (indicating its dominant non-CME origin), the intercept grows quadratically with initial axial charge as expected from the CME signal.Finally we report the AVFD predictions for the isobaric collisions planned at RHIC [14], which provide a unique means to unambiguously decipher CME from backgrounds.In such "contrast" colliding systems (specifically for 96 40 Zr-96 40 Zr versus 96 44 Ru-96 44 Ru), the bulk particle production and background correlations are expected to be basically identical while the CME signals are expected to differ by about 20% due to their different magnetic field strength by about 10%.This will be a crucial test for the search of CME, and quantitative predictions are important.In Fig. 3 we show the AVFD results for the correlations in isobar collisions.The left panel shows our prediction for pure CME signal in these two systems.By assuming the background correlation strength is mainly dependent on multiplicity, we extrapolate the departed background strength from Au-Au to the isobaric collisions (shown in middle panel).Folding CME signal with such background estimates we obtain the predictions for γ-correlations, with a visible ∼ 15% shift in peripheral region differentiating the two systems.Provided current experimental uncertainty and limitations [15], the difference at such a level should be measurable with adequate number of events from isobaric collisions.

Summary
In summary, a new simulation tool -the Anomalous-Viscous Fluid Dynamics (AVFD) framework has been developed for quantifying the charge separation signal induced by Chiral Magnetic Effect in relativistic heavy ion collisions.We find that, subject to current theoretical and experimental uncertainties, the AVFDpredicted CME signal with realistic initial conditions and magnetic field lifetime is quantitatively consistent with measurements from 200AGeV AuAu collisions at RHIC.Also, by event-by-event simulation, we find that the intercept of the v 2 dependence of the correlator γ S S −OS is directly sensitive to the CME while the slope is likely to be dominated by background contributions.Finally, we make predictions for the upcoming isobaric collisions that would be a critical test for the search of CME in heavy ion collisions.
Fig. 1.Quantitative predictions from Anomalous-Viscous Fluid Dynamics simulations for the CMEinduced H-correlations, in comparison with STAR measurements[4].The green bands reflect current theoretical uncertainty in the estimates of initial axial charge generated by gluonic field topological fluctuations.

14 + 2 Fig. 2 .
Fig. 2. (color online) (left) The v 2 dependence of γ OS −S S ; (right) The slope and intercept extracted from the v 2 dependence of γ OS −S S .

Fig. 3 .
Fig. 3. (color online) (left) AVFD predictions for CME-induced H-correlations in isobar collisions; (middle) Extrapolation of background correlations from AuAu to ZrZr and RuRu systems; (right) Predicted γ-correlations in ZrZr and RuRu collisions by folding together computed CME signals and estimated background correlations.