Quasinormal modes of charged magnetic black branes&chiral magnetic transport

We compute quasinormal modes (QNMs) of the metric and gauge field perturbations about black branes electrically and magnetically charged in the Einstein-Maxwell-Chern-Simons theory. By the gauge/gravity correspondence, this theory is dual to a particular class of field theories with a chiral anomaly, in a thermal charged plasma state subjected to a constant external magnetic field, $B$. The QNMs are dual to the poles of the two-point functions of the energy-momentum and axial current operators, and they encode information about the dissipation and transport of charges in the plasma. Complementary to the gravity calculation, we work out the hydrodynamic description of the dual field theory in the presence of a chiral anomaly, and a constant external $B$. We find good agreement with the weak field hydrodynamics, which can extend beyond the weak $B$ regime into intermediate regimes. Furthermore, we provide results that can be tested against thermodynamics and hydrodynamics in the strong $B$ regime. We find QNMs exhibiting Landau level behavior, which become long-lived at large $B$ if the anomaly coefficient exceeds a critical magnitude. Chiral transport is analyzed beyond the hydrodynamic approximation for the five (formerly) hydrodynamic modes, including a chiral magnetic wave.

understanding of the physics governing the QGP in these HICs is likely to lead to prediction of more distinct and identifiable experimental signatures. This motivates our study of strongly coupled charged thermal states in presence of magnetic fields and anomalies in this paper.
On the field theory side, some results were previously obtained for plasmas which are anisotropic due to a magnetic field. Anomaly-driven effects such as CME and CVE reveal themselves in the transport properties of a system outside of equilibrium. The hydrodynamic description of chiral transport in weak magnetic fields 2 in the case of a single axial current has been considered recently in [26,27]. In these works, part of the hydrodynamic modes were given explicitely. In this paper, we perform an independent systematic and complete calculation of the relevant hydrodynamic modes at nonzero chemical potential, charge density, and magnetic field. Whenever the relevant quantities are provided, we find complete agreement with [27], and we agree at vanishing chemical potential and charge density with the results shown in [26] for ω 1 , ω 2 , ω 3 (see also Appendix A for a discussion of the different hydrodynamic frames). A hydrodynamic analysis similar to [26,27] for a system containing both the axial and a conserved vector current has also appeared very recently [28]. There (as well as in [26]), excitations propagating perpendicular to the magnetic field were also considered, while in [27] and the present work, the focus is on propagation parallel to the magnetic field. In the weak field case the chiral magnetic wave [29] appears in a certain form, as is evident from [27,28] and from our result in Appendix A. Fewer reliable hydrodynamic results are available for plasmas in strong external magnetic fields 3 . One exception is a work considering polarization effects in equilibria with strong external (electro)magnetic fields [30], generalizing [31]. An attempt at hydrodynamic constitutive relations in presence of a strong external magnetic field has been performed in [29]. The striking finding of that latter paper are the five "shear viscosities" and the two "bulk viscosities" which appear because of the broken rotational symmetry due to the magnetic field. A strong claim was conjectured in [32], where the authors argue for the chiral magnetic wave velocity to be of a particular form for any strength of the magnetic field. The authors also find that this velocity approaches the speed of light in the limit of large magnetic field at weak (and at strong coupling in the probe limit), and relate this limit behavior to Landau level physics using a weak coupling argument. In our work, we find evidence for the relation between Landau level physics and the chiral magnetic wave velocity at strong coupling, not assuming any probe limit.
Despite previous field theory achievements, the study of non-equilibrium dynamics, especially for strongly interacting systems, has been difficult using traditional methods. However, the problem becomes particularly amenable in gauge/gravity duality [33][34][35][36], because not only does it provide new tools to study strong interactions, but also the complicated quantum evolution problem can be turned into a conceptually simple one of solving differential equations with given initial conditions. By now, it is evident that QGP is a strongly interacting system from the many experimental and theoretical studies 2 That is, magnetic fields which are of first order in derivatives. See Appendix A for details on the definition of strong and weak magnetic fields and about the derivative counting. 3 Strong magnetic fields are of zeroth order in a derivative expansion, i.e. they are O (1). See Appendix A.
supporting this [37], and is thus naturally a candidate for the application of holography. For Weyl/Dirac semimetals, this is not yet clear. But then an exploration from the stronglyinteracting perspective is even more interesting. In that vein, a host of papers explored the interplay between external magnetic fields and a chiral anomaly within holographic models of Weyl semimetals [38][39][40][41][42][43].
Within holography transport is encoded in quasinormal modes (QNMs), see for example the reviews [44,45]. 4 Numerous studies have shown that QNMs play an important role in the evolution dynamics of strongly interacting non-equilibrium systems, particularly in the relaxation process [48][49][50][51][52][53]. On the gravity side, QNM frequencies are (quasi)eigenvalues of the linearized Einstein equations describing fluctuations around a black hole or black brane gravity background. Since the operator corresponding to the eigenvalue problem is not self-adjoint, the QNM frequencies are complex in general. These QNM frequencies correspond to the locations of poles of the retarded correlators in the dual field theory [54,55]. They characterize the near-equilibrium behavior, encoding information about transport and dissipation. For correlators of conserved quantities such as the energy-momentum tensor, the QNM spectrum typically contains an infinite tower of gapped modes that are strongly damped, as well as a set of hydrodynamic modes ω = ω(k) such that ω(k) → 0 as k → 0 [55,56]. At vanishing magnetic field, quasinormal modes of charged black branes (also called Reissner-Nordström black branes) have been calculated, see [51], references and the discussion therein. Those references worked either at vanishing momentum for the fluctuations or did not consider a Chern-Simons term, i.e. in those cases effects of the anomaly remained unknown. Fluctuations about a charged black brane in presence of a Chern-Simons term but at vanishing magnetic field were considered holographically in [57,58] in the hydrodynamic limit. Previously, that system was studied holographically in the hydrodynamic limit without the Chern-Simons term [59,60]. In the present work, we are closing the gaps in the gravitational calculations and present a systematic study of all metric and gauge field fluctuations with nonzero momentum along the external magnetic field in a background of charged magnetic black branes in presence of a Chern-Simons term within a particular theory which includes the aformentioned charged black brane setups as special cases.
In this paper, we consider (3+1)-dimensional systems at nonzero temperature, with a nonzero charge density, and in the background of a constant magnetic field. Most importantly, we consider such systems in the presence of an anomaly, and analyze perturbations with nonzero momentum and frequency in order to derive transport properties. A holographic model describing such a class of systems are charged magnetic black brane solutions with Einstein-Maxwell-Chern-Simons (EMCS) theory [61]. This model is "top-down" at least for a particular value of the Chern-Simons coupling, as it can then be derived from a consistent truncation of supergravity [62][63][64][65], and it includes full backreaction from all the matter fields in its gravitational dynamics. Previously, the QNM spectrum of this model with only the magnetic field turned on has been studied [51], while [66,67] computed two shear viscosities in this anisotropic system (see also [29]). The thermodynamics and the phase structure of this system at nonzero charge density and magnetic field has also been studied [68]. Here, we shall analyze the QNM spectrum of the latter system, and how it is influenced by the strength of electromagnetic background, as well as the anomaly. We solve for the background numerically using spectral methods, which we then use as input to compute the QNMs. We investigate the effects of changing external parameters: the temperature, T , chemical potential, µ, and the strength of the magnetic field, B. We find interesting phenomena such as the appearance of Landau levels in some cases and various other gapped modes. In the hydrodynamic regime, we compare to our field theoretic calculation provided in Appendix A. In particular, we provide dispersion relations for five modes in (A.20) and (A.17). Following the modes (also beyond the hydrodynamic regime), we describe the evolution of the shear modes, see Fig. 2, of two formerly hydrodynamic modes developing a (complex) gap, see Fig. 3, and of three hydrodynamic modes, see Fig. 7, all for increasing B. Our results are summarized in Sec. 5.
The paper is organized as follows. In Sec. 2, we discuss the setup of the holographic model dual to the (3+1)-dimensional system including a chiral U (1)-anomaly. The procedure of setting up the numerics by linearizing the field equation is presented in Sec. 3. In Sec. 4 we compute the QNMs. We present how they behave as momentum and the strength of the magnetic field as well as the anomaly is varied, and we discuss their implications for transport and hydrodynamics with anomalies, and report on the Landau level behavior found in some cases. We conclude in Sec. 5. The construction of the hydrodynamic description is to be found in Appendix A. Details of the numerical method employed are recorded in Appendix B.

