Lifshitz Anomalies, Ward Identities and Split Dimensional Regularization

We analyze the structure of the stress-energy tensor correlation functions in Lifshitz field theories and construct the corresponding anomalous Ward identities. We develop a framework for calculating the anomaly coefficients that employs a split dimensional regularization and the pole residues. We demonstrate the procedure by calculating the free scalar Lifshitz scale anomalies in $2+1$ spacetime dimensions. We find that the analysis of the regularization dependent trivial terms requires a curved spacetime description without a foliation structure. We discuss potential ambiguities in Lifshitz scale anomaly definitions.


Introduction
Lifshitz scale symmetry is a symmetry under which space and time scale differently, where z is the dynamical critical exponents which represents the anisotropy between space and time, and d is the number of spatial dimensions. Non-relativistic field theories that exhibit a Lifshitz scale symmetry have attracted much attention in recent years. In nature Lifshitz scaling is a property of certain low energy systems of condensed matter, which exhibit quantum criticality (see for example [1]). The study of holographic Lifshitz systems has been initiated e.g. in [2,3] Similary to the relativistic trace anomalies (for brief reviews see e.g. [4,5]), nonrelativistic quantum field theories may exhibit Lifshitz scale anomalies, where the classical Lifshitz symmetries are broken by the quantum corrections. A cohomological analysis provides a general framework to determine the possible structures of Lifshitz scale anomalies [6][7][8][9] (see also [10,11] for a z = 2 example). Naturally, it is of importance to calculate the anomaly coefficients, which carry information about the quantum field theory. A heat kernel calculation of anomaly coefficients of a Lifshitz scalar in 2 + 1 spacetime dimensions has been performed in [11], and in 3 + 1 dimensions using a different regularization scheme in [12]. The aim of this work is to develop a general scheme for such field theory calculations. In order to do that we will first analyze the general structure of correlation functions of the stress-energy tensor in Lifshitz field theories and analyze the corresponding anomalous Ward identities. We will encounter a subtle ambiguity in the definition of the anomaly coefficients and will clarify its meaning. Next, we will develop a general framework for calculating the anomaly coefficients. It consists of two elements: a split dimensional regularization [13,14] where space (momentum) and time (frequency) integrals are regulated separately, and pole residue calculations [15] that allow to extract the anomaly coefficients without a full calculation of the correlation functions. In order to demonstrate the power of the latter we will calculate using the pole residues the trace anomaly coefficients of a relativistic scalar in two and four spacetime dimensions and obtain the known results. We will then apply the complete framework to calculate the anomaly coefficients of a Lifshitz scalar in 2 + 1 spacetime dimensions. The results agree with the heat kernel calculation [11]. Following the discussion in [7], we will show that the analysis of the regularization dependent trivial terms arising from this calculation requires a curved spacetime description that violates the Frobenius condition (and therefore has no foliation structure). This paper is organized as follows. In section 2 we present the Ward identities corresponding to the symmetries of the field theories we consider, in terms of the expectation value of the conserved currents in curved spacetime, and in terms of the flat space correlation functions of the stress-energy tensor. We also discuss a possible ambiguity in the Lifshitz anomaly coefficients. In section 3 we explain the method by which we use split dimensional regularization to extract the Lifshitz anomaly coefficients from the flat space field theory correlation functions. In both sections we start by reviewing the case of a conformal field theory and then discuss the Lifshitz case. In section 4 we apply the aforementioned method of extracting scale anomaly coefficients from the pole residues of correlation functions in dimensional regularization to calculate the trace anomaly coefficients for a free relativistic scalar, both in two and four spacetime dimensions. In section 5 we use this method in its split dimensional regularization version to calculate the anomaly coefficients for a z = 2 free Lifshitz scalar in 2 + 1 dimensions. We conclude in section 6.

Ward Identities and Correlation Functions
In this section we present the form of the Ward identities corresponding to the symmetries of the field theories we consider, both in terms of the expectation values of the conserved currents over curved spacetime, and of the flat space correlation functions of these currents. These include the anomalous scaling symmetry as well as the other symmetries that are assumed to be non-anomalous. The anomalous Ward identities will later be used to extract the anomaly coefficients from the flat space correlation functions. We start by reviewing the conformal case, and then discuss the Lifshitz (non-relativistic) case.

