Spontaneous breaking of non-relativistic scale symmetry

We analyze the mechanism of spontaneous symmetry breaking of scale invariance in Galilean invariant field theories. We show that the existence of a dynamic gapless dilaton mode depends on whether the U(1) particle number or the Galilean boost symmetry are spontaneously broken. When both scale and particle number symmetries are spontaneously broken there is one propagating gapless Nambu-Goldstone mode. Its dispersion relation is linear if the chemical potential is nonzero and quadratic otherwise. We discuss the reversibility of RG flows in such theories.


Introduction
Spontaneous symmetry breaking (SSB) underlies a large number of physical phenomena such as superconductivity, superfluidity and the generation of elementary particle masses. A spontaneous breaking of a continuous global symmetry implies via the Nambu-Goldstone (NG) theorem the existence of a gapless Nambu-Goldstone mode. In a relativistic field theory there is one NG mode for each broken symmetry generator [1][2][3]. This one to one correspondence between broken generators and gapless NG modes does not hold when spacetime symmetries are spontaneously broken [4,5]. It is also not generally the case for spontaneously broken global symmetries in non-relativistic field theories [6][7][8].
When scale symmetry is spontaneously broken in a relativistic field theory, there is a corresponding NG mode, the dilaton. The dilaton effective action encodes, for instance, the information about the A-type conformal anomaly and has been valuable in proving the a-theorem in [9]. This leads naturally to inquire about the mechanism of spontaneously broken scale symmetry in non-relativistic field theories. Such field theories have much importance in the study of low energy condensed matter systems, as well as in non-relativistic holography (e.g. [10]).
The aim of this paper is to analyze SSB of scale invariance in Galilean invariant field theories. In non-relativistic field theories space and time scale differently: x → e σ x, t → JHEP10(2017)063 e zσ t, where σ is a real parameter and x is a d-dimensional vector. z is called a dynamical exponent and we will consider the Galilean case z = 2. We will show that there is no gapless dilaton mode in a Galilean field theory unless the U(1) particle number symmetry or the Galilean boost symmetry are also spontaneously broken. Particle number symmetry is an additional symmetry of Galilean field theories that does not exist in the relativistic case. We will see that when both scale and particle number symmetries are spontaneously broken there is just one propagating gapless mode. 1 In the Galilean algebra the generator M of the U(1) particle number symmetry is a central extension. It appears in the commutator of translations P i and Galilean boosts K j [P i , K j ] = −iδ ij M, i = 1, . . . , d . (1.1) The broken symmetry generators create NG modes from the vacuum and their commutator algebra may impose relations between these modes. Such an example is the relation among the NG mode related to SSB of the U(1) particle number symmetry and the modes related to the SSB of boosts. While commutator (1.1) predicts that U(1) SSB implies that Galilean boosts are spontaneously broken, it also predicts that there is only one independent NG mode. This relationship is an example of a general structure called Inverse Higgs Constraints (IHC) [11][12][13]. Such an argument, however, does not explain why when scale and U(1) are spontaneously broken there is only one NG mode. We will derive this result by an explicit calculation of the spectrum.
Based on z = 2 dimensional arguments we may anticipate the dispersion relation of the NG mode. If we denote by v the symmetry breaking length scale, then the dispersion relation takes the general form We'll see that this is indeed the behavior, where with zero chemical potential F (0) is finite and we get to leading order a quadratic dispersion relation while with a nonzero chemical potential we get to leading order a linear dispersion relation instead of a quadratic one. There are various methods to construct the effective action of NG bosons. One way is to write all possible terms that respect all the symmetries. Another way is to couple the theory to curved external background sources. We will consider both methods in the study of the non-relativistic dilaton. 2 The paper is organized as follows. In section 2 we will outline the setup of the problem and the various cases that we consider. In the next section we will begin by considering the NG effective action based on symmetry arguments. We will see that spontaneously breaking scale invariance while maintaining the U(1) and Galilean boost symmetries does not allow for a propagating dilaton mode. We will then spontaneously break the U(1) and boost symmetries as well. We will construct the NG effective action at leading orders in the derivatives expansion. Next, we will derive the same effective action using a coupling to the Newton-Cartan curved geometry. We will analyze the spectrum of the NG bosons and 1 The case where the particle number symmetry is not broken but the Galilean boosts are broken seems to be of less physical relevance and will be discussed in the appendix. 2 Another possible method utilizes the coset construction, however we will not pursue it in this work.

JHEP10(2017)063
find one gapless propagating mode. Finally, we will discuss the possible relevance of the results to RG flow theorems in Galilean field theories and conclude with a brief summary and outlook. In appendix A we discuss the case of unbroken particle number and broken boost symmetries.

