The generation of matter-antimatter asymmetries and hypermagnetic fields by the chiral vortical effect of transient fluctuations

We study the contribution of temperature-dependent chiral vortical effect to the generation and evolution of the hypermagnetic fields and the matter-antimatter asymmetries, in the symmetric phase of the early Universe, in the temperature range $100\mbox{GeV} \le T\le 10\mbox{TeV}$. Our most important result is that, due to the chiral vortical effect, small overlapping transient fluctuations in the vorticity field in the plasma and temperature of matter degrees of freedom can lead to the generation of strong hypermagnetic fields and matter-antimatter asymmetries, all starting from zero initial values. We show that, either an increase in the amplitudes of the fluctuations of vorticity or temperature, or a decrease in their widths, leads to the production of stronger hypermagnetic fields, and therefore, larger matter-antimatter asymmetries. We have the interesting result that fluctuating vorticity fields are more productive, by many orders of magnitude, as compared to vorticities that are constant in time.


Introduction
Anomalous transport effects play important roles in particle physics and cosmology, particularly in the early Universe [1]. One important effect of this kind is the so-called Chiral Vortical Effect (CVE), which refers to the generation of an electric current parallel to the vorticity field in the chiral plasma [2]. This effect was discovered by Vilenkin who showed that a neutrino current density can result from a rotating black hole [2]. He obtained the neutrino current density in the direction of the rotation axis as where Ω is the angular velocity, µ is the chiral chemical potential of the neutrino, and T is its temperature. Thirty years after its discovery, the CVE appeared in the relativistic hydrodynamic equations as an interesting manifestation of anomalies in quantum field theory [3]. This effect has attracted much attention and has been investigated extensively in recent years, leading to a deeper understanding of the subject [4][5][6][7][8][9][10]. In particular, it has been established that in single species chiral plasma in the broken phase at high temperatures, the CVE shows up in the vector current as J cv = 1 4π 2 (µ 2 R − µ 2 L ) Ω, and in the axial current as J 5 cv = T 2 6 + 1 4π 2 (µ 2 R + µ 2 L ) Ω, where µ R and µ L are the right-handed and the left-handed chemical potentials of the species, respectively [2][3][4][5][6][7][8][9][10] 1 . Interestingly, the term proportional to T 2 indicates that there can be an axial current, even if µ R = µ L = 0.
In this study, we present the correct form of the chiral vortical current in the symmetric phase. Then, we show the prominent effects of the temperature-dependent part of this current in the symmetric phase of the early Universe close to the electroweak phase transition (EWPT). In particular, we show that even very small, but overlapping, transient fluctuations in the vorticity field and temperature of matter degrees of freedom can have important consequences, including the generation of hypermagnetic field in the absence of initial matter asymmetries. The vorticity fluctuations that we consider are about the zero background value, while the temperature fluctuations are about the finite equilibrium temperature of the plasma. The most important role of the CVE in this context is to produce the magnetic fields, either through the chiralities, or through the temperature fluctuation, the latter of which is the main focus of this work. Henceforth, we shall refer to transient fluctuations, which we take to be in the form of short pulses, simply as fluctuations.
Another anomalous transport effect is the chiral magnetic effect (CME), which refers to the generation of an electric current parallel to the magnetic field in the imbalanced chiral plasma [11][12][13][14]. It is known that, in single species chiral plasma in the broken phase at high temperatures, the CME appears in the vector current as J cm = Q 4π 2 (µ R − µ L ) B, and in the axial current as J 5 cm = Q 4π 2 (µ R + µ L ) B, where Q is the electric charge [10,11,15,16] of the species. 1 Later, we will present the correct form of the chiral magnetic current, in the symmetric phase. The chiral magnetic current originating from the electroweak Abelian anomaly, and the chiral vortical current are both non-dissipative currents which can strongly affect the generation and the evolution of the magnetic fields and the matter-antimatter asymmetries in the early Universe [17][18][19][20][21].
Observations clearly show that our Universe is magnetized on all scales [22,23]. Various models have been proposed to explain the origin of these magnetic fields [20,[24][25][26][27][28][29][30], among which the one relying on the electroweak Abelian anomaly has attracted much attention and has been considerably investigated [17][18][19]31]. There exists a relationship between the generation and the evolution of the hypermagnetic fields and the fermion number densities in this model, which is due to the chiral coupling of the hypercharge gauge fields to the fermions before the EWPT [27]. In case there is a preexisting asymmetry of the right-handed electrons, their number density is almost conserved far from the EWPT, i.e. T > 10TeV, due to their tiny Yukawa coupling. For lower temperatures, this asymmetry can be converted to the hypermagnetic helicity according to the Abelian anomaly equation, 17-19, 27, 31, 32]. The anomaly equation shows that, in a reverse process, a strong helical hypermagnetic field can generate the matter-antimatter asymmetries in the Universe, as well [33][34][35].
Another challenge in particle physics and cosmology is the excess of matter over antimatter, with the measured baryon asymmetry of the Universe being of the order of η B ∼ 10 −10 [36][37][38].
The three Sakharov conditions 2 should be satisfied in any CPT invariant model used to explain this asymmetry from an initially symmetric Universe [39].
In previous studies based on the electroweak Abelian anomalous model, it has been assumed that there is either a significant amount of matter-antimatter asymmetries to produce the hypermagnetic field, or a strong hypermagnetic field to produce the matter-antimatter asymmetries. The most important result of this study is that the matter-antimatter asymmetries and the hypermagnetic field can all be generated simultaneously from zero initial values, by considering the temperature-dependent CVE before the EWPT. To obtain this interesting result, we take into account the effects of the temperature-dependent term of the chiral vortical current on the evolution of the hypermagnetic fields and the matter-antimatter asymmetries, by considering simultaneous small fluctuations, about the background values, in temperature of the right-handed electrons and the vorticity field, close to the EWPT. We also show that fluctuations in the vorticity field are much more productive than vorticity fields that are constant in time.
To be more precise, sharp fluctuations yield results comparable to constant vorticities whose amplitudes are many orders of magnitude larger.
As mentioned above, an underlying assumption of this work is the presence of fluctuations in the plasma of the early Universe, containing all elementary particles and gauge fields. A general description of this plasma in the context of hydrodynamics, as an effective field theory, has been presented [40,41]. In such a physical system, hydrodynamic variables naturally fluctuate around their statistical averages. Such stochastic and frequent fluctuations are a commonplace in any plasma, including that of the early Universe. These fluctuations, which may stem from a sum of weakly-correlated random local events, can occur for physical observables such as the number density, velocity, and temperature [40,[42][43][44]. Even though the gauge interactions are strong and their rates are fast in this epoch, the extreme temperatures of the primordial plasma imply that fluctuations could still occur for all species. Moreover, owing to the different Yukawa and gauge couplings, and in particular the chiral nature of the electroweak sector, different chiral fermions could, in principle, experience different fluctuations. All that is required here is that the fluctuations for at least one of the species be non-identical to the rest, at least once.
Another justification for our hypothesis is the following. Realistic fluctuations of density, temperature and vorticity are usually local in space and time. Meanwhile, in a multi-component plasma local and independent density fluctuations can occur for different matter components. The concurrent occurrence of density, temperature and vorticity fluctuations leads to the possibility of occurrence of non-identical fluctuations for different species of particles. To be concrete, consider a local density fluctuation leading to a local excess of a particle species, e.g. e R , in the central region and rarefaction in the surrounding region, followed immediately by a temperature fluctuation which increases the temperature of the central region. The concurrent occurrence of these two fluctuations leads to a temporary local, and to a much lesser degree global average, increase of temperature of that particle species relative to quasi-equilibrium temperature of the plasma. All that is required now is an overlapping vorticity fluctuation to have the necessary conditions for our model. In this work, we make the simplifying assumption, as is usually done, of considering these fluctuations to be global rather than local. It should be noted that the amplitude and duration of these fluctuations are usually small enough such that the equilibrium energy density, pressure and entropy of the system can be defined. Nevertheless, as we shall show, even very small and brief fluctuations can trigger the mechanism that we propose here and affect significantly the dynamical evolution of the system, as is formulated within the equations of anomalous magnetohydrodynamics (AMHD). This paper is organized as follows: In Sec. 2, the anomalous magnetohydrodynamics equations and the evolution equations for the matter-antimatter asymmetries are derived in the ex-panding Universe. In Sec. 3, the set of coupled differential equations are solved numerically. In Sec. 4, the results are summarized and the conclusion is presented.

Anomalous Magnetohydrodynamics Equations
In this section, the anomalous magnetohydrodynamics (AMHD) equations are obtained in the Landau-Lifshitz frame in the symmetric phase of the expanding Universe. Taking the CVE and the CME into account, the Maxwell's equations for the hypercharge-neutral plasma in the expanding Universe are given as [3,20,[45][46][47][48][49][50][51][52] where R is the scale factor, H =Ṙ/R is the Hubble parameter, σ is the electrical hyperconductivity of the plasma, and v and ω = 1 R ∇ × v are the bulk velocity and vorticity of the plasma, respectively. Furthermore, the chiral vorticity and helicity coefficients c v and c B are as follows (see Appendix A) 3 , where n G is the number of generations, and N c = 3 and N w = 2 are the ranks of the non-Abelian SU(3) and SU(2) gauge groups, respectively. Moreover, µ L i (µ R i ), µ Q i , and µ u Ri (µ d Ri ) are the common chemical potentials of left-handed (right-handed) leptons, left-handed quarks with different colors, and up (down) right-handed quarks with different colors, respectively. Furthermore, 'i' is the generation index, and the relevant hypercharges are After substituting the hypercharges in Eqs. (2.7) and (2.8), we obtain Let us make the same assumptions as in our previous studies, and simplify c v and c B accordingly [18,20]. We assume that all quark Yukawa processes are in equilibrium and, because of the flavor mixing in the quark sector, all up or down quarks belonging to different generations with distinct handedness have the same chemical potential [18,53]. For simplicity, we also assume that the Higgs asymmetry is zero and obtain [18,54] (2.12) Furthermore, we assume that only the contributions of the baryonic and the first-generation leptonic chemical potentials to c v and c B are significant. As for the temperature fluctuations, it suffices to consider fluctuations in only one of the matter components, which we take to be e R . Using Eq. (2.12) and the aforementioned assumptions, we simplify Eqs. (2.10) and (2.11) to obtain where µ B = 12µ Q , and ∆T 2 = T 2 e R − T 2 is the temperature fluctuation, and T is the equilibrium temperature of the thermal bath, which includes all other components of the plasma. We set is an arbitrary profile function to be specified later, and x(T ) = t(T )/t EW = (T EW /T ) 2 is given by the Friedmann law.
In the following, we choose a simple monochromatic Chern-Simons configuration for the hypermagnetic field B Y = (1/R) ∇ × A Y , and the velocity field v = (1/R) ∇ × S [18,20]. To do this, we choose A Y = γ(t) (cos kz, sin kz, 0), and S = r(t) (cos kz, sin kz, 0), for their corresponding vector potentials [56][57][58]. Note that we have chosen a fully helical form with the negative helicity for both, the reason for which will be stated later. Let us now obtain the evolution equation for the velocity field. Neglecting the displacement current in the lab frame in Eq.
where the kinematic viscosity ν ≃ 1/(5α 2 Y T ) [60,61]. Neglecting the displacement current in the lab frame and using the aforementioned configurations, the hyperelectric field and the evolution equation for the hypermagnetic field are obtained, as follows: where ω = −k ′ v, σ = 100T , angle brackets denote spatial average, and k ′ = k/R = kT . The latter shows the increase of the hypermagnetic length scale, due to the expansion of the Universe. Note that with the choice of vector potentials for v(t) and B Y (t), the advection term v × B Y has been set to zero in the above, and we have the following simplification: In the next subsection we obtain the evolution equations for the matterantimatter asymmetries.

Evolution equations for the matter-antimatter asymmetries
Before the EWPT, the gauge fields of U Y (1) couple to the fermions chirally, in contrast to those of U em (1) in the broken phase, leading to the non-conservation of the matter currents. This shows up in the Abelian anomaly equations [62], which, for the first-generation leptons, are where ∇ µ is the covariant derivative with respect to Friedmann-Robertson-Walker (FRW) metric ds 2 = dt 2 − R 2 (t)δ ij dx i dx j , t is the physical time, and x i s are the comoving coordinates. Integrating the above equations over all space and considering the perturbative chirality flip reactions for the leptons, we obtain (see Refs. [17,18,20,63] and Appendix B for details), In the equations above, η f = (n f /s) with f = e R , e L , ν L e is the fermion asymmetry and n f is the charge density of the f th species of fermion, s = 2π 2 g * T 3 /45 is the entropy density and g * = 106.75 is the effective number of relativistic degrees of freedom, where M Pl is the Plank mass. Furthermore, the term Γ 0 appearing in the equations is the chirality flip rate of the right-handed electrons. In a similar manner, the evolution equation for the baryon asymmetry can also be obtained as (see Appendix B) Using 1Gauss ≃ 2 × 10 −20 GeV 2 , and setting the kinematic viscosity ν to zero for simplicity, we obtain the complete set of evolution equations for the matter-antimatter asymmetries and the amplitudes of the hypermagnetic and velocity fields as where M = 2π 2 g * /45, and the coefficients C i , i = 1, ..., 8, are 29) and α Y = g ′2 /4π ≃ 0.01 is the fine-structure constant for the U Y (1). We now choose the profile of temperature fluctuation β[x(T )] = ∆T 2 /T 2 , as defined in Eq. (2.13) and the paragraph below it, to be a Gaussian function of x: where β 0 is the amplitude multiplying the normalized Gaussian distribution, and x = (t/t EW ) = (T EW /T ) 2 , as defined before. The profile of the vorticity fluctuation must have an overlap with that of temperature fluctuation, in order to produce any effect. For simplicity, we choose the two profiles to be identical. That is, where v 0 is the amplitude of the velocity fluctuation. We should mention that the occurrence of any fluctuation in a plasma in a quasi-equilibrium state would normally trigger a restoring response originating from dissipative effects, such as viscous effects. Here, for simplicity we assume that the combined results of the original fluctuations and the ensuing dissipative effects have the profiles given by Eqs. (2.30,2.31).
The majority of analysis presented in the following section is for a single pulse for temperature and vorticity, as stated above, which, as we shall show, produce matter-antimatter asymmetries and helical hypermagnetic field. However, at the end of the next section we present the results for two sets of successive pulses, the latter one having the same temperature profile but with negative amplitude, showing that the first pulse is the main determinant of the outcome.

