Linear-T resistivity from low to high temperature: axion-dilaton theories

The linear-T resistivity is one of the hallmarks of various strange metals regardless of their microscopic details. Towards understanding this universal property, the holographic method or gauge/gravity duality has made much progress. Most holographic models have focused on the low temperature limit, where the linear-T resistivity has been explained by the infrared geometry. We extend this analysis to high temperature and identify the conditions for a robust linear-T resistivity up to high temperature. This extension is important because, in experiment, the linear-T resistivity is observed in a large range of temperatures, up to room temperature. In the axion-dilaton theories we find that, to have a robust linear-T resistivity, the strong momentum relaxation is a necessary condition, which agrees with the previous results for the Guber-Rocha model. However, it is not sufficient in the sense that, among large range of parameters giving a linear-T resistivity in low temperature limit, only very limited parameters can support the linear-T resistivity up to high temperature even in strong momentum relaxation. We also show that the incoherent term in the general holographic conductivity formula or the coupling between the dilaton and Maxwell term is responsible for a robust linear-T resistivity up to high temperature.

A UV completion 25 1 Introduction It has been shown that strongly correlated electron systems may have interesting universalities regardless of microscopic details.For example, the resistivity (ρ) is linear in temperature (T ) [1] ρ ∼ T , in various strange metals such as cuprates, heavy fermions and pnictides with a remarkable degree of universality. 1 This is a distinctive property compared with the ordinary metal case, in which ρ ∼ T 2 , a different universal property explained by the Fermi liquid theory.However, this well known "linear-T resistivity" puzzle has not been completely resolved because of the theoretical difficulty in dealing with strong correlation.As a novel and effective tool to address strong correlation in general, the holographic methods or gauge/gravity duality [1,14,15] has been widely used.The basic idea is to understand strongly correlated systems by mapping them to the dual classical gravity systems.
There have been many researches to understand the linear-T resistivity by holographic methods.Most researches, for example [3][4][5][6][16][17][18][19][20][21][22][23][24][25], have focused on the relation between the linear-T resistivity and the infrared (IR) geometry, which can be supported by matter fields and couplings.The IR geometry is characterized by three critical exponent: the dynamical critical exponent (z), hyperscaling violating exponent (θ) and charge anomalous parameter (ζ).In particular, Goutéraux [18], for the first time, systematically studied these geometries with momentum relaxation and characterized their scalings in terms of z, θ and ζ as well as derived the temperature scaling of the resistivity in low temperature limit.Here, x is some function of critical exponents.It is an interesting result since it gives an understanding of the linear-T resistivity based on the scaling properties of condensed matter systems, which are nicely geometrized in dual gravity models.However, at the same time, it has an important limitation.The result (1.2) is valid only in low temperature limit in the sense that T has to be very low compared with any other scales in given models.For example, with a chemical potential (µ) and momentum relaxation (of which strength is denoted by k), T /µ 1 and T /k 1 etc.Indeed, such a condition enables us to compute the power x(z, θ, ζ) analytically.However, noting that the linear-T resistivity has been observed in a 'large' range of temperatures, up to room temperature ∼ 300K in experiments, the condition of 'low' temperature limit for (1.2) will be too restrictive.
To deal with this problem, we first have to quantify 'low' or 'high' temperature compared with 'what'.We may choose the chemical potential µ as our reference, T /µ.However, it is not clear the relation between 'µ' in holography and the chemical potential, µ, in real world.Most holographic models including ours belong to bottom-up models.Without a top-down construction with a precise field theory dual, the meaning of the 'chemical potential' ('µ') is ambiguous.It is just a 'chemical potential' ('µ') for 'some' conserved U (1) charge.Furthermore, even in the case they are exactly the same quantity, it is still possible that there is some difference in numerical values, for example, µ ∼ 10'µ' or µ ∼ 0.01'µ'.Thus, for a reference scale, it will be better to choose an intrinsic scale in the model.For example, because experimental results show that the resistivity is linear in T up to T > T c from zero T , where T c is the superconducting transition (critical) temperature, T c can play a role as a reference scale.In other words, the 'low' or 'high' temperature can be defined by the temperature below T c or above T c . 2 Of course, it is possible to have a different reference scale other than T c depending on the context.Our purpose here is to motivate why a simple notion of small temperature based on T /µ < 1 might be misleading in some cases.For example, in the context of superconductor, T c is a better reference scale than µ because we have to check if the linear-T resistivity persists up to high temperature T > T c .We refer to [26] for more details.
Ref. [26] addressed this issue for the first time, to our knowledge.In this work, the Gubser-Rocha model [27] with momentum relaxation [17,28,29] have been considered with a complex scalar field to trigger superconducting instability.It has been shown that these models can exhibit the linear-T resistivity up to 'high' temperature, i.e. above T c only if momentum relaxation is strong.The importance of strong momentum relaxation was first emphasized in [30]: it was argued that, if the momentum relaxation, which is extrinsic and non-universal, is strong (quick), transport can be governed by diffusion of energy and charge, which is intrinsic and universal.Thus, the universality of the linear-T resistivity emerges with strong momentum relaxation 3 .
The work in Ref. [26] is important because it is the work studying the linear-T resistivity at 'high' temperature, while most of holographic works [3][4][5][6][16][17][18][19][20][21] have been focused in the low temperature limit.However, Ref. [26] deals with only one class of models based on the Gubser-Rocha model, so it is not clear if strong momentum relaxation is really necessary and/or sufficient to have the linear-T resistivity in general.The goal of this paper is to investigate this issue in a more general setup.
We start with a most general scaling geometry studied in [3,5,6,[17][18][19][20][21], so called axion-dilaton theories or the Einstein-Maxwell-Dilaton with Axion model (EMD-Axion model).The axion field is introduced to realize the momentum relaxation effect.The dilaton field is introduced with some potentials and couplings charaterized by three parameters (α, β, γ), in order to support the IR geometry parametrized by three scaling exponents (z, θ, ζ) as explained above (1.2), 4 which is rich enough to explore various possibilities.The IR geometry with an emblackening factor is valid only at low temperature limit.For the arbitrary temperature solutions, we need to introduce potentials and couplings supporting asymptotically ultraviolate (UV) AdS geometry [36][37][38].In general, there are many possibilities for potentials and couplings, and in this paper we consider a minimal (oneparameter) UV completion for potentials without changing couplings.This UV completion was introduced in [39] for the purpose of studying the shear viscosity to entropy ratio and includes the Gubser-Rocha model in [26] as a special case.
In short, our strategy is i) start with the various IR geometry giving the linear-T resistivity such that x(z, θ, ζ) = 1 in (1.2), which is valid only in low T limit; ii) after UV completing the geometry, obtain the arbitrary finite T dependence of resistivity; iii) change the momentum relaxation parameter to see how it affects the robustness of linear-T resistivity at high temperature.As a result, we have found that, in general, the strong momentum relaxation is still necessary to have a robust linear-T resistivity up to higher temperature, but not sufficient: the parameter range for the robust linear-T resistivity is quite limited compared with the high possibility in the low temperature limit [3,4,[16][17][18][19][20]. We have identified this parameter range which is different from the Gubser-Rocha model.In addition, we have also clarified the term which is responsible for the linear-T resistivity in axion-dilaton theories: it is the incoherent term 5 the first term in the conductivity formula (2.33) or the coupling between the Maxwell and dilaton fields (for spatial dimension 2).
We organize the paper as follows : In section 2, we review axion-dilaton theories and its low T limit properties focusing on the linear-T resistivity.In section 3, we classify all possible parameter range in the action to obtain the linear-T resistivity in IR and explain our UV-completion.In section 4, we report our results showing a robust linear-T resistivity up to high temperature and demonstrate the importance of strong momentum relaxation.In section 5, we interpret our results in more detail.We identify the term responsible for the linear-T resistivity.We also describe the effects of the parameters introduced in axion-dilaton theories on the temperature dependence of the resistivity.In section 6 we conclude.
2 Axion-dilaton theories: a quick review In this section, we make a very brief review on the axion-dilaton theory investigated in [18]. 6 The purpose of this section is to set the stage and collect the results that will be useful for our discussion later.We refer to the original work [18] or a review in [21] for more detailed explanation and derivations.

Action and equations of motion
An action of generic EMD-Axion models (or axion-dilaton theories) can be expressed as follows, where φ and χ i are scalar fields which are called dilaton and axion respectively.The terms denoted by J, Z, and V are the coupling functions and potential function.
The action yields the following Einstein equations: 5 The first term in (2.33) can be properly called 'incoherent' which means 'no momentum dragging' only for strong momentum relaxation, which is the very regime we are interested in.For weak momentum relaxation, there is an incoherent contribution from the second term too [40].The first term is sometimes called a pair-creation term, which is proper if there is no net charge.The first term also has an interpretation as the conductivity in the absence of heat flows [23]. 6See also [16] for the original work without axion. where and T = g µν T µν .The Maxwell equation, scalar equation, and axion equation are By considering the following homogeneous (meaning all functions are only functions of r) ansatz we obtain the Einstein equations ) which come from the equations corresponding to R tt , R rr , and R xx in (2.2) respectively.The prime denotes the derivative with respect to r.The Maxwell equation and scalar equation are reduced to and the axion equations are satisfied trivially.

