Flow of shear response functions in hyperscaling violating Lifshitz theories

We study the ﬂow equations of the shear response functions for hyperscaling violating Lifshitz (hvLif) theories, with Lifshitz and hyperscaling violating exponents z and θ . Adapting the membrane paradigm approach of analysing response functions as developed by Iqbal and Liu, we focus speciﬁcally on the shear gravitational modes which now are coupled to the perturbations of the background gauge ﬁeld. Restricting to the zero momenta sector, we make further simplistic assumptions regarding the hydrodynamic expansion of the perturbations. Analysing the ﬂow equations shows that the shear viscosity at leading order saturates the Kovtun– Son–Starinets (KSS) bound of 14 π . When z = d i − θ , ( d i being the number of spatial dimension in the dual ﬁeld theory) the ﬁrst-order correction to shear viscosity exhibits log-arithmic scaling, signalling the emergence of a scale in the UV regime for this class of hvLif theories. We further show that the response function associated to the gauge ﬁeld perturbations diverge near the boundary when z > d i + 2 − θ . This provides a holographic understanding of the origin of such a constraint and further vindicates results obtained in previous works that were obtained through near horizon and quasinormal mode analysis


Introduction
The framework of gauge/gravity duality [1][2][3][4] has been generalized and applied to understand strongly coupled non-relativistic field theories.In particular, a certain class of nonrelativistic field theories, dubbed as hyperscaling violating Lifshitz (hvLif) theories (which are conformal to Lifshitz theories) has been extensively explored in previous works .In fact, there are concrete examples of realizable condensed matter systems where certain correlators exhibit similar scaling behaviour as that of hvLif theories [35].Interested readers can see [25,35] for a comprehensive review of these class of non-relativistic field theories.
The gravity dual of hvLif theories can be realized as solutions to effective Einstein-Maxwell-dilaton theories [5][6][7][8][9][10][11][12][13][14][15][16][17][18][19].hvLif solutions may be embedded in string theory as null reductions of boosted black branes [38,39] (Lifshitz spacetimes which are conformal to hvLif spacetimes also admit gauge/string realizationss [40][41][42][43][44][45]). For a better understanding of this class of non-relativistic field theories, it is crucial to understand their infrared (IR) behaviour, in particular, hydrodynamics and various response functions that emerges in the low-energy limit.In previous works, the shear diffusion constant and the shear viscosity bound for hvLif theories were analysed using the membrane paradigm approach [26] as well as quasi-normal modes of the dual gravity theory [34].It was found that for a d i + 1-dimensional hvLif theory with Lifshitz exponent z and hyperscaling violating exponent θ, one must have z ≤ d i + 2 − θ for a consistent hydrodynamic expansion.When z = d i + 2 − θ, the shear diffusion constant exhibits a novel logarithmic scaling while the Kovtun-Starinets-Son (KSS) shear viscosity bound is saturated [46].For z > d i + 2 − θ, the first order solution diverges at the boundary presumably hinting towards a breakdown of the hydrodynamic expansion for this parameter regime.
In this paper, we take the approach as pioneered by Iqbal and Liu [47].The gauge/gravity duality maps the strongly coupled field theory on the boundary to the weakly coupled black hole spacetime in the bulk.However, the membrane paradigm approach to black holes endows hydrodynamic properties such as viscosity, entropy, conductivity etc. to a fictitious stretched horizon which is hovering very close to the real event horizon.Using the UV/IR point of view, Iqbal and Liu essentially attempted to relate this horizon fluid to the hydrodynamic regime of the strongly coupled field theory living on the boundary in the context of AdS gravity.It turned out that in the low-frequency, long wavelength limit (i.e.hydrodynamic limit) the evolution of retarded Green's function of the boundary with respect to energy scale is trivial.To be more precise, one can think of the radial direction of the bulk gravity theory as the energy scale of the boundary theory.Thus, the perturbed bulk Einstein's equations at linearized order can be thought of as a RG flow equation for a certain generalized response function which turns out to be independent of the radial direction at leading order.The triviality of flow of the response function implies that the corresponding transport coefficient can be expressed in terms of geometric quantities over any constant r hypersurface of the bulk theory and hence can be shown to be universal.
The aim of this work is to adapt the above approach and study the RG flow of response functions in the context of hvLif theories.The analysis is significantly more complicated due to nontrivial coupling between the shear perturbative modes with the gauge field perturbations.This is to be contrasted with previous works such as [22,48] where such flow equations were studied in the context of anisotropic gravity duals or the background resulted from higher derivative corrected action.In such cases, the holographic duals interpolate between Lifshitz or hvLif in the deep IR while it asymptotes to pure AdS near the boundary.
The starting point of our analysis is a (d + 1)-dimensional gravity dual of hvLif theory.Turning on perturbations of the form e −iωt+iqx h µν (r) and e −iωt+iqx a µ (r) respectively for the metric and gauge field, we notice the shear sector modes h xi , h ti and a i (where i runs over all boundary direction except t and x) forms a coupled set of differential equations.We associate a conjugate momenta to each of these perturbation modes and correspondingly define appropriate response functions.As one would expect, the radial flow equations for each of these response functions also follow complicated coupled non-linear differential equations.However, one must note that our principal aim is to extract the transport coefficient out of these response function in the hydrodynamic limit which one does in the language of linear response theory, adapted to the context of gauge/gravity duality.Consider a generic field theory containing an operator O which is coupled to a source ϕ.At the level of linear response they are related as where ω and q are very small frequency and momenta respectively, while G R denotes the retarded correlator for the operator O.The corresponding transport coefficient is defined as which is known as Kubo's formula.In particular, when O ≡ T xy , the corresponding transport coefficient is the shear viscosity η while for a charge current i.e.O ≡ J x , the analogous transport coefficient is the DC conductivity.Since in the above we essentially require to find the response function at zero momenta, we focus on that regime and analyse the flow equations.Interestingly, we see that indeed for q = 0, the flow equation for χ xi i.e. the response function corresponding to h xi follows a Riccati equation which leads to a constant χ xi at leading order for all values of z and θ.This behaviour is identical to that encountered in pure AdS gravity [47].However, when z = d i −θ, the first order correction to the response function has a logarithmic scaling which diverges at the boundary r → 0. This necessitates the introduction of a cut-off presumably signifying the UV scale beyond which the hydrodynamic expansion breaks down.The analysis for the response function associated with a i i.e. χ a i is more involved due to the complicated nature of the flow equation.In fact at q = 0, it turns out the variable ζ a i = ωχ a i seems to admit a hydrodynamic expansion.In order to analyse the behaviour of ζ a i , we focus on the near-horizon region and the near-boundary region separately which somewhat simplifies the analysis.At leading order itself, we see the solutions for ζ a i are different for the two different regimes.This is different qualitatively from the behaviour of χ xi which followed a trivial flow equation allowing one to write the response function at any point along the radial direction.Interestingly, we see close to the boundary, the leading piece of χ a i diverges when z > d i + 2 − θ which is identical to results obtained in earlier works [26,28,34].To further vindicate our result, one can look at the Markovianity index of the fluctuating modes in the spirit of [49].Interestingly, we observe that for z ≤ d i + 2 − θ, the fluctuations starts to behave like a non-Markovian probe.
The paper is organized as follows: In section 2, we describe our setup and define appropriate response functions corresponding to the shear gravitational modes.The general flow equations are worked out which describes the non-perturbative evolution of response function for arbitrary frequency and momenta.Section 3 focuses on the zero momenta sector and look at the transport coefficient associated with the modes h xi , h ti and a i .Finally, keeping some details of the calculations in three appendices, we end with a discussions of our main results with possible future directions along with a simple analysis of the Markovianity index results in section 4.