Numerical Solution
In this section, we obtain the numerical solutions of the evolution equations. As mentioned earlier, we investigate the effects of the temperature fluctuations of right-handed electrons, in the presence of vorticity, on the generation and evolution of the hypermagnetic field and the matter-antimatter asymmetries, in the temperature range 100 GeV ≤ T ≤ 10 TeV. We consider the temperature fluctuations as small Gaussian distributions in x, as shown in Eq. (2.30), that occur close to the EWPT. As for the vorticity field, we consider small fluctuations, whose profiles, as shown in Eq. (2.31), coincide with those of the temperature fluctuations. We also investigate the cases with constant vorticity fields for comparison. As we shall show, the former is much more interesting and will be the focus of our work, since it is not only physically more realistic, but also could yield orders of magnitude larger results for the asymmetries and the hypermagnetic field. In the following, we solve the evolution equations by considering the comoving wave number as k = 10 −7 , and setting the initial values of the hypermagnetic field amplitude and the matter-antimatter asymmetries to zero, i.e. B For our first case, we solve the coupled differential equations with the initial conditions, v 0 = 10 −5 , b = 2 × 10 −4 , and x 0 = 45 × 10 −5 , for various values of the amplitude of temperature fluctuations β 0 , and present the results in Fig. 1. As can be seen in the figure, the simultaneous occurrence of small vorticity fluctuation and temperature fluctuation for the right handed electrons leads to the generation of strong hypermagnetic fields which then produce the matter-antimatter asymmetries, all starting from zero initial values. It can be seen that by increasing the amplitude of the temperature fluctuation, the maximum and the final values of the hypermagnetic field amplitude, as well as the matter-antimatter asymmetries, increase. We have found that in our model signs of the matter-antimatter asymmetries produced and the helicity of hypermagnetic and vorticity fields, are always opposite. This is a manifestation of the generalized charge conservation as stated in Eq. (6.9) of Appendix B. 5 The Chern-Simons configuration that we have chosen for the hypermagnetic and vorticity fields has negative helicity. Figure 1 also shows that the hypermagnetic field amplitude grows to its maximum value of about 10 21 G, then decreases due to the expansion of the Universe. We have also investigated the effects of changing the amplitude of the vorticity fluctuation, and have found similar results. For our second case we solve the set of evolution equations with the initial conditions, v 0 = 10 −5 , β 0 = 5 × 10 −4 , and x 0 = 45 × 10 −5 , for various values of the width or duration of both fluctuations b, and show the results in Fig. 2. As can be seen, by decreasing the width of the Gaussian function, the maximum and the final values of the hypermagnetic field amplitude, and the baryon asymmetry increase. 6 For our third case, we solve the coupled equations with the initial conditions, v 0 = 10 −5 , b = 2 × 10 −4 , and β 0 = 5 × 10 −4 , for various values of center time of the fluctuations x 0 , and present the results in Fig. 3. As can be seen, when the fluctuations occur at an earlier time or higher temperature, the maxima and the final amplitudes of the hypermagnetic fields increase, and as a result, the matter-antimatter asymmetries increase as well. Second configuration is a constant vorticity with amplitude v 0 = 10 −2 . The results are presented in Fig. 4. As can be seen from the figure, the general trends of the evolution curves are similar. The prominent feature of this comparison is the surprising result that a fluctuation with amplitude smaller by three orders of magnitude produces results comparable with the constant vorticity configuration. Now we can address our assertion that µ/T ≪ 1 within our model, i.e., our initial conditions and results. First, note that upon using the relations µ f = (6s/T 2 )η f and s = (2π 2 g ⋆ /45)T 3 , we obtain (µ f /T ) = (12π 2 g ⋆ /45)η f = 280.95η f . The largest asymmetry in our results is obtained for η B and is shown in the Fig. 2(a). This figure shows that (η B ) max ∼ 5 × 10 −9 , leading to (µ/T ) max ∼ 10 −6 . The condition µ/T ≪ 1 has been used, for example, to justify neglect- ing terms of O(µ/T ) in Eqs. (5.5-5.8) in Appendix A, leading to the expressions for c v and c B shown in Eqs. (2.7,2.8), respectively. Moreover, as shown in Appendix B, this condition, along with the assumption of non-relativistic velocity of the plasma, imply that the CME and CVE contributions to the temporal components of the four-currents in the AMHD equations are negligible within our model.
We finally address the issue of two successive pulses with opposite temperature profiles to see by how much can the second pulse negate the results of the first. For this purpose it suffices to assume that the profile of the vorticities are unchanged. To be specific, we assume and We consider three cases in which the time separation of the pulses ∆x 0 = x 0,+ − x 0,− are 5b, b and 0.1b, where b denotes the width of the pulses. The results are shown in Fig. 5, where we also show our results for a single pulse for comparison. As can be seen, the final values of the asymmetries generated are reduced, as compared to the single pulse case, by a factor of about 5, 50 and 1000, respectively. The final values of the hypermagnetic field generated are reduced by square root of values stated above. It is interesting to note that even in the case ∆x 0 = 0.1b, the model produced η B ≃ 10 −13 , a value which can be increased easily by increasing β 0 , v 0 , and/or decreasing b.