IR analysis of the axion-dilaton theories
Our task is to find the functions B(r), C(r), D(r), A t (r), and φ(r) in the ansatz (2.4) satisfying the equations of motion (2.5)-(2.9)for a given the couplings and potential functions in J, Z and V in (2.1).
In order to have analytic scaling solutions in IR (Infrared: far from the AdS boundary), we give the following constrains to J, Z and V in IR: where the constant parameters (α, β, γ, V 0 ) are introduced.We consider the coefficient V 0 only for V without loss of generality because the overall factors for J and Z can be absorbed into the field χ i and A t .Plugging the "IR scaling coupling" (2.10) into the equations of motion (2.5)-(2.9)with the following scaling solution ansatz, we obtain the solution (z, θ, ζ, L, Q, κ) in terms of (α, β, γ, V 0 , k).We will call the exponents (z, θ, ζ) and the coefficients (L, Q, κ) "solution parameters" or "output parameters".We will call the parameters (α, β, γ, V 0 , k) "action parameters" or "input parameters" 7 .Physically, Q is proportional to the charge density and k means the strength of the momentum relaxation.Note that φ → ±∞ in IR8 as r → ∞ where J ∼ r κβ , Z ∼ r κγ and V ∼ r κα .After plugging scaling solution ansatz (2.11) into equations of motion (2.5)∼(2.9),we find that there are four possibilities to satisfy the equations of motion.We may classify these solutions according to the "relevance" of the axion and/or charge, following [18].By "marginally relevant axion" we mean the axion parameter k appear explicitly in the leading solutions and by "marginally relevant charge" we mean Q appears explicitly in the leading solutions.By "irrelevant axion (charge)" we mean k (Q) do not appear explicitly in the leading solutions but they can appear in the sub-leading solutions.Therefore, we may consider four classes as follows.
• class I: marginally relevant axion & charge (k = 0, Q = 0) Notice that the classification is based on the property of the leading solutions.To have a more complete picture, we also should consider the deformation by the sub-leading solutions: where Φ i denotes every leading order solution collectively, i is a small parameter and β i denotes the exponent of the sub-leading order and β i < 0 (β i > 0) when the IR is located at r → ∞ (r → 0).Therefore, Q = 0 in the leading solution does not mean zero density and k = 0 in the leading solution does not mean no momentum relaxation because these parameters can appear in the sub-leading solutions.In particular, if the axion is relevant, we may expect the momentum relaxation affects IR physics more strongly than the irrelevant axion cases.We present the explicit solutions for every classes in the following.
Class I: charge and axions are marginally relevant The leading order solutions read as well as a constraint between input parameters: or output parameters: Class II: charge is irrelevant; axions are marginally relevant The leading order solutions read as well as a constraint between input parameters: (2.17) Note that the Q and γ does not appear in the solutions (2.16) and the constrain (2.17) can be understood as the condition Q = 0 in (2.13) so ζ is also undetermined.The value of Q and ζ and their γ dependences will be determined in the subleading order As explained in (2.12), this subleading order solution will back-react to metric and dilaton field, which behave as ∼ r d+ζ−θ . (

2.19)
To have a stable IR geometry, we impose the constraint Class III: charge is marginally relevant; axions are irrelevant The leading order solutions read Note that the k and β does not appear in the solutions (2.22) so, compared with class I and II solutions, κβ = −2 does not hold anymore.The solutions (2.22) may be understood from (2.13) by setting k = 0 and replacing β by using (2.14).
Similarly to class II, in the subleading order, the axion (χ i = kx i ) starts playing a role and its back-reaction to metric and dilaton field behaves as ∼ r 2+κβ . ( To have a stable IR, we impose the constraint in the IR, r → ∞.In terms of α, β, γ this constraint means

25)
Class IV: charge and axions are irrelevant The leading order solutions read Note that the Q, γ, k, and β do not appear in the solutions.The solutions (2.26) can be understood from (2.22) by setting Q = 0 or z = 1.Similarly to case II the subleading gauge field yields