Holographic setup
Our goal is to study charged fluids in a particular class of quantum field theories with a chiral anomaly at strong coupling, and in the presence of a (strong) magnetic background field. Our treatment below is generally applicable, but let us first illustrate it using the N = 4 Super-Yang-Mills (SYM) theory as a concrete example. The N = 4 SYM is coupled to a U (1) symmetry, which is a subgroup of the global SU (4) R R-charge symmetry. Its matter content consists of one vector, four left-handed Weyl fermions, and six real scalars, all of which transform in the adjoint representation of the SU (N c ) color group. The fermions and scalars are charged under the U (1), while the vector is uncharged.
Consider the case in which an axial current, J α , of the U (1) is coupled to a constant external magnetic field B, see e.g. [50]. The action is then In this theory, the energy-momentum tensor, T µν , and the axial current satisfy the following conservation equations: where αβγδ is the totally antisymmetric Levi-Civita tensor in four-dimensional (flat) spacetime with 0123 = 1, and C is the chiral anomaly coefficient [3,69,70]. The Greek indices here run from 0 to 3, and they indicate field theory coordinates . However, for the more general theories which we consider below, C is arbitrary.
In order to obtain results for a charged fluid within this strongly coupled theory, we use the gauge/gravity correspondence, and we perform the relevant calculations in the dual gravitational description, which is given by an Einstein-Maxwell-Chern-Simons (EMCS) theory with an external magnetic field. The relevant fully backreacted gravitational solutions dual to the desired plasma states are the charged magnetic black branes [68,71]. Below, we first review the EMCS theory and its charged magnetic black brane solutions. We then discuss the thermodynamics in the gravity description, and match it to the dual field theory.

Einstein-Maxwell-Chern-Simons theory
The EMCS theory in five dimensions is defined by the action where L is the AdS 5 radius, 2κ 2 ≡ 16πG 5 , and g ≡ det g mn , with G 5 and g mn the fivedimensional Newton's constant and metric respectively. In the gauge sector, F = dA is the five-dimensional Maxwell field strength, and γ is the Chern-Simons coupling. We denote the five-dimensional gravitational bulk by M, and ∂M its boundary. The coordinates of M are indicated by lower case Roman indices, which run from 0 to 4. The action given in (2.4) has to be amended by boundary terms [72][73][74] of the form Here,ĝ µν is the metric induced by g mn on the conformal boundary of AdS 5 , and K is the trace of the extrinsic curvature with respect toĝ µν . The extrinsic curvature K mn is given by where ∇ is the covariant derivative, and n m is the outward pointing normal vector of ∂M. The equations of motion following from the action (2.4) read where˜ mnopq is the totally antisymmetric Levi-Civita symbol in five spacetime dimensions with˜ 01234 = 1. The Chern-Simons coupling γ on the gravity side is related to the chiral anomaly coefficient C on the field theory side via [75,76] C = γ 6 .
(2.9) This identification will be crucial for the analysis of the thermodynamics below, and for the subsequent analysis of the chiral transport effects. For simplicity, we will work in units where L = 1 and 2 κ 2 = 1 from now on.

Charged magnetic black brane solutions
The charged magnetic black brane solutions of the EMCS equations, (2.7) and (2.8), are dual to the charged plasma of interest in a particular class of quantum field theories with a chiral anomaly at strong coupling, and in the presence of a magnetic background field of any strength. Due to the presence of various scales including the temperature, the chemical potential, and the magnetic field, conformal symmetry is broken by the state of such field theories.
As we shall see below when studying the fluctuations, it is convenient to use the Eddington-Finkelstein coordinates. Working in units where the horizon of the black brane is located at z = 1 and the conformal boundary at z = 0, the metric for the charged magnetic black brane can be expressed as The field strength tensor (2.11) may be obtained from a gauge field A of the form which is symmetric in x 1 and x 2 , and points to the SO(2) rotational symmetry in the (x 1 , x 2 )-plane. Near the conformal boundary at z = 0, the functions appearing in (2.10) and (2.11) may be expanded as where u 4 , w 4 , c 4 , ρ, p 1 are undetermined coefficients, while B and µ are parameters (corresponding to sources in the field theory) determining the background solution. We see from this that the charged magnetic black brane is asymptotically AdS. Note that the leading term in P (z) has been fixed such that no explicit source for a persistent current in the x 3 -direction appears.
Near the horizon at z = 1, expansion of the same functions yields whereū 1 ,c 1 ,w 0 ,v 0 ,Ā v0 andP 0 are undetermined coefficients analogous to those in the near boundary expansion (2.13). Note that regularity at the horizon fixes the leading coefficients in u(z), c(z), and A v (z). The full solution of the functions u, c, w, v, A v , and P has to be obtained numerically. The main method used in this paper is the spectral method employed in [68]. For a detailed discussion of the numerical method, see Appendix B. We stress that these magnetic black brane solutions are valid solutions within EMCS theory for any strength of the external magnetic field B.

Thermodynamics from gravity
After a Wick-rotation of the time direction in the metric (2.10), requiring regularity at the horizon leads to the field theory temperature, T , being given by The field theory entropy density, s, can be found from the Bekenstein-Hawking entropy formula, which gives (in units where 2κ 2 = 1 which we work in) Making use of standard recipes of gauge/gavity correspondence, in particular the relation [73] T µν = lim we can extract the energy-momentum tensor of the dual conformal field theory. Similarly, from the relation [77] J µ = lim we can obtain the expectation value of the current in the dual field theory. Given our ansatz, the energy momentum tensor and the current are given by Note that the trace of the energy momentum tensor is T µ µ = −B 2 /2.

