Simplest phonons and pseudo-phonons in field theory

We study the emergence of Nambu-Goldstone modes due to broken translation symmetry in field theory. Purely spontaneous breaking yields a massless phonon which develops a mass upon introducing a perturbative explicit breaking. The pseudo-phonon mass agrees with Gell Mann-Oakes-Renner relations. We analyze the simplest possible theories featuring gradient Mexican hats and describing space-dependent order parameters. We comment on homogeneous translation breaking and the connections with holographic Q-lattices.


Introduction and motivation
The breaking of spatial translations is fundamental to many physical systems, especially but not only in condensed matter. The lack of momentum conservation and the emergence of phonons are two important consequences of different kinds of translation symmetry breaking, explicit and spontaneous respectively. Typically the realization of translation symmetry breaking is technically complicated due to the spatial dependence of the fields. Although not being a conceptual obstruction, such complication usually requires non-analytic tools.
An important exception to this state of affairs has been explored in recent years within the holographic context, with the phenomenological aim of modeling stronglycorrelated systems that break translations and lack a standard quasi-particle description. Holographic models geometrize the renormalization group flow of a conjectured quantum field theory by means of dual gravitational fields, often having a radial dependence in an otherwise symmetric and homogeneous space. Spatial features can be added, yet the computations get much more difficult since the equations of motion become partial differential equations.
Wide families of holographic models avoid the leap in complexity at the price of introducing conceptual subtleties. Such models are characterized by a homogeneous breaking of translation symmetry, where an internal global symmetry renders all the spacetime points of the broken phase equivalent [1,2,3]. 1 Basic concepts like those of unit cell or commensurability [11] are thereby absent.
The present paper exports the idea of homogeneous translation symmetry breaking to standard, non-holographic field theory. Specifically, it studies the "Mexican hat potential" for translation symmetry breaking, both in a lita daniele.musso@usc.es 1 Models featuring massive gravity [4,5,6] have been considered too, they fall into the same category because they can be realized through a Stueckelberg mechanism in terms of fields that enjoy a global symmetry [7,8]. Massive gravity has also been employed to define effective field theories for phonons [9,10]. eral and a loose sense: literally, the mechanism of the Mexican hat is realized in the gradient sector; loosely, the aim at stake is to analyze simple prototypical field theory examples of translation symmetry breaking.
Phonons have been long considered in field theory, yet they are introduced ad hoc. This is the case for Fröhlich Hamiltonians describing electron-phonon interactions, and for effective field theories of elastic media [15]. Many phenomenological questions indeed are insensitive to the microscopic details. Nevertheless, there are reasons why elucidating the origin of phonons is important, both conceptually and phenomenologically: 2 .
-Low-energy effective field theories relies on information about the symmetries: the number of low-energy degrees of freedom is in general expected to be given by counting theorems descending from the symmetry breaking pattern [16]. -Spontaneous pattern formation is at the basis of the density waves physics [19,20]. Dynamical density waves are the best candidates to explain the anomalous transport properties of many strongly-correlated electron systems (e.g. high-T c superconductors both in the "normal" and "condensed" phases [20]). -Implementation and check of general expectation about the consequences of the symmetry breaking.
In this latter sense, we prove the validity of Gell Mann-Oakes-Renner relations for pseudo-phonons: 3 their squared mass is, at leading order, linear in the perturbation which breaks translations explicitly. A perturbative explicit component of translation symmetry breaking can model the effects of weak disorder. In particular, it pins a modulated order parameter preventing it from sliding freely. Weakly pinned density waves are the best candidate mechanism to explain the bad metal phenomenology [22,23,24]. The argument goes as follows: despite having little disorder, bad metals have exceptionally low dc conductivity due to a significant spectral weight transfer to higher frequencies. In other words, the optical conductivity has a marked peak at finite frequency, which subtracts spectral weight from null frequency. The soft finite-frequency peak can have a direct connection to a pseudo-phonon associated to a weakly-pinned, modulated order parameter. This hypothesis is further corroborated by the ubiquitous presence of spatial patterns throughout the phase diagram of strongly-correlated electron systems [25,26,27].
The models studied in the present paper shed light on possible field theory duals to holographic models which break translation symmetry homogeneously. 4