Resistivity in low temperature limit
Furthermore, it turns out that the emblackening factor f (r) can be turned on (dt 2 → f dt 2 and dr 2 → dr 2 /f in (2.11)) for all classes.Then, the Hawking temperature from this emblackening factor allows us to study the axion-dilaton theories at low temperature.However, it should be emphasized that the solution with this emblackening factor is valid only for very low temperature compared with any other scales, for example, chemical potential and momentum relaxation.
In this low temperature limit, the Hawking temperature T and charge density q can be expressed in terms of solution parameters (z, θ, ζ, Q): [18,21] where the subscript H denotes horizon.The transport coefficients, for instance electric conductivity or thermal conductivity, can be calculated [18,23].Especially the electric DC conductivity formula for the EMD-Axion model is given by (2.34) Note that only the behavior of (Z, J, C) at the horizon play a role in determining the power of T via the relation between r H and T (2.31).More variables (A t , B, C, D, Z) enter for the charge density q as shown in (2.32).By using the constraints between ζ, β, κ, (2.14), (2.20), (2.24), (2.28), we can find which term in (2.33) is dominant in low temperature limit.In class I both two terms are of the same order, in class II the first term is dominant, in class III the second term is dominant, and in class IV the dominant term depends on the parameters.i.e. in low temperature limit the resistivity ρ ∼ σ −1 DC behaves as: (2.35) Note that the power x in ρ ∼ T x does not depend on the momentum relaxation k.This is because this formula is valid only in low temperature limit.At finite temperature, this is not the case as we will show.
In fact, in addition to the constrains, (2.14), (2.20), (2.24), (2.28), there are more constraints 9 for the exponents z and θ, which will be necessary to constrain our physical model later.First, since we are considering the case that IR is located at r → ∞, 10 the exponent of r of each metric component in (2.11) should be negative, which implies θ < dz , θ < d . (2.36) Second, the emblackness factor (2.30) f (r) → 1 at the UV (r → 0), which implies Third, the condition that the specific heat is positive implies where we used the relation between the entropy and temperature (2.31): s ∼ T . Fourth, the κ, Q, and L in the (2.13), (2.16) and (2.22) should be real.The κ can be rewritten by : The conditions (2.36) -(2.37) and 3 Resistivity up to high temperature From (2.35), we may obtain the condition for the linear-T resistivity.However, as we emphasized in the previous section, the temperature dependence of the resistivity (2.35) is valid only in low temperature limit, which means that the temperature is very low compared with any other scales such as the chemical potential, the momentum relaxation strength, or the superconducting phase transition scale.However, from experimental results, it is important to have a robust linear-T resistivity also at intermediate and high temperature, which we will call finite temperature to distinguish it from low temperature limit.