Thermodynamics from field theory and matching to gravity data
Eqs. (2.19) and (2.20) relate expectation values (or one-point functions) of the dual field theory energy momentum tensor and axial current on the left-hand side (LHS), to combinations of coefficients (u 4 , w 4 , c 4 , ρ, p 1 ) in the boundary expansion (2.13) of the metric and gauge field in the gravitational EMCS theory on the righ-hand side (RHS). However, each of these coefficients can also be related via gauge/gravity correspondence to physical quantities in the dual field theory by matching to the purely field theory expression of the RHS. 5 From a completely field-theoretic analysis in thermodynamic equilibrium, taking into account a strong background magnetic field of O(∂ 0 ) in the derivative expansion [30], and a chiral anomaly [78] (see Appendix A for more details), the energy-momentum tensor and the axial current are given by where a subscript "EFT" (for effective field theory) indicates the quantity is obtained from the field theory directly and not from the dual gravitational EMCS theory. A subscript "0" on the energy density, , the pressure, P , and the charge density, n, indicates the quantity is evaluated in the thermodynamic equilibrium state. Similarly, the superscript "(0)" on the chiral transport coefficients, ξ V,B = ξ V,B (µ, T ), indicates all thermodynamic quantities involved are evaluated in the thermodynamic equilibrium. Note that, at strong magnetic fields, the equilibrium partition function can depend on T , µ and on B 2 [30]. Therefore, also the equilibrium quantities 0 , P 0 , and n 0 are in general functions of T , µ, and B 2 . Note that in a particular hydrodynamic frame (the thermodynamic frame [79,80]), ξ V,B are analytically known to be related to the anomaly coefficients: Here,C is an undetermined coefficient related to the mixed gauge-gravitational anomaly [81][82][83][84], which is suppressed in the large-N limit when derived from string theory setups. Note also that the trace of the energy momentum tensor is in general given by T EFTµ µ = α 0 B 2 for any conformal field theory (such as ours) in an external B field in flat spacetime, and the contribution on the RHS is the standard trace anomaly arising from the external field, with a coefficient α 0 , see e.g. [50]. In principle, our expression for the energy momentum tensor (2.21) receives derivative corrections. However, in our case T , µ, and B, as well as the field theory metric are all constant in space and also time-independent. Hence, derivative corrections vanish [30] in this particular equilibrium state. The trace of the energy momentum tensor is given by T EFTµ which can be understood as a defining relation for the magnetic susceptibility coefficient χ BB . Further, we are lead to identify those components arising entirely from the anomaly: For a detailed study of the identifications (2.25), see Sec. 2 and Appendix C, both in [68]. The remaining equilibrium quantities can also be related to the holographic near boundary data: where the right hand side is exact, and the left hand side receives no derivative corrections because our equilibrium state has T , µ, and B which are constant in space and time, i.e. derivatives on these quantities vanish.

Fluctuations
We consider fluctuations of the metric, g mn (z, x µ ), and the gauge field, A m (z, x µ ), on a fixed background solution, which we denote byḡ mn (z) andĀ m (z) respectively. To derive the fluctuation equations, we write the metric and the gauge field as sums of a background and a fluctuation part: and we expand the equations of motion to first order in ε. It is more convenient to work in momentum space, and so we perform a Fourier transformation on the fluctuations, h mn and a m , along the spacetime coordinates of the dual field theory: where k µ x µ = −ωt + k · x. For notation simplicity, we will drop the tilde on the Fourier transformed fields from now on. The presence of the magnetic field B defines a preferred direction. For simplicity, we shall consider only the case of the momentum k being aligned with B. We can then choose a coordinate system in which k is orthogonal to the (x 1 , x 2 )-plane, and the SO(2) symmetry of the background is preserved. If k were not aligned with B, the SO(2) would be broken by the fluctuations, and they would all be coupled together in the fluctuation equations.
Given that the SO(2) stays unbroken, fluctuations may be classfied in accordance to how they transform under it. The metric and gauge fluctuations transform as helicity-2, helicity-1, and helicity-0 modes under rotations about the x 3 -axis: The equation of motions for modes of different helicities decouple, thus each helicity sector can be treated independently.
In order to consider only the physical modes of the system, we have to fix the gauge freedom. To do so, we choose a gauge where a z = 0 and h mz = 0. The equations of motion for these fields then correspond to constraints, which are first order ordinary differential equations (ODEs). We then expect constraints to arise from the following modes: Consider the helicity-2 sector. The two fluctuations h 12 and h 11 − h 22 are decoupled from each other, and their equations of motion are identical. Next, the fluctuations within the helicity-1 sector can be further decoupled. In particular, the physical modes split into helicity-1 ± subsectors consisting of modes 3) The ± subsectors are decoupled from each other, and fluctuations of each subsector satisfy three second-order differential equations and one constraint equation. Lastly, there are six physical modes, h tt , h t3 , h 33 , h 11 + h 22 , a t , a 3 , in the helicity-0 sector, and they satisfy six second order differential equations and four constraint equations.
In each of the three helicity sectors, the set of fluctuations is invariant under the following two transformations: Under the transformation R 1 , all other parameters such as µ, T , and k, as well as the background fields are remain unchanged. The equations of the helicity-2 and helicity-0 sectors are invariant under R 1 , while the equations of the helicity-1 ± subsectors are mapped into each other. Under R 2 , the sign flips in the background functions c(z) and P (z) as can be seen from their near boundary expansions in (2.13) given (2.25); all fluctuations with one leg in the x 3 -direction also flip their signs. The equations of the helicity-2 and helicity-0 sectors then stay invariant under R 2 (up to an overall sign flip), while the helicity-1 ± subsectors are again mapped into each other.

Quasinormal modes and numerical details
QNMs are solutions to linearized fluctuation equations about the charged magnetic brane background here, subject to specific boundary conditions. Since QNMs correspond to poles of the retarded Greens functions in the dual field theory [55], incoming boundary conditions are imposed at the horizon of the brane. With the black brane solution written in Eddington-Finkelstein coordinates, incoming boundary conditions are imposed by requiring regularity at the horizon. Because QNMs do not source any dual operators, at the conformal boundary we have to set the non-normalizable modes to zero. In order to find the QNMs, we have to find the QNM frequency ω for which there is a non-trivial regular solution to the fluctuations equations, subject to the boundary conditions mentioned above. The problem can be recast into a generalized eigenvalue problem for ω. 6,7 In particular, in each helicity sector the resulting second order differential equations may be schematically written as where A and B are differential operators involving the background fields,ḡ mn andā m , derivatives with respect to z, and the momentum k.
The generalized eigenvalue problem can be solved numerically by using spectral methods, in which the differential operators, A and B, are represented by matrices,Â andB. The QNM frequencies are then the generalized eigenvalues associated with these differential matrices, and the QNM functions the corresponding eigenvectors. Below, we will use the QNM functions to check explicitly whether the constraint equations are satisfied in the helicity-1 and helicity-0 sectors. For more details, we refer to Appendix B.

Quasinormal mode results
We determined the frequencies ω(k) for the QNMs in all three helicity sectors. Of particular interest are QNMs corresponding to hydrodynamic modes which satisfy ω(k) → 0 for k → 0 and modes corresponding to long-lived quasi-particles with Im(ω) being small.
As discussed above, the sets of fluctuation equations are invariant under the transformations R 1,2 given in Eq. (3.4). We can use these two transformations to restrict ourselves to only positive values of γ and B. The momentum k of the QNM is taken to be real.
Moreover, since we consider a conformal field theory, instead of considering T , µ, k and B separately, we use dimensionless quantities by normalizing µ, k and B to appropriate powers of temperature, i.e.μ = µ/T andB = B/T 2 as well ask = k/T. Sometimes, it is also convenient to normalize the dimensionful quantities to the chemical potential µ by introducingT = T /µ = 1/μ,B = B/µ 2 andk = k/µ.

