Collective diffusion and quantum chaos in holography

We define a particular combination of charge and heat currents that is decoupled with the heat current. This `heat-decoupled' (HD) current can be transported by diffusion at long distances, when some thermo-electric conductivities and susceptibilities satisfy a simple condition. Using the diffusion condition together with the Kelvin formula, we show that the HD diffusivity can be same as the charge diffusivity and also the heat diffusivity. We illustrate that such mechanism is implemented in a strongly coupled field theory, which is dual to a Lifshitz gravity with the dynamical critical index z=2. In particular, it is exhibited that both charge and heat diffusivities build the relationship to the quantum chaos. Moreover, we study the HD diffusivity without imposing the diffusion condition. In some homogeneous holographic lattices, it is found that the diffusivity/chaos relation holds independently of any parameters, including the strength of momentum relaxation, chemical potential, or temperature. We also show a counter example of the relation and discuss its limited universality.


Introduction
One of the most mysterious phenomenon in condensed matters is the ubiquitous appearance of linear in temperature resistivity. In the materials with such 'strange-metal' hallmark, the quasiparticle picture is not applicable because the resistivity can cross the Mott-Ioffe-Regel limit [1,2]. It has been long suggested [3,4] that the transport in strange metals is controlled by the 'Planckian dissipation'-the temperature in units of time through Planck's constant: τ P ∼ /(k B T ), but the explicit framework has not been built up.
In [5], Hartnoll has pointed out that when the momentum decays quickly, the collective diffusion of charge and energy controls the strange-metal transport. Inspired by the putative bound on the ratio of shear viscosity to entropy density [6], he proposed that the eigenvalues of diffusive matrix are lower bounded by the Planckian dissipation and Fermi velocity.
Translating by the Einstein relation and neglecting some thermoelectric effect, he argued that the saturation of the diffusion bound may be responsible for the ubiquity of hightemperature regimes in metals with T -linear resistivity.
The theory of incoherent metals is insightful, but there are several issues which deserve to be explored. First, the diffusion bound can be violated in various situations [7][8][9][10].
Second, some strange metals are relatively clean and their thermoelectric effect may not be especially small. Third, the characteristic velocity v P is taken as the Fermi velocity v F , which may not be well defined in strongly coupled systems. To address the last issue, the butterfly velocity v B that quantifies the speed of chaos propagation has been proposed as a natural replacement [11,12]. Actually, an interesting relation has been found by holography between the charge diffusivity and the quantum chaos which is characterized by v B and τ L . Here τ L is the Lyapunov time that is expected to indicate the Planckian dissipation in non-quasiparticle systems [13][14][15][16][17]. However, in addition to the non-universal prefactors that occur in special cases [18,19], the application of charge diffusivity/chaos relation is limited since the particle-hole symmetry must be imposed. As a result, the relation does not appear in the normal state of cuprates. On the other hand, the relation between the thermal diffusivity and chaos seems to be more robust [12,[20][21][22][23][24], but it cannot be directly translated into the statement of resistivity.
Nevertheless, by a recent local optical measurements of thermal diffusivity on underdoped YBCO (an ortho-II YBa 2 Cu 3 O 6.60 and an ortho-III YBa 2 Cu 3 O 6.75 ) crystals, it was found that the thermal anisotropy is almost identical to the value of the electrical resistivity anisotropy and starts to decrease sharply below the charge order transition [25]. The experiment was interpreted that the non-quasiparticle transport is dominated by the collective diffusion of electron-phonon 'soup'. In particular, it led to the conjecture that both charge and heat diffusivities saturate the proposed bounds [25,26].
To understand how the collective diffusion could be relevant to the intrinsic mechanism for robust metallic transport that is complicated by the extrinsic processes for momentum relaxation, it would be promising to search hints in clean systems. Indeed, in a conformal field theory (CFT) with charge doping, Davison, Goutéraux and Hartnoll (DGH) have isolated a diffusive mode by hydrodynamics [27,28], which is carried by a particular combination of the charge and heat currents. The DGH mode can likely be considered intrinsic, because the DGH current is decoupled with momentum in the conformal fluid and the DGH conductivity is universal for some clean holographic theories [29]. However, the DGH mode has not been studied when the translation symmetry is broken, partially because a simple proposal to include slow momentum relaxation in hydrodynamics [30] is inconsistent with the holographic models [28,31] and the fast momentum relaxation invalidates the hydrodynamics essentially.
In this paper, we will explore the collective diffusion of charge and energy and its relationship to quantum chaos in terms of the gauge/gravity duality with momentum relaxation.
At the beginning, we will show that there is a universal bound for thermo-electric transport.
The bound is trivial, except that its saturation indicates the decoupling between the heat current and a particular combination of charge and heat currents. The 'heat-decoupled' (HD) current is nothing but the DGH current in most of homogeneous holographic lattices. Without relying on the hydrodynamics, we will verify that the HD current can be transported by diffusion at long distances, when some thermo-electric conductivities and susceptibilities satisfy a simple condition. Then we will study the HD mode by combining the diffusion condition and the Kelvin formula [32]. Note that the Kelvin formula arises if the static limit is taken before the thermodynamic limit in the evaluation of Onsager coefficients [32]. It provides a good approximate (sometimes even exact) expression of the thermopower in various contexts including strongly correlated systems, such as the fractional quantum Hall states [32], high temperature superconductors [32][33][34], and the homogeneous Sachdev-Ye-Kitaev model [19]. As a result, we will point out that not only the heat diffusivity but also the charge diffusivity can be identified with the HD diffusivity.
Furthermore, in a Lifshitz gravity with the dynamical critical index z = 2, we will illustrate that the diffusion condition and the Kelvin formula can be realized. Meanwhile, the HD diffusivity exhibits the relationship to the quantum chaos and is equal to both charge and heat diffusivities.
In the thermo-electric systems that we care about, the fluctuations of charge and energy are coupled and hence they do not satisfy separate diffusion equations in general. Accordingly, the charge and energy diffusivities do not correspond to any eigenvalues of the diffusivity matrix. However, they still can be well defined (by their Einstein relation) and would imply interesting physics. This inspires us to study the diffusion constant of HD mode, but without imposing any diffusion condition. In other words, the HD mode is no longer purely diffusive. We will exhibit that the diffusivity/chaos relation is respected exactly in various homogeneous holographic lattices. This is different with previous results in references, where the diffusivity/chaos relation always requires certain limits on the parameters, including the particle-hole symmetry [9,11], strong momentum relaxation [12,21,22], or low temperature [20,23]. However, the relation is not universal and we will show a counter example. The limited universality and implication will be discussed.

Heat-decoupled current
The thermo-electric transport is characterized by three conductivities: electrical σ, thermal κ, and thermoelectric α 1 . They reflect the response of the charge and heat currents to small gradients of temperature and chemical potential, The DGH current is defined by a particular combination of charge and heat currents [27] where ρ is the charge density and µ is the chemical potential. Now we construct a general linear combination but write it as the DGH-like form for comparison In the present,ρ is an arbitrary quantity with the dimension of charge density. The currentcurrent correlation can be calculated bȳ We reorganize the numerator of Eq. (4): In the dc limit, the last term in Eq. (5) is non-negative, so we immediately obtain a universal boundσ where σ 0 ≡ σ − T α 2 /κ denotes the electric conductivity at zero heat current. The bound is almost trivial. The only interesting point is that, due to the bound is saturated when the general current is decoupled with the heat current. Hereafter, we will focus on the HD current 2 One can find that the HD current is equal to the DGH current in most of homogeneous holographic lattices which have 3 However, this relation does not hold when the momentum relaxation is inhomogeneous [36] or involves some non-minimal coupling [10,22,37]. In these cases, the HD current is different with the DGH current and might be viewed as an extension.

Diffusive mode
Suppose that the energy and the charge are conserved. The HD current respects the con- and carries a particular combination of charge and energy where we denote as the energy density. Using the definition of conductivities (1), the continuous equation (10), and the thermodynamic identities ∇ρ = χ∇µ + ζ∇T, with the susceptibilities defined by derivatives of pressure one can prove at long distances where we have defined Following Ref. [38], the diffusion equation can be constructed, provided that C = 0, which means At this time, the diffusion constant can be written as Some remarks are in order. First, the diffusion condition (18) can be understood as a special balance between the thermoelectric conductivity α and susceptibility ζ. The thermoelectric balance is reminiscent of the Kelvin formula which can be written as σ/χ = α/ζ. The Kelvin formula can be derived by requesting that the density gradient vanishes [34]. At that time, Eq. However, they can be cancelled in D HD due to ζ/χ = (∂s/∂ρ) T , which is indeed related to the thermodynamics of black holes [20]. At last, it is instructive to clarify the relation between the HD current with the diffusion condition and the decoupled thermo-electric currents which are transported by diffusion [28,39], see Appendix A. In the following, we will reveal a general and exact relation between D HD and v 2 B τ L in various holographic models. In particular, the diffusion condition (18) and the Kelvin formula both hold in a model of Lifshitz gravity with z = 2, leading to 4 Holographic models

Einstein-Maxwell-axion
A simple holographic framework with momentum relaxation was presented in [40]. The model contains linear axions χ i along spatial directions. We consider the four-dimensional Einstein-Maxwell-axion (EMA) theory, Here the AdS radius L and the Newton constant 16πG N are set to unity. The equations of motion (EOM) derived from the action has an isotropic background solution, in which χ i = βx i with the disorder parameter β [40]. Suppose that the horizon locates at r + . The Hawking temperature and entropy density are from which two thermodynamic response functions can be derived By building up the conserved currents and invoking the sources that are linear in time, Donos and Gauntlett derived the thermo-electric conductivities analytically [35]: Using Eq. (22), Eq. (23), and ζ/χ = (∂s/∂ρ) T , the diffusion constant (15) can be calculated.
The result is simple: We need to study the quantum chaos. The onset of chaos is characterized by the Lyapunov time τ L and the butterfly velocity v B . They can be calculated by constructing a shock wave near the black hole [41][42][43][44][45]. For the EMA model, they are [11,12] Remarkably, we have found an exact diffusivity/chaos relation As a comparison, we write down the ratio between

Non-relativistic scaling
We have interest on an Einstein-Maxwell-Axion-Dilaton (EMAD) theory studied in [46,47] where all the functions of the dilaton are assumed to have the exponential form with several parameters Λ and λ i (i = 0, 1, 2, 3). This theory admits Lifshitz-like, hyperscaling violating, analytical black-brane solutions. The line element is given by where the blackness factor is We have denoted the dynamical critical index as z and the hyperscaling violating factor θ. Besides the usual electromagnetic field, there is an additional Maxwell field which is necessary for z = 1. Two Maxwell fields have the charges (q 1 , q 2 ), respectively. We will impose the physical condition that the first U(1) current is vanishing. Then the charge density is ρ = q 2 . The parameters Λ, λ i , and q 1 are all determined by z and θ. For instance, . In [47], it has been pointed out that the black-brane solution is divergent when z = 2, indicating the logarithmic behavior. The blackness factor then becomes Using the black-brane solutions, the temperature and entropy density can be written as which is effective even for z = 2. They can lead to two response functions . (34) From Ref. [46], the dc thermo-electric conductivities can be read off, To obtain the quantities in chaos, we study a shock-wave metric which is generated by releasing a particle from the boundary at x = 0. From the EOM, one can obtain the solution of h(t, x) at later time and large distance. Then we can extract the Lyapunov time τ L = 1/ (2πT ) and the butterfly velocity v B = 2πT /m, with the screening length m given by Equation (37) is same as the one in [11], indicating that the additional Maxwell field does not change the form of the screening length. Translating the Kruskal coordinates into Thus, we have the robust relation D HD = v 2 B τ L for the generic z and θ. For comparison, the ratio between D T and v 2 B τ L depends on z and θ: Let's focus on the special case with z = 2. From Eq. (34) and Eq. (35), one can see which implies that the Kelvin formula σ/χ = α/ζ holds and D HD = D T . Furthermore, we have interest on the case with θ = 0. To be clear, we rewrite the temperature and entropy density, and also invoke the chemical potential: The last equation is necessary to calculate the quantities χ, ζ, and c µ . One can find that the diffusion condition (18) is obeyed. Collecting all the results in the exceptional case with z = 2 and θ = 0, we have a series of equalities

High curvature
The higher derivative corrections appear generally in any quantum gravity theory from quantum or stringy effects. These corrections may be holographic dual to 1/N or 1/λ corrections in some gauge theories, allowing independent values of two central charges a and c. This is in contrast to the standard N =4 super Yang-Mills theory where a = c, hence likely underpinning the violation of viscosity bound. Actually, the Gauss-Bonnet (GB) correction has been treated as a dangerous source of violation for the feature that is universal in Einstein gravity [48,49].
Consider the GB correction to the EMA theory, that is whereα is the GB coupling constant. The black-brane solution with χ i = βx i has been derived in [50]. The temperature and entropy density can be written as Here is the square of the effective AdS radius. Taking the derivatives of the entropy density, the two response functions read: We also need the thermo-electric conductivities [50] Using Eq. (46) and Eq. (47), the diffusion constant (15) can be calculated as The butterfly effect for GB gravity has been studied in [42]. We will repeat the calculation but involve the linear axions. Release a particle from the boundary and consider the five-dimensional shock-wave metric Using the EOM, one can derive the solution of h(t, x) at later time and large distance, which in turn gives the Lyapunov time τ L = 1/ (2πT ) and the butterfly velocity v B = 2πT /m, with the screening length m given by Under the Schwarzschild coordinates, one can read Then we build the diffusivity/chaos relation D HD = v 2 B τ L again. The ratio between D T and v 2 B τ L can be given by

Anisotropy
Spatially anisotropic black-brane is the first example that yields the violation of viscosity bound in Einstein gravity [51]. A careful examination of the diffusivity/chaos relation in such anisotropic model is worthwhile. In [52,53], a R-charged version of the spatially anisotropic black-brane solution was derived via nonlinear Kaluza-Klein reduction of type-IIB supergravity to five dimensions, which leads to the presence of an Abelian field in the action. The introduction of the U (1) gauge field breaks the SO(6) symmetry and thus leads to the excitations of the Kaluza-Klein modes. Consider the reduced theory that is described by a five-dimensional EMAD action There is an anisotropic solution which has been found numerically in [52]. Note that the strength of anisotropy and disorder are both characterized by the single parameter β. The analytical solution with small β has also been obtained, see Appendix D in [52]. We write them as follows: where u ≡ 1/r. The three functions (F, B, φ) are given by where (F 2 , B 2 , φ 2 ) can be read off from Eq. (135) in [52] with the replacement for our convention a, u H , q 2 → β, u + , ρ 2 u 6 To be simple, we will use the analytical solution to study the diffusion and chaos, instead of the numerical solution. The two thermodynamic quantities are We need the thermo-electric conductivities, which can be found in [53]: We will study the anisotropic chaos. The present model is the five-dimensional anisotropic model with an axion. Note that the butterfly velocity in the four-dimensional model with anisotropic lattices has been studied recently in [45]. Consider a shock-wave metric One can find that the previous form of Lyapunov time τ L = 1/ (2πT ) and butterfly velocity v 2 B = 2πT /m is not changed, but the screening length is affected by anisotropy: Using the analytical expression of the background solution and conductivities, we display in Figure 1 the temperature dependence of the butterfly velocity and two diffusivity/chaos relation for comparison. Typically, we fix the parameter β/µ = 0.1 and cut off the regime T /µ 0.1 in which β T and hence the analytical solution is not reliable. Other values of the parameter do not lead to the qualitative difference. As a result, one can see D HD = v 2 B τ L .

Non-minimal coupling
In all aforementioned models, the translation-symmetry breaking sector is minimally coupled to the gravitational and electromagnetic sectors. There are novel models which involve the non-minimal coupling between the Maxwell term and the axions [10,37]. We will focus on one of these models, i.e. the model 1 in [37]. Compared with others, its conductivities are more trivial, leading to the situation with J HD = J DGH . Actually, this model is so distinctive that it breaks various bounds on the viscosity [6], electric conductivity [54] and charge diffusivity [5]. The action is given by where J is the coupling constant and The background solution is same as that in the EMA model. Thus, the thermodynamical quantities (21) are still applicable. From Ref. [10], one can read off the thermoelectric .
Using Eq. (22) and Eq. (65), we find that the D HD is not changed by the non-minimal coupling, that is, D HD = 1/ (2r + ). Since the butterfly velocity is also same as the one of EMA theory [10], the relation D HD = v 2 B τ L remains valid. The non-minimal coupling also does not affect the ratio (27).

Nonlinear electromagnetism
Born-Infeld theory is the simplest nonlinear generalization of Maxwell electromagnetism [55]. A recent review on Born-Infeld gravity can be found in [56]. Adding two linear axions in the Born-Infeld theory gives rise to the action where b is the Born-Infeld parameter. Without the axions, the black-brane solution that is asymptotically AdS can be found in [57,58]. In Appendix B, we extend the solution to involve the axions. Using the metric ansatz (78) and the blackness factor (82), the Hawking temperature and the entropy density can be obtained where Z ef f denotes the effective electromagnetic coupling on the horizon, see Eq. (80).
From Eq. (67), one can get the derivatives In Appendix B, we have calculated the dc thermo-electric conductivities by the Donos-Gauntlett method [35]: Using the derivatives of the entropy and the thermo-electric conductivities, we obtain D HD = 1/ (2r + ) again. Performing a similar calculation as in the EMA model, we find that the Lyapunov time τ L = 1/ (2πT ) and the butterfly velocity v 2 B = πT /r + are not changed by the nonlinear electromagnetism. Then the relation D HD = v 2 B τ L is not changed either. On the contrary, the ratio between D T and v 2 B τ L is changed:

Massive gravity
The massive gravity with a reference metric is a well-known gravitational model which breaks the diffeomorphism symmetry explicitly [59][60][61]. Recently, it is applied to the gauge/gravity duality, where the reference metric imitates the mean-field disorder [62].
The action of the model is where the matrix K is defined by a matrix square root K µ ν = g µλ b λν , with b µν = diag(0, 0, 1, 1). It has been found that this model admits an analytical black-brane solution, from which the thermodynamic quantities can be obtained [62,63]. The analytical dc thermoelectric conductivities have been calculated in [64]. One can find that they are all same to the ones of EMA models, see Eq. (21) and Eq. (23). The chaos in the massive gravity has not been studied before. But we have checked that the Lyapunov time and the butterfly velocity are the same to Eq. (25). As a result, the massive gravity and the EMA model share the same diffusivity/chaos relation.

Discussion
We have studied the collective diffusion of charge and energy in the strongly coupled systems with finite density. Our focus is a particular combination of charge and heat currents that decouples with the heat current and can be reduced to the DGH current in most of homogeneous holographic lattices. We derived the condition under which the HD current can be transported by diffusion. Using the diffusion condition, Kelvin formula, and holography, we have provided a mechanism for the relation between the diffusion and chaos This might be interesting because the evidence of the relation D C ∼ D T in the normal state of the cuprates was presented by a recent experiment and it was conjectured that both charge and heat diffusivities saturate the Planckian bound [25]. Note that the relation D C ∼ D T ∼ v 2 B τ L was previously obtained in the EMA model when the momentum relaxation is very strong and in the generic holographic scaling geometries with a particle-hole symmetry [12]. Our mechanism does not require the incoherent limit and zero-density limit. Instead, the balance between the thermoelectric effect α and ζ is important, which is reflected in the diffusion conditionκ/c µ = α/ζ and the Kelvin formula Moreover, we have found a new diffusivity/chaos relation D HD ∼ v 2 B τ L that is universal in some holographic models. Most of these models are based on the EMA theory, involving various corrections: the non-relativistic scaling, the high curvature, the anisotropy, the nonminimal coupling, and the nonlinear electromagnetism. The theory of massive gravity is also taken into account. Nevertheless, it is not very difficult to find a counter example, see Appendix C. We attempt to understand the limited universality and implication as follows.
(i) There are two reasons at least that the heat diffusivity can be connected to the quantum chaos in a wide class of theories 4 . On the one hand, when the chemical potential is nonvanishing, the open-circuit heat conductivity κ is finite in the translationally invariant limit [23]. Thus, it may be intrinsic in the sense that it is not sensitive to irrelevant deformations that dissipate the momentum [65]. On the other hand, from the perspective of holography, D T = κ/c ρ can be determined solely by the near-horizon physics [12,20,23]. Thus, the universal features may be emergent due to the similarity of all horizons.
Interestingly, D HD = D C D T ζ/α is also finite in the translationally invariant limit. This can be verified readily in terms of the hydrodynamics for clean systems with charge doping: using Eq. (2) in Ref. [27], one can read σ/α = ρ/s in the dc limit and thereby D HD = D T ζρ/ (sχ) is finite. Moreover, as we pointed out in Sec. 3, the quantities χ and ζ that depend on the full bulk geometry can be cancelled in D HD . As a result, the HD diffusivity is also determined solely by the near-horizon physics.
(ii) Based on the EMAD theory, it was found that the relation D T ∼ v 2 B τ L and the Kelvin formula are both respected in the holographic models that flow to the AdS 2 ×R d fixed points in the IR [20]. Immediately, this leads to D HD = D T ∼ v 2 B τ L at low temperatures. One can find that all of the isotropic holographic models that we have studied have the AdS 2 horizon 5 , including the counter example of the relation D HD ∼ v 2 B τ L given in Appendix C. This implies that the AdS 2 horizon might be conducive to the existence of the relation D HD ∼ v 2 B τ L , but it is not a sufficient condition. (iii) In the various models studied in the main text, the relation D HD ∼ v 2 B τ L exactly holds independently of the temperature, chemical potential, and strength of momentum relaxation. In particular, the diffusion condition is not satisfied in these holographic models, except the Lifshitz gravity with z = 2. This indicates that even when the HD mode is not purely diffusive, it still has the relation to the chaos.
(iv) Although the relation D HD ∼ v 2 B τ L is not accidental, it is not universal for all holographic models. The limited universality that we have exhibited could be attributed to both the existence of IR/UV fixed points and the simplicity of those homogeneous holographic lattices. In the future, it is worth exploring whether there is an explicit physical criterion that determines when D HD ∼ v 2 B τ L holds. (v) For most holographic models in references, the bound D T v 2 B τ L is saturated at low temperatures or strong momentum relaxation. We have checked that this is true for all the isotropic models that we have studied. However, the violation of the bound D T v 2 B τ L has been found in the inhomogeneous Sachdev-Ye-Kitaev chains [67] and quasi-topological Ricci polynomial gravities [68]. Recently, by reconciling the conflict between the diffusive behavior and operator growth lightcone, Hartman, Hartnoll and Mahajan proposed a new bound D v 2 τ eq that is constituted by the diffusivity D, equilibration timescale τ eq and lightcone velocity v [69]. For the holographic theories, τ eq can be determined by the leading non-hydrodynamic quasinormal mode of black holes. The bound is obeyed in various weakly and strongly interacting theories and can be relevant to the various transport. As an example of this bound, it has been found that D T ∼ v 2 B τ eq in the EMA model with strong or weak momentum relaxation. Note that the bound would be violated if τ eq is replaced by τ L , since τ eq τ L when the momentum relaxation is weak. Let's compare D T ∼ v 2 B τ eq and D HD ∼ v 2 B τ L . It can be understood that they indicate two methods which might be useful to remedy the non-universal relation D T ∼ v 2 B τ L . One is to change the characteristic timescale and the other the diffusivity.

A Decoupled thermo-electric currents
Davison and Goutéraux have constructed two decoupled thermo-electric currents in the EMA model [28], which are correlated to the decoupling of perturbation equations in the bulk. One of the currents is expected to be transported by diffusion at all distance scales and the other carries a gapped sound mode at short distances but it is diffusive at long distances. Similar phenomenon has been found in the EMAD model [39] but has not been studied in more general holographic states. Here we will present a general expression of the decoupled currents which are transported by diffusion at least at long distances.
Define two currents by the general linear combination of the electric current and momentum These currents are decoupled J + J − = 0 provided that the coefficients γ ± have the relation Simultaneously, the matrix of the conductivities can be diagonalized, Next, we can require J ± transported by diffusion at long distances, that means Rigorously, this leads to and the diffusion constants Note that Eq. (77) denotes the eigenvalues of the diffusion matrix which have been studied in [5]. As a consistent check of the above derivation, one can find that the currents J ± with special γ ± are nothing but the decoupled currents found in the EMA model [28] and the EMAD model [39].

B Born-Infeld gravity with axions
Here we will show a black-brane solution in the Born-Infeld gravity with the axions and then calculate the dc thermo-electric conductivities. Consider the background fields taking the form Variation of the action (66) with respect to A t (r) generates the Maxwell equation Its integration leads to Inserting Eq. (78) and Eq. (80) into the Einstein equations, one can find the only nontrivial The solution is where 2 F 1 is a hypergeometric function and m is the mass parameter.
We will turn to calculate the dc thermo-electric conductivities using the Donos-Gauntlett method [35]. The consistent ansatz of the perturbation along x = x 1 reads Note that the former two modes have linear terms depending on time. This ansatz corresponds to applying an external electric field E and temperature gradient ∇T = ξT to the boundary theory. The other metric components and fields have the perturbation depending only on the radial coordinate.
In order to determine the dc conductivities, one needs to impose the regularity of fluctuation modes on the horizon: The rx component of the linearized Einstein equations gives It is important to construct two conserved currents which are independent of the radial coordinate. The first one can be obtained from the Maxwell equations: The second one is built up through introducing a two-form associated with the Killing vector k = ∂ t : where We only concern its rx component Consider the renormalized on-shell action [70] S os = (S bulk + S GH + S ct ) on−shell , where the Gibbs-Hawking term is supplied to implement a well-defined variational principle and the counterterms are invoked to cancel the UV divergence. Hereγ ab is the induced metric on the boundary, γ is its determinant and K the trace of the extrinsic curvature. The covariant currentsT µν andJ µ can then be calculated bỹ where n b is the normal vector of the boundary. One can prove that J x matches the electric current √ −γJ x on the AdS boundary exactly. Q x also matches the thermal current √ −γ γ xxT tx − µJ x , up to a term depending on the time linearly. But this term does not contribute to the dc conductivity [35]. Thus, by evaluating J x and Q x on the horizon, one can extract the dc conductivities

C A counter example
We will study an EMAD theory with the action where the gauge field coupling is Z (φ) = e −δφ and the potential of the dilaton involves three different exponential functions with the parameters δ and Setting χ i as two linear axions, an analytical black-brane solution has been found in this theory [71]. We write down the line element where f (r) = r 2 h(r) The parameter Q is related to the charge density ρ by The Hawking temperature and entropy density can be read: The thermo-electric conductivities have been obtained in Ref. [35] σ = Z + 4πρ 2 sβ 2 , α = 4πρ In terms of Ref. [72], the chaos quantities can be expressed as , v 2 B = Q + r + r + 1− 2 δ 2 +1 δ 2 + 1 2 [r + + (Q + r + ) δ 2 ] .