The problem setup
We consider non-relativistic systems in d + 1 dimensions, which are invariant under the (centrally extended) Galilean group with the addition of non-relativistic (Lifshitz) scale invariance with a dynamical exponent of z = 2. The symmetry generators satisfy the following algebra: where H, P i , L ij , K i and D are the generators of time translations, space translations, space rotations, Galilean boosts and Lifshitz scaling respectively. M is a generator of the U(1) internal symmetry that corresponds to the conserved particle number, and all other commutators are zero. We will also discuss systems which are invariant under the full Schrödinger group, which contains, in addition to these symmetries, the special conformal transformations, the generator of which we denote by C. In addition to the above algebra (2.1), the generators then also satisfy: We assume either that the system is in its vacuum state that minimizes the Hamiltonian H, or alternatively that the state of the system minimizes the modified HamiltonianH ≡ H − µM , where µ is the chemical potential of the system. 3 Such Galilean systems are usually observed in nature simply as low velocity limits of relativistic systems. In these cases, the corresponding relativistic systems have a relativistic Hamiltonian H rel and relativistic chemical potential µ rel associated with them, and one defines the non-relativistic Hamiltonian and chemical potential as: and takes the limit of v c (or pc E). It is important to note that, even in such cases, the Hamiltonian H and chemical potential µ we refer to in the following sections are the non-relativistic ones, and it is the non-relativistic Hamiltonian H that satisfies the above JHEP10(2017)063 algebra. It is also worth mentioning that, while such a Galilean symmetry often appears as the small velocity limit of a full Lorentzian relativistic symmetry, it may also appear as an emergent symmetry at low energy.
In physical systems, the particle number M generally has a non-negative spectrum, where the only state with M = 0 is the empty state that contains no particles. The case of M = 0 (with zero particle number density ρ) is therefore not physically interesting for the study of SSB, as we don't expect any symmetry to be spontaneously broken in such cases. However, in the interest of mathematical completeness (and since, as we point out later, the study of the dilaton effective action is useful for analyzing other QFT properties outside of SSB), we will briefly consider this case. The more physically interesting case, however, is that for which M > 0. In this case, the Galilean boost symmetries are necessarily spontaneously broken by the particle number current In both cases, for the purposes of this work, we will assume that space translations and space rotation symmetries (as well as the modified time translation generated byH) are not broken.
We will consider systems for which the non-relativistic scale symmetry D is spontaneously broken, that is, there is some dimensionful operator O such that [D, O] ∝ O = 0. There are then two possible cases depending on whether the operator O is invariant under the particle number U(1) symmetry.
If the operator O is U(1) invariant (and assuming all other operators with non vanishing expectation values are as well), the U(1) symmetry is not spontaneously broken. Assuming non-vanishing particle number M = 0, 4 this is a case where only the scale and Galilean boost symmetries are broken. Such a system is similar to the "framids" discussed in [15], in the non-relativistic limit. One might expect the appearance of d + 1 NG modes, corresponding to the broken scale and Galilean boosts. Since it is not clear whether such a scenario indeed appears in any physical system, we relegate this case to appendix A.
If the operator O is not U(1) invariant, and [M, O] = 0, the U(1) symmetry is spontaneously broken as well. Examples of such systems may be found in superfluids, superconductors and other types of non-relativistic condensates. In these systems, since both the scale, the boosts and the U(1) symmetries are spontaneously broken, one might initially expect to find separate NG modes for each of them. However, as explained in e.g. [12], this is not the case: due to the well-known "Inverse Higgs" mechanism, the commutation relation between K i and P j in (2.1) imposes a relation between the modes v i that would otherwise correspond to the broken Galilean boosts and the mode θ corresponding to the broken U(1) symmetry, of the form: v i ∝ ∂ i θ. Therefore, in this case, the broken Galilean boosts play no role when building the possible effective action for the NG modes -one need only consider the modes corresponding to the broken scale and U(1) symmetries (while the effective action is, of course, required to be invariant under all of these symmetries, including Galilean boosts).

JHEP10(2017)063
For systems which are Schrödinger invariant, it is clear from the relations (2.2) that, in all the cases considered, the special conformal symmetry C is necessarily also broken. However, again due to the "Inverse Higgs" constraints derived from (2.2) (and similarly to the relativistic conformal case), we do not expect a corresponding independent NG mode to appear as the algebra imposes a relation between this mode and the other ones.
When considering the case of broken U(1), it is important to distinguish between two possible sub-cases, corresponding to the existence of a non-vanishing chemical potential µ. When µ = 0, the expectation of the operator O is constant in time, since in the vacuum state: However, as pointed out in [14], when µ = 0, this is not the case and one instead obtains: assuming that O has a charge of m under the U(1) symmetry. This is accounted for in our spectrum analysis in subsection 3.1.2 by perturbing the U(1) NG mode θ around −µt.
Note that generally, the condition of non-vanishing particle number M = 0 doesn't necessarily imply a non-vanishing chemical potential. As a simple example, consider a model for a mixture of two Bose condensates in 2 + 1 dimensions with contact interactions, given by the following action (this is a special case of the model considered e.g. in [16]): where m is the Galilean mass of the condensates, and g a (dimensionless) coupling constant. This model is indeed Galilean invariant, and classically it is Lifshitz scale invariant as well.
With vanishing (non-relativisitic) chemical potential, this model still classically allows for a moduli space of (inequivalent) vacua with a finite particle number density, where both φ 1 and φ 2 gain a non-vanishing expectation value such that |φ 1 | 2 = |φ 2 | 2 = 0. Each of these vacua breaks the U(1), boosts and scale symmetry spontaneously. While this is only true without considering any quantum corrections, this example serves to demonstrate the possibility of having a finite particle density configuration with vanishing chemical potential.
To summarize, in this work we consider the following cases: • The case where only scale symmetry is spontaneously broken (the non-physical M = 0 case), • The case where both scale and Galilean boost symmetries are broken but the U(1) particle number symmetry remains unbroken (considered in appendix A), • The case where scale, Galilean boosts and the U(1) symmetry are all spontaneously broken, either with vanishing or finite chemical potential.