Helicity-2 sector under SO(2) rotations
First, let us consider the helicity-2 sector. In order to find the QNM frequencies, we rewrite the ordinary differential equation for h 12 in terms of a generalized eigenvalue problem and represent the derivatives and background fields using spectral methods. It turned out that we had to use of order N = 100 Chebyshev polynomials in order to get good convergence 8 of the QNM frequencies. How many of the QNMs we can trust depends highly on the valuesB andμ characterizing the background as well ask which specifies the momentum of the QNM.
For example, in Fig. 1  are shown in Appendix B. In particular, we conclude that we can trust the values of the QNM frequencies up to 10 −7 . For large magnetic fields, the QNMs are more sensitive tok, and some of them approach the real frequency axis. Moreover, we see that in the extremal limit,μ → ∞ (i.e. T → 0), the QNMs coalesce along the imaginary axis, presumably forming a branch cut for the exact extremal case, as was discussed already in [51]. This would be analogous to the near-extremal AdS 4 case studied thoroughly in [86,87].
Moreover, we observe that some of the QNMs approach the imaginary axis for very large magnetic fieldsB. This holds at zero as well as at finite chemical potential µ. In Fig. 2 we show two of these modes for magnetic branes, 9 i.e. with µ = 0. 8 The convergence is discussed in Appendix B. In order to improve accuracy, we use a mapping z → z 2 .
This mapping improves the convergence of the QNM frequencies and avoids the observed oscillatory behavior of the QNM frequency as a function of N . See Appendix B for more details. 9 This result is in agreement with [51], in particular Fig. 12, where their numerical results were inconclusive at those large magnetic field values.

Helicity-1 sector under SO(2) rotations
Next we consider the helicity-1 sector which is of particular interest since the Chern-Simons coupling constant enters there explicitly in the fluctuation equations. As explained above, the helicity-1 sector contains two decoupled subsectors, denoted by helicity-1 + and helicity-1 − . Both sectors consists of three second order ordinary differential equations as well as one constraint. We reformulate the three second order equations in terms of a generalized eigenvalue problem.
We show in Fig. 3 the lowest lying QNMs forμ = 10 andB = 65 (corresponding toT = 0.1 andB = 0.65). The momentum varies overk ∈ [0, 20] and we have chosen equidistant values fork. In Fig. 3 we find in each sector one QNM which does not move withk and is located at ω = −2πinT with n being an integer. These are fake QNMs satisfying the three coupled second order differential equations, but not the constraint. See Appendix B for a discussion of fake QNMs. From now on, we will discard these modes.

Hydrodynamic modes
Let us first discuss the hydrodynamic modes in the helicity-1 sector. For vanishing magnetic field, B = 0, we find two ungapped hydrodynamic QNMs, which can be identified as the momentum diffusion modes. In presence of the anomaly a term of third order in momentum k affects the momentum diffusion as was already discovered in holographic models [57,58], whereas we clarify the field theoretic origin of this term in Appendix A. In the case of B = 0, the QNMs are gapped, i.e. their dispersion relations satisfy ω(k) = 0 for k → 0. However, it is remarkable that we are still able to find agreement with hydrodynamics, at least for small magnetic fields B and counting B as first order in the derivatives, i.e.
as we explain below. In particular, analyzing the numerical QNMs for small k 1, we can fit our results to the hydrodynamic predictions.
For the numerical data, we use backgrounds with fixed temperatureT = 0.2 and magnetic field B and scan the momentum k from zero to one. We then perform a polynomial fit to the real-and imaginary part of the QNMs respectively. The upper limit for the fit range of k, denoted by k max , is chosen such that the coefficients of the fit polynomials do not change when lowering k max .  We take the helicity-1 − sector to be exemplary for helicity-1. We split our analysis in real and imaginary part of the hydrodynamics and numerical results for QNMs, respectively. In Fig. 4 we show an example for a polynomial fit to the imaginary part of the dispersion relation of the lowest QNM in the helicity-1 − sector forT = 0.2 and B = 1. For small k, the imaginary part of the dispersion relation is fitted by the polynomial Im(ω) = −0.054 − 0.0034k − 0.18k 2 − 0.031k 3 . The coefficients of the fit are also displayed in Table 1. For larger k, the third-order polynomial fit highly deviates from the imaginary part of the QNM dispersion relation, see Fig. 4. We compare the coefficients of the fit polynomial to the hydrodynamic prediction of the dispersion relation at weak as discussed in Appendix A. In (4.1) we have included some terms of third order, but not all of them, as indicated by · · · + O(∂ 3 ). Note that the energy density 0 , the pressure P 0 , the charge density n 0 as well as transport coefficients such as the shear viscosity η and the conductivity σ enter the prediction from hydrodynamics. While the energy density and the pressure can be determined from the numerical background solutions, the transport coefficients are not known explicity for our setup. Hence, we approximate the transport coefficients for the charged magnetic brane by its values for zero magnetic field, i.e. for Reissner-Nordström black brane. In particular, we use η/s = 1/(4π) and σ = (µ 2 − 6) 2 /[4(µ 2 +3) 2 ] [57,60]. This additional approximation is justified for small enough magnetic field B.
With increasing magnetic field B we can clearly track the disagreement between hydrodynamic prediction and numerics. The breakdown for increasing B shows the limits of both of our assumptions: weak field hydrodynamics and the RN approximation for the transport coefficients.
First, we show the imaginary part of the dispersion relation in Table 1. We display the fit coefficients and the corresponding hydrodynamic predictions for different B, comparing order by order in k. In particular, we find good agreement for the zeroth, second and third order in k up to surprisingly large B. Also at nonzero B, the aforementioned k 3 -term arises in presence of the anomaly and is approximating the numerical data surprisingly well. The first order in k does not fit as expected, because we have not taken into account all contributions of O(∂ 3 ) in our hydrodynamic expansion.  Table 1. Imaginary part of the dispersion relation of the hydrodynamic mode in the helicity-1 − sector forT = 0.2 and γ = 3/2. Here we display the coefficients expected from hydrodynamics and from the polynomial fit to the imaginary part of the dispersion relation of the lowest QNM.
Second, we consider the real part of the dispersion relation of the lowest QNM in Table 2. The zeroth and first order coefficients in k are well described by our prediction from hydrodynamics. Note that these coefficients can be calculated exactly with the numeri-cal background without approximating transport coefficients by its RN values, and hence tests only the weak field approximation within our hydrodynamics. Both give very good agreement with the data up to relatively large B. The second order in k does not fit the data, again due to unknown O(∂ 3 ) contributions in the hydrodynamic approximation. The third order in k is predicted to be zero in O(∂ 2 ) hydrodynamics, while our numerical data indicates a k 3 -term increasing monotonically with the magnetic field B.

Long-lived modes
Besides the hydrodynamic modes, it is important to characterize other long-lived modes with a small imaginary part in the dispersion relation for large magnetic fields.
For γ > γ c ≈ 4, we find long-lived modes in the limit of large magnetic fields. 10 In particular, we find that the real part of the QNM dispersion relation is proportional to √ B for large magnetic fields, see Fig. 5. Moreover, for γ > γ c , the imaginary part of these QNMs stay small for large B. Hence these QNMs are alongside the hydrodynamic modes the longest-lived modes. For example, for γ = 5 > γ c , the imaginary part of this QNM at B ≈ 5 · 10 3 is smaller than |Im (ω)| < 10, while the imaginary part of all other QNMs is much larger.
In contrast to that, for γ < γ c , the imaginary part of the QNMs for these modes diverges for large B and they are not the longest lived modes in our system. However, the real part of the dispersion relation is still proportional to √ B for large magnetic fields B.