Flow equations of response functions
We are considering a hvLif theory living in d = d i + 1 spacetime dimensions with Lifshitz exponent z and hyperscaling violating exponent θ.This field theory has a (d + 1)-dimensional gravity dual given by (2.1) The above metric is a solution to Einstein-Maxwell-dilaton theory (details of background solution in appendix A).The temperature of the field theory dual to the hvLif theory (2.1) is the Hawking temperature of the black brane where the event horizon is located at r = 1 r 0 .As per the holographic dictionary, the radial coordinate r can be thought of as the energy scale in the bulk theory.Our central goal in this section is to essentially set up the RG flow equations governing the response functions that we want to study.In order to obtain the RG flow equations, we turn on linearized perturbations in the bulk theory, which in general is given as, where quantities ḡµν , Āµ and φ denote background fields as given in appendix A. We turn on perturbations of the form e −iωt+iqx h µν (r) and e −iωt+iqx a µ (r) and restrict ourselves to the radial gauge (h µr = a r = 0).The shear gravitation modes h xi now couples to h ti and a i where the index i runs over all boundary coordinates except t and x.For convenience, we define the following field variables In terms of these modes the equations of motion take the form ) where k = (d i + z − θ)α.The above linearized equations of motion follow from the perturbed second order action, the details of which are provided in appendix B. In terms of the variables defined in (2.4), the action (B.1) can be recast in a more 'canonical' form as which yields (2.5)-(2.8)as the equations of motion.For completion, we state the boundary action i.e. S bdy is given by (2.10) Motivated by the equations of motion appearing in (2.5)-(2.8),we observe that the coupling term between H ti and a i appearing in the above action (2.9), namely +kH ti a ′ i can be rewritten as −kH ′ ti a i along with a boundary term.Thus, the effective Lagrangian reads as bdy . (2.11) The conjugate momenta for the modes H xi , H ti and a i are defined respectively as, (2.12) The above definitions immediately yield Corresponding to each of the modes H ti , H xi and a i , we associate a response function given by, In terms of the response functions, the constraint equation (2.7) takes the form (2.17) Using the equations (2.12)-(2.17),we can eventually write down the generalized flow equations for the response functions χ xi , χ ti and χ a i which takes the form Note that the above set of coupled differential equations are exact i.e. they describe the flow of response functions for generic values of frequency and momentum.Although, they are complicated and analytically intractable, we are however interested in the hydrodynamic regime which is essentially the limit where the frequency and momenta ω and q are much smaller than the temperature scale i.e. q ≪ T 1/z ∼ r 0 and ω ≪ T ∼ r z 0 .Further, it is evident from (1.2) that the q = 0 sector is relevant for evaluating shear transport coefficient.Thus, we will focus exclusively on the q = 0 sector of the flow equations (2.18)-(2.20).