JHEP10(2017)063
Our main focus is on the third case, which seems relevant to various condensed matter systems. For each case, we will derive the possible form of the effective action for the NG modes, and study their spectrum. Finally, it is important to note that the analysis of effective actions is also useful for the study of other QFT properties other than SSB. In relativistic theories, for example, it has been used to study quantum anomalies and RG flow properties of field theories (see [9]). As such, even the cases which are less physically relevant in the context of SSB may be helpful for drawing conclusions on other properties of non-relativistic QFTs. In section 4 we indeed comment on the possibility of applying some of our conclusions to the study of RG flows in Galilean field theories.

Symmetries
In this subsection we derive the first few terms in a derivative expansion of the effective action of the non-relativistic NG boson following a spontaneous breaking of scale invariance ( x, t) → (e σ x, e zσ t), using the symmetries of the theory.

Unbroken U(1) particle number and Galilean boost symmetries
The dilaton is a real field which we will denote by τ . It carries zero U(1) charge. Under scale transformation the dilaton transforms as where L is a general polynomial in these differential operators. Under z = 2 non-relativistic boost transformation with boost parameters u, the dilaton and its derivatives transform as The ∂ i term is invariant while the ∂ t term is not. This basically forbids any time derivative term in the effective Lagrangian, hence there is no dynamical dilaton.
Consider this more explicitly. Zeroth order derivative terms can't by themselves generate dynamics, however, they are allowed provided the scale transformation is fixed correctly as follows: where Λ is a dimensionless constant. First order derivatives terms are total derivatives, so next we consider the second order derivatives terms. By rotational symmetry, the only JHEP10(2017)063 possible two derivatives terms are (∂ t τ ) 2 and (∂ i τ ) 2 . Thus, up to two derivatives the effective dilaton Lagrangian should read The first term is not invariant under z = 2 non-relativistic boost transformation. One can repeat the same analysis at higher order in derivatives arriving to the same conclusion that there is no dynamical dilaton. Note, that if one works in Lifshitz field theory but without imposing Galilean boost invariance, one can have a dynamical dilaton. For example, the z = 2 scale invariant terms This is the Lagrangian for a free Lifshitz real scalar, which is not boost invariant.

Broken U(1) particle number symmetry
In the following both scale and U(1) particle number are spontaneously broken and we will consider a complex scalar field φ = e ∆τ +imθ that carries a U(1) charge m and scaling dimension ∆ (where τ corresponds to the dilaton and θ to the U(1) NG mode). As explained above, Galilean boosts must also be broken but the corresponding NG modes are directly related to the U(1) NG mode, and are therefore not independent modes. 5 Under z = 2 boost transformation one has or τ invariant and The zero derivatives term φφ * is invariant and gives (3.4) with ∆ = d+2 2 . At leading order in derivatives we can write a scale, rotation and boost invariant Lagrangian with ∆ = d 2 . More generally, we can build a boost invariant structure (up to a phase) −i∂ t φ− 1 2m ∂ 2 i φ. It can be used for the construction of a boost and scale invariant Lagrangian, so to next order in derivatives one has also

JHEP10(2017)063
where φ has U(1) charge m and scaling dimension ∆, and ∆ = d 2 − 1. In terms of τ and θ Of course, we may similarly construct other invariant terms of this order in derivatives (or higher) using the field φ, the covariant derivative operator −i∂ t − 1 2m ∂ 2 i and appropriate exponents of τ to compensate for the dimension. Note that in subsection 3.3 we use a term of the form (3.10) as an example of an higher order term for the purpose of analyzing the spectrum of the theory, however using other higher order terms constructed this way does not change the main results of that subsection.
We can also consider the case of theories which are invariant under the full Schrödinger group, which contains in addition the special conformal transformation, given by: 12) or in terms of τ and θ: .
It is easily verified that the previously mentioned structure −i∂ t φ − 1 2m ∂ 2 i φ is covariant under these special conformal transformations (that is, invariant up to a phase and a scale factor), only when ∆ = d 2 (or when m = 0, which corresponds to terms that depend only on τ and its spatial derivatives). Therefore the Lagrangians (3.4) (with ∆ = d+2 2 ) and (3.9) (with ∆ = d 2 ) are indeed Schrödinger invariant. However, the requirement for Schrödinger invariance restricts the higher order Lagrangian (3.10) to the case of ∆ = d 2 (or alternatively m = 0).