Helicity-0 sector under SO(2) rotations
Finally, we discuss the helicity-0 QNMs and their frequencies. Note that as in the helicity-1 sector, the Chern-Simons coupling γ also appears explicitly in the equations of motion of the helicity-0 sector.
In Fig. 6 we display the four lowest QNMs in the helicity-0 sector forμ = 10 and B = 65. The momentum varies fromk = 0 tok = 20, while the Chern-Simons coupling is γ = 3/2. In particular, we identify three hydrodynamic modes and one fake mode. Two of these modes originate from the well-known sound modes and the third one originates from the charge diffusion mode. Note that these three modes stay hydrodynamic modes at B = 0, i.e. they remain gapless.

Hydrodynamic modes
For the three hydrodynamic QNMs identified, we perform again a fit in k to the numerical results for the QNM dispersion relation of the form In particular, we find that one of the modes, e.g. ω 0 , is purely diffusive for B = 0, in agreement with analytical results for the Reissner-Nordström case and with the hydrodynamic predictions summarized in Appendix A. This agreement with hydrodynamics holds even for small magnetic fields B as we discuss now. Taking the hydrodynamic prediction for the diffusion constant D 0 and sound velocity v 0 , see eq. (A.23), together with the Reisser-Nordström approximation for the thermodynamic quantities such as entropy density and the enthalpy, as well as for the conductivity we get v 0 = − Bγ 4 6 + µ 2 3 + µ 2 , D 0 = 6 + µ 2 12 + 4µ 2 . The same computation can be done for the hydrodynamic modes ω ± . Using the hydrodynamic prediction (A.24), and again approximating by thermodynamic expressions of Reisser-Nordström black branes, we get which compares to real part of the first order in k fit coefficient: The holographic calculation reveals also some interesting properties for hydrodynamics at strong magnetic fields. For any value of the magnetic field, we find three hydrodynamic modes ω 0 and ω ± . Let us first consider the imaginary part of each dispersion relation, i.e. Γ ± and D 0 . While counting magnetic fields of order derivative, i.e. at weak B, we conclude in (A.24) that the attenuations are the same for the sound modes ω ± , i.e. Γ + = Γ − . This, however, is not true for strong B, as can be seen in Fig. 7. For B 1, we see that Γ + = Γ − . In particular, the figure suggests that Γ + approaches zero for large enough magnetic fields, while Γ − first approaches D 0 , before both approach zero together. From the real part of the dispersion relation we can read off the sound velocities v ± and v 0 . If we treat the magnetic field of order derivative, the velocities v + and v − satisfy v + − v − = 2c s = 2/ √ 3. This is indeed the case as we can see from Fig. 7. However, for magnetic fields B 1, we clearly see deviations from it. It appears that v − approaches v 0 in the limit of large magnetic fields. While performing the fit, we noticed that a term linear in k is sufficient to fit the real part of the frequency, and a term quadratic in k is sufficient to obtain a reliable fit for the imaginary part up to fairly large values of B ≈ 8. While we show the evolution of velocities and attenuations for γ = 3/2 < γ c in Fig. 7, we have checked that their behavior does not change significantly when we perform the same calculation for γ = 5 > γ c . One of the helicity-0 modes (the former charge diffusion mode) turns into a chiral magnetic wave mode [32] 11 with velocity v 0 and attenuation D 0 displayed from zero to large B in Fig. 7. At large magnetic field the chiral magnetic wave velocity v 0 asymptotes to -1, i.e. to the magnitude of the speed of light. This is particularly intriguing as that same chiral magnetic wave QNM shows Landau level behavior at large magnetic field, i.e. the real part of the frequency is proportional to √ B and the imaginary part asymptotes to zero. This provides evidence at strong coupling for the relation between Landau level occupation and the large B behavior of v 0 . That relation was proposed by Kharzeev and Yee in [32] based on weak coupling reasoning. This can also be understood as supporting Kharzeev and Yee's conjecture that the form of the chiral magnetic velocity is valid at arbitrary B.

Long-lived modes: Landau Levels in the helicity-0 sector
Besides the hydrodynamic modes, we find other modes which are long-lived in particular for large magnetic fields. These modes behave similarly to Landau levels, as discussed in 11 [32] works with an axial and a vector current, defining the chiral magnetic wave as an excitation that involves a coupling between axial and electric charge. In this present work we have only one current and thus only one charge, namely the axial one. Nevertheless, the hydrodynamic dispersion relations of the relevant excitations are identical as can be seen from Appendix A.
[85] for a similar model. As in the helicity-1 sector, we identify QNMs whose real part scales as √ B for large magnetic fields as can be seen from the left panel of Fig. 8. In particular, we find that for µ = 0.01, the real part of the dispersion relation of the two next-to-lowest Landau levels may be fitted by Re(ω) ∼ ±1.8 √ B and Re(ω) ∼ ±2.2 √ B for large magnetic fields B. In the left panel of Fig. 8 we investigated the k-dependence of the next-to-lowest Landau level. In particular, we fixed the magnetic field B = 7.4 and the chemical potential µ = 0.01, and plotted the real part of the dispersion relation as a function of k, the momentum along the magnetic field. The numerical data is best fitted by Re(ω) = 0.98 (5.1) 2 + k 2 , the expected k-dependence for Landau levels if we extrapolate the dispersion-relation at weak coupling.
Moreover, we investigated the influence of finite charge density on the dispersion relation of the Landau levels. Again, the characteristic √ B-behaviour of the real part of the frequency for large magnetic fields is unchanged. This is to be expected since for µ √ B the behaviour due to the magnetic field dominates. However, as displayed in Fig. 9, the chemical potential modifies the real and imaginary part of the dispersion relation for small magnetic fields.  T Figure 10. Imaginary part of next-to-lowest Landau level as a function ofB for γ = 5 > γ c .
As in the case of the helicity-1 sector, the anomaly coefficient γ plays an important role whether the modes are long-lived. For γ < γ c ≈ 4, the imaginary part seems to diverge, and hence the QNM is not long-lived. In particular, the imaginary part of the dispersion relation is fitted best by Im(ω) ∼ − √ B. However, as displayed in Fig. 10 for γ > γ c ≈ 4, the QNMs, corresponding to the next-to-lowest Landau levels, approach the real axis and hence are long-lived excitations for large magnetic fields.