Zero momentum response functions
Before we proceed to study the q → 0 limit of the flow equations we derived in the preceding section, it is imperative to talk about solutions of the field equations in the q → 0 limit.An earlier work [34] analysed the field equations assuming a hydrodynamic expansion in the dimensionless parameters Ω = ω 2πT and Q = q (2πT ) 1/z .One can however reabsorb the constant temperature factor in each term of the hydrodynamic expansion and simply write the fields as an expansion in ω and q.
Starting with the equations of motion (2.5)-(2.8),a gauge invariant combination H i was defined as The fields H i and a i formed a system of coupled differential equations which were solved up to first order in the hydrodynamic expansion.For the redefined field H, it was observed that for z < d i + 2 − θ the terms in the hydrodynamic expansion of the field variables can be solved order-by-order.When z = d i + 2 − θ, the first order correction to H i scales logarithmically and seems to diverge close to the boundary.The logarithmic scaling is suggestive of the emergence of a new scale in the UV limit.In the regime when z > d i + 2 − θ, the first order correction to H i diverges suggesting a breakdown of the methodology for parameters in this regime.The solution to the combination H i up to first order in the hydrodynamic expansion is given by The gauge field fluctuations, a i satisfy a second order non-homogeneous differential equation.Upon imposing regularity on a i , the leading solution takes the form where C 0 is an arbitrary non-zero constant.The first order piece does not have a closed form solution but can be written as an integral and thus cannot give us further insight into its behaviour.The reader can find details and methodology of solving for a i up to first order in [34].Since in the current context, our interest is to explore the flow equations, we will not further concern ourselves with solutions to the fields H xi , H ti .We will however make certain assumptions about them which will help us in dealing with the complicated flow equations we derived in the preceding section.Motivated from the form of H i and a i as given in (3.2) and (3.3), we will assume that the perturbations H xi , H ti and a i for ω = 0 and q = 0 admit a hydrodynamic expansion of the form φ(r, ω, q) = φ (−1) (r, ω, q) + φ (0) (r, ω, q) + φ (1)  (3.4) where φ(r, ω, q) represents any one of the perturbative modes H xi , H ti or a i .The leading term φ (−1) (r, ω, q) is parametrically an O( 1 ω ) quantity while φ (n) (r, ω, q) ∼ O(ω n ) ∼ O(q n ).The first term in the above expression can be generically of the form while the O(1) term and the n-th order term in the hydrodynamic expansion will take the general schematic form respectively.In (3.5) and (3.6), all the exponents a, b and k are strictly positive.The functions b a (r), g a (r) and h (a,b,k) (r) are all regular in the interval 0 < r < 1 r 0 .The sum in the n-th order term has a 'prime' to denote that it is a constrained sum such that a + b = n ≥ 1.The above schematic forms of each term in the hydrodynamic expansion of the field variables is well behaved in the limit q → 0. A comparison of (3.3) with the schematic forms as given in (3.5) and (3.6) tells us that for a i , the O( 1 ω ) term is identically zero; the O(1) term consists of a single term with a = 1 while g 1 (r) ∼ f (r) − iΩ 2 .A comparison of the above schematic expansion with H i as given in (3.2) is difficult, since it appears as as a linear combination of H xi , H ti and an integral over a i .We can still comment on the heuristic behaviour of the response functions that follows from the above assumptions regarding the hydrodynamic expansion of the field variables.
The structure of (2.13)-(2.16)along with (3.5)-(3.6)implies that the response function can be written schematically as where χ (n) (r, ω, q) denotes a term which is O(q n ) ∼ O(ω n ) in the hydrodynamic expansion but is determined by the explicit forms of b a (r), g a (r) and h (a,b,k) (r) while F (r) is some specific function depending on which mode is under consideration.More specifically, where W[f, g] = f g ′ − f ′ g denotes the Wronskian for the pair of functions f and g.In the case, neither of these are linear combinations of various powers of the ratio q ω , we simply recover Armed with the above heuristic analysis, we further closely look at the following terms appearing in (2.18)-(2.20).
• The last term appearing in (2.18) can be written as where H(r) is a function of r whose details we are not concerned with for the purpose of this analysis.For the sake of simplicity, if we assume (3.5) is not in fact a linear combination of various powers of the ratio q ω , the leading behavior of this term for non-zero q and ω is given by where the "• • • " represents terms that are higher order in q or ω while Hn (r) represents various functions of r.Now, every expression Hn (r) involve ratios of derivatives of the family of functions b a (r), g a (r) and h (a,b,k) (r) that appear in the hydrodynamic expansion of the field variables.To be more explicit, . (3.12) At this point, we further make the assumption that for non-zero q and ω, each of these functions i.e. b a (r), g a (r) and h (a,b,k) (r) appearing in the field expansion of H ti and H xi are non-constant, non-trivial functions of r.Clearly, under the above assumption, this term vanishes when q → 0. We will subsequently infer from the equations of motion at q = 0 that χ ti = 0 for this sector however, the term that we just discussed does not have any singularity or does not go to any constant as q → 0. Physically speaking, χ ti presumably contains a leading O(q) piece which ensures that the ratio q 2 χ ti → 0 as q → 0 while presence of higher powers of the momenta q in subsequent higher order terms ensure it vanishes identically as q → 0.
• In (2. 19), we see the last two terms can be schematically written as Although this term will indeed vanish at q = 0, since χ ti = 0 in this sector, our assumptions up to this point dictates a possible O(ω) contribution as to the flow equations as q → 0. Hence, we keep this term in the limit of vanishing momenta.
• Finally, we need to study the final term in (2.20) which we schematically write as q ω 2 P (r) Note that from (3.3) and comparing with the expansion (3.4), it follows that for the field a i , the family of functions b a (r) = 0 identically.Again, for simplicity, assuming H xi does not have have a linear combination of terms at O(1/ω), we get the leading behaviour of the last term as q ω 2 P (r) Thus, in the q → 0 limit, we recover the simplified flow equations as ) The above equations can also be derived by turning on perturbations of the form e −iωt h µν (r) and choosing the radial gauge h µr = 0. Before we proceed with the detailed analysis of the flow of response functions, the equations of motion in the q = 0 sector simplifies significantly to give Thus, in the q = 0 sector, the mode H xi further decouples from H ti and a i .The constraint equation (3.21) clearly implies By the assumptions we made in (3.4)-(3.6),we see that This indeed is consistent with (3.17) and renders the equation trivial.Thus the q = 0 sector requires us to analyse two independent equations governing the flow of H xi and a i given by (3.16) and (3.18) respectively.