Geometry
One can construct the NG boson effective action by geometrical considerations: coupling the field theory to a curved background, promoting the symmetries to local ones, looking at all possible actions in this framework and taking the flat background limit. This was done for the relativistic case [17] where the curved background is a Riemannian geometry with a metric tensor G µν . One constructs a Weyl invariant metricĜ µν = e −2τ G µν and writes the effective action in terms of scalar terms constructed from it.
In the non-relativisitic case one has to use the Newton-Cartan (NC) geometry instead (see e.g. [18,20]). The NC geometry is built from a time direction described by a 1-form n µ , a spatial metric h µν orthogonal to n µ and a U(1) gauge field A µ which couples to the conserved particle number current. Further, one defines a vector v µ that satisfies v µ n µ = 1 and induces a metric h µν satisfying (3.14)

JHEP10(2017)063
These definitions are not unique because we can redefine v µ using an arbitrary vector ψ ν v µ → v µ + h µν ψ ν , (3.15) and redefine h µν correspondingly: together with the following redefinition of A µ : These transformations are the Milne boosts and one can define Milne boost invariant objects, as follows: Since we are interested in the (z = 2 anisotropic-)Weyl (scale) and U(1) symmetries, we need the Weyl and U(1) transformations of the NC objects. The U(1) gauge symmetry transforms only A µ out of the basic structures: while the Weyl symmetry transforms n µ and h µν : There are two equivalent ways to construct the NG boson effective action for the scale and U(1) spontaneous breaking. One is to introduce spectator fields, τ for the scale symmetry and θ for the U(1) symmetry, with τ transforming under scale like a dilaton τ → τ + σ, and θ transforming under U(1) as θ → θ + Λ where A µ → A µ + ∂ µ Λ (and both invariant under Milne boosts). Using these spectator fields and the NC structures h, n, A, we can then define Weyl and gauge invariant geometric quantitiesĥ,n,Â as follows: From these structures, we can build boost invariant scalars. This procedure is made straightforward by defining the boost invariant metricĝ and vectorv µ A as follows: Finally, we take the limit where the geometry is flat: and thus remain only with the spectator fields. This proposal is a generalization of the relativistic case discussed above, where the dilaton factor was added to the metric G µν JHEP10(2017)063 to compensate for the Weyl variation. In our case we have the boost invariant structures g µν and v µ A given in (3.18), as well as n µ and h µν . Here again, to enforce Weyl invariance we add τ , and to enforce U(1) invariance we add θ. The structuresĝ µν ,v µ A ,n µ ,ĥ µν are therefore Milne boost invariant (by the construction of g µν , v µ A , n µ , h µν ), Weyl invariant (by the use of the compensator τ ), and U(1) invariant (by the use of the compensator θ).
The other equivalent way to construct the effective action, which is the one we will pursue, is to take Milne boost invariants and perform Weyl and U(1) transformations with parameters τ and θ, respectively. This will evidently give the same answer (up to the sign of the fields τ and θ, which will be opposite to the one in the first method, and therefore to the convention used in [9]), since the Weyl and U(1) transformations will force the appearance of τ and θ in exactly the right form such that if they themselves transform under Weyl and U(1), the whole expression would be invariant.
Using the boost invariant scalars, we will get the effective action for τ and θ after restricting to flat geometry (3.23). Note, that the flat geometry restriction is respected by a combination of a Milne boost and a U(1) transformation. Under this combination, θ transforms in the same way as it transformed under non-relativistic boost transformations in the previous section (3.8). The reason is that under the appropriate Milne boost transformation parametrized by ψ µ = (0, u) at flat geometry, A µ transforms as in (3.17), which in this case takes the form In order to compensate for this transformation, the U(1) transformation should be (again, at flat geometry) Λ = 1 2 u 2 t − u · x, which produces the correct boost transformation of θ. Note also that the flat geometry restriction is respected by an additional combination of Milne boosts, U(1) transformations and anisotropic Weyl transformations, which corresponds to the special conformal transformation (3.13) (see [18]). Therefore the effective action obtained using this geometric method necessarily corresponds to the full Schrödinger invariant case.
We wish to list the invariant scalars to leading order in derivatives. Note that derivative counting should be done after restricting to flat geometry. Consider the Milne boost invariant metric and vector (3.18). The simplest geometrical term is the cosmological constant term (constant up to the det(γ µν ) factor, where γ µν ≡ h µν + n µ n ν , that contributes the τ dependence and ensures Weyl invariance), which matches the non-dynamical term e (z+d)τ (3.4) discussed previously.
The next simplest scalar one can build is g µν v µ A v ν A . We perform Weyl transformation and U(1) transformation and then restrict to flat geometry to get an expression in terms of τ and θ. Counting derivatives naively in this expression, g may contribute one derivative and v A may contribute one space derivative, so we might be lead to think that in total we have 3 derivatives. However, since h µν n µ = 0, the full expression, when restricted to flat geometry, has only one time derivative or 2 space derivatives. Written explicitly, we get: The other term at this order in derivatives is the spatial Ricci scalarR corresponding to the standard Levi-Civita connection of the metric induced on the space foliation. 6 Since JHEP10(2017)063 the metric induced on the foliation is invariant under Milne boosts (i.e. g µν u µ w ν is boost invariant for any space tangent vectors u µ , w µ ),R is boost invariant as well. It is also gauge invariant, and therefore it depends only on τ . By Weyl transformation we get There are more terms, e.g. a µ a µ where a µ ≡ −L v n µ (see e.g. [20]), but when we restrict to conformally flat background and consider integration by parts, they are all equivalent to the terms above. Therefore, these two expressions complete the list of terms up to one time derivative or two space derivatives (i.e. the same order in z = 2 Lifshitz scaling counting). From the last two expressions, (3.24) and (3.25), we construct the leading order Lagrangian for the complex field φ (3.9). For higher orders in derivatives, one can define a boost invariant affine connection from the structures n µ , h µν , v µ A , g µν , as well as the corresponding Riemann tensor R λ µσν and Ricci tensor R µν ≡ R σ µσν (see [18]). One can then obtain various higher derivative boost invariant scalars, such as R µν h µν and R µν v µ A v ν A , as well as purely spatial ones such asR 2 and a 4 . As before, by performing Weyl and gauge transformations on a linear combination of these scalars and restricting to flat geometry, one can construct the higher order Lagrangian (3.10) (as well as other higher derivatives terms, which do not affect the conclusions of the analysis in the next subsection).