Discussion
In this paper, we have investigated transport and dissipation in a thermal plasma in a magnetic field of arbitrary strength at strong coupling. This plasma state is defined within a (3+1)-dimensional field theory with chiral anomaly. In the dual gravitational calculation we evaluate and analyze quasinormal modes (QNMs) of Einstein-Maxwell-Chern-Simons theory in the asymptotically AdS charged magnetic black brane background. For a particular value of the Chern-Simons coupling, the action is a consistent truncation of type IIB supergravity, and the dual field theories are well-known, including N = 4 SYM. We have investigated the effects of having a finite temperature T , chemical potential µ, a non-zero external magnetic field B, on the QNMs. We have chosen the momentum of the QNMs, k, to be (anti-)parallel to B.
In parallel to the holographic calculation, we have computed the weak B hydrodynamic description of our system. The five (formerly) hydrodynamic poles were computed at nonzero B. We extend previous treatments keeping our derivation general and providing expressions for the various velocities and attenuations. At leading order in derivatives, we have also worked out an expression for the energy momentun tensor and axial current containing polarization effects (2.21). These hydrodynamic results are collected in Appendix A.
In a systematic study of all metric and gauge field fluctuations propagating along the magnetic field, we investigated QNMs both within and outside the hydrodynamic regime. We find that some of the QNMs are long-lived and show characteristic behavior known from Landau levels, see Figs. 8, 9 and 10. The helicity-2 QNMs do not exhibit any hydrodynamic modes. However, some of these QNMs appear to have vanishing real part in the limit of large B, see Fig. 2. There are five hydrodynamic modes in total. Two of those modes at van-ishing B are momentum diffusion modes appearing in the helicity-1 sector. At nonzero B these modes acquire a complex gap and can not be considered hydrodynamic modes according to our previous definition. However, they are still well described by the hydrodynamic dispersion relation  Fig. 7. These are hydrodynamic modes (ungapped) even at B = 0, as can be seen from the dispersion relations (A.20). One of the helicity-0 modes (the former charge diffusion mode) turns into a chiral magnetic wave mode 12 with velocity v 0 and attenuation D 0 displayed from zero to large B in Fig. 7. At large magnetic field the chiral magnetic wave velocity v 0 asymptotes to -1, i.e. to the magnitude of the speed of light. This is particularly intriguing as that same chiral magnetic wave QNM shows Landau level behavior at large magnetic field, i.e. the real part of the frequency is proportional to √ B and the imaginary part asymptotes to zero. This provides evidence at strong coupling for the relation between Landau level occupation and the large B behavior of v 0 . That relation was proposed by Kharzeev and Yee in [32] based on weak coupling reasoning. This can also be understood as supporting Kharzeev and Yee's conjecture that the form of the chiral magnetic velocity is valid at arbitrary B. In [32] a holographic and a field theory result for the chiral magnetic wave velocity, and a holographic result for its attenuation are provided. That holographic calculation was based on a probe-brane approach, which neglects backreaction by definition. In that sense, the present work extends [32] to fully backreacted geometry and a solution valid at arbitrary B. We further provide a hydrodynamic and holographic calculation of the chiral magnetic wave attenuation D 0 , as well as the velocities v ± and attenuations Γ ± of the former sound modes. Remarkably, all of these velocities obtain contributions from the anomaly, as seen explicitly e.g. in (A.23) and (A.24). This is also true for the two helicity-1 modes, see e.g. (A.17). Dispersion relations for all five hydrodynamic modes are provided in the Landau frame and in the thermodynamic frame [79,80], see Appendix A.
For the future, it would be interesting to relax the restriction of having the QNM momentum k parallel to the external magnetic field B; see [28] for a field theory discussion of a related system. This would allow many more modes to arise, and thus even richer phenomenology. It would also be interesting to have two U(1) gauge fields in the gravity theory [88], i.e. introducing an axial and a conserved vector current in the dual field theory. This is relevant for testing some predictions for chiral magnetic waves and for Weyl semimetals and their surface states, see [38][39][40][41][42]. Moreover, it would be interesting to include mixed gauge-gravity anomalies and study their effects on the QNMs, see for example [89] in which a mixed anomaly appears to lead to a phase transition. In this paper we have considered hydrodynamics in the regime of weak external gauge fields. It would be very interesting to work out the hydrodynamics in the case where the gauge fields are strong, i.e. zeroth order in a derivative expansion, building on [30] and [29]. Lastly, it would be desirable to extend our setup systematically to magnetohydrodynamics as viewed from a modern perspective [90].

A Hydrodynamics
In this appendix, we derive expressions for the location of the poles of retarded correlation functions of the energy momentum tensor, T αβ , and the axial current, J µ , in the hydrodynamic approximation. These poles are identified with the lowest lying QNMs of the corresponding gravitational fields via the gauge/gravity correspondence. Hence, this section will serve to predict the QNM frequencies at small frequency values and small momenta.

A.1 General framework and definitions
Hydrodynamics can be considered as an effective field theory, or more precisely, as an expansion of n-point functions in terms of (small) gradients of the hydrodynamic variables. The defining relations for hydrodynamics are the constitutive equations (or 1-point functions) for the energy momentum tensor and the conserved current(s), as well as the corresponding conservation equations already defined in the main text in (2.2), see for example [91,92]. The constitutive equations express the conserved and hence long-lived quantities, i.e. the energy momentum tensor and the conserved current(s), in terms of temperature T , chemical µ, and the fluid velocity, u µ . These hydrodynamic variables can be considered as fields in hydrodynamics when considered as an effective field theory. As such, they are only defined modulo field redefinitions leaving the physical quantities invariant. Fixing this freedom of field redefinitions is fixing a "hydrodynamic frame". A well known example is the Landau frame in which the heat current vanishes, while we will be mostly working in the thermodynamic frame [79,80] defined below.
We consider here a system in the presence of external sources, viz. the external gauge field, A µ , and the external metric, g µν . We are interested here in a non-trivial gauge field background, A µ = (µ, −x 2 B/2, x 1 B/2, 0), with a non-vanishing chemical potential µ, and a magnetic field B in the x 3 -direction, plus a first order correction a µ . However, for the metric background we take it to be Minkowskian, η µν = diag (−1, 1, 1, 1), plus a first order correction h µν .
The system will respond to the metric and gauge sources with corrections to the equilibrium values for all the hydrodynamic variables , (A.1) with the expansion parameter ε (not to be confused with the energy density ). Here we have choosen the particular case in which the momentum of the response is parallel to the magnetic field, i.e. k || B, where the vector arrow denotes the spatial part of the four-vectors. 13 Those corrections in the hydrodynamic variables will cause corrections in the thermodynamic quantities , P , n, which we expand as where here and in the remainder of this appendix partial derivatives with respect to T are evaluated at fixed µ and vice versa, unless otherwise stated. Below, we will make frequent use of the following thermodynamic relations: Quantities of relevance to the constitutive equations used in the following subsections are the projector, electric field, the magnetic field and the vorticity: While the chiral transport coefficients in the thermodynamic frame are defined by where the T 2 terms can be neglected in the large N limit as they originate from supressed gauge-gravitational anomalies [81][82][83][84].

A.2 Strong field thermodynamics (B ∼ O(1))
In this subsection, we consider strong magnetic field thermodynamics, more precisely zeroth order in derivatives hydrodynamics for time-independent quantities. A full hydrodynamic treatment of the strong magnetic field case would be interesting but is beyond the scope of this work. 14 By "strong field" we mean a background magnetic field of zeroth order in derivatives, i.e. the field strength F ∼ B ∼ O(1), with B being the magnitude of the magnetic field. When a system is placed in a strong external field, strong field effects need to be taken into account, for example the polarization or magnetization of the medium [30,31]. Such polarization effects to zeroth order in a derivative expansion can be parametrized utilizing a polarization tensor M µν = 2∂P/∂F µν . Furthermore, effects of the chiral anomaly can contribute to thermodynamic quantities at zeroth order in derivatives [78]. Together, 13 Another interesting choice would be to consider the momentum of the response to be perpendicular to the magnetic field, see [28] for such a choice in a system with an axial and a vector current. 14 For an attempt at this, see [29].
these effects at a strong external magnetic field B can be expressed in terms of equilibrium constitutive equations for the energy momentum tensor and the axial current Note that all quantities above have to be evaluated at zeroth order in derivatives, e.g. u µ = (1, 0, 0, 0). When the magnetic field is chosen to point in the x 3 -direction, then these relations directly reduce to the equilibrium energy momentum tensor (2.21) and the equilibrium axial current (2.22). Note that in the Landau frame the heat current vanishes, q µ = 0, hence the off-diagonal components in the energy momentum tensor vanish in Landau frame. Physically this simply means that in that frame the observer is traveling with the heat within the fluid. The interesting point is actually that such a heat current exists in equilibrium at strong magnetic fields and in the presence of an anomaly.

