Abstract

In this paper, based on the basic principles of thermodynamics, we explore the hydrodynamic regime of interacting Lifshitz field theories in the presence of broken rotational invariance. We compute the entropy current and discover new dissipative effects which are consistent with the principle of local entropy production in the fluid. In our analysis, we consider both the parity even and the parity odd sector upto first order in the derivative expansion. Finally, we argue that the present construction of the paper could be systematically identified as that of the hydrodynamic description associated with spin waves (away from the domain of quantum criticality) under certain limiting conditions.

1. Overview and Motivation

Low temperature phases of various strange metal systems have been a mysterious issue since their discovery. One of the key reasons for this lies behind the fact that unlike the ordinary metallic systems, strange metals do not possess any quasiparticle description close to its Fermi surface. These systems are therefore known as Non-Fermi liquids(the enthusiastic reader is also referred to the review[1] and references therein) [25].

Matter in this quantum critical domain (so called quantum liquids) exhibits certain fascinating features as the temperature is lowered substantially as follows.

The ratio between the specific heat and the temperature turns out to be very large indicating the fact that the effective mass for electrons could be many orders larger in magnitude than the free electron mass near its Fermi surface.

The resistivity grows linearly with temperature, namely, , which is also strange from the point of view of a Fermi liquid theory.

All the above features therefore strongly suggest the breakdown of the usual Fermi liquid theory in the domain of quantum criticality and as a result the low temperature physics of strongly correlated many body systems has remained as a puzzle for the past couple of decades.

However, for various reasons, it is widely believed that the physics within this quantum critical domain is mostly governed by the obvious scaling symmetries of the underlying Quantum Critical Point (QCP). It is by now quite evident that the isometry group associated with QCPs is that of Lifshitz type; namely, the scaling symmetries are anisotropic at the fixed point [6]:where is the dynamic critical exponent [7, 8].

The purpose of the present paper is however not to explore physics within this domain of quantum criticality; it is rather slightly away from it by explicitly breaking symmetries associated with the fixed point. In the following, we illustrate this in a bit detail.

It turns out that QCPs with Lifshitz isometry group explicitly break Lorentz boost invariance and on the other hand preserve rotational symmetry as well as the translational invariance. The long wavelength dynamics associated with these fixed points has been explored extensively from both the perspective of a fluid description and the AdS/CFT duality [918]. Keeping this literature in mind, in the present paper, we are however interested in the question whether there exists any notion of hydrodynamics beyond the domain of quantum criticality and a possible systematic formulation of it based on the fundamental principles of Thermodynamics. In the following, we systematically quote the purpose of the present analysis:

We would like to explore the validity of the hydrodynamic description in a situation when the full Lifshitz algebra is broken down to some reduced symmetry group of it. It would be indeed quite interesting to explore what happens if we relax the condition of rotational symmetry and in particular how does it affect the parity odd transports of the theory. Naturally, it relaxes many of the previous constraints of the theory and allows additional tensor structures in the constitutive relations which were previously incompatible with the underlying symmetry of the system.

The ultimate goal of our analysis is to write down an entropy current [1924] in the presence of broken symmetry generators and explore how does it constraint the additional transports of the system which have emerged due to the symmetry breaking effect. In other words, we would like to explore whether there exist new dissipative effects (which are consistent with the principle of local entropy production) away from the domain of quantum criticality. In this regard, we consider both the parity even and the parity odd sector of the hydrodynamic constitutive relations and in particular explore how the parity odd transports are affected due to the presence of anisotropy in the system. In this sense, the present analysis is valid for all of those fluid systems whose underlying isometry group comprises that of the broken generators of the Lifshitz algebra.

Two questions are quite obvious at this stage: How do we introduce anisotropy in our system and What are the low lying physical d.o.f. of the system that one would like to interpret in spite of the fact that hydrodynamics is an effective description and does not rely on the microscopic details of the system.

The anisotropy that we introduce in our analysis is in fact quite similar in spirit to that with the earlier analysis [25, 26]. However, as we shall see shortly, the interpretation as well as the motivation for introducing such anisotropy is completely different in spirit from that of the relativistic scenario.