Spectrum analysis
In the following we will analyze the spectrum of the low energy theory of τ and θ. We use the notation φ ∆,m = e ∆τ +imθ for a complex scalar field of dimension ∆ and mass m, and φ k as a shorthand for φ ∆ k ,m k . The leading order boost, U(1) and scale invariant Lagrangian which we derived in the previous sections reads where the dimensions ∆ i are fixed by scale invariance to the following values: Note that this is not the most general expression one could build up to this order in derivatives -we could take a linear combination of similar terms using fields φ ∆,m with various values of ∆ and m, while compensating for the dimension by multiplying each by an appropriate exponent of τ (as in (3.10)). However, for the cosmological constant term (the Λ term) and the leading term in derivatives (the A term), one can always rewrite these terms in the form given in (3.26) (using just φ 0 and φ 1 with no extra exponents of safely assume that nµ satisfies the Frobenius condition and therefore induces a foliation of the spacetime manifold into equal time slices. See [20] for further discussion.

JHEP10(2017)063
τ ), with an appropriate choice of the parameters Λ, A and m 1 . For the subleading term in derivatives (the B term), this is not true in general, but since we will be mainly interested in the leading contributions of the Λ and A terms, we will assume the form given in (3.26) as an example for the contribution of subleading terms in derivatives. We therefore use the Lagrangian (3.26), for which all τ dependence is through the φ 0 , φ 1 and φ 2 fields. It is also important to note that the subleading term in (3.26) is not Schrödinger invariant, as mentioned in subsection 3.1.2. Full Schrödinger invariance restricts this subleading term to the form (3.10) with ∆ 2 = d 2 . However, using such a term instead in (3.26) does not change the main results of this subsection.
We derive the equations of motion by varying the action defined by (3.26) with respect to τ and θ, through φ 0 , φ 1 and φ 2 . We have δS = δS δφ i δφ i + c.c., where: and their complex conjugates.

(3.29)
Separating the variation by τ and θ we have where and denote the real and imaginary parts respectively.
We will perturb the equations of motion around fixed values whereμ is the chemical potential. 7 The particular value ofτ can be absorbed by a redefinition of A,B and Λ, so we will take it to be zero. Note however that v ≡ e −τ quantifies the scale of the state we are perturbing about.

JHEP10(2017)063
We express δS δφ 1 φ 1 and δS δφ 2 φ 2 using τ and θ (see (3.28)): where here denotes keeping only terms that contribute up to linear order in the τ and θ perturbation. Similarly, again up to linear order, we have (3.33) Combining the last two equations, the variation of the action with respect to τ and θ is given by the following expressions up to linear order in the perturbation: (3.34) The zeroth order equation requires and we see that the static potential value (cosmological constant) Λ corresponds to a nonzero chemical potentialμ (although a non-zero chemical potential is possible even when Λ vanishes). Consider first the case of a non-zero chemical potentialμ. We get the following linearized equations of motion to leading order in momentum k and energy ω,

JHEP10(2017)063
where δτ denotes the Fourier transform of the τ perturbation around 0, and δθ denotes the Fourier transform of the θ perturbation around −μt. From these equations we obtain the following dispersion relation: We can see that in this case, in the limit k → 0, the dispersion relation is linear ω ∼ k and δτ ∼ ωδθ. This is consistent with the analysis in [19]. Note, that since δτ ∼ ωδθ, the perturbation is mainly in θ, which may also justify ignoring the τ contribution in the leading order superfluid effective field theory as was implicitly done in [19]. Also note that the stability of the modes in the k → 0 limit requiresμ > 0.
The result (3.37) implies the speed of sound v s ≡ 1 ∆ 1μ = 2 dμ . This can be easily understood from dimensional analysis considerations: in the presence of a z = 2 Lifshitz scale invariance, we can expect the relation between the chemical potentialμ and the conserved particle number (or mass) density ρ to be of the formμ = Cρ 2/d , where C is some dimensionless parameter. From standard thermodynamic relations, the speed of sound will be given by: Consider next the case Λ =μ = 0. The linearized equations are given by: When A = 0 we have at leading order for small values of k (and therefore ω) the nonrelativistic dispersion relation Thus, we find one gapless mode at large length scales compared to the breaking scales. The corrections can be computed to give ω 2 = αk 4 + βk 6 + . . . , where the coefficients α, β, . . . are determined from expanding the determinant of (3.39) to growing orders in k 2 . Note that the limit m 1 → 0 takes us from the broken U(1) to the unbroken U(1) case. The dispersion relation (3.40) blows up and we are left with no propagating mode. In addition to the gapless mode, we have also a gapped mode as can be seen by setting k = 0 in (3.39) and we get Finally, note that we can obtain the case of U(1) SSB without scale invariance from the above equation (3.34). Take τ to be a constant rather than a dynamical field, and to first order the variation with respect to θ gives