UV-completion of potentials
At finite temperature, the couplings and potential in (2.10)only valid in IR have to be UV-completed.In order to have asymptotically AdS space in UV, we impose the following conditions [16]: near boundary (r ∼ 0).In other words, for V (φ), where ∆ is the conformal dimension of the dual operator of φ and φ → 0 in UV.For simplicity we set AdS = 1 from here.
In principle, there will be many possibilities to construct V, Z, J satisfying (2.10) in IR and (3.1) and (3.2) in UV.See for example [35,36,39,41].Here, for simplicity, we choose one minimal way studied in [39]: and three cases for V (φ) where c is a constant and we will set c = 0. Thus, our choice of the dilaton mass makes (3.4) satisfy the UV condition (3.2).The particular choice of V 0 we make is not necessary but for simplicity.Thanks to this choice, our V (φ) becomes a function of only α. 11We consider three kinds of V (φ) depending on the sign of θ in (3.5).To understand the necessity of these three different forms let us first find the relation between the signs of α and θ.Near IR r → ∞, we choose φ → +∞, which means κ > 0 from (2.11).When κ > 0, the relations between α and θ in (2.13), (2.16), (2.22) imply that the signs of θ and α are opposite, and if θ = 0 then α = 0. Thus, we find that the asymptotic potentials near IR read which precisely match to our IR condition V ∼ V 0 e αφ in (2.10).For example, the first potential will take a form of V 0 e −αφ for α < 0, which is not consistent with our IR condition V ∼ V 0 e αφ in (2.10).