Response function χ xi at q = 0
As argued in the previous section, in the q = 0 sector, the χ xi flow equation decouples from the χ ti flow equation and (2.18) simplifies to If we demand regularity of χ xi at the horizon, we clearly see the RHS of the above is singular at r = 1 r 0 .This forces us to choose leading to the boundary condition In the hydrodynamic regime, we are allowed to write a perturbative expansion for χ xi (r, ω) where O(ω 2 ) represents higher order terms beyond the linear one.Plugging in the above expansion, in (3.25), the leading order piece follows Physically, the above equation tells us that the RG flow of the χ xi is trivial at leading order remaining unchanged as we go along the radial direction.Along with the boundary condition (3.27) that we just derived, we have The O(ω) equation which gives the flow of χ xi , is given by, The solution to the above equation is ) where C and C ′ are integration constants for the two cases of the Lifshitz exponent z while 2 F 1 [a, b, c; r] represents the hypergeometric function.We then come across the following two cases, Case I • z = d i − θ: Using the boundary condition (3.27), we can fix the constant of integration to be where γ is the Euler-Mascheroni constant and ψ(x) is the polygamma function which is singular over the set non-positive definite integers.Taking into account the null energy condition (A.5), we focus when d i − θ > 0 and z > 1. Since, this solution is true when z = d i − θ, the argument in the polygamma function cannot be 0.However, z−d i +θ d i +z−θ = −1 gives z = 0 which violates our the assumption of z ≥ 1.For all other parameter values of (z, θ) the null energy condition ensures that ψ( z−d i +θ d i +z−θ ) is non-singular.
Case II • z = d i − θ: In this case too, plugging in the boundary condition (3.27), we get, which then gives the full solution Clearly the divergent nature of the solutions as r → 0, hints at a possible breakdown of the analysis when z = d i − θ near the boundary.
Earlier works [28,34] used perturbative techniques to evaluate 2-point correlator of the stress-energy tensor which seemingly broke down when z > d i + 2 − θ.However, an analysis of the response function corresponding to H xi i.e. χ xi seems to carry through for all values of the Lifshitz exponent.As mentioned earlier in (1.2), shear viscosity up to leading order is thus given by which inturn saturates the KSS bound of η s = 1 4π .Also, note that the first order correction for either cases, namely z = d i − θ and z = d i − θ is positive since r 0 r < 1 thus following the bound.However, when z = d i − θ, we see the first order correction to be logarithmic and is actually divergent at the boundary when r → 0. This enforces us to put a cut-off suggesting the emergence of a new scale.
Interestingly, earlier works [40,41] constructed families of Lifshitz geometries as dimensional reduction of AdS null deformations.Specifically, starting with AdS 5 null deformation, one can perform a reduction along one of the light-cone coordinates, namely x + which results in a 4-dimensional metric of the form (2.1) with z = d i = 2 and θ = 0. Thus, dimensional reduction of null deformed AdS 5 results in a metric which falls in the family of hvLif solutions constrained by z = d i − θ.In light of this observation, it will be interesting to understand the logarithmic scaling of the first order contribution to χ xi from the perspective of the deformed higher dimensional theory.