JHEP10(2017)063 4 On a non-Relativistic a-theorem and the Frobenius condition
The relativistic dilaton effective action was valuable for the proof of the a-theorem [9] in 3 + 1 dimensions, i.e. the coefficient of the A-type conformal anomaly a satisfies a IR < a UV . In the following we will make a few comments on the non-relativistic Galilean case and the feasibility of using similar arguments to prove an RG flow theorem in case such a theorem indeed holds.
In [9], the RG flow of a generic relativistic field theory from a UV to an IR fixed point was studied by weakly coupling the theory to a dilaton as a conformal compensator, and matching the conformal anomalies between the UV and IR theories. In particular, the A-type anomaly of the theory contributes to the effective action of the dilaton in the IR, and therefore the S-matrix of dilaton scattering. The a-theorem then follows from unitarity of the theory.
In the case of non-relativistic field theories, invariance under a Lifshitz scale symmetry implies the following Ward identity for the stress-energy tensor: which is just a generalized version of the conformal tracelessness condition. However, similarly to the relativistic case, the scale symmetry can be violated due to quantum anomalies (analogous to the conformal anomalies) [20][21][22][23][24][25][26]. The expectation value of the stress-energy tensor on a curved spacetime manifold then no longer satisfies identity (4.1). It instead acquires an anomalous contribution: 8 where A is a local functional of the backgrounds fields, and the infinitesimal (anisotropic) Weyl transformation of the effective action is given by: For Galilean invariant theories in d + 1 dimensions, it has been suggested in [18] that these Lifshitz anomalies correspond to conformal anomalies of relativistic field theories in d + 2 dimensions defined on a manifold with a null isometry, via a null reduction procedure. In particular, this suggests a possible A-type anomaly in these Galilean theories, which corresponds to the Euler density anomaly term of the d + 2 dimensional relativistic theory. In [20], this possibility was confirmed for a 2 + 1 dimensional Galilean theory via a cohomological analysis of the Wess-Zumino consistency condition (an explicit expression for this A-type anomaly is given in equation (5.13) of [20]). However, it was observed that this A-type anomaly exists only when one assumes the Frobenius condition is violated by the curved spacetime NC structure, i.e. when the 1-form n µ does not satisfy:

JHEP10(2017)063
and therefore does not induce a foliation of the spacetime manifold into equal-time slices. When such a foliation structure exists, this A-type anomaly term becomes cohomologically trivial, and can be removed by adding an appropriate local counter-term to the effective action, of the form: where a µ ≡ −L v n µ is the acceleration associated with n µ (see equation (5.14) in [20]). The existence of the A-type scale anomaly in the Galilean case suggests it may be possible to follow a similar argument to the one given in [9] to prove an a-theorem for Galilean field theories. Consider a theory which is invariant under the Galilean group, that flows between a UV and an IR z = 2 Lifshitz fixed points. This theory can be coupled to a non-relativistic dilaton field τ by multiplying any dimensionful parameter by an appropriate exponent of τ compensating for its dimension, thereby rendering the theory scale invariant. We can also add a kinetic term for the dilaton, however as we have seen in previous sections in order to have a boost invariant, dynamic term one has to involve the U(1) particle number Goldstone mode θ. 9 We can then choose this kinetic term to be of the form (3.9) 10 (with an arbitrary value of m). The coupling of the matter to the dilaton can be made arbitrarily weak by taking the coefficient of this term A to be much larger than all other dimensionful parameters.
Similarly to the way equation (3.2) in [9] was derived, we can write an IR effective action for the non-relativistic Galilean theory coupled to the dilaton. As in the relativistic case, the IR effective action will have a contribution NR IR from the non-relativistic Lifshitz invariant Galilean field theory (that replaces the relativistic CFT) in the infrared. It will have a contribution from the invariant terms, L dilaton , corresponding to the local effective action of the dilaton as discussed in the previous sections, which is of the general form (3.26) (with possibly more terms of the 2 time or 4 space derivative order or higher). And finally, it will have a contribution from the A-type anomaly of the theory. This contribution can be calculated by replacing the Weyl parameter σ in (4.3) by τ , substituting into A the expression for the A-type anomaly (given in [20] for 2 + 1 dimensions) evaluated on a NC background which is Weyl transformed 11 with τ as the parameter, solving the resulting equation and restricting to flat space (see e.g. [17]). Alternatively, it can be obtained from the conformal A-type anomaly contribution in d + 2 dimensions via a null reduction. The result for 2 + 1 dimensions is similar to the 3 + 1 dimensional relativistic case, and given by: where a UV and a IR are the coefficients of the non-relativistic A-type anomaly in the UV and the IR theories respectively (including the dilaton contribution). Note that the A-type 9 Alternatively one can use the Galilean boosts modes vi (as considered in appendix A). We will focus here on the U(1) mode, but the other option isn't significantly different. 10 Alternatively we could discuss an RG flow triggered by some operator acquiring a VEV that spontaneously breaks both scale invariance and the U(1) particle number symmetry, leading to a dynamic NGB effective action of the form discussed in previous sections. 11 Since we are assuming the U(1) symmetry is not anomalous, the anomalous contribution is gauge invariant, and so will not depend on θ if a gauge transformation is performed with θ as the parameter.