Main results
The main results of the paper are: 1. The emergence of phonons in a generic class of field theories (1), and the characterization of their dispersion relation (21). 2. The addition of a perturbative term that breaks translation symmetry explicitly and the appearance of a mass for the phonon (31). 3. Realization of homogeneous translation symmetry breaking in a purely field theoretic model (Section 4). 4. Construction of toy-models for the formation of concomitant density waves at an angle (Section 6).

Spontaneous breaking of translations
Consider the following action for a scalar complex field φ where A, B, F and G are real positive numbers and i, j are spatial indexes, i, j = 1, 2. The metric signature is (−1, 1, 1). The terms controlled by the couplings B, F and G are fourth-order both in the spatial derivatives and in the field. It is convenient to parametrize the complex field in terms of the modulus and the phase fields φ(t, x, y) = ρ(t, x, y) e iϕ(t,x,y) .
Model (1) is spatially isotropic but breaks Lorentz invariance explicitly, it enjoys spacetime translation invariance and global U (1) phase rotations of φ (namely shift symmetry for the phase field ϕ). The Euler-Lagrange equations for φ are given by and its complex conjugate equation EOM[φ * ]. Recall that the variational problem leading to the equations of motion assumes δφ = δφ * = 0 at infinity. Passing to the "polar representation" (2), one can obtain the equations of motion for the modulus and phase fields by means of the combinations Consider the following static but x-dependent ansatz: and plug it into the modulus equation of motion (4), thus obtaining while the equation of motion for ϕ (5) is automatically satisfied. The solutions to (7) are 5 To have real solutions we demand B + F > G. The energy density of a static configuration of the form (6) is given by and it is spatially homogeneous. In particular, on the solutions (8) and (9) the energy attains the following values Notice that E(k (3,4) ) represent degenerate global minima because A > 0. In Figure 1 we plot a specific example.
Having specified the background of interest (6), let us consider now the fluctuations and parametrize the complex fluctuation field δη by its real and imaginary parts The quadratic action for the fluctuations in Fourier space is given by where q = (ω, q x , q y ) and The entries of the matrix M defined in (15) are given by The eigenvalues e 1,2 of M are given by complicated expressions. To highlight the physical characteristics of the two modes, it is convenient to expand the dispersion relations e 1,2 = 0 for low momenta, thus obtaining where k is given by (9), and In (21), the positivity of the quartic term in momenta requires G ≥ F . We also require F = 0. Indeed, for F = 0 there would be neither pure longitudinal (q y = 0) nor pure transverse (q x = 0) propagation in (21). Note that F = 0 implies that we are keeping in the action (1) a term with more than one spatial derivative acting on a single field. Since a similar term in time derivatives would lead to Ostrogradsky instabilities [32,33], the condition F = 0 constitutes an obstruction to building simple relativistic generalizations of model (1). Combining the requirements above we have Equation (21) describes a gapless mode which is the phonon. 6 It represents the Nambu-Goldstone mode arising from the spontaneous breaking of the product of phase shifts and translations to its diagonal subgroup. The phonon dispersion relation is not standard, because ω is not linear in the momentum 7 , and the propagation is in general not isotropic (despite being x ↔ y symmetric). The dispersion relation given in (20) corresponds instead to a gapped mode. Both its mass and the leading q x term are not essentially related to the couplings F and G, but the transverse propagation at quadratic order would vanish for G = 0. In general such gapped mode does not propagate isotropically and its dispersion relation is not symmetric under x ↔ y exchange.
As a final comment, studying the fluctuations about a background (8) with k = 0, one finds that it is locally unstable.

Adding a small explicit translation breaking
We add to the action (1) a perturbation which breaks translations along x explicitly: where n is a perturbative coupling, and f κ is an explicit symmetry-breaking term controlled by the parameter κ.
Concretely we take: which vanishes for κ → 0. We consider again the ansatz (6), this time the equation of motion (5) for the phase field ϕ is not automatically satisfied, but takes the form To solve (25) we fix Upon considering (26), the equation of motion for the modulus field ρ reduces again to (4). The quadratic action for the fluctuations gets modified by the perturbation, in particular the entries of the matrix The only difference with respect to the spontaneous case is given by a new −2nκ 2 term in M τ τ . Expanding the two dispersion relations e 1,2 = 0 in low momenta and in small n, one obtains 8 and where k and κ are given by (26) and (9). The squared mass 2nκ 2 in (31) is linear in the perturbative coupling n which controls the explicit breaking, this agrees with Gell Mann-Oakes-Renner expectation for a pseudo Nambu-Goldstone mode, thereby (31) describes a pseudo-phonon. 9 Note also that the explicit breaking affects both the dispersion relations of the gapped and the Nambu-Goldstone modes at the quadratic level in q x . This is particularly relevant for the pseudo-phonon as it lowers the leading order at which the longitudinal momentum enters the dispersion relation. As far as the leading transverse propagation is concerned, the explicit breaking term does not alter the qualitative picture.