Symmetries and Ward Identities In Curved Spacetime
Consider a conformal field theory in d dimensions. The theory can be coupled to a curved spacetime manifold equipped with a background metric g µν (or alternatively vielbein structure e a µ ). Suppose this theory is described by the classical action S(g µν , {φ}), where {φ} stands for the dynamic fields of the theory. This action reduces to the flat space action of the theory when g µν → η µν , and is invariant under the following symmetries: 1. Diffeomorphisms, given by: 1) or in vielbein notation (here we also include local Lorentz transformations): δ D ξ e a µ = ξ ν ∇ ν e a µ + ∇ µ ξ ν e a ν , δ L α e a µ = −α a b e b µ .
(2.2) 2. Weyl transformations, given by: 3) The stress-energy tensor is defined as the variation of the action with respect to the metric: T µν = 2 |g| δS δg µν , (2.4) where g ≡ det g µν . The stress-energy tensor satisfies classical Ward identities corresponding to the above symmetries. From diffeomorphism invariance (and local Lorentz invariance) we have the conservation and symmetry of the stress-energy tensor: and from Weyl invariance we get the traceless property: When conformal anomalies are present, the expectation value of the stress-energy tensor no longer satisfies identity (2.6). Instead it is modified: where A(g µν ) is a local scalar function of the metric. The possible form of A is restricted by the Wess-Zumino consistency condition to a linear combination of possible expressions (see e.g. [5,[16][17][18][19] for details): the Euler density of the background manifold E d (Atype anomaly), the various Weyl invariant densities of the manifold (B-type anomalies) and other terms that can be cancelled by adding local counterterms to the curved space effective action of the theory (trivial terms). In d = 2 dimensions, the only such expression is the Euler density E 2 = R (where R is the Ricci scalar of the manifold), so that: (2.8) In d = 4 dimensions, these expressions consist of the Euler density E 4 (A-type anomaly), the Weyl tensor squared W 2 (B-type anomaly) and R (trivial term): While these expressions are universal, the value of their coefficients depends on the content of the theory. Note that the coefficients of the trivial terms may depend on the regular-ization scheme, but those of the anomalies are regularization independent. One method to obtain the value of these coefficients is to extract them from the flat space n-point correlation functions of the stress-energy tensor. For later reference, for a free scalar in d = 2 dimensions, the anomaly coefficient has been determined to be (see e.g. [20]): . (2.10) For a free scalar in d = 4 dimensions, the coefficients have been determined to be 1 (see [4,[21][22][23][24][25] and references therein): (2.11)

Ward Identities For Correlation Functions
We denote by W (g µν ) the effective action of the conformal field theory on curved spacetime, defined by: e i W (gµν ) ≡ Dφ e i S(gµν ) , (2.12) where φ stands for the dynamic fields in the theory. Two types of correlation functions can be defined for the stress-energy tensor: First, the connected Feynman n-point correlation functions, given by the expectation value 2 : . . . O n (x n ) e i S − Non-connected terms, (2.13) where O 1 , . . . , O n are local functions of the dynamic and the background fields, and Z 0 is the flat space partition function. Alternatively, one may define the correlation functions as variations of the effective action W with respect to the metric: T µ 1 ν 1 (x 1 ) . . . T µnνn (x n ) W ≡ (−i) n−1 2 n |g(x 1 )| . . . |g(x n )| δ n W δg µ 1 ν 1 (x 1 ) ...δg µnνn (x n ) . (2.14) We will refer to these as the variational correlation functions. Since the stress-energy tensor itself depends on the metric, these two types of correlation function do not coincide. One obtains the following relations between the two types of two point functions and three point functions in flat spacetime 3 (see [24]): T µν (x)T ρσ (y) W = T µν (x)T ρσ (y) F , (2.15) 1 The coefficient βc here corresponds to a dimensional regularization scheme. The values cited here of the coefficients β, βa, β b , βc correspond to our conventions for the definitions of the stress-energy tensor and the Riemann tensor. 2 In the following sections we may omit the F subscript when referring to Feynman correlation functions. 3 In deriving these formulas, one point functions have been dropped as they correspond to massless "tadpole" diagrams and therefore vanish in the flat space limit in a dimensional regularization scheme.
The Ward identities satisfied by the variational correlation functions can be obtained by taking derivatives of the curved spacetime Ward identities (2.5)-(2.6) (or (2.7) in the presence of conformal anomalies) with respect to the background metric. These identities have been derived in [24][25][26] for Euclidean signature. We repeat them here modified to our use of Lorentzian signature. For diffeomorphism invariance, one obtains the following identities for two and three point correlation functions: For the (anomalous) Weyl invariance, one obtains the following identities: (2.20) Using identities (2.19)-(2.20), one can extract the conformal anomaly coefficients from the flat space correlation functions.

Symmetries and Ward Identities In Curved Spacetime
Consider a non-relativistic field theory in d + 1 spacetime dimensions with a Lifshitz scale symmetry of the form (1.1). Assume this theory can be coupled to a curved spacetime manifold equipped with a metric g µν (or alternatively vielbein structure e a µ ) and a 1-form t α corresponding to the time direction at each point (or the normalized n α : n α n α = −1), as described in [6,7] 4 . Suppose this theory is described by the classical action S(g µν , t α , {φ}), where {φ} are the dynamic fields of the theory (or alternatively S(e a µ , t α , {φ}) in vielbein formalism). Further assume that this action can be defined such that it reduces to the flat space action of the theory when 5 g µν → δ µν = diag(−1, 1, 1, . . .), t α → (1, 0, 0 . . .), and is invariant under the following symmetries: 1. Time-direction preserving diffeomorphisms (TPD). These are the curved spacetime generalization of space rotations. As explained in [6,7], in our covariant notation these take the form of standard diffeomorphisms, given by: 21) or in vielbein notation (here we also include local Lorentz transformations): (2.22) 2. Anisotropic Weyl transformations. These are the local generalization of Lifshitz scaling, given by: where P µν = g µν + n µ n ν is the spatial projector, or alternatively using the vielbeins: (2.24) We define various field theory currents as the variation of the action with respect to the background fields. The stress-energy tensor can be defined in two possible ways, either using the metric or vielbein descriptions: (2.25) We also define the variation of the action with respect to the time direction 1-form: not coincide. Instead they are related by the following identity: From the symmetries and these definitions, one can obtain the classical Ward identities these currents satisfy over curved spacetime (see [6,7]). From TPD invariance follow the identities: (2.29) or equivalently in terms of T µν (e) : (2.30) These identities imply that in the flat spacetime limit, T µν (g) is symmetric but not conserved, while T µν (e) is conserved but not symmetric. The latter then corresponds to the conserved energy and momentum currents of the flat space theory, and we therefore choose to use it and the vielbein description in the following sections 6 .
From anisotropic Weyl symmetry, one obtains the following Ward identity: where D µν ≡ P µν − z n µ n ν . Note that similarly to the relativistic case, assuming that a Weyl-invariant coupling of the theory to curved spacetime exists is equivalent to assuming that the flat space stress-energy tensor can be improved to satisfy equation (2.31) (and still remain conserved).
In a field theory for which TPD invariance is not anomalous, the expectation value of the stress-energy tensor satisfies the Ward identities (2.30). However, when the Lifshitz scaling symmetry is anomalous, its expectation value no longer satisfies identity (2.31). Instead it is modified: where A(e a µ , t a ) is a TPD invariant and local function of the background fields. The possible form of A is restricted by the Wess-Zumino consistency condition to a linear combination of possible expressions (analogous to E 4 , W 2 and R of the (3+1)-dimensional conformal case), and was obtained for several cases in [6-8, 10, 11]. While these expressions are universal, the value of their coefficients depends on the content of the theory. Like the conformal case, the coefficients of the trivial terms may depend on the regularization scheme, but those of the anomalies are regularization independent. Our goal in this work is to extract these coefficients from the flat space n-point correlation functions of the stressenergy tensor.
For the case of a Lifshitz field theory in 2 + 1 dimensions with z = 2, which is the case we study in section 5, it was shown in [7] that generally (when the Frobenius condition is not assumed) the possible expressions in A consist of an infinite set of linearly independent anomalies and trivial terms, with increasing number of derivatives. However, the only anomalies that possibly contribute to the three point correlation functions of the stressenergy tensor (and therefore relevant to the calculations we perform here) are given by the following expressions: in the two derivatives sector, and: in the four derivatives sector. The various geometrical structures in these expressions are defined as follows: • (K S ) µν ≡ 1 2 L n P µν is a generalization of the extrinsic curvature of the foliation induced by n µ to the non-Frobenius case. We also define is antisymmetric and vanishes when n µ satisfies the Frobenius condition. We also define K A ≡ 1 2 K A µν¯ µν .
• a µ ≡ L n n µ is the acceleration vector associated with n µ .
•∇ µ is a space projected covariant derivative.
• R is a generalization of the Ricci scalar of the foliation induced by n µ to the non-Frobenius case.
For further discussion and more detailed definitions see [7] and appendix B. The possible trivial terms relevant to our calculations are listed in appendix E.

Ward Identities For Correlation Functions
Let us denote by W (e a µ , t a ) the effective action of the Lifshitz field theory on curved spacetime, defined by: where φ again stands for the dynamic fields in the theory. We define two types of correlation functions for the stress-energy tensor. First, the connected Feynman n-point correlation functions, defined similarly to the relativistic case (2.13), where O 1 , . . . , O n are local functions of the dynamic and the background fields, and Z 0 is the flat space partition function. Alternatively, we define the correlation functions given by variations of the effective action W with respect to the background fields: . (2.36) We will refer to these as the variational correlation functions. Like the relativistic case, these two definitions for the correlation functions do not coincide, since the curved space field theory currents depend on the background fields. Instead these two types of correlation functions are related via relations, that are obtained by differentiating equation (2.35). For the two and three point functions, one obtains: (2.38) In the flat space limit, these relations reduce to the following: 7 These relations make it clear that, while the flat space Feynman correlation functions depend only on the field theory currents (given as the first derivative of the action with respect to the background fields), the variational correlation functions depend on higher derivatives of the action as well. The implication is that, if there is more than one way to couple the field theory to curved spacetime while preserving all symmetries discussed in subsection 2.2.1 and leaving the flat space currents the same, then any such coupling would produce the same Feynman correlation functions but different variational correlation functions. As we will show, this leads to a possible ambiguity in the relation between the flat space Feynman correlation function and the consistent anomalies obtained from the Wess-Zumino analysis, which does not occur in the relativistic case. One method to derive Ward identities directly for the flat space Feynman correlation functions is via a change of variables in the path integral. Suppose that the operator δ implements an infinitesimal transformation that corresponds to a symmetry of the theory, and transforms both the dynamic and the background fields so that δS = 0. Let δ dyn denote an operation that transforms only the dynamic fields, leaving the background fields unchanged, and similarly let δ bg transform only the background fields, so that the following is satisfied: δ dyn φ = δφ, δ dyn e a µ = 0, δ dyn t a = 0, and δ = δ dyn + δ bg . These operators are explicitly given by: (2.42) By performing a change of variables φ →φ = φ + δφ in the path integral (2.13), and assuming for now that that the symmetry corresponding to δ is non-anomalous (so that the change of integration measure does not contribute in the context of the dimensional regularization scheme we are using here), we obtain: Noting that δ dyn S = −δ bg S, the following identity can then be derived: Using TPD (as given in (2.22)) in identity (2.44) and choosing ξ σ (w) = δ σ a δ(w − x), we get the following Ward identities for the flat space two and three point Feynman correlation functions of the stress-energy tensor: 8 where for a scalar field φ, for example, we have (δ D dyn φ)(x, w) = δ(w − x)∂ a φ(w). Similarly, when applying identity (2.44) to the anisotropic Weyl transformation (2.24) and setting σ(w) = δ(w − x), we obtain the following identities for the flat space correlation functions: where for a Lifshitz scalar 9 In the presence 8 See appendix A.3 for details on our notations for the Ward identities. 9 Here and in the following sections, "Lifshitz scalar" will refer to a scalar φ with a second order time derivative kinetic term in the action as used in section 5, so that it has a Lifshitz dimension of [φ] = d−z 2 .
of Lifshitz scaling anomalies, expressions (2.47)-(2.48) will not vanish, and in general will be equal to some linear combinations of contact terms instead.
An alternative method to derive the Ward identities for the variational correlation functions is by taking derivatives of the curved spacetime Ward identities (2.30) and (2.31) (or (2.32) in the presence of Lifshitz anomalies) with respect to the background fields. This method allows one to directly relate the flat space correlation functions to derivatives of the consistent anomalies as derived in [7], and extract the anomaly coefficients. Taking the first and second derivatives of identity (2.32) in the flat space limit, we obtain the following Ward identities for the two and three point variational correlation functions: (2.50) Using the relations (2.39)-(2.40) in these expressions, W W can be written as linear combinations of the expressions in the Ward identities (2.47)-(2.48) and other expressions that vanish in the absence of Lifshitz anomalies: bc (x, y, z), where J ρα bc (x, y, z) and K ρα bc (x, y, z) are given by: When no anomalies are present, the expressions I W and K are expected to vanish due to identity (2.44). The expression J can also be shown to vanish in this case, by noting that: 10 ( and using the following expression for (δ W bg T ρ b )(x, y) (where again σ(w) = δ(w − x)): (2.56) and the anisotropic Weyl scaling properties of T ρ b : This scaling property can be derived by using the definition of the stress-energy tensor (2.25), expressing the operator δ W in terms of the anisotropic Weyl transformations of the background and dynamic fields as in (2.42), exchanging the order of derivatives and using the anisotropic Weyl invariance of the classical action S. 11

Ambiguity of the Anomaly Coefficients?
As mentioned earlier, unlike the relativistic case, the relation between the flat space Feynman correlation functions and the Lifshitz anomaly coefficients (that is, the coefficients on the RHS of equation (2.32) of the various anomalous terms as obtained from the Wess-Zumino consistency conditions) may be ambiguous in some cases. This is due to the Weyl invariant coupling of the theory to curved spacetime being non-unique.
As an example, consider any Lifshitz invariant theory for which d = z, that contains a real Lifshitz scalar φ, so that the field φ is dimensionless. In the absence of anomalies, one expects the two point function D(x)φ 2 (y) F to vanish. When Lifshitz anomalies are present, it will instead be equal to some contact term, and from dimensional analysis we have: where c is some constant. Suppose that A 0 (e a µ , t a ) is any local functional of the background fields with the following 2 properties:

2.
A 0 is a Weyl invariant density of dimension d + z, that is: Next, consider adding to the curved spacetime classical action S(e a µ , t a ) a term of the form: S 0 ≡ d d+1 w eA 0 (w)φ 2 (w). (2.60) The new actionS = S + S 0 is still invariant under TPD and anisotropic Weyl transformations, and it coincides with S in flat spacetime. Moreover, the flat spacetime conserved currents derived fromS (including the stress-energy tensor) are the same as those derived from S. As a result, the anomalous Ward identity expressions I W (as defined in (2.47)-(2.48)) remain unchanged in flat spacetime. However, expression W W does change. Due to the assumed properties of A 0 , the term S 0 satisfies: and therefore does not contribute to J (as defined in (2.53)) either. The only contribution of S 0 to the anomalous identity W W is from the expression K (as defined in (2.54)). It follows that the change in W W due to the S 0 term is given by: . (2.62) If we choose A 0 = aÃ, whereÃ is a possible B-type anomaly of the theory and a is a constant, we conclude that S 0 contributes an additional −ac to the coefficient of theÃ anomaly, without changing the action or the current operators of the flat space theory. 12 The freedom to add such a term to the action and thereby change the corresponding anomaly coefficient was previously pointed out in [27] for a purely spatial anomaly of a free Lifshitz scalar in d = z = 2, but it is in fact more general. It could be applied, for example, to any of the consistent anomalies in the d = z = 2 case (which are all B-type, see [6,7,10,11]), for any theory that contains a Lifshitz scalar. This freedom may suggest that, unlike the relativistic case, knowing the flat space action, the flat space currents and their Feynman correlation functions is not enough in these cases to determine the consistent anomaly coefficients -one needs to specify the full curved spacetime action, and by choosing different couplings of the theory to curved spacetime one may obtain any value for any of the B-type anomaly coefficients.
Stated differently, knowing all the Feynman correlation functions of the flat space field theory is not sufficient in order to determine the curved spacetime action. In order to construct the full curved spacetime action we have to know all the variational correlations functions and the ambiguity is in the relation between these two types of correlation functions, the variational and Feynman. This ambiguity can be avoided if we add another ingredient to the discussion. The φ field in (2.60) is a Log-correlated field and is therefore ill defined when we take the large volume limit. Consistency of the quantum field theory at infinite volume forbids such an operator in the correlation functions. 13 The ambiguity may still have consequences for field theories on a finite volume spacetime, or with other modifications of the IR physics that take care of the Log divergence. We leave this for future studies.
In the following sections, we use split dimensional regularization to calculate the Lifshitz anomaly coefficients for a free Lifshitz scalar in 2 + 1 dimensions and z = 2. The above discussion implies that a specific coupling of the theory to curved spacetime needs to be specified. However, since we are performing the calculation for a free theory in an infinite volume and no physical IR regulator, we will use the minimal coupling (which is Weyl invariant in this case) and will not allow for Weyl invariant couplings of the form (2.60).

The Regularization and Renormalization Method
In this section we present the split dimensional regularization scheme we employ and explain the method by which we use it to extract the Lifshitz anomaly coefficients from the field theory correlation functions. We start by reviewing the relativistic conformal case for reference, and then explain the non-relativistic Lifshitz case.

Conformal Anomaly From Dimensional Regularization
In the standard relativistic dimensional regularization scheme, one defines the theory and calculates various quantities in a general dimension d, and then analytically continues the obtained expressions to dimension d = d phys − ε, where d phys is the physical dimension. Suppose that I (n) (d, p i , m) is some n-point correlation function written in momentum space and calculated to one-loop order in perturbation theory using the corresponding 1PI Feynman diagrams, where p µ i (i = 1, . . . , n) are external momenta and m is an IR mass regulator. Generally after analytic continuation of the dimension, I (n) will take the form: 14 is an expression which is regular around ε = 0, and is a linear combination of terms of the form: 2) 13 We thank Z. Komargodski for this comment. 14 Beyond one-loop order, one has to first cancel the possible subdivergences using the appropriate counterterms. The pole in ε can then be of higher order. In the following sections we focus on a free theory, and therefore on the one-loop case.
where g(ε, p i , m) is a scalar expression. Expanding around ε = 0 we have: where I (res) (n) is the residue of the ε pole, and I is the renormalized correlation function. 15 It is a well known property of relativistic field theories that the residue I (res) (n) is always a polynomial in the external momenta and the mass regulator (see for example [28][29][30]). This can be shown by taking derivatives of I (n) with respect to the external momenta (and mass regulator) enough times so that the corresponding Feynman diagram no longer diverges. As long as there are no IR divergences in the physical dimension, one can safely take the limit m → 0 in I (ren) (n) to obtain the physical renormalized correlation function 16 (the correlation functions of the stress-energy tensor in the cases studied here are indeed free of IR divergences, even in d = 2 dimensions, as will be explained in subsection 4.1.2).
A useful property of scale anomalies is that, in some cases, one can calculate them from the ε pole residue alone: Suppose the theory has a symmetry (such as a scaling symmetry) that is not explicitly broken by the dimensional regularization scheme itself, with a corresponding Ward identity of the form: where {I k } is a set of correlation functions and T (ε) is a linear operator that takes expressions of the form (3.2) to expressions of the same form, and may or may not depend on the dimension. Since the symmetry is not broken by dimensional regularization, the unrenormalized correlation functions satisfy identity (3.4). Therefore we can deduce: The anomalous Ward identity is then given by: Since the LHS of equation (3.6) is finite, we can immediately draw two conclusions from it: 2. If T does not depend on ε then T I (res) k = 0 and there is no anomaly. For example, the operator T corresponding to diffeomorphism invariance does not explicitly introduce new factors that depend on d, and therefore as long as the dimensional regularization itself does not break this symmetry, it will not be anomalous.

Expansion in the External Momenta
Equation (3.6) allows one to calculate the anomalous Ward identity from the ε pole residue. This is useful, since there is no need to calculate the full correlation functions in order to extract their divergent part -it can be obtained simply by expanding the Feynman diagram integrand in powers of the external momenta. Suppose the integrand is h(p i , q, m), where q is the internal loop momentum and h has a mass dimension l and therefore satisfies: Rescaling q and m by a factor of 1/λ where λ → 0, we get: We next expand in powers of λ around λ = 0 to obtain: 17 Thus in order to calculate the pole residue, the only integrals left to evaluate are the ones over the expressions (h (k) ) i 1 ...i k µ 1 ...µ k (q, m). Each of these expressions is a linear combina-tion of terms of the form: When performing the integration over q we may use the Lorentz symmetry of the integral to make the standard replacement: ..µ 2n is a sum of all possible unique products of n metric factors, i.e.: (3.14) Using (3.13), the integral overh(p i , q, m) can be written in terms of the following known integral (see e.g. [30][31][32]): where k = l − 2r + 2s ≤ d phys + l (so that d phys + 2r − 2s ≥ 0). Note that the integral (3.15) contributes to the pole residue a term proportional to m 2r−2s+d phys . Since terms with 2r − 2s + d phys > 0 vanish in the m → 0 limit, one only needs to compute terms in the expansion with 2r − 2s + d phys = 0 (those of order k = d phys + l) in order to obtain the ε pole residue. For these terms, the integral (3.15) has the following pole: Using this procedure of expanding the Feynman diagram integrand in external momenta, computing the pole residue via equations (3.13)-(3.16) and applying formula (3.6), one can calculate the anomalous Ward identities corresponding to scale symmetry without calculating the full correlation functions. This will be especially useful in the Lifshitz case, where the Feynman diagrams are more difficult to fully evaluate.

Lifshitz Anomalies From Split Dimensional Regularization
In the non-relativistic case, in analogy to the relativistic case, we use a split dimensional regularization scheme (first suggested in [13,14] for regularizing gauge theories in the Coulomb gauge). We follow the general scheme defined in [15], in the context of nonrelativistic field theories. We start by defining the theory in a general number of time dimensions d t and space dimensions d s (while keeping the critical dynamical exponent z constant).
In flat spacetime, the theory is defined on a manifold M = M t × M s , where M t is a d t -dimensional time manifold and M s is a d s -dimensional space manifold, such that is invariant both under rotations in the time manifold and in the space manifold separately. Spacetime coordinates will be denoted by x µ = (xμ, xμ) whereμ = 1, . . . , d t are time indexes,μ = 1, . . . , d s are space indexes and µ = 1, . . . , d t + d s are spacetime indexes. We define a flat metricδμν = diag(−1, . . . , −1) on M t , andδμν = diag(1, . . . , 1) on M s . We also define the time projector on M asδ µν = diag(δμν, 0) and similarly the space projector asδ µν = diag(0,δμν), so that δ µ ν =δ µ ν +δ µ ν . Given a vector v µ on M, we denote its time projection byv µ ≡δ µ ν v ν , and its space projection byv µ ≡δ µ ν v ν . 19 Similarly to the relativistic case, we calculate various expressions for general d t and d s values, and then analytically continue them to non-integer is some n-point correlation function calculated to one-loop order in perturbation theory using the corresponding 1PI Feynman diagrams, where p µ i (i = 1, . . . n) are external momenta and m is an IR mass regulator. Generally after analytic continuation of the dimensions, I (n) will take the form (see [15]): is an expression which is regular around ε t = ε s = 0 and is a linear combination of terms of the form: where g(ε t , ε s , p i , m) is a scalar expression (with respect to time and space rotations).
Since there are two different dimensional regularization parameters in this case, unlike the relativistic case, in order to renormalize the expression one must choose a particular way in which one takes the limit (ε t , ε s ) → (0, 0), with each choice leading to a different renormalized expression. Defineε(ε t , ε s ) to be some coordinate on the two dimensional regularization parameter space such thatε(0, 0) = 0 and such that the transformation in ε lif while keepingε constant, we have: where I (res) (n) is the residue of the ε lif pole, and I (ren) is the renormalized correlation function. Like the relativistic case, the residue I (res) (n) (ε) is a polynomial in the external momenta and the mass regulator for any value ofε (see [15]), and therefore represents contact terms in coordinate space. This can be shown by taking derivatives of I (n) with respect to the external momenta (and mass regulator) enough times so that the corresponding Feynman diagram integral no longer diverges. As long as there are no IR divergences in the physical dimension, one can safely take the limit m → 0 in I (ren) (n) to obtain the physical renormalized correlation function (the correlation functions of the stress-energy tensor in the case studied here are indeed free of IR divergences, as will be explained in section 5).
As mentioned earlier, the renormalized function I (ren) (n) depends on the choice of the parameterε that is kept constant as we take the limit (ε t , ε s ) → (0, 0). Suppose we instead choose a different parameterε (ε t , ε s ), that corresponds to the renormalized correlation function I (ren) (n) . Then we have the following relation between I : differ from each other by contact terms, as expected from a change in the renormalization scheme. Also note that the possible change in the renormalized expressions as a result of the choice of ε is completely described by the single parameter α. In order to account for all possible choices, we leave α as a free parameter in our calculations, and use the following choice of ε:ε which is chosen such that moving fromε toε = ε t we get the factor ∂ε The inverse transformation from (ε lif ,ε) to (ε t , ε s ) is given by: As in the relativistic case, we can calculate the Lifshitz scale anomaly coefficients from the ε lif pole residue. Suppose the theory has a symmetry (such as a Lifshitz scale symmetry) that is not explicitly broken by the split dimensional regularization scheme itself, with a corresponding Ward identity of the form: where {I k } is a set of correlation functions and T (ε lif ,ε) is a linear operator that takes expressions of the form (3.18) to expressions of the same form, and may or may not depend explicitly on the time and space dimensions. Since the symmetry is not broken by the regularization, the unrenormalized correlation functions satisfy identity (3.23). We therefore have from (3.19): The anomalous Ward identity is then given by: Since the LHS of equation (3.25) is finite, we can again draw the following conclusions: 2. If T does not depend on ε lif andε then T I (res) k (ε) = 0, and there is no anomaly. Therefore, as long as TPD invariance is not explicitly broken by the split dimensional regularization scheme (as is the case with the free scalar we consider in the following sections), we don't expect it to be anomalous. 20 It is important to consider the consequences of changing the choice of the parameterε (and thereby the renormalization) on the anomalous Ward identity (3.25). By changing our choice fromε toε, we know from (3.20) that the change in the anomalous Ward identity is given by: Since ∂f ∂ε (0, 0) is a local expression (a polynomial in the external momenta), the term T (0, 0) ∂f ∂ε (0, 0) represents a trivial solution of the WZ consistency condition (one that can be cancelled by a local counterterm). We therefore expect only the coefficients of trivial terms to depend on α. This is consistent with the general expectation that only coefficients of trivial terms can be regularization dependent.

Expansion in the External Momenta
Equation (3.25) allows us to calculate the anomalous Ward identity from the ε lif pole residue. This is especially useful in the Lifshitz case, since the denominators of the propagators are generally polynomials of degree 2z, and the Feynman diagram integrals are therefore more difficult to fully evaluate than they are in the relativistic case. Their divergent parts, however, can again be obtained simply by expanding the Feynman diagram integrands in powers of the external momenta. Suppose the integrand is h(p i ,p i ,q,q, m), where q is the internal loop momentum and h has a Lifshitz dimension 21 l and therefore satisfies: Rescalingq by a factor of 1 λ , andq and m by a factor of 1 λ z (where λ → 0) we get: We next expand h(λ zp i , λp i ,q,q, m) in powers of λ to obtain: is a polynomial in the external momenta, given by: we conclude that the integral over h(p i ,p i ,q,q, m)−h(p i ,p i ,q,q, m) has a divergence degree of zd t + d s + l − k 0 − 1. If we choose k 0 = d phys s + z + l , the integral over h −h converges in d phys s + 1 dimensions, and therefore the integrals over h(p i ,p i ,q,q, m) andh(p i ,p i ,q,q, m) have the same ε lif pole.
Thus in order to calculate the pole residue, the only integrals left to evaluate are the ones over the polynomial coefficients in the expressionsh(p i ,p i ,q,q, m). Each of these is a linear combination of terms of the form: When performing the integrations overq andq, we may use the time rotation and space rotation symmetries to make the replacements given in equations (3.13)-(3.14) separately forq products (using theδ µν metric) and forq products (using theδ µν metric). The integral overh(p i ,p i ,q,q, m) can then be written in terms of the following known integral (see [15]): . Note that the integral (3.33) contributes to the pole residue a term proportional to m 2r−2I+2s/z+1+d phys s /z . Since terms with 2r−2I +2s/z+1+d phys s /z > 0 vanish in the m → 0 limit, one only needs to compute terms in the expansion with 2r−2I +2s/z+1+d phys s /z = 0 (those of order k = d phys s + z + l) in order to obtain the ε lif pole residue. For these terms, the integral (3.33) has the following pole: Note that since r, s and I appear in the correlation functions as non-negative integers, it is possible to relate the different poles X(d t , d s , r, s, I, z) appearing in the correlation function of a given order for a fixed value of z and with various values of r, s, I using the recursive property of Gamma functions Γ(n + 1) = nΓ(n).
In conclusion, this procedure of expanding the Feynman diagram integrand in external momenta, computing the pole residue via equations (3.13)-(3.14), (3.33)-(3.34) and applying formula (3.25) enables us to calculate the anomalous Ward identities corresponding to Lifshitz scale symmetry without calculating the full correlation functions. In the following sections we use this procedure to calculate the anomaly coefficients for the case of a free z = 2 Lifshitz scalar in 2 + 1 dimensions.

The Conformal Scalar Field and its Weyl Anomalies
In this section we review the calculation of the conformal anomaly coefficients for the relativistic free and massless scalar field using dimensional regularization and the procedure described in section 3. We include this example as a reference for the calculation of Lifshitz anomaly coefficients for the non-relativistic Lifshitz scalar given in section 5 using a similar procedure. The calculation was performed both for two and four spacetime dimensions. In both cases the Weyl anomalies agree with the known results found in literature [4,[20][21][22][23][24][25].
In order to calculate the conformal anomalies for the free scalar, one must first define the flat space stress-energy tensor such that it satisfies the Ward identities (2.5)-(2.6) (i.e. it is conserved, symmetric and traceless). One way to do this is to define the theory over a curved spacetime manifold such that the action is invariant under both diffeomorphisms and Weyl transformations (see section 2), and derive the stress-energy tensor from the curved spacetime action using the definition (2.4).
The conformal coupling of a free relativistic scalar field φ to curved spacetime is given by the following action (see e.g. [23]): where d is the number of spacetime dimensions, and R is the Ricci scalar of the background manifold. This action in indeed diffeomorphism and Weyl-invariant. The improved stressenergy tensor calculated from this action using the definition in equation (2.4) is (see [5]): Using the equations of motion, one can check that it indeed satisfies the conservation and tracelessness Ward identities: for any number of spacetime dimensions d. 22

The Relativistic Free Scalar in Two Dimensions
We start with the calculation of the Weyl anomaly of a relativistic free scalar field in two spacetime dimensions from the two point correlation function of the stress-energy tensor. We demonstrate two ways of performing the calculation: First, by preforming the full calculation of the one-loop diagram using dimensional regularization. Second, using the procedure of extracting only the divergent part of the diagram and using it to compute the anomaly coefficients, as described in section 3. The results agree with the known ones from the literature.

The Full Calculation
In the first way of calculating the anomaly, the two point correlation function of the stressenergy tensor is fully calculated using standard dimensional regularization. The Feynman rules and diagram used for the calculation are given in appendix C.1.
The full evaluation of the expression that corresponds to the diagram (C.2) was performed using the massless integral formulas given in appendix C.2. The final result for the two point correlation function is given by the following expression: where p is the external momentum of the diagram, ε is defined by d = 2 − ε and the basic integral I 1 is given by: (4.6) One can expand the integral I 1 around the physical dimension d phys = 2 to obtain: It is easy to verify that the result in equation (4.5) indeed satisfies the conservation Ward identity (2.17), which in Fourier space takes the form: 23 As expected, this identity holds separately on the finite part and on the pole part of (4.5) and is not anomalous. In addition one can verify that, when tracing over the full unrenormalized expression (4.5), the Ward identity corresponding to Weyl invariance is satisfied: However, this identity is only satisfied on the pole and finite parts of expression (4.5) together. Thus after performing renormalization of the correlation function we have from equation (3.6): The ε pole residue of expression (4.5) is given by: (4.11) Tracing over (4.11) (in d = 2 − ε dimensions) and using equation (4.10) we get: which is the well known anomalous Ward identity in two dimensions (see e.g. [16,20]).

Poles Calculation
Using the procedure described in subsection 3.1, the pole residue of the two point correlation function of the stress-energy tensor can also be computed without evaluating the full expression. Starting from the expression for the two point function given in (C.2), we expand the integrand in powers of the external momentum and extract the terms proportional to m 0 (where m is an IR mass regulator). The expression for the two point function contains only five different types of integrals in this case. Their relevant ε poles around two spacetime dimensions, obtained using expansion in the external momentum p, are given by the following expressions: (4.14) where: Note that although some of these expressions (4.13)-(4.17) are individually IR divergent at d = 2, the total expression for the two point function of the stress-energy tensor (and indeed any correlation functions of the stress-energy tensor) is not IR divergent in this case. This follows from the form of the action (4.1) and the improved stress-energy tensor (4.2): The action at d = 2 contains only derivatives of the field φ, and only the O(ε) part of it contains φ with no derivatives. It follows that the stress-energy tensor (and any other operator constructed by taking variations of the action with respect to the metric) has the same structure -at d = 2 it contains only derivatives of φ. These operators are therefore not Log correlated at d = 2, and their correlation functions do not diverge in the IR. In terms of dimensionally regulated Feynman diagrams, for any propagator in the 1-loop diagram with momentum q µ , the d = 2 part of each of its adjacent vertexes is of order O(q). Thus the terms in the integrand of the diagram that contribute to the IR divergence when q → 0 are of order O(ε 2 ), and therefore vanish as ε → 0 after integration is performed. This justifies the assumption made in subsection 3.1, that the correlation functions of the stress-energy tensor are free of IR divergences. A similar argument will be valid for the Lifshitz case, as will be explained in the next section.
Substituting these poles in the general expression for the two point function (C.2), we find exactly the same pole residue for the two point correlation function as in equation (4.11), and therefore, by the use of (4.10), the same anomalous Ward identity.

The Relativistic Free Scalar in Four Dimensions
We next turn to the example of a relativistic free scalar field in four spacetime dimensions. Again using dimensional regularization and the procedure described in subsection 3.1 with d = 4 − ε, the anomalous contributions to the Ward identities (2.19) and (2.20) can be computed from the ε pole residues of the two and three point correlation functions of the stress-energy tensor.
The details needed for the calculation of the correlation functions as defined in equation (2.14), including the Feynman rules for the vertexes and the expressions for the diagrams, are given in appendix C.1. The pole residues were extracted via a power expansion in the external momenta as described in subsection 3.1. Since these calculations involve a very large number of terms, they were performed using a computer script that was written for this purpose.
The pole residue of the two point correlation function, as obtained from the power expansion in the external momentum p, is given by: (4.20) It can be easily checked that this expression satisfies the conservation identity (4.8). However, tracing over the indexes µ, ν in (4.20) yields an expression that does not vanish for a general dimension d. Using equation (3.6), one obtains: This is in agreement with the anomalous Ward identity (2.19), due to the trivial R term in the Weyl cohomology of the relativistic theory in four dimensions (see e.g. [24,25]).
The pole residue of the three point correlation function can also be calculated using expansion in the external momenta. However, the final result is too long to be shown here. As expected, this pole residue and the pole residue of the two point function together satisfy the following conservation identity: where δ 2 |g|A (q, p) ρσαβ represents the Fourier transformed second variation of the anomalous contribution to the Ward identity (2.7) with respect to the metric (see appendix A.2 for notations). After calculating the second variation of the anomaly and trivial term densities listed in equation (2.9) with respect to the background metric, the results can be substituted into the RHS of the anomalous Ward identity (4.23).
The coefficients β a , β b and β c can then be extracted by comparing the two sides of the identity. The values of the coefficients obtained using this procedure are the same as the ones given in (2.11) and therefore agree with those found in the literature.

The Free z = 2 Lifshitz Scalar Field and its Scale Anomalies
In this section we study the Lifshitz anomaly of a Lifshitz z = 2 scalar field in 2 + 1 dimensions, using the method of split dimensional regularization and renomarlization described in subsection 3.2, applied to the two and three point correlation functions of the stressenergy tensor. Up to second order in the background fields, we find one anomaly in the two derivatives sector, and no anomalies in the four derivatives sector. This is in agreement with the results previously found in [11] using a heat kernel calculation. The value of the anomaly coefficient also agrees with the result in [11]. As expected, only the coefficients of the trivial terms which appear in the two and three point correlation functions depend on the regularization parameter α, defined in subsection 3.2 -they are all regularization dependent, and can be removed by adding the appropriate counterterms to the effective action. We also show that trivial terms that correspond to a curved background structure that violates the Frobenius condition (that is, terms that vanish when the Frobenius condition is assumed) appear in the case we study with non-vanishing coefficients. Therefore, the curved spacetime description of these terms, and their cancellation via a counterterm, requires giving up the foliation structure of the background manifold (see [7] for further discussion).

The z = 2 Free Lifshitz Scalar Field in General Spacetime Dimensions
The Lifshitz anomalies of a free Lifshitz scalar field with a dynamical critical exponent z = 2 in 2 + 1 dimensions have been considered in several previous works [11,27]. The flat space action of the theory is given by: In order to apply the method of split dimensional regularization and renormalization described in subsection 3.2, one must first couple the theory described by the action (5.1) to a curved background manifold with a general number of space dimensions d s and time dimensions d t . This coupling must be done in a way that preserves the curved spacetime symmetries detailed in subsection 2.2.1, as well as local rotations of the time directions, and smoothly reduces back to (5.1) in the limit of flat spacetime and d s → 2, d t → 1.
This can be accomplished in several ways (see appendix B for details). Here we choose to use the following curved spacetime action: where d t and d s are the numbers of time and space dimensions respectively and {n (i) µ } (i = 1, . . . , d t ) is an orthonormal set of 1-forms that represents the local d t time directions on the manifold. The coefficients ξ 1 , ξ 2 , ξ 3 , ξ 4 are given by: is any orthogonal matrix in d t dimensions that depends on the spacetime coordinates x.
The action (5.4) is invariant under time and space rotations, as well as the following Lifshitz scaling transformation: where λ is a parameter of the scaling transformation, and κ is dimensionless under the scaling transformation. The flat space stress-energy tensor, as derived from the action (5.2) by taking the variation of the action with respect to the vielbeins (according to the definition (2.25)), is given by: Note that this expression is symmetric in its spatial components (Tμν = Tνμ), and its temporal components (Tμν = Tνμ), indicating both time and space rotations invariance, but not in its combined space/time components (Tμν = Tνμ), since there is no Lorentz inveariance. It is also regular in the limit d s → 2, d t → 1. These two properties are crucial for the assumptions underlying the procedure outlined in subsection 3.2 to be satisfied. 24 Using the equations of motion in flat space, given by: it is easily verified that the stress-energy tensor (5.6) satisfies the following Ward identities for any values of d t and d s (these are just the flat space versions of identities (2.30)-(2.31), generalized to d t > 1): In the particular case of d s = 2, d t = 1, the stress-energy tensor (5.6) reduces to the following expression: (5.10)

The Lifshitz Anomaly of the Free Scalar
In this subsection, we apply the procedure outlined in subsection 3.2 to calculate the Lifshitz anomaly coefficients of the free Lifshitz z = 2 scalar field in 2 + 1 dimensions, using a split dimensional regularization scheme. The general steps of this calculation are as follows: 1. We consider the Feynman diagrams contributing to the flat space variational two and three point correlation functions of the stress-energy tensor, as given in (2.39)-(2.40), in a general number of time and space dimensions.
2. For each of these diagrams, we extract the ε lif pole residue of the diagram by expanding the integrand in powers of the external momenta as explained in subsection 3.2, utilizing the formulas (3.13)-(3.14) and (3.33)-(3.34) (Note that we use the choice of ε given in (3.21)).
3. We substitute the pole residues into the LHS of Fourier transformed versions of the Ward identities (2.49)-(2.50), and then use formula (3.25) (taking the limit (ε lif ,ε) → 0) to obtain the anomalous contribution to these Ward identities. 4. We calculate the first and second order variation of each of the possible independent anomaly and trivial term densities corresponding to this case 25 , as listed in [7], with respect to the vielbeins in the flat space limit. Note that we do not consider any n-point functions with n > 3 in this work, and therefore only the coefficients of terms which are at most second order in the background fields can be calculated.
5. We substitute these expressions into the RHS of the Ward identities (2.49)-(2.50), compare to the results of step 3, and extract the various coefficients. 26 For the purpose of calculating the pole residues of the correlation functions of the stressenergy tensor, it is convenient to drop terms in (5.6) that carry coefficients proportional to the parameter ε lif . These terms will only contribute expressions which are regular in the limit (ε lif ,ε) → 0 to the correlation functions. Therefore, when calculating these pole residues, one can ignore such terms. After dropping these terms, the stress-energy tensor takes the form: The expression for the Feynman diagram vertex that corresponds to the stress-energy tensor (5.11) is given in (D.3). All other Feynman rules needed for the calculation of the relevant Feynman correlation functions can be found in the appendix D. Note that although the correlation functions of the stress-energy tensor in this case seem like they might be IR divergent, similarly to the d = 2 relativistic case as discussed in subsection 4.1.2 this is in fact not the case. The terms in the action (5.2) of order O(ε 0 lif ) contain only derivatives of the field φ, and only the terms of order O(ε lif ) and higher contain φ with no derivatives. The stress-energy tensor (5.11) and other operators defined as variations of the action with respect to the vielbeins (such as δ 2 S δe b αδe c γ ) therefore have the same structure. These operators are therefore not Log correlated at 2 + 1 dimensions, and their correlation functions are not expected to diverge in the IR. In terms of the corresponding Feynman diagrams regulated using split dimensional regularization, for any propagator in the 1-loop diagram with spacetime momentum Q µ , the O(ε 0 lif ) part of each of its adjacent vertexes is of order O(Q) or O(Q). The terms in the integrand of the diagram that contribute to the IR divergence when Q → 0 are of order O(ε 2 lif ), and therefore vanish as ε lif → 0 after integration is performed. This justifies the assumption made in subsection 3.2, that the correlation functions of the stress-energy tensor are free of IR divergences.
Since these calculations involve a very large number of terms, they were performed using a computer script that implements the previously described steps. For the purpose of presenting the results of the calculations, we introduce the following notations for the Fourier transformed versions of the Ward identities of the flat space correlation functions: 27 (2π) dt+ds δ (P + Q + K) I W ) ρα bc (x, y, z) are divided into two separate sectors according to the total number of derivatives: a two derivatives sector (n D = 2), and a four derivatives sector (n D = 4). This is in accordance with the definitions and the discussion in [6,7] (higher derivative sectors do not appear in the two and three point functions).
At this point, we would like to reiterate that since we study the correlation functions only up to the three point function level, we are able to study only the anomalies and trivial terms (as found in [7]) that appear in this level, that is, only the ones which are at most second order in the background fields.

Results for the Two Point Function
The pole residues of the two point correlation function of the stress-energy tensor were extracted using the method described in subsection 3.2. Taking the weighted Lifshitz trace over these pole residues and following the previously mentioned steps of the calculation yields the following result for the anomalous contribution to the two point Ward identity: 28 Note that although the expression on the LHS of (5.16) seems to depend on two different momenta P and Q, due to the delta function δ(P +Q) in the definition (5.13) these momenta are not independent, and the RHS of (5.16) involves only one independent momentum. The same comment holds for other similar identities in this section. When comparing this result to the first order variation (in the flat space limit) of the cohomologically trivial terms 29 (or coboundaries in cohomological terminology) in the two derivatives sector (E.9) and in the four derivatives sector (E.31) with respect to the vielbeins, one finds the following results: In the two derivatives sector we have: where [δX] αβ (P ) represents the Fourier transform of the first order variation of the expression X with respect to the vielbeins in flat space, as defined in (E.1), and α is the regularization parameter defined in subsection 3.2.
In the four derivatives sector we get: where C 0 is a free parameter, whose value cannot be extracted from the two point correlation function. This is due to the fact that, in the basis of first order trivial terms we are using here, the first order variations of the trivial term densities in the four derivatives sector are linearly dependent in the flat space limit (that is, there is a linear dependence between the first order variations of F 1 , F 3 and F 6 with respect to the vielbeins in flat space). When looking at higher point correlation functions, this dependence is removed and the coefficient C 0 can be extracted, as we indeed show in the next subsection studying the three point function.
We have also confirmed the pole residue of the two point function satisfies the conser- 28 The results for the anomalous Ward identities in this section are in dt = 1 time dimension. In this case, we useP andK to denote the time components of the momenta P and K respectively (see appendix A.3). 29 Only trivial terms are first order in the background fields.
vation Ward identity: This is expected from the argument made in subsection 3.2 that the conservation Ward identity (2.45) holds separately on the pole part and on the regular finite part of the correlation functions.

Results for the Three Point Function and the Anomaly
Using the previously mentioned calculation steps, we obtained the following result for the anomalous contribution to the three point Ward identity in the two derivative sector 30 Writing this expression as a linear combination of the second order variations with respect to the vielbeins of the anomaly and trivial term densities in the two derivative sector (given in (E.11)-(E.12)), we get: where [δ 2 X] αγ ef (P, K) represents the Fourier transform of the second order variation of the expression X with respect to the vielbeins in flat space, as defined in (E.2). Note that, as expected, the coefficient of the anomaly A (2,0,0) is independent of the regularization parameter α and agrees with the result found in [11], while the coefficient of the trivial term F 1 depends on α, and vanishes when α = 3 2 . The full result for the anomalous contribution to the three point Ward identity (2.50) in the four derivative sector is given in (E.33) in the appendix. Comparing this result to the second order variations of the anomaly and trivial term densities in this sector we get the following linear combination: αγ ef (P, K) + C δ 2 F 12 αγ ef (P, K), where the trivial terms (the F i -s) are defined in (E.31) and the anomalies are defined in (2.34) or (E.32). C is again a free parameter that cannot be extracted from the three point correlation function, similar to the appearance of the free parameter C 0 in the two point level (5.18). Note that the three point function fully determines the coefficients of F 1 , F 3 and F 6 , and therefore the value of the free parameter C 0 from equation (5.18). Alternatively, the result in (5.22) can be written in terms of the second order variations of the scalars (φ's) defined in (E.13): It is apparent from these results that the coefficients of the 3 possible anomalies in the four derivative sector that contribute to the three point correlation function vanish. This is consistent with previous results (see [11,27]), that considered only the case where the Frobenius condition is satisfied (and therefore only the anomaly A 4 (0,4,0) = R +∇ α a α 2 ). All non-vanishing terms are proportional to α, and therefore represent only trivial terms that can be removed by adding local counterterms to the effective action. It is also important to note that these trivial terms contain contributions from φ 28 , φ 29 , φ 32 , φ 33 , φ 38 and φ 39 . These are terms that vanish in the Frobenius case (for which K A = 0, see [7] for details), and therefore their coefficients cannot be extracted from a curved spacetime coupling that assumes the existence of a foliation structure. We conclude that for α = 0, violating the Frobenius condition is essential for describing the obtained trivial terms in curved spacetime, and for constructing the appropriate counterterms to cancel them.
Finally, we have again verified that the pole residues of the various two and three point functions satisfy the conservation Ward identity: ρα af e (Q, P, K) = 0, (5.24) which is consistent with the argument made in subsection 3.2 that the conservation Ward identity (2.46) holds separately on the pole part and on the regular finite part of the correlation functions (so that TPD invariance is not anomalous).

Summary and Outlook
In this work we developed a general scheme for field theory calculations of Lifshitz scale anomalies. We analyzed the general structure of correlation functions of the stress-energy tensor in Lifshitz field theories and constructed the corresponding anomalous Ward identities. We presented a subtle ambiguity in the definition of the anomaly coefficients and clarified it. Our framework for calculating the anomaly coefficients was based on a split dimensional regularization where space (momentum) and time (frequency) integrals are regulated separately, and pole residue calculations that allowed to extract the anomaly coefficients without a full calculation of the correlation functions. In order to implement the calculational scheme we had to analyze the coupling of a nonrelativistic, Lifshitz invariant field theory in d t time dimensions and d s space dimensions to a curved spacetime manifold. This generalized the curved spacetime structure introduced in [6,7] to the case of multiple time directions.
We considered as a particular example the z = 2 free scalar field theory in 2 + 1 spacetime dimensions. We showed that the only non-zero anomaly coefficient (of those appearing in the three point functions) is in the two derivatives sector, which agrees with the heat kernel calculation in [11]. In order to account for some of the trivial terms arising from this calculation, we found it necessary to give up the Frobenius condition (and the corresponding foliation structure) of the curved spacetime description of the theory. This is because these terms are not in the span of possible anomalous contributions one obtains when assuming the Frobenius condition.
There are many directions for further studies that follow from our analysis. It would be interesting to use the general scheme developed here to calculate the anomaly coefficients of other Lifshitz field theories. One can also generalize the discussion and consider nonrelativistic field theories that exhibit non-relativistic boost invariance. In such cases one has in addition to the B-type scale anomalies also A-type ones [7] (see also [33]). In the relativistic CFT case, scale anomaly coefficients multiply universal terms in entanglement entropies. It would be of interest to analyze the entanglement entropy structure in Lifshitz field theories and the role of the scale anomaly coefficients. In the CFT case, A-type anomaly charges exhibit RG properties, i.e. a decrease from the UV to the IR. The nonrelativistic versions of these are still lacking. Finally, it would be interesting to ask whether there are experimental observables of the anomaly charges in non-relativistic systems such as low energy condensed matter ones. One potential path to consider is the hydrodymanics of such systems [34,35] and the role of scale anomalies in such descriptions.

Acknowledgments
We would like to thank Itamar Hason, Carlos Hoyos, Zohar Komargodski, Adam Schwimmer and Stefan Theisen for valuable discussions and comments. Many of the calculations in this paper were performed using xAct [36] and xTras [37], tensor computer algebra packages for Mathematica. This work is supported in part by the I-CORE program of Planning and Budgeting Committee (grant number 1937/12), the US-Israel Binational Science Foundation, GIF and the ISF Center of Excellence. A.R.M gratefully acknowledges the support of the Adams Fellowship Program of the Israel Academy of Sciences and Humanities. I.A is thankful for the support of the Alexander Zaks fellowship program.

A Notations and Conventions
In this appendix we describe the notations and conventions used in this paper.

A.1 Notations and Conventions in Curved Spacetime
Throughout this paper we use Greek letters (µ, ν, ρ, . . .) to denote spacetime indexes, both in the relativistic and non-relativistic cases. For the relativistic case, we use a spacetime metric g µν with signature diag(1, −1, −1, . . .) and denote the flat space metric by η µν . For the non-relativistic case, we use a signature of diag (−1, −1, . . . , 1, 1, . . .) (with negative signs for the time dimensions and positive for the space dimensions), and denote the flat spacetime metric by δ µν .
In both cases we use the standard torsionless Levi-Civita connection associated with the spacetime metric g µν . That is, the covariant derivative of a vector A µ is given by: where the Christoffel symbols are given by: We use the following convention for the Riemann tensor: while the Ricci tensor and scalar are given by: When using the vielbein formalism, the vielbeins e a µ are defined in the relativistic case such that: and in the non-relativistic case: We use the following notations for the determinants of the metric and the vielbeins: We use the following formulas for variations of the metric: δg σρ δe c α = −g σλ g τ ρ δg λτ δe c α , (A.10) δg ρν δe c α = η ac (δ α ρ e a ν + δ α ν e a ρ ) , (A.11) Finally, for reference we give here the expressions for the Weyl tensor squared and the Euler density of the background manifold in the (3 + 1)-dimensional case: (A.14)

A.2 Notations and Conventions for the Relativistic Scalar Case
As mentioned in appendix A.1, we use Greek letters (µ, ν, ρ, . . .) to denote spacetime indexes. In the relativistic case, we use a flat spacetime metric of the form: η µν = diag(1, −1, −1, . . .). We use the following conventions for the Fourier transforms of two and three point correlation functions in the relativistic case: where the lower-case letters p, k, q denote spacetime momenta. We also use the following notations for the Fourier transformed two and three point correlation functions of the stress-energy tensor and the variations of the action in flat space: .
Finally, we use the following notation for the Fourier transformed second variation of the expression X with respect to the background metric, evaluated in flat spacetime:
We use capital letters P, K, Q to denote spacetime momenta. The notationsP andP then refer to the spatial and temporal projections of the momentum P , respectively. We use the following conventions for the Fourier transforms of two and three point correlation functions in the non-relativistic case: FT I (2) (x = (x,x) , y = (ŷ,ȳ)) ≡ d dt+ds x d dt+ds y I (2)  We also use the following notations for the Fourier transformed two and three point correlation functions of the stress-energy tensor and the variations of the action with respect to the vielbeins in flat space: . The subscript refers to the relevant symmetry: D corresponds to TPD symmetry, whereas W corresponds to anisotropic Weyl symmetry. The superscript (n) refers to the number of points in the correlation function. The notation for the Fourier transforms of these Ward identities is given in (5.12)-(5.15).

B Non-Relativistic Curved Spacetime with Multiple Time Directions
In this appendix we discuss the coupling of a non-relativistic, Lifshitz invariant field theory in d t time dimensions and d s space dimensions to a curved spacetime manifold. We first generalize the curved spacetime structure introduced in [6,7] to the case of multiple time directions. We then discuss the local symmetries of the theory over curved spacetime. Finally we construct the curved spacetime action that corresponds to the free z = 2 Lifshitz scalar.

B.1 Curved Spacetime Definitions
Consider a non-relativistic field theory defined over a spacetime manifold with d t time dimensions and d s space dimensions, such that it is invariant both under time rotations and space rotations. In order to define the theory over a curved spacetime manifold, we generalize the structure introduced in [6,7]. We require the background manifold to be equipped with a metric g µν , or alternatively vielbeins e a µ , as well as a distribution of dimension d s , corresponding to the space directions at each point of the manifold.
This distribution can be represented by the cotangent subbundle of 1-forms that annihilate space tangent vectors. Suppose this subbundle is spanned by a basis of d t linearly independent 1-forms t Then physical quantities, such as the curved spacetime action of the theory, will depend on the subbundle spanned by {t (i) µ }, but not on the choice of basis. We therefore expect them to be invariant under transformations of the form: t where L ij is an invertible matrix that depends on the spacetime coordinate x. Alternatively, we can choose an orthonormal basis of 1-forms n (i) µ , satisfying: and require invariance under local rotations of the time directions -transformations of the form n (i) where Λ ij is an orthogonal matrix that depends on the coordinate x.
Using the set of 1-forms n (i) µ we can make several definitions on the background manifold. A tensorT αβ... is space tangent if it satisfies: for any 1 ≤ i ≤ d t . Any tensor can be rendered space tangent by projecting it on the space directions using the space projector (or spatial metric), defined as: 31 The covariant derivative of the 1-form n (i) µ can be decomposed as follows: αβ is symmetric, and given by: where the space projected Lie derivativeL n (i) of a space tangent tensorT αβ... is defined as follows:L We denote its trace by K (i) αβ is antisymmetric and given by: and we define the generalized acceleration vector a β as trace over the time indexes: a β ≡ a (ii) and from (B.1) it is easy to show that b α . One can also obtain the following expression for the space projected Lie derivative of the 1-form n (j) µ in the direction of n µ (i) : so that the generalized acceleration vector is also given by: Finally, c (ijk) is a scalar given by: that satisfies c (kji) = −c (ijk) (again from (B.1)).
Next we define the space tangent covariant derivative of a space tangent tensorT αβ... as follows:∇ Note that the spatial metric P αβ is covariantly constant under this derivative: Operating with the commutation of two space tangent derivatives on a space tangent tensor T αβ... , one obtains the following expression: where R αρµν is a space tangent tensor defined by: αβ is the total space tangent component of ∇ α n (i) β and R αρµν is the standard Riemann curvature associated with the covariant derivative ∇ µ (see appendix A for our conventions). Similarly to the one time direction case (see [7]), the tensor R αρµν does not have all of the standard symmetries of the Riemann tensor. It is therefore useful to define a modified Riemann tensor: which satisfies the usual Riemann tensor symmetries except for the second Bianchi identity. We then define the equivalents of the Ricci tensor and scalar for this modified Riemann tensor R αρµν as follows: Note that from the above definitions, one gets the following identity for the divergence of a space tangent vectorV µ :

B.2 Symmetries Over Curved Spacetime
Next we turn to discuss the symmetries of the curved spacetime field theory. Like in the one time direction case, the symmetries of the flat space Lifshitz theory translate to local symmetries over curved spacetime: First, as mentioned in subsection 2.2.1 for the one time direction case, we require timedirection preserving diffeomorphism (TPD) invariance that corresponds to space rotation symmetry in flat space. In this case, these are diffeomorphisms with a parameter ξ that satisfies: Similarly to the one time direction case (see [6,7]), we can extend these to the full diffeomorphism group by having the 1-forms t (i) µ transform appropriately: Second, as previously mentioned, we require invariance under local time rotations of the form n (i) where Λ ij is a spacetime dependent orthogonal matrix. In infinitesimal form, these transformations are given by: where ω ij is a transformation parameter that satisfies ω ji = −ω ij . From this transformation and the definition of the derivativesL n (i) ,∇ µ we obtain the following for a space tangent tensorT αβ... :

(B.22)
Using these formulas and the various definitions from appendix B.1, the following time rotation transformation rules can be derived: Finally we require anisotropic Weyl invariance, which is the local version of Lifshitz scale invariance in flat space (5.5). In the case of multiple time directions, the infinitesimal anisotropic Weyl transformation is given by: From this transformation and the definition of the derivativesL n (i) ,∇ µ we get the following formulas for a space tangent tensorT αβ... : where I[T ] is the rank of the tensorT αβ... . Using these formulas and the various definitions from appendix B.1, the following anisotropic Weyl transformation rules can be derived: (B.26)

B.3 Action of the Free z = 2 Scalar
Consider the free Lifshitz scalar in d s + d t dimensions with a dynamical critical exponent of z = 2. Its flat space action is given in (5.4). In order to couple it to a curved spacetime manifold, we have to define its curved spacetime action S(g µν , t (i) α , φ) such that it is invariant under TPD, local time rotation and anisotropic Weyl transformations, and it reduces to the action (5.4) in the flat space limit. We also require it to be regular in the physical limit d t → 1, for the sake of using the split dimensional regularization procedure as described in subsection 3.2.
Suppose that φ has a scaling dimension s under anisotropic Weyl transformations, that is: For the temporal part of the action, using the transformations given in appendix B.2 we can find a linear combination of L n (i) φ and K For the spatial part of the action, we can find a linear combination of∇ 2 φ, a µ∇ µ φ, a 2 φ and∇ µ a µ φ which is anisotropic-Weyl-covariant, given by: 29) or alternatively, we can use∇ 2 φ, a µ∇ µ φ, a 2 φ and Rφ to obtain the following anisotropic-Weyl-covariant linear combination: 32 Finally, for the action to be anisotropic-Weyl-invariant with z = 2, the dimension of the scalar field is required to satisfy: Combining these expressions, we arrive at two possible options for the curved spacetime action of the free z = 2 scalar (corresponding to the two options for the spatial part). The first option is given by: where: and d lif ≡ 2d t + d s . The second option is given by: (B.34) 32 Of course, any linear combination of these two options could also be used. where: (B.35) These actions are indeed invariant under TPDs, local time rotations and anisotropic Weyl transformations. They are also regular in the d t → 1 limit, as required.
In this work we have chosen to use the first option. Note, however, that for the purpose of the Lifshitz anomalies calculation in 2 + 1 dimensions (as done in section 5), the results would be the same for both options, since the difference between them is proportional to ε lif and therefore does not contribute to the ε lif pole residues of the flat space correlation functions (see subsection 3.2 for details).

C Relativistic Scalar Field -Feynman Rules, Vertexes and Integrals
In this appendix we give some details for the calculations of the two and three point correlation functions which are required for computing the conformal anomaly coefficients of the relativistic scalar field in two and four spacetime dimensions, as explained in section 4.

C.1 Feynman Rules and Diagrams
The propagator of the relativistic scalar field is given by: where, as explained in subsection 3.1, m is an IR mass regulator later taken to be zero. The (Fourier transformed) two and three point Feynman correlation functions of the stressenergy tensor are given by the expressions: and These expressions correspond to Feynman diagrams of the form given in figure 1 for the two point function of the stress-energy tensor, and figure 2 for the three point function of the stress-energy tensor. The vertexes V µν T are given by: where A T , B T , C T and D T are defined by: Note that terms which carry coefficients of order O(ε) do not contribute to the pole residues of the correlation functions, and can therefore be ignored for the purpose of our calculations. For example, the terms C T and D T in equation (C.4) can be neglected when calculating the pole residues around two spacetime dimensions. However, these terms cannot be neglected in four dimensions since their coefficients are no longer proportional to ε.
The expression for the correlation function (A.19) (which corresponds to a Feynman diagram of the form given in figure 3) is the following: where the vertex V µνρσ is defined by: where:

D Lifshitz z = 2 Scalar Field -Feynman Rules, Vertexes and Integrals
In this appendix we give some details for the calculation of the two and three point correlation functions which are required for computing the Lifshitz anomaly coefficients of a z = 2 free scalar field in 2 + 1 dimensions, as explained in section 5. These include the expressions for the Feynman diagrams, the Feynman rules for the propagator and all the vertexes needed.
The propagator of the Lifshitz z = 2 scalar is given by: We denote the external momentum of the two point Feynman diagram by P µ = P µ ,P µ , and the "running" loop momentum by Q µ = Q µ ,Q µ . The expression for the two point function Feynman diagram, as illustrated in figure 1, is given by: where (as in the relativistic calculation) m is an IR regulator later taken to be zero.   It is straightforward to calculate the expression corresponding to the vertex V µν from the stress-energy tensor given in (5.11). The result is: The expression (D.3) corresponds to the vertexes in figures 1 and 2, and to the right vertex in figure 3.
The expression for the Feynman diagram corresponding to the three point function of the stress-energy tensor, as illustrated in figure 2, is given by: ] .