JHEP10(2017)063
anomaly contribution has no time derivatives. This is a consequence of the fact that it is a U(1) singlet.
Similarly to the relativistic case, the dilaton action L dilaton has additional higher order scale invariant terms, whose couplings are non-universal and cannot be fixed. However, there is an important difference: in the relativistic case, the 4 derivatives terms in L dilaton are distinguishable from the contribution of the A-type anomaly, that is the anomaly contribution term cannot be reproduced by a linear combination of the allowed invariant terms in L dilaton . This is not true for the Galilean case in 2 + 1 dimensions, since the anomaly contribution term in (4.6) can also be obtained from a local contribution of the form (4.5) to L dilaton . This observation can be understood as a consequence of the Frobenius condition: in order to obtain the dilaton effective action in flat space, one naturally works with a conformally flat background, which necessarily satisfies the Frobenius condition (4.4) and therefore has a foliation structure. On such a background, as discussed in [20], there is no A-type anomaly, as it becomes cohomologically trivial. In order to be able to distinguish the anomaly contribution from non-universal contributions to L dilaton , one would have to instead look at a field theory defined on a curved NC background that violates the Frobenius condition.
It is also important to note here the role played by the special conformal transformation. In the relativistic case, the contribution of the A-type anomaly to the dilaton effective action is invariant under global scaling transformations. It is not, however, invariant under special conformal transformations, and it is this property that prevents it from being included in L dilaton from the point of view of the flat space theory. In the nonrelativistic case, imposing full Schrödinger invariance does not have the same consequence, as the (purely spatial) contribution of the A-type anomaly is in fact invariant under the Schrödinger special conformal transformations.
Another difference compared to the relativistic case is that these higher order terms in L dilaton may contribute to the dilaton 2 to 2 scattering, in contrast to the relativistic CFT case where the higher order terms didn't contribute at leading order. There are two other notable differences that have already been mentioned: the first is that, while the non-relativistic anomaly term is a U(1) singlet (4.6), the dilaton effective action L dilaton in the non-relativistic case involves the U(1) particle number Goldstone mode. Second, in the relativistic CFT case where the RG flow is triggered by conformal SSB, a cosmological constant leads to a gapped mode in contradiction to the Goldstone theorem, and is not allowed in the SSB effective action. On the other hand, in the non-relativistic case as we showed above, a cosmological constant is allowed and simply leads to a chemical potential.
These differences, and especially the first one, namely the fact that in a conformally flat background the A-type anomaly contribution is indistinguishable from that of a trivial term, seem to suggest that it is not straightforward to generalize the proof of the a-theorem to the Galilean case.

Summary and outlook
We studied the mechanism of spontaneous symmetry breaking of scale invariance in nonrelativistic field theories that possess Galilean boost invariance. We showed that there is a JHEP10(2017)063 dynamical gapless mode only if the U(1) particle number symmetry or the Galilean boost symmetries are spontaneously broken too. The dispersion relation of the gapless mode depends on the breaking pattern and on the chemical potential. It is quadratic in the case of spontaneously broken particle number symmetry and zero chemical potential, and linear in all other cases. We constructed the leading terms in the dilaton effective action in two ways: first by using symmetry arguments and second by employing the Newton-Cartan geometrical structure.
The effective action of the dilaton in relativistic field theories encodes the information of trace anomalies. We considered the question whether and how the non-relativistic scale anomalies are encoded in the non-relativistic dilaton effective action. We found that there is a major difference between the relativistic and non-relativistic cases. The construction of the dilaton effective action in flat space requires working with a conformally flat background. However, such a background satisfies the Frobenius condition and therefore implies that the A-type anomaly is cohomologically trivial [20]. Thus, in contrast to the relativistic case, in order to distinguish the anomaly contribution from non-universal contributions to the dilaton action, one has to consider a curved NC background that violates the Frobenius condition. The study of such field theories and their consistency is an important challenge that we leave for the future. This is likely to shed light on the question whether non-relativistic RG theorems analogous to the relativistic ones exist. Another interesting topic left for the future, which is potentially linked to the structure of the RG flows, is the role of the special conformal generator in the symmetry algebra of non-relativistic field theories. This is the non-relativistic version of scale versus conformal invariance of relativistic field theories.

A The case of broken boosts with unbroken U(1) symmetry
In this appendix we derive the first few terms in a derivative expansion of the effective action of the non-relativistic NG bosons for the case of broken Lifshitz scale and Galilean boost symmetries, but unbroken U(1) symmetry. We also analyze the corresponding NG modes spectrum.