In our analysis, we introduce anisotropy by means of four vectors, pointing along a specific direction. We associate a term with those four vectors which corresponds to the fact that the pressure along two different directions is different even in the rest frame of the fluid. With the introduction of this anisotropy, our next task turns out to be to write down all possible tensor structures which are possible to construct out of the basic variables of the fluid system under consideration, namely, the fluid velocity (), temperature (), chemical potential (), and the anisotropy (). More specifically, one should be able to write down a number of tensor structures within the constitutive relations with arbitrary transports which could be fixed later by demanding the positivity of the so called entropy current. Moreover, it is to be noted that as the present theory is considered to be a slow deformation () of the original Lifshitz hydrodynamics, therefore one must recover the known results of the Lifshitz theory as the fixed point is approached ().

As far as the d.o.f. are concerned, the answer is rather tricky. However, we claim that the present analysis could be used to model the hydrodynamic description associated with spin waves [2730] in certain appropriate limit. The reason for this claim is the following. First of all, one should note that spin waves could be thought of as being low lying excitation of certain material like dimer antiferromagnet above its ground state. In other words, spin waves could be thought of as being the gapless(this renders that the underlying theory is scale invariant) excitation (Goldstone modes) associated with the broken symmetry (at the level of hydrodynamics such symmetry breaking effects should be manifested in terms of the emergence of the pressure anisotropy in the constitutive relations of the conserved quantities). These excitations therefore clearly break the underlying rotational invariance of the system and they vanish as one approaches the domain of quantum crticality [2730]. In the following, we illustrate these points in detail. In particular, we try to explain what do we exactly mean by spin wave hydrodynamics.

It turns out that the classical hydrodynamic description of spin waves is valid as long as spin wave fluctuations are smaller than the correlation length [27, 28]:where is the value of the chemical potetial exactly at the critical point and the energy of the spin wave excitation could be thought of . In other words, we are still in the region where the effective correlation is still much greater than the low temperature thermal fluctuations of the system. Therefore the hydrodynamic description is still valid. Our ultimate goal would be to investigate whether it is at all possible to write down an entropy current under such circumstances.

The organization of the rest of the paper is the following. In Section 2, we provide the basic anisotropic construction for the most general charged dissipative fluids in the presence of a global anomaly. In Section 3, we compute the entropy current for both the parity even and the parity odd sector of the theory. Finally, we conclude in Section 4.

Note Added. The entropy production in Lifshitz fluids might have direct or indirect consequences on gravitational particle production through the remarkable identification of gravitational dynamics as an entropic force [31] whose relativistic generalization leads to Einstein’s general theory of relativity.

2. Basics

The purpose of the present calculation is to explore the consequences of broken rotational invariance on chiral Lifshitz hydrodynamics [14]. However, to start with we first introduce the notion of chiral anomaly in the context of Lifshitz hydrodynamics as (the reader should not get confused with the index as the index associated with internal symmetry (non-abelian) generators. This index simply corresponds to the fact that to start with we have number of anomalous charges in our system. However, in the subsequent analysis we would perform our computation for the simplest case with . The corresponding result for multiple abelian charges could be generalized easily) [14]where the anomaly coefficients are given by the rank 3 tensor . Here, s are the anomalous currents and are the abelian field strength tensors of the theory. The electric and the magnetic field strengths are, respectively, given by

We start our analysis by writing down the constitutive relations (compatible with the reduced symmetries of the system) for both the stress tensor and the charge currents associated with the underlying hydrodynamic description. We choose to work in the so called Landau frame [14].

We first focus on the stress tensor part. It typically consists of the following two pieces (the reader should be aware of the fact that to start with Lifshitz field theories are QFTs for which . As a result such theories by no means could be connected to relativistic QFTs with Lorentz invariance (for which, ). In other words, Lifshitz QFTs do not possess a smooth Lorentzian limit and therefore it is not correct to interpret Lifshitz QFTs as nonrelativistic QFTs simply because they lack Lorentz boost generators. On the other hand, because of the presence of the gapless excitation in the theory, one might still hope to find some upper bound for the speed of particles in the theory which could be identified as the speed of light () as well. However, one could switch over to the corresponding nonrelativistic regime by taking appropriate () limit which would render nonrelativistic QFTs with broken Galilean boost invariance [9]):where is the nondissipative (ideal) part of the stess tensor. On the other hand, is the dissipative piece of the stress tensor that includes all the first order derivative expansions such as that in the Landau frame [9]:

Since the full Lifshitz isometry group is now broken down to some reduced subset of it, therefore, to start with and in principle we are allowed to incorporate more generic tensor structures in the constitutive relation which were previously disallowed due to the presence of the underlying rotational invariance. Therefore in some sense the analysis of the present paper generalizes all the previous studies in the literature [914]. On top of it, due to the presence of additional transports in the theory, the present analysis also turns out to be a bit more intricate compared to its relativistic counterpart. In the following we elaborate these issues into much more detail.