Conclusion
In this study, we have investigated the contribution of the temperature-dependent CVE to the generation and evolution of the hypermagnetic fields and the matter-antimatter asymmetries, in the symmetric phase of the early Universe and in the temperature range 100 GeV ≤ T ≤ 10 TeV. The CVE has two possible sources in a vortical plasma, one from the chiralities and the other from the temperature of the particles. The former has been investigated in the literature much more than the latter. Here, we have focused on the latter in the form of transient temperature fluctuations, and have shown its important role in the production and evolution of the hypermagnetic fields and the matter-antimatter asymmetries. The transient fluctuations that we have considered are in the form of sharp Gaussian shaped pulses. In particular, we have shown that small simultaneous and transient fluctuations of vorticity about zero background value and temperature of some matter degrees of freedom about the equilibrium temperature of the plasma, close to the EWPT, can generate strong hypermagnetic fields and large matterantimatter asymmetries, even in the absence of any initial seed for the hypermagnetic field or any initial matter-antimatter asymmetries. Furthermore, we have shown that, an increase in the amplitude of temperature or vorticity fluctuations leads to the production of stronger hypermagnetic fields, and therefore, larger matter-antimatter asymmetries. This outcome has not been observed in any of the previous studies. In some studies which only take the CME into account, either an initially strong hypermagnetic field produces matter-antimatter asymmetries, or initial large matter-antimatter asymmetries strengthen a preexisting seed of hypermagnetic field [17][18][19]. In some other studies which also include the CVE and assume large initial chiralities in the vortical plasma, a seed of the hypermagnetic field is produced which then grows due to the CME [20].
In this work, we have considered a simple monochromatic helical configuration for the vorticity and hypermagnetic fields with a negative helicity, which ensures the production of the desired positive matter-antimatter asymmetries. Furthermore, we have shown that, either an increase in the amplitude of the temperature or vorticity fluctuations, or a decrease in their widths leads to the production of stronger hypermagnetic fields, and therefore, larger matterantimatter asymmetries. We have also shown that fluctuations in vorticity are several orders of magnitude more productive than constant vorticity. Within our model, the temperaturedependent CVE is the dominant effect as compared to the CME. We have also shown that when there are two sets of successive pulses with opposite temperature profiles, the first pulse is the main determinant of the outcome. A generalization of this work would be a stochastic analysis of fluctuations of both temperature and vorticity.