Numerical methods
Because the potential is changing, the scaling solution (2.11) is not valid anymore so we start with the following ansatz : where u := r/r h .In this coordinate, the horizon and the boundary are located at u = 1 and u = 0 respectively.Since we want the geometry to be asymptotically AdS d+2 near boundary, we impose the conditions U (0) = 1 and S(0) = 0.
In our set-up, there are three dimensionful parameters: the chemical potential µ = A t (0), the Hawking temperature T , and the momentum relaxation strength k.For numerical analysis, we will take d = 2.We choose the chemical potential as our scale so the dimensionless parameters are T /µ and k/µ.
Our controlling parameters are (α, β, γ), which determine the whole action as well as the IR behavior of the model via (3.4) and (3.5).We want to find the parameter range (α, β, γ) which yield the linear-T resistivity from low temperature to high temperature.Therefore, we start with the restricted range of (α, β, γ), which gives the linear-T resistivity in low temperature regime.Fig. 1 shows such a range in (α, β, γ) space.It can be understood as follows.
• class IV: There is no range.See the explanation below (2.40).
The potential we consider is valid only for β < 0 because it corresponds to κ > 0 due to (2.13), (2.16) and (2.22) so correctly match the IR potential (3.7).If we want to consider the case β > 0 we need to reconstruct the potential (3.5) accordingly.Note also that the sign of α is opposite to θ so i) for α > 0 region, the first potential in (3.5) should  be used ii) for α = 0, the second potential should be used.iii) for α < 0 region, the third potential in (3.5) should be used.
4 Linear-T resistivity from low temperature to high temperature We have fine-gridded the surface in Fig. 1.By trying out all gridded data set (α, β, γ), we have found that, only for a small range of parameters near the parameter set yields the linear-T resistivity from low temperature to high temperature when the momentum relaxation is strong.See Fig. 2(a) for the numerical results for this parameter set with a strong momentum relaxation k/µ = 20.It is a plot for the resistivity (ρ) vs temperature T , showing ρ ∼ T x with x ∼ 1.To show the value of the exponent x more clearly, we make another plot for ∂ ln ρ ∂ ln T in Fig. 2(b), where we see x ∼ 1 for T /µ 5.As we explain in section 5.2, a small neighborhood of the point (−1/ √ 3, −2/ √ 3, √ 3) also exhibits the linear-T resistivity.For a purpose of comparison, we also show typical plots for non-linear-T resistivity in Fig. 3.The parameters used in Fig. 3 were chosen as the same as the one in Fig. 9(b) and Fig. 11.Note that Fig. 2 corresponds to the third potential (α < 0) in (3.5) while Fig. 3(a) and Fig. 3(b) correspond to the first (α > 0) and second (α = 0) potential in (3.5) respectively.It turns out that the third potential or α < 0 is more advantageous to have a linear-T resistivity than the others.We will discuss about it more in sec.5.2.

Momentum relaxation effect
In this subsection, we will show how the momentum relaxation affects the linear-T resistivity behavior in the finite temperature region.In Fig. 4, we find that if the momentum relaxation becomes smaller, the temperature range of the linear-T resistivity becomes shorter.The different colors in the figures represent different momentum relaxation strength: {red, orange, yellow, green, blue, purple} means k/µ = {1/10, 1/2, 1, 5, 10, 20}.From Fig. 4(a), one may think, for all k/µ we considered, the resistivity looks linear in T for certain range of T .However, in fact it is not.By carefully reading off the exponent x in ρ ∼ T x in Fig. 4(b), we find that in most cases x becomes quickly deviated from 1 as T increases from zero.Nevertheless, we stress all the curves in Fig. 4(b) go to 1 as T goes to zero.This confirms our numerics are consistent with the analytic formula (2.35).It turns out that the strong momentum relaxation is important to have robust linear-T resistivity up to high temperature.For example, the range for the linear-T resistivity is around up to T /µ ∼ 5 for k/µ = 20 and T /µ ∼ 2 for k/µ = 10 as shown in Fig. 4(b).

Interpretations of numerical results
-15 - In the previous section, we have shown the importance of large momentum relaxation to have a robust linear-T resistivity at finite temperature.To have a better understanding on the mechanism for this observation, in this section, we want to answer the following two questions.
1. (section 5.1) To have a robust linear-T resistivity, which term is important among two terms in (2.33)?We will show it is the first term, so called the pair creation term.

First term or second term?
In general, the conductivity (2.33) consists of two terms.The first term is called the pair creation term (σ DC,pc ) and the second term is called the dissipation term (σ DC,diss ) [18].For d = 2, (2.33) reads Thus, we may ask which term is responsible for the linear-T resistivity at finite temperature.
To answer this question we compare σ ,DC,pc σ DC and σ DC,diss σ DC in Fig. 5 for three representa-tive cases 12   Class I : where (z, θ, ζ) can be calculated by using (2.13), (2.16), (2.22).We can determine which (α, β, γ) corresponds to which class by checking the equality and inequalities (2.14), (2.21), and (2.25).For all cases, the resistivity is linear in T in low temperature limit as shown in (2.35).Class I (5.2),Class II (5.3), and Class III (5.4) are shown in Fig. 5(a), Fiq.5(b) and Fig. 5(c) respectively.The green and blue curves represent the momentum relaxations, k/µ = 5 and k/µ = 10.We find that, in general, at high temperature, the pair creation term contribute more dominantly.It can be understood by more active pair creation at high temperature.At strong momentum relaxation, naively, one may think that the second term is always suppressed because there is k 2 factor in the denominator.However, this is not obvious because all the other factors q, C H , and J H are also implicitly functions of k.Indeed, we find that the dissipation term is dominant at low temperature in Class III in Fig. 5(c).
In particular, Fig. 5(a) is for Eq. ( 5.2), which exhibits a robust linear-T resistivity.In this case, the larger momentum relaxation is, the more pair creation term is dominant.Thus, we find that the pair creation term, the horizon value of Z (Z ), is responsible for the linear-T resistivity at finite temperature.
Furthermore, from Fig. 5(c), we may understand why class III case is hard to exhibit the robust linear-T resistivity at finite temperature.As temperature increases, the dominant mechanism for conductivity is changed: at low temperature, the dissipation term dominates while at high temperature the pair creation term dominates.It will be more difficult to have a universal physics from two different mechanism. 13