A.1 Symmetries
We consider the case where Galilean boost symmetries are spontaneously broken but particle number symmetry is conserved. In this case, we expect to find NG fields v i corresponding to the broken Galilean boosts generators, and a NG field τ corresponding to the JHEP10(2017)063 broken scale generator as usual. The transformations of these fields under non-relativistic scale transformation are clearly given by: whereas their transformations under Galilean boosts are given by: Using the field v i , one can build the covariant differential operator which is dimensionful but boost invariant. By acting with this operator on various covariant objects, we can obtain the possible effective action terms which include time derivatives. Furthermore, we may use simple space derivatives at will, as long as we keep rotation invariance. Finally, we compensate for the dimension using an exponent of τ of the corresponding dimension.
Listing the allowed terms up to two derivatives, including only derivatives of v i , we have: including only derivatives of τ we have: and finally including both τ and v i derivatives we also have: A cosmological-constant-like term Λ is also allowed by symmetries, but is forbidden in the NGB effective action in this case since it induces a gap. In parity violating theories in 2 + 1 dimensions one may also have terms such as ij ∂ i v j , but here we restrict ourselves to parity invariant theories (or to d > 2).
Overall, the Lagrangian density takes the following form up to second order in derivatives: (A.7) For the purpose of studying the linearized equations of motion, it is sufficient to expand this effective action up to second order in the field perturbations of τ around some valuê τ and v i around 0. The particular value ofτ can be absorbed by a redefinition of the coefficients in (A.7), and will be set to 0. Note that, since all terms in the action (A.7) are of second order in the field perturbations (except for the one time derivative term e −2τ ∂ t τ JHEP10(2017)063 which is a total derivative), the exponentials do not contribute at this order. Up to second order in derivatives and field perturbations (and taking into account integration by parts), one is left with the following independent terms: Note that out of these leading order terms, only v i ∂ i τ and ∂ t τ ∂ i v i mix the v i and τ perturbations. We may also consider the case of fully Schrödinger invariant theories, which requires in addition invariance under the special conformal transformation. For the field τ , this transformation is given in (3.13), whereas for v it is given by: Not all of the terms in (A.7) are invariant under this transformation -this added requirement restricts the possible terms to the following independent combinations: (A.10) In particular, out of the two mixing terms only the two-derivatives one is allowed. As we will see below, the other term can be ruled out by stability considerations as well.

A.2 Geometry
Similarly to subsection 3.2, we can also construct the NG boson effective action from geometric considerations. Following a similar path, we consider all possible gauge invariant expressions (since the U(1) particle number symmetry is unbroken in this case), perform Milne boost and Weyl transformations (see (3.15)-(3.17), (3.20)) and restrict to flat geometry (3.23).
To that end, we first define the following space tangent tensors (see e.g. [20]): where F µν is the field strength tensor associated with the gauge field A µ , as well as the previously mentioned a µ ≡ −L v n µ and the spatial Ricci scalarR. The zero derivative invariant scalar that can be constructed is the cosmological constant term Λ mentioned in subsection A.1 (and ruled out since it creates a gap). At first JHEP10(2017)063 order in derivatives, the only non-vanishing scalar that can be constructed is K S , which is a total derivative and therefore does not contribute to the effective action (again, we consider only parity invariant terms). The possible independent two derivative scalars, and their contributions to the effective action (after performing Milne boost and Weyl transformations and restricting to flat geometry) are as follows: 12 as well as theR given in (3.25) (as before, a 2 is related toR via integration by parts after restricting to flat geometry). These terms indeed agree with the terms found in subsection A.1 for the full Schrödinger invariant case (which, as mentioned in subsection 3.2, corresponds to this geometric construction).

A.3 Spectrum analysis
In the following, we analyze the spectrum of the low energy theory of the τ and v i perturbations. Starting from the scale and boost invariant Lagrangian density (A.7) derived in previous subsections, we expand it up to second order in the field perturbations and keep only the independent terms (A.8) to obtain: We can see that the transverse v mode follows a linear dispersion relation of the form: which is independent of the longitudinal and τ modes, while the longitudinal and the τ modes are mixed. For these mixed modes, one obtains the following equation: Solving for ω 2 , we get the dispersion relation: It can be clearly seen from this relation that when δ = 0, ω(k) becomes complex as k → 0. The stability of the modes therefore requires the coefficient E of the one-derivative mixing term v i ∂ i τ to vanish (as mentioned in subsection A.1, the same condition is required by Schrödinger invariance). Imposing this requirement, we get a linear dispersion relation ω ∝ k for the longitudinal v and τ modes. 13 Also note that when H = 0, the τ and v modes are mixed. We thus conclude that, making the physically reasonable assumption that unstable modes are forbidden, one finds d + 1 NG modes for the d broken Galilean boost symmetries and one broken scale symmetry as one would expect, all with linear dispersion relations, and with the longitudinal boost and dilaton modes being possibly mixed. 13 Note that positivity of the Hamiltonian requires that sign(B) = sign(C) = sign(G) = − sign(A) = − sign(F ). Given that E = 0, these conditions also ensure that there are no unstable modes. Also, note that we assume here α, β, γ = 0, otherwise analysis of the k → 0 limit requires expanding the action to higher order in derivatives.