The starting point of our analysis is quite straightforward. We write all possible distinct tensor structures which are allowed by the (reduced) symmetry of the system (consistent with the Landau frame criterion) and we express our constitutive relations as a linear combination of each of these structures. The coefficients associated with these tensor structures are finally constrained by the principle of positivity of entropy production. Those transports which do not seem to be consistent with this physical constraint are set to be zero although they are allowed by the symmetry of the system. Keeping these facts in mind, the ideal part of the stress tensor could be formally expressed as (this form of the stress tensor is unique and is indeed consistent with the Landau frame criterion, namely, [9]. Moreover, here is the background metric with signature )where is the difference between the transverse and the longitudinal pressures (note that the above choice (7) also stems from the fact that at the ideal level the fluid stress tensor must be diagonal) [25, 26]. Here, is the usual fluid velocity and is the spacelike vector that defines the anisotropy in the system and is pointing along the longitudinal axis (the longitudinal axis defines the anisotropy in the system along the direction perpendicular to the reaction plane. is the pressure along this longitudinal axis, whereas on the other hand is the pressure along the direction perpendicular to this axis) such that they are mutually orthogonal to each other; namely (it turns out that the Landau frame criterion (6) and the velocity normalization (8) are intimately related to each other. Note that the Landau frame criterion (6) looks quite natural in the sense that it simply implies the fact that in the rest frame of the fluid the zeroth component of the stress tensor must be equal to the rest energy (density) of the fluid and this is the universal claim that should be valid for all classes of QFTs irrespective of their underlying symmetry group. Now it turns out that, given the stress tensor (5), the above criterion could be satisfied if we choose the specific normalization (8) for the velocity field () which therefore turns out to be a natural choice even for Lifshitz systems) [25, 26],

Considering the local rest frame of the fluid, It is indeed quite straightforward to showwhich is clearly consistent with the fact that the momentum density vanishes in the local rest frame of the ideal fluid.

Before we proceed further, one nontrivial yet simple check is inevitable. We consider the Lifshitz Ward identity (note that the Ward identity is always satisfied if the underlying scaling symmetry is preserved which is precisely the case for our analysis as the spin wave excitation is gapless [27, 28]) [9]: which for the present case yieldswhere is the number of spatial dimensions. This clearly suggests that the corresponding equation of state is modified due to the presence of anisotropy (spin wave excitations) in the system. Suppose for the time being that we are dealing with an uncharged fluid where the temperature () is the only scale in our theory. Therefore, it turns out that the above equation of state (12) could be obtained from the Euler relation if we demand the following nontrivial scaling relations: where the exponent clearly vanishes near the fixed point of the theory and as a result one recovers the standard scaling relations [9].

We now focus on the dissipative part of the stress tensor (5). After a straightforward computation, one could in fact express the dissipative contributions to the stress tensor (5) into the following two pieces; namely (note that here we have assumed as there do not exist anysources for anisotropy [25, 26]),where the first piece (here, is the socalled vorticity factor [32, 33])contains the usual contributions to the stress tensor that do not include any effects of anisotropy [14]. On the other hand, the second piece (the transport associated with the term should be considered in a more general sense. When we discuss the parity even sector it should be thought of as the parity even transport while for the parity odd sector we would replace it accordingly by parity odd transports)contains all the first order dissipative corrections to the stress tensor (5) that explicitly include the effects of anisotropy. At this stage it is customary to note that in order to arrive at the above relations (16), (17), we have used the following definitions and/or the identifications; namely,with as the projection operator. At this stage it is also noteworthy to mention that the entity contains all the first order dissipative corrections in the basic hydrodynamic variables, namely, the fluid velocity (), chemical potential (), and the temperature (). Moreover, here the entity in the Landau frame [14].

Finally, the charged current could be formally expressed as (in principle one could have added a term proportional to at the zeroth order level in the derivative expansion. However, as far as the present analysis isconcerned we switch off all the components of the electric current along the direction of the anisotropic vector field [25, 26]) where the function contains all possible dissipative corrections in the first order derivative expansion; namely (one could formally express ),such that in the Landau frame, .

3. Entropy Current

The purpose of this section is to construct a systematic hydrodynamic description away from the quantum critical region with Lifshitz isometry group. In order to do that, in this paper we illustrate the simplest case with .