(D.4)
The expression for the Feynman diagram corresponding to the correlation function (A.25), as illustrated in figure 3, is given by: , where the expression δ 2 S αβγδ (P, K, Q, L) ≡ δ αα δ γγ δ b β δ c δ δ 2 S α γ bc (P, K, P + K − Q, Q) corresponds to the insertion of the second variation of the action δ 2 S δe b αδe c γ in flat space. P and K in this expression refer to the external momenta, whereas the two momenta Q and L correspond to the two scalar propagator lines attached to the vertex, as illustrated in figure 3.
The expression for the vertex δ 2 S αβγδ (P, K, Q, L) is quite long, and the remainder of this appendix is dedicated to the details of its calculation. For convenience, we split it into several parts as follows. Ignoring the terms of order O(ε lif ), the action (5.2) takes the form: Let us define: and so that: The contribution of the first part of the action (D.7) to the vertex δ 2 S αβγδ (P, K, Q, L) can be easily derived: (D.10) For the contribution of the second part of the action (D.8), it is convenient to use the following identity (which follows from (B.19)): S II can then be written as follows: where I 1 , I 2 and I 3 are defined by: and (D. 15) The second variation of the expression S II with respect to the vielbeins, evaluated in flat space, can be decomposed as follows: , where all variations with respect to the vielbeins in this expression are evaluated at flat space. Each instance of the first order variations of I 3 that appear in the second, third and fourth terms in (D. 16) contributes an expression of the following form to the vertex: where the momentum P here refers to one of the external momenta associated with the vertex, and Q to the momentum of one of the scalar propagator lines attached to it. The contribution of the expression d dt+ds x∂ µ∂ µ φ(x) δ 2 I 1 (x) δe b α(y)δe c γ (z) (contained in the last term of (D.16)) to the vertex is given by: Finally, the contribution of the expression d dt+ds x∂ µ∂ µ φ(x) δ 2 I 2 (x) δe b α(y)δe c γ (z) (also contained in the last term of (D.16)) to the vertex is given by: Using the expressions (D.17)-(D. 19) and the standard relations: δe (x) δe a µ (y) flat space = δ a µ δ (x − y) , (D.20) the total contribution of S II to the vertex δ 2 S αβγδ (P, K, Q, L) can be assembled as follows: 2 αβγδ (P, K, Q, L) .

E Variations of the Anomaly and Trivial Term Densities
In [7], the possible forms of Lifshitz scale anomalies (and the cohomologically trivial terms) were derived for the case of 2 + 1 dimensions with a dynamical critical exponent of z = 2, using a cohomological formulation of the Wess-Zumino consistency condition in curved spacetime. In this appendix we detail the first and second order variations of these anomaly and trivial term densities 33 with respect to the background vielbeins, evaluated in flat spacetime. These variations are used in conjunction with the anomalous Ward identities for the flat space correlation functions given in (2.49)-(2.50) to extract the anomaly and trivial term coefficients in section 5. The appendix is divided into two parts, for the two and four derivatives sectors respectively. Throughout this appendix we use the following notation for the Fourier transformed variation of some scalar X with respect to the vielbeins: where the variations are evaluated in flat space, and the Fourier transforms use the conventions (A.21)-(A. 22). Note that the expressions in this appendix are all in d t = 1 time dimension, and therefore the time components of the momenta are given by: whereP ≡ P µ n µ andK ≡ K µ n µ , as explained in appendix A.3.