Response function χ a i at q = 0
Recall from our earlier definition (2.16), that the response function χ a i associated to a i is defined as To reiterate, the flow equation for the response function χ a i decouples from that of χ xi and χ ti in the limit q → 0 to yield, The structure of (3.38) is significantly different from the flow equation of χ xi .Assuming a Laurent expansion in ω for the function χ a i (r, ω), we see that in general it must have a term which goes as 1 ω along with regular terms.Thus, like the earlier case of χ xi , it does not make sense to naively perform a hydrodynamic expansion of containing only positive powers of ω.However, we define the new field in terms of which (3.38) becomes Imposing regularity for ζ a i along the radial direction demands us to write the boundary condition as which eventually yields, One must note that (3.40) is exact in ω and consistent with a hydrodynamic expansion of the form 16πG N .Unlike the earlier case of χ xi , we see here that at leading order ∂ r ζ a i follows a nontrivial flow equation given by Thus, we see for this response function, the RG flow is not trivial and it actually changes along the radial direction.Solving the above equation yields complicated solutions which one cannot use easily to construct further subleading contributions that are higher order in ω.
To circumvent the issue, we follow a different strategy.We will analyse the flow equation successively in the near horizon and the near boundary region.

Near horizon region:
In order to analyse the flow near the horizon, we define a new radial coordinate ρ given by In turn, the blackening factor can be written as Thus in the near horizon region, the flow equation can be approximated as (3.47) An ansatz consistent with a hydrodynamics description may be written as nh (0) = 0.The resulting equation for Which on solving naively yields a solution of the form − i c 1 +Aρ+B log ρ where A and B are constants depending on r 0 , d i , z and θ while c 1 is an arbitrary constant which remains unfixed even after imposing the relevant boundary condition for ζ (0) nh (ρ).This is because the very boundary condition (3.42) is specified at a singular point of the equation.We can however choose a cutoff surface at ρ = ǫ (which can be thought of as a stretched membrane) hovering at a distance ǫ outside the real horizon at 1 r 0 where ζ identically in the near horizon region.This in turn leads to the simple equation at O(ω) i.e.
The above along with (3.42) implies Eventually, the equation at O(ω 2 ) is given by The solution to the above equation consistent with the boundary condition (3.42) is Thus in the near horizon region, we have, Using (3.39) and (3.45), we see that in the near horizon region, we can write, Clearly, χ a i being a constant at leading order exhibits trivial RG flow and is thus qualitatively similar to χ xi .However, one must note that this is true only in the near horizon region.
Near boundary region: In this regime, we can approximate the blackening factor f (r) ≈ 1 which simplifies (3.40) to Assuming a series expansion in ω of the form we see that the leading order satisfies an equation of the form whose solution is given by Assuming reality of the gauge field (A.3) i.e. z > 1 the null energy condition (A.5) implies Thus, near the boundary, Clearly, when bdy (r) → 0 as r → 0, however for d i +2−z−θ < 0, we see a divergent solution as r → 0 while it goes to a constant as r → 0 when z = d i +2−θ.In fact, due to the functional form of the solution (3.60), its limit as r → 0 will be independent of the constant c 1 which will remain unfixed for any Dirichlet condition imposed at the boundary.Hence, for z > d i + 2 − θ it seems such a hydrodynamic description for the gauge field response function will simply breakdown near the boundary.
Starting with a AdS d i +3 dimensional boosted black brane, performing a boost and taking an appropriate double scaling limit involving the boost parameter and horizon radius yields the so-called AdS d i +3 plane wave.Subsequently reducing along x + and identifying x − ≡ t yields (2.1) where the Lifshitz exponent z and hyperscaling violating exponent θ are related by [38] Clearly, from the above expressions it follows that z = d i +2−θ.This is precisely the point in the (z, θ) parameter space where we see the leading behaviour of ζ a i near the boundary is a constant.From the viewpoint of the AdS d i +3 boosted black brane, this is suggestive that the hydrodynamic analysis for such effective theories obtained as null reductions break down.However, a concrete understanding of this breakdown would require further detailed analysis concerning the stability of such spacetimes which we plan to carry out in subsequent works.
It is interesting to notice that this condition was recovered in earlier works [28,34].In particular, [34] studied QNM modes in the black brane background given by (2.1).As described earlier in section 3, the gauge invariant combination (3.1) has a solution given by (3.2) up to first order in the hydrodynamic expansion for z < d i + 2 −θ.For z = d i + 2 −θ, the first order term develops a logarithmic scaling while it diverges near the boundary when z > d i + 2 − θ.Our current analysis suggests it is the behaviour of the perturbations in the background gauge field i.e. a i near the boundary which is presumably the cause of this divergence.Thus, the RG analysis seems to be suggestive of the fact that it is the hydrodynamic expansion of a i which breaks down causing its response function to yield an unphysical answer when z > d i + 2 − θ.Further, we should contrast this with [26,28] which were near-horizon analysis, also led to the same restriction on the Lifshitz exponent z.In our current analysis, the divergence seem to occur in the boundary theory as r → 0. Earlier work [22] studied hvLif solutions as solutions to theories with higher derivative corrections.Null energy conditions and stability criteria led to certain regions in the (z, θ) parameter space that were identified as physically allowed.The criteria that we obtain above i.e. z < d i + 2 − θ seems to be an independent bound which cannot be obtained by NECs or stability criteria.
The first order equation is given by (3.64) Owing to the null energy condition (A.5) and reality of the gauge fields, which implies z > 1, we see that lim which leaves the constant c 2 which remains unfixed.Finally, the equation governing the second order contribution is given by whose solutions are listed in appendix C. From (3.39), it follows that the response function associated to a i near the boundary is given by The response function χ a i at leading order exhibits non-trivial dependence on the radial coordinate and thus shows a very distinct behaviour compared to the response function χ xi .