α, β, γ dependence
In the previous subsection, we have investigated the effect of the momentum relaxation on the resistivity in finite temperature regime.In this subsection we want to investigate the effect of {α, β, γ} or {z, θ, ζ} on the resistivity.
Because we found that, in general, the large momentum relaxation is important to have a robust linear-T resistivity, here we fix the momentum relaxation to be large, say k/µ = 10.For a systematic study we first need to choose a potential in (3.5).Thus, we have three cases depending on the sign of θ θ > 0 , θ = 0 , θ < 0 .
(5.5) 12 Here, we changed z and ζ by one for comparison.In terms of α, β, γ our choice looks more complicated. 13If the momentum relaxation is weak, the dissipation term is always dominant so there is no crossing in Fig. 5(c).We thank Blaise Goutéraux for pointing this out.

This classification is equivalent to
respectively because ακ = − 2θ d with a positive dilaton φ (κ > 0).For a given potential we have three classes Class I,II, and III explained in sec.2.2.The classes are determined by the parameter range.This parameter range and the effect of α, β, γ on resistivity are best explained by figures, from Fig. 6 to Fig. 12.
Because there are common features in a class of figures (Fig. 6, Fig. 9, Fig. 11) and another class of figures (Fig. 7, Fig. 10, Fig. 12) we explain them here for all of them.1. (Fig. 6, Fig. 9, Fig. 11) show allowed region of (β, γ) for a given α.The red region, blue region, and the black line between them correspond to class II, class III, and class I respectively.
2. The dotted line in blue and red corresponds to the parameters giving the linear-T resistivity in low temperature limit.Let us imagine that we file up the dashed lines for every α.Then, the blue(red) dashed lines form a blue(red) surface in (α, β, γ) space.This surface is nothing but the surface in Fig. 1.In other words, the dashed lines in figures (Fig. 6, Fig. 9, Fig. 11) are cross-section of Fig. 1 at a given α.The region above (below) the dashed line corresponds to x > 1 (x < 1), where x is defined in the relation ρ ∼ T x in low temperature limit.
3. Therefore, the black dot in (Fig. 6, Fig. 9, Fig. 11) gives the linear-T resistivity as a class I case.Because all the other color dots do not belong to the dotted line they do not give the linear-T resistivity in low temperature limit.We chose these deviated points to investigate the effect of α, β, γ.
4. α increase as the purple dot → the black dot → the red dot in (Fig. 6 and Fig. 9).This was shown as an arrow in (Fig. 7(a), Fig. 10(a)).
5. β increase as the blue dot → the black dot → the orange dot in (Fig. 6, Fig. 9, Fig. 11).This was shown as an arrow in (Fig. 6. γ increase as the dark yellow dot → the black dot → the green dot in (Fig. 6, Fig. 9, Fig. 11).This was shown as an arrow in (Fig. 7. (Fig. 7, Fig. 10, Fig. 12) shows the exponent x in ρ ∼ T x .The colors of the curves are chosen to be the same as the colors of the dots in (Fig. 6, Fig. 9, Fig. 11).Recall that the region above (below) the dashed line corresponds to x > 1 (x < 1), where x is defined in the relation ρ ∼ T x in low temperature limit.Our results show that this tendency remain at finite temperature in general.
8. As a consistency check of our numerics, we have confirmed that the numerical values of x in the limit T /µ → 0 for every curve in (Fig. 7, Fig. 10, Fig. 12) agree with the analytic expressions (2.35) i.e. in the low temperature limit x → 1 for class I, x → 3 ) for the red curve (class III) we find that x ∼ 0.71 and x ∼ 1.31 respectively. 14These agree with our numerical results.One may think that x ∼ 0.71 is fine but x ∼ 1.31 looks quite different from the value in our numerics (the red curve in the limit of T /µ → 0).However, this is because the red curve belongs to class III.As we showed in Fig. 5(c) in class III the second term of (5.1) is dominant but the first term's contribution is still not negligible at low temperature.This contamination by the first term is reflected on our numerics.If we go to extremely low T /µ we will find a good agreement, which we have checked. 15he third potential in (3.5) (α < 0 and θ > 0) The reference point is the black dot in Fig. 6(b), which is (5.7) Fig. 6 shows the allowed region of (β, γ) for a given α: α = 1.4α 3 for Fig. 6(a), α = α 3 for Fig. 6(b), and α = 0.6α 3 for Fig. 6(c).
For fixed (β, γ) = (β 3 , γ 3 ), the α decreases from the purple dot (α = 1.4α 3 ) to the black dot (α = α 3 ) and to the red dot (α = 0.6α 3 ).For these three points the change of x in ρ ∼ T x is shown in Fig. 7(a).The value of x increases as α increases at low temperature T /µ while it does not change at high temperature T /µ.
For fixed (α, γ) = (α 3 , γ 3 ), the β increases from the blue dot (β = 1.4β 3 ) to the black dot (β = β 3 ) and to the orange dot (β = 0.6β 3 ).For these three points the change of x in ρ ∼ T x is shown in Fig. 7(b).The value of x decreases as β increases in general, but it does not change much in the intermediate temperature range 1 < T /µ < 3.
In brief, these can be summrized as follow.The region above (below) the dashed line corresponds to x > 1 (x < 1) in the relation ρ ∼ T x in low temperature limit.Our results show that this tendency remain at finite temperature in general.
From Fig. 7 one might wonder if we decrease α and increase γ the shift-up effect and shift-down effect cancel each other and the linear-T may be obtained again in another point different from the black dot.Indeed, it is true as shown in the purple curve in Fig. 8.By this way, we have found that there is some range of parameters yielding the linear-T resistivity.However, this does not work in the other way, i.e. first increase α and then decrease γ.See the red curve in Fig. 8.The difference between two ways is the class that the ending point arrives at.The former belongs to class II and the latter belongs to   class III.As we showed in Fig. 5(c) the dominant mechanism for conductivity changes as temperature increases in class III, so it is not easy to keep the linear-T resistivity across this change.
Because of very complicated coupled dynamics at finite temperature, it is not easy to have some analytic understanding of our observation.However, we speculate that it may have something to do with the vanishing potential near IR, V ∼ e αφ , for α < 0. For α 0 the potential will diverge or be constant.The first potential in (3.5) (α > 0 and θ < 0) The reference point is the black dot in Fig. 9 Fig. 6 shows the allowed region of (β, γ) for a given α: α = 0.6α 1 for Fig. 6(a), α = α 1 for Fig. 6(b), and α = 1.4α 1 for Fig. 6(c).

Conclusion
In this paper, we studied the linear-T resistivity up to finite temperature in more general cases: axion-dilaton theories or the EMD-Axion models.In order to study resistivity from low to high temperature, we start with the low temperature limit or IR limit.In this limit, resistivity can be analyzed by the scaling geometries supported by the asymptotic potential and couplings in IR (2.10) They are characterized by three parameters α, β, γ, in terms of which, the necessary conditions for the linear-T resistivity has been well studied in [18] and summarized in (2.35).In addition, there are many constraints for α, β, γ coming from physical conditions.For instance, the specific heat should be positive and our geometry should be stable under small perturbation [18,21].The constraints are expressed in (2.14), (2.21), (2.25), (2.29), and (2.40).Considering all, we have explicitly identified a parameter region which yields the linear-T resistivity in low temperature limit.This region is displayed as a two dimensional surface in three dimensional α, β, γ space.See Fig. 1.
To study resistivity at finite temperature, not only in the limit of low temperature, we have UV-completed V (φ) as with the same J(φ) and Z(φ) in (6.1).Contrary to the low temperature limit, no analytic solution is available so we need to resort to numerical analysis.By fine-gridding the surface of the linear-T resistivity in low temperature limit, in Fig. 1, we have systematically searched the parameters yielding the linear-T resistivity up to high temperature.We found that the point and its small neighborhood give the linear-T resistivity from low temperature to high temperature if momentum relaxation is strong.We have found (6.3) by systematic and hard work by brute force.However, unfortunately, we do not have a deep understanding on the precise physical mechanism for (6.3) yet.The values (and the size of its neighborhood) will be changed with different UV completions.Thus, our main point is not the values (6.3) but some qualitative discoveries we have obtained, which will be summarized as follows.
1. Large momentum relaxation is a necessary condition to obtain robust linear-T resistivity from low temperature up to high temperature.See Fig. 4.
2. Among two terms in the conductivity formula (2.33), the pair-creation term (the first term) is responsible for the linear-T resistivity.In terms of geometry, this is nothing but the horizon value of the coupling Z(φ) with the Maxwell term for d = 2 as shown in (2.1) and (5.1), i.e. e γφ(r h ) = r γκ h .(6.4) It is a very simple formula.However, it is not easy to have a good intuition, because r h is a functions of (T, µ, k) determined by complicated coupled dynamics of various fields.Only in low temperature limit, r h ∼ T so things are simplified.
3. In class III the dominant mechanism for conductivity is switched from the second (dissipation) term to first (pair-creation) term in the conductivity formula, as temperature increases.See Fig. 5. Thus, we may expect that it is not easy to have a universal property, the linear-T resistivity, in class III because of this mechanism change.Indeed, the parameters we found for linear-T resistivity belong to class I and II.
4. In class I, since both deformations (axion and charge) are marginal, the condition T Max(µ, k) is enough for the geometry to be captured by the IR scaling geometry.Thus, by increasing momentum relaxation k, it is more possible to have a larger range of linear-T resistivity than the other classes. 16 17  5.We also discussed how the conductivity behavior changes as α, β, γ changes for every potential in section 5.2.We speculate that the third potential in (6.2) may be favorable for the robust T dependence of resistivity because the potential V ∼ e αφ vanishes near IR for α < 0. For α 0 the potential will diverge or be constant.
It has been shown in [26] that the Gubser-Rocha model with linear axion exhibits the linear-T resistivity up to high temperature.This model is also included in our set-up as shown in (3.8).It corresponds to the parameter set (α, β, γ) = (1/ √ 3, 0, −1/ √ 3), which is the green dot in Fig. 13.It belongs to the first potential in (6.2) and class I. 18 In this case, the DC conductivity (σ DC ) is where Q is a complicated function of (T, µ, k).Thus, the first (pair creation) and second (dissipation) term contribute in the same way through Q.However, it has been shown that the resistivity is linear in T up to high temperature only for large momentum relaxation, k/µ 1, so the first term dominates.This agrees with what we have found in this paper.In order to identify the power x in ρ ∼ T x in a more precise way we made a plot of ∂ ln ρ/∂ ln T .This method is good enough to find a linear-T behavior all the way from zero T to high T .However, if there is a residual resistivity at zero T (i.e.T = constant + T x ) or if there is the linear-T resistivity after some temperature T 1 (i.e. for T > T 1 > 0, T = constant + T x ), our measure ∂ ln ρ/∂ ln T may not be able to capture it.Thus, in such

. 40 )
Note that Eqs.(2.28) and(2.40)imply that both 2 − ζ and d − θ − κβ can not be 1 in class IV.Thus, the class IV solution does not allow the linear-T resistivity in the low temperature limit.Being interested in the linear-T resistivity, we study only class I, II and III solutions in the following section.