A.3 Weak field hydrodynamics (B ∼ O(∂))
By "weak field" we refer to a background magnetic field which is of first order in derivatives, i.e. the field strength is F ∼ B ∼ O(∂). In other words, here we consider a background magnetic field which is of the same order as the spatial derivatives, and hence B ∼ k 3 = k. So, for example, any term that contains one power of the magnetic field and one spatial derivative is of second order: B∂ 3 (. . . ) ∼ O(∂ 2 ). We start from the constitutive equations to first order in derivatives T µν = u µ u ν + P ∆ µν + u µ q ν + u ν q µ + τ µν , (A.12) with u µ τ µν = 0, u µ ν µ = 0, and u µ q µ = 0. In order to determine the location of the poles of the hydrodynamic correlation functions, we have inserted the constitutive equations (A.12) into the conservation equations (2.2) and linearized in ε. We have then solved that system of linear equations for the corrections to the hydrodynamic variables, namely for T 1 , µ 1 , and for the three spatial components of u µ 1 , while the time component of u µ 1 is fixed by normalization u µ u µ = −1. The result of this exercise is provided in the following subsections sorted according to the helicity of the corrections and sources with respect to the SO(2) rotation symmetry in the (x 1 , x 2 )-plane. 15

A.3.1 Poles in helicity-1 sector
There are two poles in the correlation functions of the helicity-1 components of the energy momentum tensor T µν and the axial current J µ , located at which originates from the two solutions when expanded in k and B. The second line in (A.17) is third order in B and k, and it will receive corrections from second order contributions to the constitutive equations. It is remarkable that in the limit of B = 0, explicit checks confirm that the k 3 -term in the pole does not appear to receive any corrections from terms of second order in the constitutive relations in the Reissner-Nordström case. In that case (A.17) reduces to which is the standard charge diffusion k 2 term plus a contribution from the anomaly. This expression is in agreement with the holographic result [57,58]. Note that this k 3 -term in equation (A.17) is not present in Landau frame at first order in the derivative expansion as it originates from the heat current contribution q µ in the energy momentum tensor in thermodynamic frame. It was shown in [93] that this k 3 term appears when working to second order in the constitutive equations in Landau frame.

A.3.2 Poles in helicity-0 sector
For convenience, we define first some useful notation: w 0 = 0 + P 0 is the enthalpy density, s 0 = s 0 /n 0 is the entropy per particle, c 2 s = (∂P/∂ ) s is the speed of sound, c n = T 0 (∂s/∂T ) n andc P = T 0 (∂s/∂T ) P are the specific heats at constant density and pressure respectively 16 , and α P = −(1/n 0 )(∂n/∂T ) P is the thermal expansivity at constant pressure.
There are three poles in the helicity-0 sector: For the generalized diffusion pole, which we refer to as a chiral magnetic wave pole from now on, we have where D 0 is the generalized diffusion coefficient. For the generalized sound poles, we have In the RN limit, B = 0 and ζ = 0 17 . We have thus 26) In particular, in the RN limit we have the relatioñ This relation actually holds more generally than just in the RN case. It is true whenever P 0 ∝ 0 , and this can be seen as follows. By using Eq. (A.6), it can be shown that Then by Eq. (A.5), we have for constant P 0 and thus 0 and thus the desired result.
In the limit of vanishing charge density, n 0 = 0, the velocity and attenuation of the generalized diffusion pole are given by For the sound poles at n 0 = 0, with the additional simplifying assumption of P 0 ∝ 0 , velocities and attenuations are given by Note that when n 0 = 0, d = T ds and dP = sdT , hence c 2 s = ∂P/∂ = (∂P/∂T )/(∂ /∂T ).

A.3.3 Comparison to Landau frame results
A field theory calculation [27] similar to ours has been carried out while this paper was in preparation. Whereas our holographic model prompted us to work in the thermodynamic frame, Landau frame was chosen in [27]. Here we repeat our calculation in the Landau frame for comparison, and we point out differences between the two hydrodynamic frames. 18 In the Landau frame, the chiral conductivities are given by Note the relation between the anomaly coefficients defined in the two frames: The chiral magnetic wave pole is given by and the generalized sound poles by Note that only the "velocities" v L 0 and v L ± are different compared to their counterparts in the thermodynamic frame, while the "diffusion coefficient" and "attenuation coefficients" are identical in the two frames. Lastly, the RN limit in the Landau frame is the same as in the thermodynamic frame, while in the n 0 = 0 limit, We point out that the expressions for v L 0 , v L ± , D 0 and Γ ± agree with the results of [27]. In particular, when n 0 = 0, v 0 and D 0 from (A.37) and (A.30) respectively are identical to v χ and D χ defined in (42) of [27] 19 . We can also compare the expressions for v χ and D χ at n 0 = 0 in [27] to our Landau frame results, (A. 35). Note that our calculation is more 18 We find full agreement between those results provided explicitly in [27] and ours. We find only partial agreement between our results and the results of [26]. In particular for the spin-0 modes at vanishing chemical potential and momentum we were able to find agreement with what the authors call ω1, ω2, ω3. The other results reported in [26] we can not reproduce unless we make further assumptions, e.g. about the relation between the energy density and the pressure. 19 Note the definition in [27] is given at vanishing chemical potential, whereas we allow for µ0 = 0.
general, as we do not impose any restrictions such as δp = 0, which was imposed in [27]. The comparison shows agreement between the velocity v L 0 in our (A. 35) and v T χ given in (45) of [27]. We also provide the attenuation coefficient, D L 0 , in (A.35) and (A.23), which is left undetermined as "D T " in [27].
The particular collective excitation with the velocity v L 0 given in (A.35) was coined "thermal chiral magnetic wave" in [27]. But the classification according to the SO(2) rotation symmetry, seems to support the claim that (A.35) is merely a generalized version of the chiral magnetic wave [32]. In addition to a background magnetic field, a background vorticity has been considered in [27], and additional effects such as a "thermal chiral vortical wave" and two modified sound modes with a vorticity-dependent velocity were reported. Certainly these three modes bear intriguing effects, we would, however, still refer to them simply as the hydrodynamic modes of the helicity-0 sector, one of them being the chiral magnetic wave mode.
The helicity-1 poles in Landau frame are given by It is important to recall that we are working with a single axial current and an energy momentum tensor here, in contrast to QCD, which has an axial and a vector current. Hence, there are more hydrodynamic modes to be expected in QCD due to the interplay between its axial current and its vector current. The weak field hydrodynamic description of such a system has been considered in [28].