Discussion and Conclusion
In this paper, we have studied and analysed the RG flow equations governing the shear response in hvLif theories from the holographic viewpoint.The presence of U(1) gauge fields along with a dilaton complicate the analysis significantly since certain gauge field perturbations i.e. a i couples to the shear tensor modes h xi and h ti .Focusing on the q = 0 sector, our central observations are: • The shear viscosity at the leading order seems to saturate the KSS bound for all values of z and θ.Earlier works failed to make any statement about shear viscosity for z > d i + 2 − θ.This analysis gets around that issue of breakdown of hydrodynamic expansion for z > d i + 2 − θ.However, for the special value of z = d i − θ, we see a very interesting logarithmic correction at the first order.This does not violate the KSS bound but, necessitates the introduction of a UV cutoff to control potential divergences at the boundary.This particular logarithmic behaviour of the subleading correction to shear viscosity for z = d i − θ seems to be a novel feature.Further, as discussed in previous works [40,41], dimensional reduction of null deformed AdS 5 results in z = 2 Lifshitz theories (they have θ = 0) and is consistent with z = d i − θ = 2. Given this observation, it is natural to ask if such logarithmic correction for z = d i − θ can be explained from the perspective of the higher dimensional null deformed AdS 5 theory.It will be interesting to further explore the hydrodynamics of theories dual to such null deformed background.
• In the response function for a i , we observe non-trivial flow even at leading order in χ a i .We have performed the analysis in the near horizon and the near boundary region with appropriate approximations.In the near horizon region, the qualitative behaviour of χ a i seems to mimic that of χ xi .However, the near boundary analysis reveals a leading behaviour which scales as χ a i ∼ r d i +2−z−θ .The response function happens to be convergent provided z < d i + 2 − θ.Thus, it seems this bound obtained in earlier works [26,28,34] can be interpreted as a regularity condition on the response function of the gauge field perturbations a i .Earlier works constructed a linear combination involving all the perturbation modes which obfuscated the source of this constraint.Our present analysis seems to suggest that it is the gauge field perturbations exclusively which are responsible for the the constraint z < d i + 2 − θ.
An aside on Markovianity index: At this point one can ask for a more physical origin for the constraints observed in this paper.In other words, we want to understand if the breakdown of the hydrodynamic expansion for a certain parameter range, namely z > d i + 2 − θ has a more deeper origin or is simply a bug of these non-relativistic gravity duals.Towards that vein one perform a Markovianity index analysis of the perturbations in the probe limit in the spirit of [49].To be more elaborate, [49] studied probes couples to conserved currents in an AdS-Schwarzschild background.The effective coupling of the probe field is characterized by a single parameter, namely the Markovianity index M. Probe fields with M > −1 exhibits short-lived memory and behave analogous to the massive scalar probes.Probes with M ≤ −1, however, retain long-term memory.In the current context, the metric petrubations we study are coupled to conserved current i.e. the stress tensor.
More precisely, [49] starts from the effective action of a probe scalar of the form describing a massless Klein-Gordon field minimally φ M minimally coupled to gravity with metric being same as (2.1) and M ∈ R being some designer parameter modulating the coupling.To reiterate more concretely, With this designer scalar the central observation of [49] is that the scalar probe field φ M is Markovian if M > −1 or else its non-Markovian.Written in terms of the usual Fourier modes the scalar field equation takes the form in the zero momentum limit 1 1 with q = 0 the equation becomes Comparing (4.2) with (2.6) in the limit ω → 0 one can check that in this case the designer parameter turns out to be Interestingly, this implies that constraining the perturbations to be Markovian also forces the probe to obey the null energy condition (A.5).In other words The situation with (2.5) or (2.8) is much more complicated due to the coupling between the fields.One can simplify the situation by considering the near boundary region for (2.5) where f (r) ∼ 1.In this regime, a comparison between (4.2) and (2.5) in the limit ω → 0 reveals 2 M = d i − z − θ + 1.
Again imposing the Markovianity condition M > −1, we interestingly have which is exactly the limit that we have obtained through the gauge field perturbations.Therefore the regularity condition of z < d i + 2 − θ., analysed explicitly in the present analysis and observed earlier in [26,28,34] can also be attributed to the fact of the probes being Markovian.The above calculations although rudimentary seems to be hinting towards a connection between Markovianity index of the fluctuations and the breakdown of hydrodynamic expansion.An elaborate investigation of this issue is beyond the scope of this paper, which we hope to address in future works.
Our strategy of analysing the near horizon and near boundary regions separately opens up some possible new directions in the hydrodynamics of hvLif theories.Our analysis is restrictive in the sense that we analysed the q = 0 sector only.A natural extension will be to understand the q = 0 sector and check if one recovers any new transport coefficient at linear order in q.Another interesting question will be to explore if regularity conditions imposed on response functions of higher order transport coefficients leads to any further constraint on the Lifshitz exponent z.We are looking forward to analysing the flow equations both analytically and numerically to comment on higher order transport coefficients which are yet unexplored in the literature.We subsequently plan on studying the RG flow of response functions that arise in the sound channel and scalar channel.
2 Near the boundary with ω → 0 and q → 0 equation (2.5) becomes where the term containing a ′ i gets dropped due to (3.3).

B Perturbations to hvLif spacetimes
The perturbed action up to second order terms is given by +2J (r, q, ω)h ti h xi + M(r)a ′ 2 i + N (r, q, ω)a 2 i , (B.1) where the various functions appearing in the action is given by: A(r) = −2r All of them vanish in the limit r → 0 (near boundary) thus leaving the constants c 3 , κ 1 and κ 2 unfixed.

. 43 )
Also, the demand of regularity gives us the ζ a i at the horizon which depends explicitly on the frequency ω.Thus, regularity in the context of the above hydrodynamic expansion implies ζ (m) a i (1/r 0 ) = 0 for all m = 1 while ζ (1)