Comments on homogeneous translation symmetry breaking
A solution of the form (6) breaks the product of translations along x and ϕ-shifts to the diagonal subgroup, so a translation along x can be compensated by a global phase shift. This is the hallmark of homogeneous translation symmetry breaking: any spacetime point is equivalent to any other up to a global internal transformation.
There is no unit cell and k should not be strictly speaking interpreted as the wave vector of a lattice. 10 8 Recall that e1,2 are the eigenvalues of the matrix M . 9 Gell Mann-Oakes-Renner relations for phonons were obtained in holographic models in [28]. For an experimental study of gapped sliding modes, see for example [21]. 10 Relatedly, holographic models with backgrounds that break translations homogeneously do not feature commensurability effects for the formation of stripes [11].
The canonical energy-momentum tensor of model (1) is given by The momentum density vanishes T tx = T ty = 0 on a solution of the form (6). On the equations of motion we have ∂ µ T µν = 0, thus the Ward-Takahashi identities for translations are satisfied. On a solution (6) the energy density is given by where L 0 is the Lagrangian density written in (1) considered on the background (6). Similarly, the spatial components of the energy-momentum tensor are: where in the second step of (34) we have used the equation of motion (7). We define the pressure Comparing (33) with (34) and (35) we obtain p = T xx = T yy . The pressure is thus isotropic and the equation of state is given by 11 The U (1) current density is given by on a solution of the form (6) it vanishes too.
Remarkably, the explicit breaking term (24) does not introduce a source into the translation 1-point Ward-Takahashi identity. Indeed, consider the φ field transformation under a diffeomorphism and take ξ µ = δ µ x ξ. One has that the explicit breaking term in (23) transforms as follows: where in the last passages we have used both the ansatz (6) and the condition (26). The triviality of the 1-point Ward-Takahashi identity for translations, in spite of the presence of a term which breaks translations explicitly, is a peculiarity of homogeneous breakings. These have already been studied in holographic models, see [28] for instance. Analogous arguments show that the explicit term (24) does not introduce a source into the Ward-Takahashi identity of the U (1) symmetry either.
5 Adding a potential V (φ * φ) Consider the action (23) with the addition of a generic potential term respecting the global U (1) symmetry, Apart from providing a more generic situation, a potential V (φ * φ) is important for a specific reasons: as shown later, the dispersion relation of the phonon becomes in general linear in q x when a potential V (φ * φ) is considered. The explicit breaking term nf κ (x, φ * , φ) in (40) still requires κ = k, which in turn solves EOM[ϕ] = 0 automatically. The equation of motion (7) for the modulus field gets instead modified into The energy is given by and its minimization with respect toρ and k returns respectively the equation of motion (41) and the extra condition The composition of the equation of motion (41) with the condition (43) implies A variation δk in (6) corresponds to δϕ = δk x, which does not vanish for asymptotically large |x|, it even diverges. In particular, it lies outside of the hypothesis assumed for the variational problem to derive the equations of motion. One can nevertheless take δk variations and derive (43) relying on a regularization of the system to a finite box 12 and focusing exclusively on the dynamical character of the symmetry breaking. Said otherwise, one can neglect the boundary term in the variation of the action or, equivalently, assume that the boundary conditions be free. Since the minimization with respect to δρ alone already reproduces the only non-trivial equation of motion (41), the study of the energy minimization with respect to both δρ and δk selects a subset among the solutions to (41) (considered with free boundary conditions). Requirement (43) is crucial to lead to a flat direction, and in turn to a Nambu-Goldstone mode.
The entries of the matrix M (defined in (15)) for the quadratic fluctuation action become where (44) has been already considered, so we are taking fluctuations over a background that satisfies the energy minimization condition (43).
The eigenvalues e 1,2 of M yield the dispersion relations of two modes through e 1,2 = 0. Let us expand these dispersion relations in low momenta and in small n: and Where again κ and k are given by (26) and (9), we have introduced the positive quantity (see comment below (9)) and we assumed V (ρ 2 ) > 0, in accordance with a positive concavity for the energy (42) with respect to δρ variations. The potential V (φ * φ) does not alter the overall qualitative picture already obtained in Section 3.2; the dispersion relations (48) and (49) still describe respectively a gapped mode and a pseudo-phonon with an n-linear squared mass. Note however that V (φ * φ) would make a qualitative difference in the purely spontaneous case: if n = 0, the leading q x term in the dispersion relation of the phonon is either quadratic or quartic depending on V (ρ 2 ) being trivial or not. The phonon dispersion relation is where the longitudinal speed of propagation c ph is proportional V (ρ 2 ). Importantly, one can repeat the analysis of Section 4 also in the presence of V (φ * φ), still getting the same equation of state = −p obtained in (36). 13 Specifically, the on-shell value of the Lagrangian density L 0 is affected by the potential, but Equations (33), (34) and (35) are still valid.