3.1. Thermodynamic Identities

As a first step of our analysis, we first construct some basic thermodynamic identities which will be required in the subsequent analysis. Considering the following identity (in order to avoid confusion, we denote chemical potential as ),and using (3), it is indeed quite trivial to note that which therefore clearly suggests that at the zeroth order level in the derivative expansion of and one must have which naturally yieldswhere is the so called entropy density such that .

We would now like to interpret the entity above in (22). To do that, let us first note that in the so called boosted frame one could formally express the four velocities aswhich clearly satisfy the above constraint in (8). Here, are the arbitrary boosts along three different spatial directions. Using (23), it is indeed quite straight forward to showwhere the function could be formally expressed as

Using (25), one could now formally express (20) as

Following the arguments of [25, 26], next we introduce the so called generalized energy density functionwhich naturally yieldsThe above equation (28) could be thought of as the modified first law of thermodynamics in the presence of anisotropies. The first two terms are precisely the contributions arising from the entropy density as well as the charge density associated with the system whereas, on the other hand, the last term is precisely the anisotropic contribution to the energy density of the system.

Finally, using (26) and (28), it is quite useful to obtain the following set of identities:which we use in the subsequent discussions.

Before we conclude this section, it is now customary to derive the so called equation for the entropy current [34] in the presence of the pressure anisotropy. Using (3), it is indeed quite trivial to note down the following conservation equation:which by virtue of the arguments as mentioned above in (21) gives rise to nontrivial contributions at the quadratic order in the derivative expansion; namely,whereis the so called entropy current of the system. Clearly, from (32), it is quite evident that in the presence of pressure anisotropies, the so called entropy current of the system receives nontrivial corrections which is naturally reflected in its divergence equation (31). The first part on the R.H.S. of (31) is the usual contribution to the entropy current in the absence of pressure anisotropy [14]. The entity , on the other hand, summaries all the anisotropic contributions to the divergence equation; namely,

3.2. A Note on the Parity Even Sector

Before we proceed further, at this stage it is noteworthy to mention that the parity breaking terms would definitely modify the entropy current as expected from the earlier analysis [34]. However, before moving to the parity odd sector, it is customary to check the parity even sector in the presence of pressure anisotropy (in the subsequent discussion of this subsection we consider only the parity even part of (31). In principle, one could do this since the two sectors do not get mixed up even in the presence of anisotropy [14]. A careful investigation would reveal the fact that this is the unique feature of first order dissipative hydrodynamics which essentially truncates constitutive relations upto leading order in ). As the reader might have noticed so far that due to the presence of the pressure anisotropy in the system even in the parity even sector of the stress tensor (17) (as well as the charge current (22)) we encounter a number of additional transports which need to be fixed first. We expect to fix these coefficients purely from the condition of the positivity of the entropy current; namely, .

We would first consider the parity even sector in (31); namely,where each of the individual entities could be formally expressed as

At this stage, it is also noteworthy to mention that some of the transports appearing in the isotropic part of the parity even sector also appear in its anisotropic counterpart. As a consequence of this, all the previously determined constraints [14] on the transport coefficients of the (isotropic) parity even sector (in the constitutive relations) might not hold good and one therefore needs to redo the calculations. In the following we enumerate each of the entities in (34) separately. A straightforward computation yields the following:

(I)

(II) where

(III)

(IV) where

(V)

(VI)

(VII)

(i) Observations. In order for the entropy current (34) to be positive definite, from the above computations (36), (38), (39), (41), (42), (43), and (44) we observe the following (in the following discussion we assume, ).

Three of the following transports need to be positive definite, namely,

The rest of the transports satisfy nontrivial algebraic identities, part of which could be simplified further in order to obtain

The rest of the constraints could be formally expressed asIn other words, once we move away from the quantum criticality (by breaking symmetries associated with it) a number of additional transports appear which are found to obey nontrivial constraints among themselves. However, all of these new transports cannot be determined uniquely as we have lesser number of constraints available. It turns out that once we determine a few of them, the rest is trivially fixed due to constraints (46)-(47) present in the system.

3.3. Some Useful Relations

We now focus solely on the parity odd sector of (31). Before we proceed further, it is now customary to derive certain useful identities that would be required in the subsequent analysis of the entropy current.

(i) Derivation of . We start our analysis by noting down the following identity [14]:

Next, using conservation equation (3) we find

Finally, using (7), (24), and (49) we obtain where the vector could be formally expressed as

