Anomalous transport from holography: part II

This is a second study of chiral anomaly-induced transport within a holographic model consisting of anomalous U(1)V×U(1)A\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$U(1)_V\times U(1)_A$$\end{document} Maxwell theory in Schwarzschild–AdS5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_5$$\end{document} spacetime. In the first part, chiral magnetic/separation effects (CME/CSE) are considered in the presence of a static spatially inhomogeneous external magnetic field. Gradient corrections to CME/CSE are analytically evaluated up to third order in the derivative expansion. Some of the third order gradient corrections lead to an anomaly-induced negative B2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^2$$\end{document}-correction to the diffusion constant. We also find modifications to the chiral magnetic wave nonlinear in B. In the second part, we focus on the experimentally interesting case of the axial chemical potential being induced dynamically by a constant magnetic and time-dependent electric fields. Constitutive relations for the vector/axial currents are computed employing two different approximations: (a) derivative expansion (up to third order) but fully nonlinear in the external fields, and (b) weak electric field limit but resuming all orders in the derivative expansion. A non-vanishing nonlinear axial current (CSE) is found in the first case. The dependence on magnetic field and frequency of linear transport coefficient functions is explored in the second.


Introduction and summary
Fluid dynamics [1,2] is an effective long-wavelength description of most classical or quantum many-body systems at nonzero temperature. It is defined in terms of constitutive relations, which relate thermal expectation values of conserved currents to thermodynamical variables and external fields. The derivative expansion in fluid-dynamic variables such as the velocity or charge densities accounts for deviations from thermal equilibrium. At each order, the derivaa e-mail: yybu@post.bgu.ac.il b e-mail: lublinm@bgu.ac.il c e-mail: sharon.amir24@gmail.com tive expansion is fixed by thermodynamic considerations and symmetries, up to a finite number of transport coefficients, such as the viscosity, diffusion constant and conductivity. The latter are not calculable from hydrodynamics itself, but have to be determined from underlying microscopic theory or experimentally.
Although fluid dynamics has a long history, the theoretical foundations of relativistic viscous hydrodynamics are not yet fully established. The Navier-Stokes hydrodynamics leads to violation of causality: the set of fluid dynamical equations makes it possible to propagate signals faster than light. To overcome this problem, simulations of relativistic hydrodynamics are usually based on phenomenological prescriptions of [3][4][5][6], which admix viscous effects from second order derivatives, so as to make the fluid dynamical equations causal. References [3][4][5][6] introduced retardation effects for irreversible currents, which, via the equations of motion, become additional degrees of freedom. In other words, one needs to include higher order gradient terms in the derivative expansion in order to obtain a causal formulation. In general, causality is violated if the derivative expansion is truncated at any fixed order. It is supposed to be restored when all-order gradient terms are included, which we refer to as all-order resummed hydrodynamics. Resummed hydrodynamics is UV complete in a sense that it has a well-defined large frequency/momenta limit. Yet it is an effective theory of hydrodynamic variables only, 1 which emerges after most of the degrees of freedom of the underlying microscopic theory are integrated out.
The most general parity-even linear in external fields and charge density off-shell constitutive relation for a vector current has the following form: where ρ is a vector charge density and the diffusion D, electric/magnetic conductivities σ e/m are functionals of spacetime derivatives. In terms of hydrodynamic expansion, the constitutive relation (1) provides all-order resummation of gradients of the fluid-dynamic variables (the charge density ρ) and external fields ( E and B). In momentum space D and σ e/m are functions of frequency ω and momentum squared q 2 (assuming isotropic medium), which we refer to as transport coefficient functions (TCFs). Via the inverse Fourier transform, TCFs appear as memory functions in the constitutive relation [8].
For a holographic charged plasma dual to U (1) Maxwell theory in Schwarzschild-AdS 5 TCFs were studied in depth in [9]. The derivative resummation in the constitutive relation was implemented via the technique of [7,[10][11][12], which was originally invented to resum all-order velocity gradients (linear in the velocity amplitude) in the energy-momentum tensor of a holographic conformal fluid. 2 It is important to stress that this linearisation procedure is a mathematically well-controlled approximation: the perturbative expansion corresponds to a formal expansion in the amplitudes of fluiddynamic variables and external fields, without any additional assumptions. In this respect, the implemented approximation is identical to that of the linear response theory based on twopoint correlators.
Our technique follows closely the original idea of [14], which relates fluid's constitutive relations for the boundary theory to solving equations of motion in the bulk. However, an important new element of our formalism is that it is not based on current conservation (i.e., "off-shell" formalism), which makes it essentially different from the "onshell" formalism of [14]. Constitutive relations and TCFs can be uniquely determined from dynamical components of the bulk equations only, while the constraint component in the bulk is equivalent to continuity equation on the boundary. Chiral anomalies emerge and play an important role in relativistic QFTs with massless fermions. The anomaly is reflected in three-point functions of currents associated with global symmetries. When the global U (1) currents are coupled to external electromagnetic fields, the triangle anomaly renders the axial current non-conserved, where J μ /J μ 5 are vector/axial currents, and κ is an anomaly coefficient. For SU (N c ) gauge theory with a massless Dirac fermion in fundamental representation, κ = eN c /(24π 2 ), and e is an electric charge which below will be set to unit.
The chiral separation effect (CSE) [54,55] is another interesting phenomenon induced by the anomalies. It is reflected in a separation of the chiral charges along an external magnetic field at finite density of vector charges. Chiral charges can also be separated along external electric field, when both vector and axial charge densities are nonzero, the so-called chiral electric separation effect (CESE) [56,57].
In [68] we went beyond [9] focussing on transport properties induced by the chiral anomaly. The holographic model was modified to be anomalous U (1) V × U (1) A Maxwell theory in the Schwarzschild-AdS 5 case. Under various approximations, off-shell constitutive relations were derived for vector/axial currents. In a weak external field approximation, all-order derivatives in the vector/axial currents were resummed into six momenta-dependent TCFs: the diffusion, the electric/magnetic conductivity, and three anomaly-induced TCFs. The latter generalise the chiral magnetic/separation effects. Beyond weak external field approximation, nonlinear transports were also revealed when constant background external fields are present. Particularly, the chiral magnetic effect, including all-order nonlinearity in magnetic field, was proven to be exact when all external fields except for a constant magnetic field are turned off. Nonlinear corrections to the currents' constitutive relations due to electric and axial external fields were computed.
In the present work we continue the study of anomalyinduced transports within the holographic model of [68]. No axial external fields will be turned on in this work. As in [9,68] we work in the probe limit so that the currents and energy-momentum tensor decouple. In dual gravity, the probe limit ignores the back-reaction of the gauge dynamics on the geometry. The holographic model under study consists of two Maxwell fields in the Schwarzschild-AdS 5 black brane geometry. The chiral anomaly is holographi-cally realised via the gauge Chern-Simons actions for both Maxwell fields. Such a holographic setup can be realised via a top-down brane construction of D4/D8/D8 [69].
Before diving into the details presented in the following sections, we summarise our main results. The paper is split into two largely independent parts. In the first part, we consider a setup in which a static but spatially varying magnetic field is the only external field that is turned on. Then the constitutive relations for the vector/axial currents are where ρ/ρ 5 are vector/axial charge densities, the underlined terms in J i /J i 5 are the chiral magnetic/separation effects. G i , H i contain derivatives of ρ, ρ 5 , B and are defined in Sect. 3. It is important to stress that, in contrast to the above discussion of linearised hydro, (3) and (4) are exact, without any approximations for ρ, ρ 5 , B. The nonlinearity of the CME/CSE in external magnetic field B is completely accounted for by the chemical potentials μ, μ 5 . The non-derivative part of (3) is consistent with the "nonrenormalisability" of CME [48,49,70,71]. However, as will be clear from (48) and (49), the derivative corrections in G i , H i introduce new effects, which do modify the original CME. Particularly, the currents along the direction of B get affected.
When ρ, ρ 5 , B vary slowly from point to point, G i , H i can be calculated order-by-order within boundary derivative expansion. Let us introduce a scaling parameter λ: Then derivative counting goes by powers of λ. Up to second order in derivative expansion, we calculated G i , H i and chemical potentials μ, μ 5 without any further assumptions. Given that these results are rather lengthy, we postpone to present them in Sect. 3; see (48) and (49) and (50). At third order O(∂ 3 ), for G i , H i we calculated only terms that are linear in ρ, ρ 5 ; see (51) and (52) for a complete listing. Among these third order terms, the diffusion constant D 0 (i.e., the DC limit of the diffusion function D) gets a negative B-dependent correction To the best of our knowledge, this is the first anomaly-induced correction to the diffusion constant, and, being negative, it happens to violate the universal form of [72].
With the third order results for J μ and J μ 5 , we also computed the dispersion relation for a free mode that can propagate in the medium: where B means a constant magnetic field. The first term in (7) represents the chiral magnetic wave (CMW) [70]. Interestingly, we see nonlinear in B corrections to both the speed of CMW and its decay rate. Note that we also expect the emergence in (7) of the following terms: ( q · B) 2 , q 2 ( q · B), ( q · B) 3 , q 2 ( q · B) 2 and ( q · B) 4 . However, our ability to determine coefficients of these terms is limited by the undertaken approximations.
In the second part of this work, we focus on a special setup which is experimentally accessible in condensed matter systems. 4 CME emerges from a nonzero axial chemical potential μ 5 , which is usually assumed to have some background profile. It is, however, possible to induce ρ 5 (and thus μ 5 ) dynamically through the interplay between the electric and magnetic fields, as is clear from the continuity equation (2). Specifically, we are ready to consider a constant magnetic field B and a time-dependent but spatially homogeneous electric field E(t). For simplicity the charge densities ρ, ρ 5 will be assumed to be spatially homogeneous too 5 . From (2), ρ could be set zero. The constitutive relations for the vector/axial currents are where V j (1), A j (1), G i and H i depend on ρ 5 , E and B nonlinearly and will be computed below. Our study is further split into two parts. In Sect. 4.1, V j (1), A j (1), G i and H i will be evaluated perturbatively within the gradient expansion (5). These perturbative results can be found in (65)- (68). In Sect. 4.2, we will consider another approximation-linearisation of the constitutive relations in the external electric field. In the linearised regime, we assume the following scaling for ρ 5

, E and B:
where σ e is a q 2 = 0 limit of the electric conductivity introduced in (1), while τ 1,2 are new TCFs. As with other TCFs, they are functionals of time derivative operator and become functions of frequency ω in Fourier space, At the linear level (in external fields and hydro variables), the transport coefficient functions in [9,68] were proved to be frame independent. Along this line of proof, we expect that σ e , τ 1,2 are also independent of the hydro frame choice. Imposing the continuity equation (2), the electric current is put on shell, where the transverse conductivity σ T is not affected by the magnetic field in contrast to the longitudinal conductivity σ L , which gets corrected by the magnetic field via the chiral anomaly. In Sect. 4.2 the TCFs τ 1 , τ 2 will be first analytically evaluated in the hydro limit and then numerically for arbitrary frequency.
While there is some overlap between our results and the literature, differences between the present study and those of [73][74][75][76] must be clarified. Utilising the weak electric field approximation (10), [73] analytically evaluated the magnetic field dependence of the longitudinal conductivity σ L in DC limit, while [74] calculated its ω-dependence. Back-reaction effects on σ L were considered in [76]. References [77] performed a similar study, focussing on time evolution of the induced vector current, given some specially chosen initial profile for the electric field. All the studies [73][74][75][76] focussed on a weak electric field, in which the axial current vanishes. Therefore, our nonlinear results and particularly the axial charge separation current (66) appears as new. As for the linearised setup (10), [73][74][75][76] imposed the continuity equation and replaced the axial charge density ρ 5 in favour of the external electric and magnetic fields, so the vector current there is on shell. This is in contrast to our off-shell formalism. As we argued in our previous publications [7,[9][10][11][12], only offshell construction reveals transport properties of the system in full. Particularly, there are three independent TCFs (σ e and τ 1,2 ) in the constitutive relation (11), all of which we are able to determine separately, compared to only two independent conductivities in (13).
Another difference worth mentioning is that we explicitly trace all the effects in the induced current that arise from the relative angle between E(t) and B fields. This is in contrast to [74,77], which limited their study to the case of parallel fields only, primarily focussing on the longitudinal electric conductivity σ L . By varying the relative angle between E(t) and B fields, one can separate the anomaly-induced effects (parametrised by τ 1 and τ 2 ) from the ones that are not related to the anomaly (σ e ).
The paper is structured as follows. In Sect. 2 we present the holographic model and outline the strategy of deriving the boundary currents from solutions of the anomalous Maxwell equations in the bulk. Section 3 presents the first part of our study: CME/CSE with static but varying in space magnetic field. In Sect. 4, CME/CSE in the presence of constant magnetic and time-varying electric fields are analysed. This study is further split into two subsections. The exploration of nonlinear phenomena in the induced vector/axial currents is done in Sect. 4.1. In Sect. 4.2 we focus on the linearised regime (10) and calculate the dependence of AC conductivity on magnetic field. Section 5 presents the conclusions. Two appendices supplement computations of Sects. 3 and 4.

The holographic model: U(1) V × U(1) A
The holographic model is the U (1) V × U (1) A theory in the Schwarzschild-AdS 5 case. The chiral anomaly of the boundary field theory is modelled via the gauge Chern-Simons terms in the bulk action, where and the counter-term action S c.t. is The field strengths M N P Q R is the Levi-Civita symbol with the convention rt x yz = +1, and the Levi-Civita tensor is M N P Q R / √ −g. Our choice for (16) is based on minimal subtraction, that is, the counter-term does not make a finite contribution to the boundary currents.
In the ingoing Eddington-Finkelstein coordinates, the spacetime metric is where f (r ) = 1 − 1/r 4 , so that the Hawking temperature (identified as temperature of the boundary theory) is normalised to π T = 1. On the constant r hypersurface , the induced metric γ μν is The equations of motion for V and A fields are as follows.
Constraint equations : where The boundary currents are defined as which, in terms of the bulk fields, are where n M is the outpointing unit normal vector on the slice , and ∇ is compatible with the induced metric γ μν . The currents (24) are defined independently of the constraint equations (21). Throughout this work, the radial gauge V r = A r = 0 will be assumed. Consequently, in order to completely determine the boundary currents (25) it is sufficient to solve the dynamical equations (20) for the bulk gauge fields V μ , A μ only, leaving the constraints aside. The constraint equations (21) give rise to the continuity equations (2). In this way, the currents' constitutive relations to be derived below are off shell.
It is useful to reexpress the currents (25) in terms of the coefficients of the near-boundary asymptotic expansion of the bulk gauge fields. Near r = ∞, where In (26) the constant term for A μ is set zero given that axial external fields are turned off in our present study. The holographic dictionary implies that V μ is a gauge potential of external electromagnetic fields E and B, When obtaining (26) and (27), only the dynamical equations (20) were utilised. The the near-boundary data V (2) μ and A (2) μ have to be determined by completely solving (20) from the horizon to the boundary. The currents (25) become The remainder of this section is to meant to outline the strategy for deriving the constitutive relations for J μ and J μ 5 . To this end, consider finite vector/axial charge densities exposed to external electromagnetic fields. Holographically, the charge densities and external fields are encoded in asymptotic behaviours of the bulk gauge fields. In the bulk, we will solve the dynamical equations (20) assuming some charge densities and external fields, but without specifying them explicitly.
Following [9] we start with the most general static and homogeneous profiles for the bulk gauge fields which solve the dynamical equations (20), where V μ , ρ, ρ 5 are all constants for the moment. The regularity requirement at r = 1 fixes one integration constant for each V i and A i . As explained below (27), the constant in A μ is set to zero. Through (29), the boundary currents are Hence, ρ and ρ 5 are identified as the vector/axial charge densities. Next, following the idea of fluid/gravity correspondence [14], we promote V μ , ρ, ρ 5 into arbitrary functions of the boundary coordinates, Then (30) ceases to be a solution of the dynamical equations (20). To have them satisfied, suitable corrections in V μ and A μ have to be introduced: where V μ , A μ will be determined from solving (20). Appropriate boundary conditions have to be specified. First, V μ and A μ have to be regular over the whole integration interval of r ∈ [1, ∞]. Second, at the conformal boundary r = ∞, we require which amounts to fixing the external gauge potentials to be V μ and zero (for the axial fields). Additional integration constants will be fixed by the Landau frame convention for the currents, The Landau frame choice can be identified as a residual gauge fixing for the bulk fields. The vector/axial chemical potentials are defined as Generically, μ, μ 5 are nonlinear functionals of densities and external fields.
In terms of V μ and A μ , the dynamical equations (20) are In the following sections we will present solutions to the dynamical equations (37)-(40) under the two independent setups discussed in Sect. 1.

CME/CSE with time-independent inhomogeneous magnetic field
In this section we consider the case in which the magnetic field is the only external field that is turned on. The magnetic field is assumed to be varying in space, but it should be time independent to avoid creating an electric field. There is no restriction on the charge densities ρ, ρ 5 . From the general results (26) and (27), In obtaining large r estimates for V t and A t , the frame convention (35) was used to fix the coefficients of 1/r 2 in the near-boundary expansion for V t , A t (thus those of V t and A t ). The dynamical equations (37)-(40) get simplified, For generic profiles of ρ, ρ 5 and B( x), nonlinearity makes it difficult to solve (42)- (45). To explore the general structure of the vector/axial currents, we rewrite the dynamical equations (42)-(45) into integral forms. In this way, the nearboundary asymptotic expansion for V μ and A μ could be extracted from the integral forms of (42)- (45). For simplicity, we address the details in Appendix A. Substituting the near-boundary behaviour (94)-(97) into (29) produces the results (3) and (4). As mentioned below (3) and (4), G i , H i are functionals of ρ, ρ 5 , B and are presented in (98) and (99). The formal analysis establishes the structure of J μ /J 5 μ , particularly the "non-renormalisation" of CME and its gradient corrections.
We proceed with the hydrodynamic gradient expansion for J μ /J 5 μ . This requires us to perturbatively solve the dynamical equations (42)-(45) within the boundary derivative expansion (5), The corrections V μ and A μ are expandable in λ, At each order in λ, V [n] μ and A [n] form a system of ordinary differential equations in r -coordinate, which can be solved via direct integration over r . The results for V [n] μ and A [n] up to n = 2 can be found in Appendix A; see (100)-(106).
Substituting the first order solutions (100)-(102) into (98) and (99) generates a hydrodynamic expansion for G i , H i up to second order in the gradient expansion (throughout this work, the electromagnetic fields are thought of as of first order in derivative counting), Meanwhile, the second order results (103) and (104) give rise to the gradient expansion of chemical potentials (36) In principle, the second order results (103)-(106) could be inserted into (98) and (99), producing derivative expansion for G i (x = ∞) and H i (x = ∞) up to third order. However, at third order O(∂ 3 ), computing G i , H i becomes quite involved. Therefore, at third order O(∂ 3 ) we decided to track only terms linear in ρ, ρ 5 . As a result, we are able to identify the first anomalous correction to the diffusion constant D 0 due to magnetic field. The final expressions are G [3] i where # 1 in (51) is given by the integral where b 2 (r ) is given in (109). The underlined term in (51) is a κ 2 B 2 -correction to the diffusion constant D 0 . Given that the lowest order anomalous correction to the diffusion constant is negative, it is interesting to explore this effect further for arbitrary magnitude of the magnetic field, which, however, goes beyond the scope of the present study.
Our results for J μ and J μ 5 can be used to explore dispersion relations for free modes propagating in the chiral medium. We consider a constant magnetic field only. Let us take a plane wave ansatz for the vector/axial charge densities Then the continuity equation (2) becomes which has a nontrivial solution when and only when where Solving (56) leads to the B-corrected dispersion relation, as summarised in (7).

CME/CSE with constant magnetic and time-dependent electric fields
Creating systems with chiral imbalance (μ 5 = 0) experimentally is problematic. In this section we consider a special setup in which the axial chemical potential μ 5 is not imposed externally but rather is induced dynamically through the chiral anomaly. This setup is of particular interest due to intriguing possibility for it to be realised experimentally in chiral condensed matter systems. Consider a constant magnetic field B and a time-dependent homogeneous electric field E(t). We also assume the charge densities to be spatially homogeneous as well. 6 The continuity equation (2) degenerates to which implies that the vector charge density is constant while the axial charge density has nontrivial time dependence inher- 6 While from the continuity equation (2) the charge densities can still have a nontrivial spatial-dependence, we found that such spatial inhomogeneity of the charge densities would make the gradient resummation out of control.
ited from E(t). The setup under consideration is Under the frame convention (35), the corrections V μ and A μ of (33) depend on r and t only. As a result, the dynamical equations (37)- (40) are reduced to

Nonlinear phenomena: general analysis and derivative expansion
The objective of this subsection is to show that beyond linearised limit (69) the setup (60) also induces a non-vanishing axial current J 5 , which has been omitted in the literature. To this end, as in Sect. 3, we first give a fully nonlinear analysis for the dynamical equations (61)-(64), followed by perturbative calculations for V μ , A μ within the derivative expansion (5). All calculational details are addressed in Appendix B. As in Sect. 3 the formal analysis are based on rewriting the dynamical equations (61)-(64) into integral form, from which one could deduce the near-boundary asymptotic behaviours for V μ , A μ . The results can be found in (112)-(115). Plugged them into (29), the near-boundary behaviour for V μ , A μ presented in (112)-(115) is translated into boundary currents (8) and (9). Generically, the quantities V i (1), (8) and (9) where C is a Catalan constant and # 1 is known numerically only, Up to second order in derivatives O(∂ 2 ), the chemical potentials (36) are 7 Evaluated on shell via (2), the axial current J i 5 is fully nonlinear in the amplitude of the electric field E(t), as clear from (66).

Linear in E phenomena
In the previous subsection we focussed on hydrodynamic regime, in which we were able to identify some nonlinear phenomena. Below, we proceed with an alternative approximation, that is, the weak electric field approximation (10): The scaling of ρ 5 follows from the continuity equation (59). Both corrections V μ and A μ are of order O( ) too. The dynamical equations (61)-(64) get further simplified, Integrating (70) and (72) over r once, we get where the frame convention (35) was used to fix the integration constant. Equation (74) makes it possible to decouple (71) and (73) become 7 While we suspect that the chemical potential μ is zero to all orders in the gradient expansion, we have not been able to prove this.
The homogeneity property of (76), combined with the regularity requirement at r = 1 and vanishing boundary condition at r = ∞ for A i , fixes A i = 0 completely. From (74), Therefore, at order O( ), the axial current J 5 = 0 as read off from (29). This is in contrast with the nonlinear analysis of Sect. 4.1.
The differential equation (76) is linear in the correction V i . Therefore, (76) can be solved via the technique developed in [7,[10][11][12]. The bulk equations reduce to linear inhomogeneous partial differential equations while the inhomogeneous terms are built from boundary derivatives of the fluiddynamic variables and external fields. The equations then can be exactly solved using the Green function formalism: the bulk fields are decomposed in terms of all possible basic vector structures constructed from the fluid-dynamic variables and external fields. These decomposition coefficients (components of the inverse Green function) are functions of the holographic radial coordinate and functionals of the boundary derivative operators. The functional dependence of the decomposition coefficients on the boundary derivative operation encodes all-order linear derivatives in the constitutive relations. Transformed into momentum space, the bulk equations give rise to ordinary differential equations for those decomposition coefficients, which are RG-like equations in AdS space. Solving the RG-like equations completely determines the fluid's constitutive relations and all transport coefficients. Below we implement these steps.
V i is decomposed as 8 where The decomposition coefficients C i satisfy the partially decoupled ordinary differential equations (ODEs), While C 1 does not feel the effect of the magnetic field, C 2,3 have a nontrivial dependence on the magnetic field via κ 2 B 2 . Near r = ∞, the pre-asymptotic expansions of the C i are where the c i are boundary data and have to be fixed through full solution of (80)-(82) from the horizon r = 1 to the conformal boundary r = ∞. From (29), the conductivities of (11) are determined by the boundary data c i , The ODE for C 1 was solved in [9]. The conductivity σ e , which is computed from C 1 , was completely determined and explored in [9], while only q = 0 limit enters into our current study (the results are quoted below). We therefore focus on the remaining two conductivities τ 1 , τ 2 , both induced by the chiral anomaly. As is obvious from (80)-(82), τ 1 , τ 2 depend on the magnetic field via κ 2 B 2 .
Using the continuity equation (59), the constitutive relations (11) are put into a linear response form, from which on-shell current-current correlators can be read off. Since the electric field is the only external perturbation that is turned on, it is possible to compute only a subset of all two-point correlators in the theory, where J i J j is split into transverse (G T ) and longitudinal (G L ) components with respect to the direction of B.
To determine the remaining current-current correlators we would have to introduce additional field perturbations, particularly an axial external field, which is beyond the scope of this paper.
To evaluate the TCFs τ 1 , τ 2 , we have to completely solve the ODEs (80)- (82). We first analytically solve them when ω = 0. As a result, the DC limits τ 0 1 (for arbitrary B) and τ 0 2 (up to leading B 2 -correction) are known analytically, where [z] is the Gamma function, and # 2 is known numerically only, When the magnetic field is very strong, τ 0 1 and τ 0 2 behave similarly The result for τ 0 1 is in agreement with [74,77]. In the DC limit ω → 0, when the magnetic field is very strong, the on-shell vector current (13) behaves as which is in agreement with [77]. When ω → 0 (DC limit), the current-current correlator is dominated by the chiral anomaly-induced effects ∼ τ 0 1 . The DC limit is of interest for experiments with electric fields turned on adiabatically, such as the ones considered in [77]. Meanwhile, when ω → 0, the longitudinal conductivity σ L in (13) is parametrised as where σ 0 e = 1, τ 1 1 is the coefficient of iω in the hydrodynamic expansion of τ 1 . For illustration, in Fig. 1 we show the κBdependence of τ 0 1 , τ 0 2 (divided by 5 to match scales), τ 1 1 and Re σ 0 L . The behaviour of Re σ 0 L agrees perfectly with that of [74].
In our calculation, Re σ 0 L acquires a negative correction due to the magnetic field and eventually vanishes when the magnetic field gets large; see Fig. 1. This is in contrast with many related studies of negative magnetoresistivity, the phenomenon of enhancement of the longitudinal DC conductivity due to the magnetic field [78][79][80][81][82][83]. However, taking a strict DC limit in σ 0 L is problematic due to the explicit 1/ω divergence. The latter is frequently regularised by introduction of axial charge dissipation effects via shifting the frequency ω → ω + i/τ 5 , where τ 5 corresponds to some relaxation time. The physics of this axial charge relaxation is beyond the scope of the present work. It was addressed within the holographic approach in [73][74][75][76]. These studies primarily rely on the Kubo formula. For arbitrary ω, we resort to numerical methods and solve ODEs (80)- (82) for representative values of κB. The numeri-cal procedure is identical to that of [68] and for all the numerical details we refer the reader to this publication. In Fig. 2 we show the ω-dependence for τ 1 and τ 2 for sample choices of κB. In Fig. 3 we plot the normalised TCFs τ 1 /τ 0 1 and τ 2 /τ 0 2 . Overall, τ 1 and τ 2 display quite similar dependences on the frequency ω. After some oscillations, both τ 1 and τ 2 approach zero asymptotically. Approach to the asymptotic regime, however, depends on strength of the magnetic field. When κB is increased, the asymptotic behaviour is delayed towards larger ω. What is more intriguing is that increasing κB renders τ 1 and τ 2 to develop a resonance-like enhancement at finite ω. This could be an interesting experimentally observable feature. For very strong magnetic fields κB → ∞, the chiral anomaly-induced effects would be pushed to the UV, corresponding to early time effects, such as in [77].
In Fig. 4 we show two-point correlators G T,L for different choices of κB. However, it is difficult to appreciate the anomaly-induced effects from Fig. 4 because in the correlators they get mixed with non-anomalous ones. To illuminate κB-correction to G L , in Fig. 5 we plot the difference δG L = G L − G T . From these plots, the effect of the chiral anomaly on the induced vector current is seen more clearly. We again notice a remarkable relative enhancement at intermediate values of ω.

Conclusions
In this paper we continued explorations of the chiral anomaly induced transport within a holographic model containing two U (1) fields interacting via Chern-Simons terms. For a finite temperature system, we computed off-shell constitutive relations for the vector/axial currents responding to external electromagnetic fields.
When a static spatially inhomogeneous magnetic field is the only external field that is turned on, we showed that both the CME and the CSE get corrected by derivative terms; see (3) and (4). Within the derivative expansion, we analytically calculated corrections up to third order in the expansion; see (48) and (49) and (51) and (52). Apart from the derivative corrections to CME and CSE, the diffusion constant D 0 was found to receive a negative anomaly-induced correction; see (6). The dispersion relation of the chiral magnetic wave was also found to be modified; see (7).
In the second part of our study, we focussed on the case of time-varying electric and constant magnetic fields without any externally enforced axial charge asymmetry, though the E(t) · B term in the continuity equation (59) generates the axial charge density ρ 5 (and thus μ 5 ) dynamically. For such a configuration of the external fields, we first analysed the most general constitutive relations for the vector/axial currents; see (8) and (9). Then, within the derivative expansion, we explicitly calculated the currents up to third order at nonlinear level; see (65) and (66). When put on shell, the axial current J 5 is fully nonlinear in the external electric field.
Employing another approximation, we linearised the constitutive relations assuming the electric field is weak (10). Within this approximation the axial current is zero, while the "off-shell" vector current is parameterised by three frequency-dependent transport coefficient functions: the electric conductivity σ e , and two chiral anomaly-induced conductivities τ 1 , τ 2 ; see (11). In the DC limit, we analytically computed these conductivities; see (88) and (89). Then, for generic ω, the numerical plots were presented in Sect. 4.2. Based on these studies, we notice that the anomaly-induced effects get enhanced at some finite frequency ω, whereas the position of the maximum and strength of the effect depends on the external magnetic field. It might be an effect worth looking for experimentally.