Comments on lattice toy models
We still consider a setup with two spatial directions, x and y, and construct vacua with two coexisting spacedependent configurations (each one similar to that of Section 3) characterized by two vectors k 1 and k 2 . To this purpose we need two independent scalar fields, φ 1 and φ 2 .
We consider doubling the model (40) and adding a crossterm whose role is to choose vacua where k 1 and k 2 form a specific angle θ; namely k 1 · k 2 = |k 1 ||k 2 | cos(α) and we want to realize α = θ dynamically. Calling L[φ * , φ] the Lagrangian in (40) with n = 0, we consider The term in λ leads in fact to the minimization of For field configurations where k 1 and k 2 form an angle θ, the term in λ vanishes and the energy of the configuration is the same as that of two non-interacting spacedependent condensates (i.e. λ = 0). Since the λ term is the square of a quantity that vanishes on the background, also at the level of linear fluctuations there are no effects due to λ. In particular, the stability analysis of the background is not affected by λ.
Model (52) enjoys a global U (1) × U (1) symmetry and leads to vacua where |k 1 | = |k 2 |. It can however be generalized to cases with a smaller symmetry. As far as the λ term is concerned, there is an interesting alternative to (52), which preserves only a global U (1) symmetry and still leads to the minimization of (53). The phases of the two fields are locked and it is no longer true that a global internal symmetry transformation can compensate for a generic translation ∆x. This is a way to partially breaking the homogeneity down to transformations satisfying analogous result was found in [13] for modulated holographic phases where, despite the presence spontaneous modulations, the equilibrium stress-energy tensor assumes the perfect fluid form.

Discussion and future directions
The models studied in this paper can be generalized to spacetimes with higher dimensionality. The addition of transverse spatial directions does not affect the essential aspects of the computations and the features of the lowenergy modes. 14