(ii) Derivation of . We start our analysis with the following identity [14]:

Following the same steps as above we note that

Finally, using (7), (24), and (52) we find

3.4. Parity Odd Sector

Following the prescription of [34], we modify our entropy current as which yields the divergence equation as (at this stage it is noteworthy to mention that here by we explicitly mean all the parity odd terms of the entropy current (31))where the function could be formally expressed as

In the following, we enumerate each of the individual terms one by one.

(I)

(II) At this stage it is noteworthy to mention that in deriving the above relation we have used the conservation equation (3) in order to obtain (in our derivation we have assumed that which corresponds to the fact that the gradient of the pressure difference is orthogonal to the anisotropy vector [25, 26]) Therefore, to be precise, we are looking for an on shell cancellation in the parity odd sector of the entropy current where the last term corresponds to the anisotropic contribution to the acceleration of the fluid particles [14].

(III)where .

(IV)

(V)

(VI)

Finally, collecting all the individual pieces in (56) we arrive at the following set of constraint equations, namely,where the entity could be formally expressed as

The above set of (65) precisely corresponds to the anisotropic generalization of the earlier results [14] obtained in the context of anomalous hydrodynamics at quantum critical (Lifshitz) point. Our ultimate goal would be to solve (65) in order to obtain parity odd transports away from quantum critical point. In order to proceed further, we however treat the above set of equations perturbatively in the anisotropy and/or the spin wave excitation . We consider the following perturbative expansion of anomalous transports, namely, where and . The zeroth order solutions are indeed trivial and have been found earlier in [14]. However, the first nontrivial correction appears at the leading order in . The corresponding set of equations at leading order turns out to be (for the sake of simplicity we do not quote the zeroth order solutions here. Interested reader may have a look at [14]. These solutions are purely fixed by the anomaly. However, we will use them explicitly in order to explore solutions at leading order in the anisotropy)

The above set of (68) could be simplified further in order to obtain the following reduced set of equations, namely, where the entity could be formally expressed as

Considering the general identity (here we define ), From (69), we arrive at the following set of equations:

Equation (72) is precisely an artifact of the symmetry breaking effect which nontrivially modifies the previously obtained constraint equations [14]. In order to proceed further, from (29) we note

Substituting (73) into (72) we find

From the first two equations in (74), we findwhere is some integration constant that does not depend on any of the thermodynamic variables. Using (75), the last two equations of (74) could be further modified asClearly, the above equation (76) is inadequate to solve all the five different transports in a unique fashion (similar conclusions have been drawn earlier by other authors in [35]). Therefore what we find from the above computation is that the entropy current cannot uniquely fix all the parity odd transports once we break the underlying symmetry of the system. However, it indeed imposes nontrivial constraints among different transports as we have also noticed earlier in the parity even sector.

4. Summary and Final Remarks

We now summarise the key findings of our analysis. The most significant outcome of our analysis is the establishment of the hydrodynamic description away from the quantum criticality. We claim that such hydrodynamic description could be identified with the corresponding hydrodynamic description of spin waves in the appropriate limit. The reason for this claim rests over the facts that the spin wave excitations are naturally diminished as the fixed point is approached and they spontaneously break the underlying rotational symmetry of the system. Both of them have been incorporated in the present analysis. Our analysis reveals that it is indeed possible to find out new dissipative effects when we explicitly break the underlying symmetry associated with QCPs. We compute the entropy current in order to fix these arbitrary coefficients in the constitutive relation. It turns out that none of these transports could be fixed uniquely. Nevertheless, they are nontrivially constrained by their respective constraint equations. We perform our analysis for the most general charged dissipative fluids where we consider both parity preserving and parity violating effects at leading order in the derivative expansion.

Finally, as mentioned earlier, for realistic models of strange metal systems one should actually consider the so called nonrelativistic limit of the Lifshitz hydrodynamics that has been discussed so far. This will eventually lead us towards a formulation of nonrelativistic QFTs with broken Galilean boost starting from a theory with broken Lorentz boost invariance. This could be done in two different ways. One could either take a limit of the theory [9] or simply go to some light cone frame [36, 37]. However, one should take care of the fact that in the nonrelativistic limit the so called scaling relations will also get modified leading to a modified Ward identity [9]. We leave all these issues for the purpose of future investigation.

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that there are no conflicts of interest.

Acknowledgments

The author would like to acknowledge the financial support from UGC (Project No. UGC/PHY/2014236).