B Numerical methods. Convergence tests
The equations of motion for the background and the equations of motion for the fluctuations on such a background are given by (2.7) and (2.8) using the ansatz (2.10),(2.12) or (3.1),(3.2) respectively. The equations obtained are solved numerically with a spectral method. In this section, we discuss the numerics for the nonlinear background equations as well as the nonlinear generalized eigenvalue problem for the QNMs in detail. The whole numerics is set up in Mathematica. This enforces to use several performance optimizations, which we will discuss for the background and the QNMs.
Spectral methods as a high accuracy method have the benefit of exponential convergence for analytic solutions. Yet the background magnetic field B introduces logarithmic terms, giving an algebraic convergence rate. Additionally spectral methods allow to include the singular points on the boundary, that appear in holographic models. Therefore the background equations and the eigenvalue problem can be solved on the whole domain, including both the conformal boundary and the horizon.

B.1 Background
The equations of motion (2.7) and (2.8) for the background are given by six nonlinear ordinary differential equations and one constraint equation for the six unknown functions, u(z), v(z), w(z), c(z) and A v (z), P (z). Expanding the unknown functions in Chebyshev Polynomials evaluated on a Chebyshev-Lobatto grid gives a nonlinear algebraic problem. This nonlinear algebraic problem is solved by means of a Newton-Raphson Method providing an initial guess. The constraint is incorporated in the boundary conditions. For this project we also subtract the leading logarithms as in [68] and calculate all background data with 100 gridpoints, which we justify in the following discussion on QNMs. A detailed description of both the nonlinear system and the subtraction of leading logarithms can be found in the more general setup of [68].
To be able to use large gridsizes and to get high precision results we use the following advanced techniques in our numerical code. In the first step we carefully make sure to keep the equations of motion in a very short form, in particular after extracting the logarithms, terms like e.g. u(z) 2 = 1 + z 4 u 4 + z 4 ln(z) B 2 6 2 should not be expanded, which can be the case when using simplify.
We than replace the functions f (z) = {u(z), v(z), w(z), c(z), A v (z), P (z)} by symbolic quantities by f = {u, v, w, c, A v , P }. Instead of replacing derivatives of f by their spectral representation, we also replace f (z) by df = {du, dv, dw, dc, dA v , dP } and their second derivatives by d 2 f . This has the advantage, that numerical values for the derivatives are calculated only once in each step using spectral differentiation matrices, such that e.g. the numerical value of the symbol du is calculated once and then replaced everywhere it appears in the equations of motion or in the Jacobian. The replacement is done with a faster Dispatch-rule in Mathematica instead of a trivial replacement-rule.
In addition we calculate the Jacobian only in the first step of the Newton-Raphson iteration symbolically and store it for the further iterations. The Jacobian is calculated using the chain rule, i.e. calculating where ∂(df ) l ∂(f ) j and ∂(d 2 f ) l ∂(f ) j are given by a block matrix of the spectral differentiation matrix. Moreover to calculate backgrounds with constantT and varyingB or constantB and varyingT , we replace the boundary condition A v (0) = µ at the conformal boundary by a condition for the temperature u (1) = 4πµT . This is equivalent to promoting µ to an extra parameter and implementing u (1) = 4πµT as an extra condition.

B.2 Quasinormal Modes
The fluctuation equations are discretized in the same way, where we now have to insert the numerical values for the previously calculated background. We use the second order Despite being in Eddington Finkelstein coordinates, the helicity-0 sector has a quadratic ω dependence, giving a nonlinear eigenvalue problem. This arises from additional ω contributions in the determinant of the metric. The helicity-1 and helicity-2 sectors are still linear in ω. The quadratic eigenvalue problem can be reduced to a linear eigenvalue problem by doubling the number of fields, introducing h new ij (z) = ω h ij (z). In addition to the numerical techniques for the background, we replace the derivatives of the background and the fluctuations again symbolically and calculate their values only once. Moreover we use the background equations of motion to eliminate all second derivatives of background fields in our QNM equations. To interpolate we use the fast and more accurate Clenshaw algorithm. The linear eigenvalue problems are solved with Mathematica Eigenvalue/Eigensystem, to get the QNM and the eigenfunctions respectively. All QNM data in this project are calculated also with 100 gridpoints and 60 digit precision, which means solving a 100 × 100 matrix for helicty-2, two 300 × 300 matrices for helicty-1 and a 800 × 800 matrix for helicty-0 sector.
By plugging back the eigenvalue and the eigenfunctions into the equation of motion, we check if we have a true solution of the full system including constraints. In our analysis QNM appear, that converge and fulfil the equations of motion, but not the constraint(s). It happens that those QNM additionally do not change with k. Accordingly these QNM are identified as numerical artefacts and are rejected for the analysis. We refer them as fake QNM which also show up in the continued fraction method [94,95]. In particular, the fake QNM may be traced back to a degeneracy of the ingoing/outgoing horizon conditions in the Poicare chart. At the horizon, ingoing and outgoing modes have a term of the form (1 − z) ±iω/u 1 , which for ω = −i n u 1 /2, with n ∈ N, is degnerate from the power series in 1 − u. Taking into account u 1 = 4πT , the fake QNMs lie at ω = −2πi T n with n ∈ N, as expected.

B.3 Helicity-2 sector
The helicity-2 sector is described by two decoupled equations of motion. The convergence in helicity-2 sector is slow due to logarithms in the background and logarithms in the QNM eigenfunctions. It further shows an even/odd gridsize oscillation, meaning that the QNM value for even/odd gridsizes approaches the correct value from above and below respectively.  Table 3. Constraint for k = 0 andk = 20,B = 65,T = 0.1 in the helicity-1 − sector Therefore we have to introduce a z → z 2 mapping. This moves the logarithms to higher orders z n log(z) → 2 z 2n log(z) , (B.6) which improves the convergence. Interestingly we find that this allows to use simplified boundary conditions at z = 0. In principle one has to specify the new boundary condition f (z) = 0 for all fields f as well as the boundary condition before the mapping. In our case the condition f (z) = 0 was sufficient to reproduce the results before the mapping, since our boundary conditions are solely given by the equations of motion. All calculations, and in particular the convergence plots shown below, are obtained with this mapping.
With the z 2 mapping we already obtain good results with 40 gridpoints, see 11.

B.4 Helicity-1 and helicity-0 sector
Helicity-1 is given by two decoupled subsectors, helicity-1 ± s. Each subsector consists of three equations of motion and one constraint and has to be solved separately. We exemplary show the lowest eigenmodes and the error in the constraint in table 3 for helicity-1 − . The error is obtained by plugging the corresponding eigenvector, which is the numerical solution for the fluctuations evaluated at the gridpoints, into the constraint equations and taking the norm. For example the QNM 0. − 2.20477 i in table 3 does not fulfil the constraint and does also not change fromk = 0 tok = 20. This mode is identified as a fake QNM. This mode also appears in helicty-1 + . Despite being a fake mode, it converges nicely, as shown in fig. 12 In fig 13 we show the slower convergence for largerB andk. Comparing the convergence at the same point in parameter spaceB = 65,T = 0.1 fork = 0, shown in fig.  12, andk = 20, shown in fig. 13, reveals the slowing down of the convergence for largerk, which requires to use large gridsizes. The helicity-0 sector is given by six equations of motion and four constraint. It is in addition quadratically in ω, therefore we introduce two additional fields,known as doubling trick, to get a linear system again. Convergence does not pose a problem in the helicity-0 sector, however calculations are time consuming compared to helicity-1 and helicity-2.