Need for higher derivatives
We pursued the spontaneous breaking of translation symmetry following a "Mexican hat" strategy applied to spatial gradients, which implies that quartic terms in the spatial gradients are a necessary ingredient. This is apparent already from the equation of motion (7) where, upon setting B = F = G = 0, only trivial solutions with either k = 0 orρ = 0 remain. Studying the fluctuations about a background (6) with n = 0 and V (φ * φ) = 0, the Nambu-Goldstone frequency does not feature independent longitudinal and transverse propagation unless the coupling F is non-trivial, see (21).
To the purpose of finding the simplest possible models, we restrained the attention to terms whose order in spatial derivatives equals the order in the fields. 15 There are three observations about this: (i) it is remarkable that in order to obtain a propagating phonon the Lagrangian in (1) needs to be already quite complicated; (ii) it would be interesting to repeat the present analysis allowing for all the possible consistent terms. 16 (iii) Higher spatial derivatives are analogous to frustration. 17 14 The models could actually be considered also in 1 + 1 dimensions, suppressing the y direction. Still, the essential computational points about the pseudo Nambu-Goldstone modes and propagation along x would remain unaltered. On the interpretation level there are however extra subtleties: spontaneous symmetry breaking in 1 + 1 dimensions is impeded by Coleman-Mermin-Wagner-Hohenberg theorem. Such obstruction is avoided by quantum field theories in the large N limit, roughly because the fluctuations which would spoil the condensation of the order parameter are suppressed. It would be interesting to consider whether, already in 1+1 dimensions, the models studied here could be regarded as effective descriptions of theories in the strict large N limit. Two relevant references on these points are: [37] where the details of symmetry breaking in (large N ) holographic field theories in 1 + 1 dimensions have been analyzed; and [38] which argues on the existence of Nambu-Goldstone bosons in 1 + 1 dimensions, even without a strict spontaneous symmetry breaking. 15 In [17] a similar claim is made in relation to generic effective field theories for phonons. Comments on power-counting schemes for theories enjoying shift symmetries are given in [18]. 16 In model (1), the order 4 for spatial derivatives and 2 for temporal ones can be argued by assuming z = 2 non-relativistic scaling for the time coordinate and time-reversal invariance. 17 In discretized theories a modulated vacuum arises typically because of the competition between different couplings, e.g. the next-neighbor coupling competing with the next-to-nextneighbor coupling. If one discretizes the model (1), something similar would happen, for instance the next-to-next-neighbor The F term in (1) features more than one derivative applied on the same field. This could not be avoided even considering spatial partial integrations. Such a term, if covariantized, would lead to multiple time derivatives applied on the same field, which in turn would produce Ostrogradsky instabilities [32,33]. As a consequence, the models studied in this paper do not admit a trivial relativistic generalization.
The possibility of breaking translations spontaneously in relativistic theories is interesting in relation to the systematic approach of [17] which classifies all theoretically possible condensed matter systems in terms of their spontaneous breaking of Poincaré invariance. 18 The models of the present paper could possibly be thought of as effective low-energy descriptions of relativistic UV theories, assuming that a spontaneous Lorentz symmetry breaking has occurred along the renormalization group flow at a scale above the cutoff of the effective theory. Studying the possibility of such completions is a future direction. The embedding of the present models into a systematic effective field theory framework is a future perspective too. To this regard one could, or perhaps should, consider topological terms in the action and couplings with generic dependence on φ.

The role of boundary conditions
In the model (40), which features a potential V (φ * φ), it proved to be essential to minimize the energy with respect to both the parameters of the ansatz (6), k and ρ. This amounts to minimizing the energy without fixing the boundary conditions for the field φ. It is important to recall that the apparent larger freedom implied by relaxing the boundary conditions is actually compensated by an additional requirement on the solutions: the extra condition descending from the energy minimization with respect to variations of k.
In order to stress the role played by the boundary conditions we can consider a somewhat complementary example: a kinetic symmetry breaking where the translations are actually broken by boundary conditions. 19 Let us take the simple example of a free real scalar field whose Lagrangian is just ∂ µ ψ∂ µ ψ. This scalar field can be intuitively related to the phase field ϕ in (6) and in fact one can consider an ansatz ψ = kx. The equation of motion is satisfied for any value of k. To avoid problems with diverging fields, one can regularize the space to a finite box and impose boundary conditions that are compatible with the ansatz for a specific value of k. Said otherwise, if the boundary conditions are fixed, they dictate the value of k regardless of energy considerations. 20 There is a transcoupling due to the F term competes with the next-neighbor coupling due to the A term. 18 The generation of spatially modulated vacua in relativistic field theories has been considered in [34]. 19 See [35] which contains a similar discussion in relation to kink solutions induced by boundary conditions. 20 If instead the boundary conditions are free, the value of k is arbitrary. lation symmetry breaking, but it is forced by a "kinetic" constraint instead of being generated dynamically. 21

Topologically non-trivial configurations
The gradient Mexican hat potential discloses various possibilities to construct non-trivial topological objects. We briefly comment some instances, which are however regarded as future prospects.
A possible topologically non-trivial configuration is the gradient kink. When the "potential" for gradients allows for degenerate absolute minima (e.g. the solutions (8)), there are sectors where the solutions feature a jump in k necessary to connect two degenerate minima. If the direction of the gradient is along the direction of the step profile, the gradient kink is longitudinal and can be mapped to a normal kink by redefining the gradient field as the fundamental field, ∂ x φ → ψ.
Another slightly more exotic possibility is to combine the topological non-triviality in the gradient sector and that resulting from possible degenerate minima in the potential V (φ * φ).
As a final remark, the structure of the the models studied in the present paper, and specifically the fact that they contain terms with different signs and different scaling properties under spatial dilatations, makes them avoid Derrick's no-go theorem for the existence of finite energy solitons [36].