APPENDIX A
In this appendix we present the expressions for the chiral vorticity and helicity coefficients, c v and c B given by Eqs. (2.7), and (2.8). First, we start with the relevant and well known expressions in the broken phase. In relativistic hydrodynamics, the energy and number currents can flow separately in the presence of dissipative processes, therefore the definition of the flow will not be trivial [64]. Some common hydrodynamic frames in which the equations of AMHD may be formulated include the Landau-Lifshitz (or energy) frame [65], the Eckart (or conserved charge/particle) frame [66], and the more recently introduced anomalous "no-drag" frame [67]. In the latter, as its name suggests, a stationary obstacle experiences no drag, even when the energy and charge currents are present. In the Landau-Lifshitz frame, the energymomentum tensor T µν and the total electric current J µ for a plasma consisting of a single species of massless fermions of both chiralities are given by where F αβ = ∇ α A β −∇ β A α is the field strength tensor, p and ρ are the pressure and the energy density of the plasma, ρ el is the electric charge density, u µ = γ (1, v/R) is the four-velocity of the plasma normalized such that u µ u µ = 1, and γ is the Lorentz factor. 7 In the above equations, ν µ and τ µν denote the electric diffusion current and viscous stress tensor, respectively [3], Q R (Q L ) denotes the electric charges of the right-handed (left-handed) fermions, B µ = (ǫ µνρσ /2R 3 )u ν F ρσ is the magnetic field four-vector, and ω µ = (ǫ µνρσ /R 3 )u ν ∇ ρ u σ is the vorticity four-vector, with the totally anti-symmetric Levi-Civita tensor density specified by ǫ 0123 = −ǫ 0123 = 1. Furthermore, in the Landau-Lifshitz frame, the CME and CVE coefficients for chiral fermions are given as [3,[45][46][47][48][49] 7 Note that the self-consistency of our calculation in which the diagonal Einstein tensor obtained from the FRW metric is used implies that not only should the electromagnetic field density be small compared to the energy density of the Universe [68], but also the bulk velocity should obey the condition | v| ≪ 1, or equivalently γ ≃ 1 and u µ ≃ (1, v/R).
where T is the temperature and n R (n L ) is the right-handed (left-handed) charge density . 8 As we shall show explicitly in Sec. 3, µ R,L /T ≪ 1 within our model, i.e., our initial conditions and results. Hence, n R,L ≃ 1 6 µ R,L T 2 and ρ = 3p ≃ π 2 30 g * T 4 , and Eqs. (5.5-5.8) may be simplified as follows Moreover, since µ/T ≪ 1, we will consider only the Ohmic part of the diffusion current given by the first term, where σ = σ R + σ L is the electrical conductivity, and E µ = F µν u ν is the electric field four-vector. In order to carry over these results to the symmetric phase of the early Universe, it suffices to replace the electromagnetic field by the hypercharge gauge field and the electric charges of different particle species by their relevant hypercharges. Taking into account all three generations of leptons and quarks, we can easily obtain the chiral vorticity and helicity coefficients c v and c B , given by Eqs. (2.7), and (2.8), using Eqs. (5.9-5.12). The four-vectors B µ , ω µ , a µ , and E µ in the limit v ≪ 1 are given below, where a µ = Ω µν u ν is acceleration four-vector, Ω µν = ∇ µ u ν − ∇ ν u µ is vorticity tensor, a i = RΩ 0i is three vector acceleration. We have also used the assumption ∂ t ∼ ∇. v in the derivative expansion of the hydrodynamics, so v × E ≃ v 2 B and we have ignored the terms of O(v 2 ).

APPENDIX B
In this appendix we start with the basic expression for the anomaly given by Eq. (2.20) and obtain its explicit form which eventually leads to Eqs. (2.21,2.24,2.25). The equations of AMHD consist of energy-momentum conservation, Maxwell's equations and the anomaly relations. These equations may be expressed covariantly as whereF µν = 1 2R 3 ǫ µνρσ F ρσ is the dual field strength tensor, C R,L are the corresponding rightand left-handed anomaly coefficients, and j µ R,L are the fermionic currents given as follows where n R,L and V µ R,L are the chiral charge density and chiral particle diffusion current, respectively [3,[47][48][49]. The anomaly equations given above can be written out as Using Eqs. (6.8, 5.9-5.12) and x = (T EW /T ) 2 = t/t EW , we obtain (6.10) where M = 2π 2 g * /45. Here the electric field E is not an independent variable and is determined by Eq. (2.18). From the approximation of a vanishing displacement current, we have ∂ x E = − E/x. As is well known, in the symmetric phase of the early Universe T > 100 GeV, and in Sec. 3 we show that the condition µ/T ≤ 10 −6 holds within our model. These conditions, in conjunction with the low-velocity limit considered in this work, highly suppress the contributions of the temporal components of the CVE, CME and diffusion currents. To ascertain this claim, we have solved Eqs. (2.24)-(2.26) both with and without the additional temporal components and present the differences in Fig. 6. Comparing the scales of this figure and Fig.1, we observe the following for the differences: ∆η/η ∼ 10 −4 and ∆B Y /B Y ∼ 10 −9 .