Figure 1 .
Figure 1.The (α, β, γ) region which gives the linear-T resistivity in low temperature limit.The (red/blue) surface corresponds to the (class II/class III) solution.The black line, which is the intersection of the red and blue surface, corresponds to the class I. Class IV does not allow the linear-T resistivity.
The power x where ρ ∼ T x .The red dashed line is a guide to the eye.

Figure 5 .
Figure 5.The relative contributions of two terms to the conductivity in (5.1).The green and blue curves represent k/µ = (5, 10) respectively.

3 Figure 6 .
Figure 6.Allowed region of (β, γ) for a given α < 0. In Fig.7we show the resistivity for the parameters corresponding to the dots in (a), (b), and (c).See the items 1-5 in sec, 5.2 for the meanings of colors, lines, and dots.In Fig.8we show the resistivity for the squares in (a) and (c).

Figure 7 .Figure 8 .
Figure 7.The exponent x in ρ ∼ T x .The colors of the curves are chosen to be the same as the colors of the dots in Fig. 6.

3 Figure 9 .
Figure 9. Allowed (β, γ) region for a given α > 0. In Fig. 10 we show the resistivity for the parameters corresponding to the dots in (a), (b), and (c).See the items 1-5 in sec, 5.2 for the meanings of colors, lines, and dots.

Figure 10 .
Figure 10.The exponent x in ρ ∼ T x .The colors of the curves are chosen to be the same as the colors of the dots in Fig. 9.

Figure 11 .
Figure 11.Allowed region of (β, γ) for α = 0.In Fig.12we show the resistivity for the parameters corresponding to the dots.See the items 1-5 in sec, 5.2 for the meanings of colors, lines, and dots.

Figure 12 .
Figure 12.The exponent x in ρ ∼ T x .The colors of the curves are chosen to be the same as the colors of the dots in Fig. 11.