Comments on phenomenology: transport, sound and helical orderings
The models studied here can be embedded in larger theories and provide sub-sectors where translations are broken spontaneously or pseudo-spontaneously. The phenomenological properties of the larger system, like for example transport, would however depend on the coupling of the translation breaking sub-sector with the larger system and, of course, on the characteristics of the latter itself.
One cannot study the compressional sound mode from the equation of state (36) because neither nor p depend on the volume. 22 The first law of thermodynamics actually coincides with the equation of state (36). Remarkably, 21 The simplicity of the free real scalar ψ just described is deceitful: being quadratic in the fields, the model presents a gapless excitation even when the breaking is kinetic, such gapless mode however is not present in more general theories that feature kinetic translation breaking; in particular this gapless mode is not a Nambu-Goldstone because the breaking is not dynamical. To stress this point, one can consider again the complex field model (40) without letting the boundary conditions be free; k would be related toρ and the mode which corresponded to the phonon would get a mass proportional to V (ρ 2 ). 22 An analogous equation of state has been studied in [50] and claimed to be consistent with a vanishing speed of sound. Quintessence models have similar equations of states too, see [51]. these comments hold independently of the presence of a potential V (φ * φ), see Section 5. The sound mode is instead described upon interpreting the real field ρ sin(ϕ − kx) whose fluctuation is given by τ (introduced in (14)) as a displacement field in a target space. Since the model contains higher derivatives terms, one cannot just adopt the standard formulae for the speed of sound in terms of the elastic modulus. One cannot either adopt hydrodynamic formulae, as the system at hand is at zero temperature and zero density. Rather, the sound mode is studied by the fluctuation analysis of the quadratic action, as described in previous sections.
As a possible future application, it would be interesting to gain intuition on the qualitative behavior of holographic low-energy modes by means of a purely field theoretical toy-model. 23 This is especially interesting in order to match the finite temperature hydrodynamic modes to those of an appropriate T = 0 theory.
Eventually, the models at hand are technically similar to helical orderings studied in ferromagnetic systems [52]. 24 There the frequency of the spin-wave modes has been claimed to depend linearly on the longitudinal momentum and quadratically on the transverse momentum [53], as in (49) above.

Acknowledgments
This work has been funded by the Spanish grants FPA2014-52218-P and FPA2017-84436-P by Xunta de Galicia The paper is dedicated to the memory of Rosetta Cervetto.

A Related literature
Some useful references related to the present paper are: -An analysis of pseudo Nambu-Goldstone bosons based on the study of Ward-Takahashi identities has been performed for a U (1) symmetry in a generic quantum field theory in [41], with also a description of its 23 For a recent review on holographic sound modes see [39]; for a description of a holographic model closer to the field theories of the present paper we refer to [40]. 24 Or helical orderings in 3 He [54].
holographic implementation. Similar analyses in a nonrelativistic context are performed in [42,43]. -The study of spatial pattern formation is intimately related to the study of phase transitions from a translationally invariant fluid to a solid. We refer to the study of Landau which opened the field [44] and to two recent papers relying on similar techniques [58,59]. The branch of research related to stripe and density wave formation is huge too, for a review paper aiming at a systematic organization in the context of the copper-oxides we refer to [20]. -The problem of counting theorems for Nambu-Goldstone bosons in general circumstances entailing spacetime symmetries and non-relativistic contexts has gathered recent interest, see for instance [16,45,46,47,48]. -The spontaneous locking of internal and external symmetries as the central mechanism for studying effective field theories through coset constructions is described for instance in [49,17]. -The interactions between phonons are dictated by the symmetry breaking, this point has been studied in [15]. It is noteworthy that [15] claims that phonon selfinteraction terms are necessarily present on the basis of consistency and the current algebra. As shown in the main text, a similar conclusion emerges also from the description of the translation symmetry breaking dynamics. A precise comparison among these two claims is left for the future. -The Q-lattice strategy has been exploited in holography also to address time-dependent systems [55,56]. It would be interesting to consider whether suitable modifications of the models described in the present paper could allow one to study time translations in analogy to spatial translations. -For a pedagogic review on effective field theories with phonons interacting with electrons see [57].