## Higgs and Goldstone Modes in Crystalline Solids

Marco Vallone 1, a)

Dipartimento di Elettronica e Telecomunicazioni, Politecnico di Torino, corso Duca degli Abruzzi 24, 10129 Torino, Italy.

In crystalline solids, the acoustic phonon can be described either as a Goldstone or as a non-Abelian gauge boson. However, the non-Abelianity of the related gauge group apparently makes the acoustic phonon a frequency-gapped mode, in contradiction with the other description. In a different perspective overcoming this contradiction, both acoustic and optical phonon - the latter never appearing following the other two approaches - emerge respectively as the gapless Goldstone (phase) and the gapped Higgs (amplitude) fluctuation mode of an order parameter arising from the spontaneous breaking of a global symmetry, without invoking the gauge principle. In addition, the Higgs mechanism describes all the phonon-phonon interactions, including a possible perturbation of the acoustic phonon's frequency dispersion relation induced by the eventual optical phonon, a peculiar behavior able to produce mini-gaps inside the phonon Brillouin zone.

## I. INTRODUCTION

The spontaneous breakdown of a continuous symmetry implies the emergence of a massless bosonic particle for each broken generator of the involved symmetry group. This is the Nambu-Goldstone theorem 1,2 in a nutshell, forming the basis of the Standard Model of the fundamental interactions. In addition, the Higgs mechanism 3 , described as a case of spontaneous symmetry breaking (SSB), plays a prominent role in the Standard Model, providing the mass for the gauge bosons of the electroweak interactions and for the fermions 4 .

These concepts are seemingly far removed from physics of conventional semiconductors and metals commonly employed in electronics and optoelectronics industry. However, the Higgs mechanism, as described by P. W. Higgs in 1964 for particle physics 3 , is the relativistic analog of the plasmon phenomenon described by P. W. Anderson one year before in superconductivity 5 , and the existence of AndersonHiggs modes in condensed matter physics - in superconductors, cold-atoms in periodic lattices, in Bose-Einstein condensates, in antiferromagnets, in charge density waves 6-8 etc. -seems pervasive. However, beside plasmons phenomenology 5 and related plasmonics 9,10 , presently the Higgs mechanism has been acknowledged in semiconductors physics only in somehow exotic materials, like e.g. in topological insulators, Weyl semimetals, cuprates 11-16 , etc., although its role could be much more diffuse even in standard semiconductors.

On the other hand, the Goldstone theorem is very general and it holds also in the non-relativistic condensed state of matter, where the massless bosons correspond to collective excitations with wavevector k and gapless frequency dispersion relation w k , that is w k → → 0 0, where k = | k | . As a few examples, spin waves in the Heisenberg model are bosons arising when the ground state of the Heisenberg Hamiltonian is magnetically ordered 17 ; collective density excitations in superconductors arise from the spontaneous breaking of the electronic phase rotational U ( 1 ) symmetry 18,19 . Finally, as a more common example, in crystalline solids, collective excitations associated to lattice vibration modes are the acoustic phonons 20 ,

a) Electronic mail: e-mail marco.vallone@polito.it

and correspond to Goldstone modes emerging from the breaking of a continuous spatial symmetry, the translational invariance, broken by the presence of the crystal lattice 21,22 . It is said that all these collective excitations originate from one of the so-called emergence principles, in this case the Goldstone theorem, since they emerge from the very beginning.

Nevertheless, in condensed matter systems as well as in high energy particle physics, interactions are mediated by gauge bosons, appearing when local gauge invariance with respect to a given symmetry group is requested for the system's Lagrangian density 23 .

In a seemingly similar way, it has been shown 24 that by gauging the spatial translational group T ( 3 ) in crystal lattice, three gauge bosons appear to provide the local gauge invariance for the Lagrangian under the action of T ( 3 . ) They can be identified with the three acoustic phonons and, as a major feature of this approach, the elastic properties of solids and the acoustic phonon's dynamical equations can be described in close analogy with General Relativity field equations. In fact, it turns out that the acoustic phonon travels in the crystal acting as a wavelike perturbation of the lattice, similarly to the graviton in vacuum, that travels as a wavelike perturbation of a locally flat differentiable manifold, both obeying very similar field equations 25-30 . Acoustic phonons arise in this case not as Goldstone, but as gauge bosons.

A major concern regards whether the two descriptions could possibly be in contradiction and to what extent. Ref. [24] describes in detail the linear limit of the acoustic phonon's gauge theory, where the two descriptions appear in agreement, in particular providing in the long wavelength limit ( k → 0) the same gapless frequency dispersion law w k = csk , where cs is the sound velocity in the given medium. However, the cited work does not explore in depth the consequences of the non-Abelianity of the involved gauge group on the ensuing dispersion relation w k .

When more than one ion is present in the lattice elementary cell, another kind of phonon - the optical phonon 20 - constitutes a further and independent vibration mode. In polar semiconductors, the electron-longitudinal optical (LO) phonon emission is the dominant intersubband scattering mechanism responsible for electrons thermalization. However, although LO-phonons are of crucial importance in semiconductors transport theory and electron dynamics 31-35 , in the gauge the-

ory of crystal lattice's interactions they are left apart, not arising as gauge bosons.

glyph[negationslash]

In order to clarify these important points, in section II we give a short recap about the emergence of acoustic phonon as Goldstone boson. In section III, after an introduction about the issues behind gauging a spatial symmetry, we go beyond the linear limit of the theory, finding that the obtained dispersion relation may result gapped (that is, w k → 0 = 0), in sharp contrast with the gapless dispersion relation characteristic of Goldstone bosons.

In section IV we present a different and more general approach, showing that both acoustic and optical phonons may arise from a SSB which sets out an order parameter f . It follows that the Higgs mechanism plays an important role: the amplitude fluctuations of f are Higgs modes, for which a mass-like term appears in the Lagrangian, that makes the mode gapped, and we identify them with the optical phonons. Conversely, the phase fluctuations of f (Goldstone modes) are the acoustic phonons, here revisited in a much more general way that also sheds light on the seeming contradiction arisen when interpreting the acoustic phonon either as Goldstone or as gauge boson. Finally, in section VI main ideas and findings are summarized.

In this work A a and B b are contravariant and covariant four-vectors, ¶m is the partial derivative ¶ / ¶ x m , where standard Greek indices a b m , , , ... = 0 , ... 3 are indices of spacetime coordinates on a four-dimensional differentiable manifold with metric g mn and connection G a mn , latin indices ( i , j , k , ... = 1 2 3) mark spatial components, whereas Greek , , indices with an 'hat' ˆ m n , ˆ , ... = 0 , ... 3 are used for indices of local four-dimensional frames (vierbeins or tetrad indices) on a flat Lorentzian space-time with Minkowski metric hmn ˆ ˆ = diag (+ 1 , -1 , -1 , -1 . ) The Einstein's summation over repeated index is always understood, and not-italicized 'i' is the imaginary unit.

## II. GOLDSTONE THEOREM AND ACOUSTIC PHONONS

Apart from superconductivity and exotic materials, things are quite complicated even in the well-known world of solidstate crystals. Let us consider the Lagrangian density

$$\mathcal { L } = \mathrm i \psi _ { k } ^ { \dagger } \partial _ { t } \psi _ { k } - \frac { 1 } { 2 m ^ { * } } \nabla \psi _ { k } ^ { \dagger } \cdot \nabla \psi _ { k } \quad \quad ( 1 ) \quad \text{The} \quad \text{(} \quad \text{(} )$$

describing the dynamics of a non-relativistic Pauli electron with wavefunction y k ( t , r ) and effective mass m ∗ , free to move in the crystal and obeying the Bloch theorem 20 . L is symmetric (i.e. invariant) under global, continuous transformations described by the Galilei group G 36,37 (the corresponding relativistic formulation of L is Lorentz-invariant). If H is the group of time-translations, the generators of the quotient group G H / are the momentum p (translations), the angular momentum J (rotations), and the boosts K . Concerning the group of spatial translations T ( 3 ) in the ordinary space (one of the subgroups of G H / ), the elements of the infinitesimal form of T ( 3 ) are the operators UT = 1 -i k · d R , where d R is an infinitesimal displacement of the crystal ions in space around a Bravais lattice 20 translation vector R . UT may be written as UT = 1 -e j pj , where e j are real parameters, and the momentum components pj = -i¯ h ¶ j are the three generators of T ( 3 , whose Lie algebra is described by ) [ pi , pj ] = 0. Similar considerations could be made for the other two subgroups.

The underlying crystal lattice can be described 38 as a potential -V ( f ) to be inserted into L and eventually depending from several fields f = { f j } . If f 0 is a the value of f that minimizes V ( f ) , the system's ground state (the system's true vacuum) is the state for which f = f 0. We can expand the potential V around its minimum, obtaining at the second order

$$V ( \phi ) = V ( \phi _ { 0 } ) + \frac { 1 } { 2 } \left ( \phi - \phi _ { 0 } \right ) ^ { k } ( \phi - \phi _ { 0 } ) ^ { j } M _ { k j }, \quad ( 2 )$$

where

$$M _ { k j } = \left ( \frac { \partial ^ { 2 } V } { \partial \phi ^ { k } \partial \phi ^ { j } } \right ) _ { \phi _ { 0 } }$$

is a symmetric matrix that in the Lagrangian plays the role of a mass term, whose eigenvalues give the eventual masses of the fields { f j } .

The crystal lattice breaks both translational and rotational invariance making the system's ground state not symmetric , and this can happen even if the Lagrangian is symmetric. To be more clear, if we again consider the translations, the fields f j transform under the action of UT as f j → f j + a f ( ) p j , where a f ( ) is an infinitesimal parameter. Nevertheless, the symmetry of the Lagrangian can remain exact, provided

$$V ( \phi ^ { j } ) = V ( \phi ^ { j } + \alpha ( \phi ) p ^ { j } ). \text{ \quad \ \ } ( 4 )$$

Taylor expanding the Eq. (4), the same condition can be written as

$$\alpha ( \phi ) \frac { \partial V ( \phi ) } { \partial \phi ^ { j } } = 0,$$

and differentiating the Eq. (5) with respect to f k around the potential minimum f 0, we obtain

$$\left ( \frac { \partial \alpha ( \phi ) } { \partial \phi ^ { k } } \frac { \partial V } { \partial \phi ^ { j } } \right ) _ { \phi _ { 0 } } + \alpha ( \phi _ { 0 } ) \left ( \frac { \partial ^ { 2 } V } { \partial \phi ^ { k } \partial \phi ^ { j } } \right ) _ { \phi _ { 0 } } = 0. \quad ( 6 )$$

glyph[negationslash]

The first term vanishes, since f 0 is a minimum of V . Regarding the second term, it can be zero if a f ( 0 ) = 0, but this would be a trivial case: the symmetry would be exact, V ( f ) would be merely an arbitrary constant, it would not be necessary to introduce a field f , and the crystal lattice simply would not exist. Instead, if a f ( 0 ) = 0, it must be

$$\left ( \frac { \partial ^ { 2 } V } { \partial \phi ^ { k } \partial \phi ^ { j } } \right ) _ { \phi _ { 0 } } = M _ { k j } = 0,$$

stating that the field is massless. In summary, the crystal lattice makes a massless field f j to arise, the Goldstone excitation associated to p j . In the end, we expect three Goldstone bosons, one for each broken translations' generators, and the ground state's symmetry is said to be spontaneously broken.

It is possible to see that apart from p , also J and K are broken generators (i.e. the system's ground state is not symmetric under the corresponding transformations they generate), but it has been shown that they do not give rise to Goldstone bosons 30,39 . Furthermore, it should be remarked that at this level nothing can be said about phonon-phonon interactions.

For uniformity, it is also possible to introduce a pseudorelativistic notation, describing the acoustic phonons emerged by the Goldstone theorem as sound-like four-momenta p m states in a locally flat manifold M with metric g mn . The tangent space to any point of M is a four-dimensional manifold T M with the same local pseudo-Minkowski metric expression hm n ˆ , ˆ as in special relativity, and known as acoustic or sound metric 24,40,41 . Space-time coordinates on M are expressed as x m =( cst , r ) , and the world ' pseudo ' means that the relevant velocity is the sound speed cs in the considered medium (that of course is not a limit velocity).

Following this convenient formalism, instead of T ( 3 ) we can consider the translations group T ( 4 ) with elements

$$U _ { T } = 1 - \epsilon ^ { \mu } p _ { \mu } & & ( 8 ) \quad \pi ^ { - } \\ & \Phi \Phi$$

on the space-time M , where the parameters e m are point dependent. The four T ( 4 ) generators are p m = -i¯ h ¶m , and the crystal lattice breaks the three spatial components p 1 ... 3, but not p 0. Since p m p m = 0 (they are sound-like four-momenta in M ), a wave-like solution for them yields the gapless dispersion relation w k = csk , valid for small k , that makes the Goldstone bosons massless quasi-particles, as required 42,43 . The periodic structure of the lattice itself is known to produce Brillouin-zone folding and the appearance of gaps in the phonon spectrum, i.e., phonon stop bands, for wave vectors satisfying the Bragg condition 20 . It must be stressed that the present approach does not address this aspect of the problem, since at this level the lattice is treated as a continuum, an approximation valid for wavevectors k glyph[lessmuch] p / a , where a is the lattice constant.

## III. CHALLENGES IN GAUGING SPATIAL SYMMETRY

The existence of Goldstone bosons is closely related to the gauge theories of interactions. The possibility to choose freely a local parameter without changing the physics of a system was declared in 1929 by H. Weyl as a general principle 23 , known as local gauge invariance . It turns out that each time a Lagrangian density L ( y k ) is requested to be invariant under the action of an element UG of a Lie group G of local transformations y k → UG y k , one or more compensating fields arise, mediating an interaction (see the description for the classic case of electromagnetism e.g. in Ref. [4], for which the compensating - or gauge - field is the vector potential A m that mediates the electromagnetic interaction, and the involved gauge group is the Lie group U ( 1 , ) whose elements induce phase rotations on y k ).

In somehow similar way, it is tempting to describe acoustic phonons as gauge fields 24 , like in electromagnetism or, more generally, in Yang-Mills theories 44,45 , considering the local

(coordinate-dependent) form of UT , that is UT = 1 -e j ( r ) pj . Although UT is unitary, the Lagrangian L in the Eq. (1) is manifestly not invariant under the action of UT on y k because of the spatial derivatives in the expression of the generators. However, the naive pathway to request the Lagrangian L to be invariant under local spatial translations promoting T ( 4 ) to the role of gauge group is not trivial. Considering its importance, before going further it is worth clarifying why gauging a spatial translation is so different from gauging e.g. a wavefunction's phase transformation. In the gauge theory of electromagnetism, since the involved symmetry transformation is a phase rotation, it is said that an internal symmetry is gauged, because the transformation does not involve spacetime coordinates. In addition, all this has a simple geometrical meaning: with reference to Figure 1, given three differentiable manifolds M (the base space), L P ( ) (the fiber) and E (the total space), and a Lie group G (the structure group), a fiber bundle is defined as a topological structure ( p : E → M , L G , ) where, for each point P ∈ M a fiber L P ( ) is associated to P , such that a neighborhood VP of P is mapped in E according to p -1 ( VP ) , and that p -1 ( VP ) is homeomorphic to the Cartesian product VP × L . The fiber bundle is defined in conjunction with a group G which acts as a transformation group on the fiber, representing the different ways the fiber can be viewed as equivalent. The covariant derivative D m is a connection on E allowing to parallel transport 46,47 a vector in E . We may also visualize the n -dimensional manifold M as a n -surface, with an internal space L associated to each point P ∈ M with given topological structure. Specifically, in electromagnetism M is the ordinary four-dimensional space-time of the special relativity, a phase transformation on y k is a 'vertical' automorphism of the bundle and a diffeomorphism of L , such that the fiber above each point in M is unchanged (Figure 2(a)). In fact, since only the phase is changing, the point in the total space E does move, but just 'vertically' along L P ( ) . This is the gauge freedom carried by the fiber bundle within its fibers.

By contrast, trying to gauge a space-time transformation group like T ( 4 , ) it is quite clear that translations in M induced by T ( 4 ) are 'horizontal' diffeomorphisms of M , since a translation by definition changes the point P in the base

FIG. 1. Scheme of a fiber bundle, a map ( p : E → M , L G , ) , with the base space M and the fiber L , on which the associated group G acts. A trajectory g in M is mapped to g ′ in the total space E .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00918v3-figure-a8bb27b7-a1c2-4ec5-b561-5db67c49fd7b.jpg)

space, involving spatial coordinates. For this reason, an invariance under the action of T ( 4 ) is said to be an external symmetry. In this context, it is not easy to identify a fiber as an internal space L for which E = M × L locally holds, preventing a straighforward definition of a fiber bundle ( p : E → M , L T , ( 4 )) (Figure 2(b)).

## A. Gauge theory of spatial translations

Gauging T ( 4 ) in the pseudo -Lorentzian, locally flat manifold, equipped with the acoustic metric and possibly suitable to describe acoustic phonon dynamics, propagation and interactions, is mathematically not too different from the gauging of T ( 4 ) on the standard Lorentzian manifold of the General Relativity 48 . Nevertheless, the latter revealed a hard goal to obtain 26,29 . Everything in this field started from the pionieering work by R. Utiyama 49 in 1956, a good starting point, but not entirely convincing, especially because Utiyama did not address the central problem: translations are a diffeomorphism in M , and not in a fiber, as a gauge theory would require, a problem not solved neither by subsequent works by T. Kibble 50 and D. Sciama 51 . Following the usual Yang-Mills formalism, they defined the covariant derivative as D ˆ n = ¶n ˆ + W ˆ n , where W ˆ n is a gauge (compensating) field that, written in terms of its components R m ˆ n , reads

$$W _ { \hat { \nu } } = R ^ { \mu } _ { \hat { \nu } } p _ { \mu }. \quad \ \ ( 9 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Employing Eq. (9), D ˆ n was written as

$$D _ { \hat { \nu } } = h _ { \hat { \nu } } ^ { \mu } \partial _ { \mu } \quad \text{(10)}$$

where

$$h _ { \hat { \nu } } ^ { \mu } = \delta _ { \hat { \nu } } ^ { \mu } - i \hbar { R } _ { \hat { \nu } } ^ { \mu }.$$

The field h m ˆ n defines a set of four orthonormal vectors h m ˆ n ¶m , the tetrad or vierbein, provided the covariant index ˆ n is a local Lorentz index. The correct understanding of h m ˆ n was given by Y. M. Cho 52 , who developed a gauge theory of translations

FIG. 2. (a) The usual picture of a gauge theory for an internal symmetry: a gauge transformation in E moves the point p -1 ( P ) to p -1 ( P ) ′ , preserving the point P on M . (b) Translations are defined on the base manifold M itself, where the gauge transformation changes the point P to P ′ .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00918v3-figure-3930b480-8221-44fe-9019-596b68396dc8.jpg)

with a Yang-Mills-type Lagrangian, where the gauge potentials were correctly interpreted as translational connections (in particular, they are the nontrivial part of the vierbein fields h m ˆ n ), and not as general coordinate transformations on the base manifold (as e.g. in Utiyama 49 ), that would have been not correct. D. and G. Grensing 25 came to similar conclusions, obtaining the gauging of the Poincaré group in a form that allowed to express General Relativity as a gauge theory of this symmetry group. More recently, the Yang-Mills theory of the affine group (the semidirect product of translations T ( 4 ) and general linear transformations GL ( 4 , R ) ) was formulated 53,54 , where tetrads have been identified with nonlinear translational connections, for which the given h m ˆ n ¶m expression is a simplified yet correct version of the general formulation.

These concepts were not developed for solid-state physics, but rather to describe gravity, now correctly formulated as a Yang-Mills theory with T ( 4 ) as gauge group. Nevertheless, very interestingly the same formulation holds and can be exploited to develop a Yang-Mills theory of elasticity in solid crystals. In fact, Eqs. (9), (10) and (11) also hold when considering, instead of the standard space-time, the crystal lattice's manifold M with metric g mn , and with its pseudo -Minkowski tangent space T M equipped with the sonic metric hmn ˆ ˆ . In this case, the field components R m ˆ n can be identified with the tensor components of the crystal elasticity field, whereas in theory of gravity they are related to the Ricci's tensor 55 . This allows to formulate the physics of lattice vibrations in crystals (acoustic phonons) in a way that results somehow similar to the physics of gravitational waves in the ordinary space-time, as described in part in Ref. [24] and developed more in detail here in section III B.

## B. Yang-Mills theory of acoustic phonons in crystal

The Lie group T ( 4 ) can be promoted to the role of gauge group 24 , defining a gauge-covariant derivative D ˆ m as a connection that makes L gauge invariant, and consistent with the meaning expressed by Eq. (10). Following the proposed pseudo-relativistic formalism and the standard gauge prescriptions, the infinitesimal displacements of the crystal lattice induce local translations y k ( x m ) → y ′ k ( x m ) = UT ( x m ) y k ( x m ) of the electronic wavefunction, where all the ordinary derivatives must be replaced by gauge-covariant derivatives D ˆ n = h m ˆ n ¶m . The field W ˆ n must transform in turn according to

$$W _ { \hat { \nu } } \rightarrow W _ { \hat { \nu } } ^ { \prime } = U _ { T } W _ { \hat { \nu } } U _ { T } ^ { \dagger } - ( \partial _ { \hat { \nu } } U _ { T } ) U _ { T } ^ { \dagger }, \quad \ \ ( 1 2 )$$

and Eq. (12) implies

$$\left ( D _ { \hat { \nu } } \psi _ { k } \right ) ^ { \prime } = U _ { T } \left ( D _ { \hat { \nu } } \psi _ { k } \right ),$$

that makes the Lagrangian L gauge invariant under the local action of T ( 4 . )

The translational symmetry manifests itself through the occurrence of conserved currents j m . In the present case, the relevant gauge charges are the three generators pk = ∫ d x jk 3 ( x ) = ∫ d xT k 3 0 ( x ) of the spatial translations, whose corresponding currents are the components of the symmetric

energy-momentum tensor T k 0 . Furthermore, the total energy (Hamiltonian) operator H = ∫ d xc s 3 2 T 00 ( x ) is written in terms of p 0, the unbroken generator of translations in the time coordinate.

A major difference with respect to electromagnetism is the fact that the elements of T ( 4 ) do not commute, even if the group generators on a flat space-time commute, [ p ˆ m , p ˆ n ] = 0. In fact, if UT ( x ) and UT ( y ) are two elements of T ( 4 , it is )

$$[ U _ { T } ( x ), U _ { T } ( y ) ] = - \hslash ^ { 2 } \left [ \varepsilon ^ { \mu } ( x ) \partial _ { \mu } \varepsilon ^ { \nu } ( y ) - \varepsilon ^ { \mu } ( y ) \partial _ { \mu } \varepsilon ^ { \nu } ( x ) \right ] \partial _ { \nu }, \quad \text{free if} \quad \\. \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

glyph[negationslash]

that in general is nonzero for x = y , making non-Abelian the local form of the gauge group. The intrinsic difference between the structure of the global and local versions of T ( 4 ) can be stated more formally writing the currents algebra

$$\left [ j _ { \mu } ( x ), j _ { \nu } ( y ) \right ] = - i \hbar { \left ( } \partial _ { \mu } j _ { \nu } ( y ) - \partial _ { \nu } j _ { \mu } ( x ) \right ) \delta ^ { 3 } \left ( x - y \right ), \ \left ( 1 5 \right ) \quad \text{comm}$$

glyph[negationslash]

where the commutator in general is nonzero for x = y . Hence the generators' algebra can be written as 56

$$[ p _ { \mu }, p _ { \nu } ] = \mathrm i g f _ { \mu \nu } ^ { \gamma } p _ { \gamma } \, \quad \ \ \ ( 1 6 ) \, \stackrel { \mathbf \Theta \, \circ } { \mathbf w } \, \stackrel { \mathbf \Theta \, \circ } { \mathbf i }$$

glyph[negationslash]

where f g mn are the structure constants for the currents' Lie algebra and g is the coupling constant of the theory. It is f g mm = 0, and f g mn = -f g nm = 1 for m = n . In a flat space-time T M of course it is [ p ˆ m , p ˆ n ] = 0, but this is not true on M . This is a central point of the theory, and the algebra expressed by Eq. (16) makes non-Abelian the local symmetry, and it is responsible for the phonon-phonon interactions. On this basis, it also follows

$$[ W _ { \hat { \mu } }, W _ { \hat { \nu } } ] = \mathrm i g f _ { \beta \gamma } ^ { \alpha } R _ { \hat { \mu } } ^ { \beta } R _ { \hat { \nu } } ^ { \gamma } p _ { \alpha } \quad \quad \quad ( 1 7 ) \quad \quad \ \.$$

Making use of Eq. (10), the commutator [ D ˆ m , D ˆ n ] can be written in the alternative forms

$$[ D _ { \hat { \mu } }, D _ { \hat { \nu } } ] = G _ { \hat { \mu } \hat { \nu } } ^ { \alpha } p _ { \alpha } = G _ { \hat { \mu } \hat { \nu } } \quad \quad ( 1 8 ) \quad \underset { \Phi \Phi } { \text{aco} }$$

where the field strength tensor G ˆ ˆ mn and its components G a ˆ ˆ mn are given by

$$G _ { \hat { \mu } \hat { \nu } } = & \partial _ { \hat { \mu } } W _ { \hat { \nu } } - \partial _ { \hat { \nu } } W _ { \hat { \mu } } + [ W _ { \hat { \mu } }, W _ { \hat { \nu } } ] \\ G _ { \hat { \mu } \hat { \nu } } ^ { \alpha } = & \partial _ { \hat { \mu } } R _ { \hat { \nu } } ^ { \alpha } - \partial _ { \hat { \nu } } R _ { \hat { \mu } } ^ { \alpha } + \text{i} g f _ { \beta \gamma } ^ { \alpha } R _ { \hat { \mu } } ^ { \beta } R _ { \hat { \nu } } ^ { \gamma }. \quad \quad ( 1 9 ) \quad \text{whe}$$

It is worth noting that Eqs. (19) resembles expressions typical of classical Yang-Mills theories, although it should not be forgotten that the symmetry group rules an external symmetry, and the gauge field W ˆ m contains spatial partial derivatives, since the infinitesimal group generators act as differential operators.

The Lagrangian density for the free gauge field is

$$\mathcal { L } _ { R } = \frac { 1 } { 4 } G _ { \hat { \mu } \hat { \nu } } ^ { \alpha } G _ { \alpha } ^ { \hat { \mu } \hat { \nu } }, \quad \quad ( 2 0 ) \quad \ m \bar { \alpha } \\ \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

that explicitly reads

$$\exp \text{entry teams} \quad & \mathcal { L } _ { R } = & \frac { 1 } { 2 } \left ( \partial _ { \hat { \mu } } R _ { \hat { \nu } } ^ { \alpha } \partial ^ { \hat { \mu } } R _ { \alpha } ^ { \hat { \nu } } - \partial _ { \hat { \mu } } R _ { \hat { \nu } } ^ { \alpha } \partial ^ { \hat { \nu } } R _ { \alpha } ^ { \hat { \mu } } \right ) \quad & \text{can} \quad & \text{the} \\ & + i g f _ { \alpha } ^ { \gamma \delta } R _ { \gamma } ^ { \hat { \mu } } R _ { \delta } ^ { \hat { \nu } } \left ( \partial _ { \hat { \mu } } R _ { \hat { \nu } } ^ { \alpha } \right ) \\ & - \frac { g ^ { 2 } } { 4 } f _ { \beta \gamma } ^ { \sigma } f _ { \sigma } ^ { \delta \varepsilon } R _ { \hat { \mu } } ^ { \beta } R _ { \hat { \nu } } ^ { \gamma } R _ { \delta } ^ { \hat { \nu } } R _ { \varepsilon } ^ { \hat { \nu } } \quad & \text{(21)} \quad & \text{into} \quad & \text{to} \quad & \text{$\lambda$}$$

and describes the dynamics of a field with cubic and quartic (self-interacting) terms. Writing the Eulero-Lagrange equation

$$\partial ^ { \hat { \mu } } \frac { \partial \mathcal { L } _ { R } } { \partial \left ( \partial ^ { \hat { \mu } } R _ { \alpha } ^ { \hat { \nu } } \right ) } = \frac { \partial \mathcal { L } _ { R } } { \partial R _ { \alpha } ^ { \hat { \nu } } }$$

and imposing the Lorenz gauge (that allows for great simplification of all the expressions), the equation of motion for the free field R a ˆ n results (defining l = g 2 and ✷ = ¶m ¶ ˆ ˆ m )

$$\Box R ^ { \alpha } _ { \hat { \nu } } + \lambda \, f ^ { \sigma } _ { \beta \gamma } f ^ { \delta \alpha } _ { \sigma } R ^ { \beta } _ { \hat { \mu } } R ^ { \gamma } _ { \hat { \nu } } R ^ { \hat { \mu } } _ { \delta } = 0. \text{ \quad \ \ } ( 2 3 )$$

This equation does not include any mass-term (it would be a term linear in R a ˆ n ), but a cubic self-interacting term is present, coming from the quartic term in the Lagrangian.

All the components are highly coupled and difficult to manage. In order to understand the underlying physics, it is particularly interesting to consider in better detail the case of a simple one-dimensional atomic chain along x , assuming without loosing generality - that the perturbation is along the longitudinal direction (longitudinal acoustic phonon), hence R 0 0 = 0. From Eq. (16), that expresses the group commutation rules, the motion equations for the nonzero field components ˆ R = R 1 0 (the longitudinal acoustic phonon field) and ˆ S = R 1 1 (the spatial strain field) decouple, yielding

$$\left ( \partial _ { t } ^ { 2 } - c _ { s } ^ { 2 } \partial _ { x } ^ { 2 } \right ) \hat { R } + \lambda c _ { s } ^ { 2 } | \hat { R } | ^ { 2 } \hat { R } = 0 \text{ \quad \ \ } \left ( 2 4 \right )$$

$$\left ( \partial _ { t } ^ { 2 } - c _ { s } ^ { 2 } \partial _ { x } ^ { 2 } \right ) \hat { S } = 0.$$

ˆ S x t ( , ) describes the fluctuations of the crystal strain field, that in higher dimensional domains are coupled to phonons, which may get scattered by them. Regarding the equation for the acoustic phonon ˆ R x t ( , ) , the parameter l leads to a nonlinear oscillatory solution for ˆ R with a dispersion relation given by 57,58

$$\omega _ { k } ^ { 2 } = \sqrt { \frac { \lambda c _ { s } ^ { 2 } } { 2 } } \rho _ { R } ^ { 2 } + c _ { s } ^ { 2 } k ^ { 2 }$$

where r R is an integration constant depending on the cell characteristics. Eq. (24) can be derived by the Eulero-Lagrange equation from a Lagrangian that is a simplified form of Eq. (21),

$$\mathcal { L } _ { R } = \frac { 1 } { 2 } \left ( \eta ^ { \hat { \mu } \hat { \nu } } \partial _ { \hat { \mu } } R \partial _ { \hat { \nu } } R - \frac { 1 } { 2 } \lambda R ^ { 4 } \right ), \quad \quad ( 2 7 )$$

the well-known Lagrangian with f 4 -potential, describing a massless self-interacting scalar field 4 . Before concluding this section, we recap the path followed so far: the local gauge invariance imposed for the electron's Lagrangian in Eq. (1) makes three gauge fields to arise, identified with acoustic phonon modes, one longitudinal and two transverse. Their Lagrangian is the Eq. (21), for which a simplified, onedimensional form is the Eq. (27), describing phonon-phonon interactions through the nonlinear (quartic) term proportional to l .

It is also worth stressing that Eq. (26) describes a possible frequency gap √ l c s 2 / 2 r 2 R at the long-wavelength side of the spectrum. This gap, not described by the classical theory of acoustic phonons, is due to the fields' selfinteraction term coming from the nonzero commutators in Eq. (17) and Eq. (19), in turn originating from the current algebra in Eq. (15) and Eq. (16). However, since the coupling constant can be very small, the frequency (energy) gap may be a very tiny mini-gap, probably beyond practical observability limit. Moreover, it is also worth remarking that the present gauge theory is a massless, non-Abelian gauge theory, not a massive one. Apart from issues in its renormalizability, a mass-term in Eq. (27) would be a term proportional to R 2 , not existent in the present formulation, beside the quartic term coming from the nonlinearities. Nevertheless, the final result in the phonon dispersion relation would be similar (a mini-gap in the spectrum), but its very origin would be very different: a massive boson, and not a self-interacting massless boson as in the present case.

glyph[negationslash]

In summary, we may say that in solid crystals a Goldstone mode is a massless excitation, which very existence is due to the broken, continuous translational symmetry caused by the lattice, and its dynamics can be described in the framework of a non-Abelian gauge theory. It is worth considering that Eq. (26) reduces to the well known Debye form 20 w k = cs | k | , provided the nonlinear term in the Lagrangian can be neglected. Nevertheless, a relevant and quite unexpected result is that, in its general form, w k results gapped (i.e. w k = 0 in the limit k → 0) owing to self-interactions, if the latter are relevant, even though the gauge field is a massless boson.

As a side feature already noticed in Ref. [24], the commutator [ D ˆ m , D ˆ n ] satisfies a cyclic identity that can be written as

$$D _ { \sigma } G _ { \hat { \mu } \hat { \nu } } + D _ { \hat { \nu } } G _ { \hat { \sigma } \hat { \mu } } + D _ { \hat { \mu } } G _ { \hat { \nu } \hat { \sigma } } = 0, \quad \ \ ( 2 8 ) \quad \frac { \infty } { \text{ilv s} }$$

an equation known as the Bianchi identity for the field strenght tensor. The tensor R n ˆ m plays the same role of the Ricci tensor in General Relativity, obeying the same dynamic equations implied by Eq. (28), and describing the local perturbation of the sonic metric hmn ˆ ˆ of the solid crystal caused by the propagation of an acoustic phonon, in strict analogy with the space-time Lorentz metric perturbation induced by a gravitational wave propagating as described in General Relativity.

## IV. OPTICAL PHONON: A HIGGS MODE

In several binary compounds of interest in electronics, as e.g. GaAs, InP, GaN and many others, the crystal can be considered as composed of two interacting sublattices, one for each atomic species. Textbook toy-models describe the relevant physics considering a one-dimensional Bravais lattice with two ions per primitive cell. When treated as a set of elementary classical oscillators, they show that two distinct oscillation modes exist 20 :

a ) the two atoms in the cell may oscillate at frequency w k around their equilibrium position with the same direction and phase (Figure 3(a)), keeping unchanged their relative dis-

FIG. 3. Scheme of a diatomic one-dimensional chain showing the inphase (a) and out-of-phase (b) oscillating modes. (c) The Mexicanhat potential V ( f ) , leading to the excitation of the amplitude (or Higgs) and phase (or Goldstone) modes, s and q .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00918v3-figure-7baa5920-6712-479b-856c-0cf4043ee562.jpg)

placement s 0. The perturbation is the acoustic mode, an elastic wave propagating in the lattice at the sound velocity cs , and the corresponding collective ions excitation is the acoustic phonon, whose dispersion relation has an expression that in the long wavelength limit can be approximated as w k = cs | k | , indicating that the acoustic mode can be excited with arbitrarily small energy (there is no gap in w k );

glyph[negationslash]

b ) conversely, the two atoms in the cell may oscillate opposite to each other (Figure 3(b)), and their relative displacement s oscillates around an equilibrium position s 0 along a well-defined direction q 0 respect to a reference, defining an oscillating cell's dipole moment, if the solid is ionic. This is the well known optical mode, and the associated collective excitation is the optical phonon , whose role in semiconductors transport theory and electron dynamics is of great importance. It is found that w k → 0 → w 0 = 0: in other words, the optical phonon dispersion w k is gapped , and the optical mode's excitation costs a fixed nonzero energy ¯ h w 0 also for k → 0.

We have found that both acoustic and optical phonons may be described as arising from a spontaneous breaking of a global symmetry, and that they are respectively a Goldstone (phase) and a Higgs (amplitude) mode 8,59,60 of an order parameter oscillation. Beside its importance as a further sign of a widespread presence of the amplitude mode in condensedmatter systems with broken continuous symmetry, the present approach provides a unified formulation of phonons, where the coupling between the acoustic and the optical phonon modes arises in a natural way.

In order to show this important point, let us define the displacement between the two atoms in the cell as a complex field f =( s 0 + s ) exp ( i q ) . We can describe the dynamics of

the collective oscillation modes by the Lagrangian of the field f

$$\mathcal { L } _ { \phi } = \frac { 1 } { 2 } \left ( \partial _ { \mu } \phi \, \partial ^ { \hat { \mu } } \phi ^ { * } - V ( \phi ) \right ),$$

where V ( f ) is the interacting potential between the two atoms in the cell. With the choice V ( f ) = ( 1 / 4 ) l f | | 4 , we would obtain the Eq. (27), able to describe acoustic phonons. We could assume for V ( f ) the form

$$V ( \phi ) = \frac { 1 } { 2 } \frac { \mu ^ { 2 } c _ { s } ^ { 2 } } { \hbar { ^ { 2 } } } | \phi | ^ { 2 } + \frac { 1 } { 4 } \lambda | \phi | ^ { 4 }, \text{ \quad \ \ } ( 3 0 ) \text{ \ \ and } \\ \text{As } \text{.}$$

typical of the f 4 theories 4 with a mass term proportional to f 2 . However, if the two parameters m 2 (where m has the dimension of a mass) and l are both positive, the expectation value of f in the physical vacuum state is &lt; f &gt; = 0 (the minimum of V ( f ) ). In this case the system is in a disordered phase, since there is not a preferred direction for the atoms' displacement oscillation. Furthermore, the system's vacuum energy is invariant under a U ( 1 ) transformation, since L f is invariant under a global U ( 1 ) transformation of the field f . Consequenty, this model cannot describe optical phonons, for which q = q 0 results a preferred direction for the system. In addition, the introduction of the mass term is arbitrary, beside making the theory non-renormalizable. Furthermore, the theory would describe massive acoustic phonons, against all evidences.

Instead, observing that the optical phonon defines an oscillating cell's dipole moment, we may inspire us to ferromagnets, assuming that the global phonon field might have a preferred direction, violating a symmetry of the Lagrangian. In this case, the related field theory has a continuous, hidden symmetry that is spontaneously broken by a term in the Lagrangian that cannot be derived, but only inserted in V ( f ) in an ad hoc manner, as in the electroweak theory 4,38 . To this end, we proceed in the standard way according to the linearsigma model 59,60 . If the crystal formation makes m 2 negative, the potential V ( f ) assumes the well known 'Mexicanhat' shape shown in Figure 3(c), and all the couples ( s q , ) on the circle with radius | f | = √ -m 2 c s 2 / ( ¯ h 2 l ) = s 0 are V ( f ) minima. They identify degenerate vacuum states, all with the same energy, and among them the system chooses the ordered state for which &lt; f &gt; =( s 0 , q 0 ) , that becomes the true vacuum state of the system. The angle q 0 identifies the preferred direction imposed by the cell structure and, without loosing generality, we can set q 0 = 0 as reference x -axis. Hence, the Lagrangian in Eq. (29) with the potential in Eq. (30) and negative m 2 has a spontaneously broken symmetry.

The next important step is to write the Lagrangian L f in terms of the fields s and q , small fluctuations of f = ( s 0 + s q , ) around the true vacuum f 0 =( s 0 , 0 : )

$$\mathcal { L } _ { \phi } = \mathcal { L } _ { \theta } + \mathcal { L } _ { \sigma } + \mathcal { L } _ { \sigma - \theta }, \text{ \quad \ \ } ( 3 1 ) \text{ \ \ lati }$$

where

$$\mathcal { L } _ { \theta } = & \frac { 1 } { 2 } \partial _ { \hat { \mu } } \theta \partial ^ { \hat { \mu } } \theta - \frac { \lambda \theta ^ { 4 } } { 4 } \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$\mathcal { L } _ { \sigma } = & \frac { 1 } { 2 } \left ( \underbrace { \partial _ { \hat { \mu } } \sigma \partial ^ { \hat { \mu } } \sigma - \frac { \omega _ { 0 } ^ { 2 } } { c _ { s } ^ { 2 } } \sigma ^ { 2 } } \right ) - \sqrt { \frac { \lambda \omega _ { 0 } ^ { 2 } } { 2 c _ { s } ^ { 2 } } } \sigma ^ { 3 } - \frac { \lambda \sigma ^ { 4 } } { 4 } ( 3 3 )$$

$$\underset { \text{We} } { \text{ns} } \quad \mathcal { L } _ { \sigma - \theta } = - \sqrt { \frac { \lambda \omega _ { 0 } ^ { 2 } } { 2 c _ { s } ^ { 2 } } } \sigma \theta ^ { 2 } - \frac { \lambda } { 2 } \theta ^ { 2 } \sigma ^ { 2 } \,.$$

Here, for later convenience, we have introduced the real and positively defined parameter w 2 0 = -2 m 2 c s 4 / ¯ h 2 , a factor playing the dynamical role of a mass-term. The Lagrangians L q and L s respectively describe the dynamics of fields q and s . As expected and as it is typical for the Higgs' mechanism, q is massless (Eq. (32) does not include a term proportional to w 2 0 ), whereas in Eq. (33) a mass-term with the correct sign is associated to the field s . The Lagrangian L s -q describes the s -q coupling terms and will be treated in section V.

It is remarkable that the vacuum state is no more U ( 1 ) invariant (it is non-symmetric in q ), although Eq. (29) and Eq. (31) describe the same system: we only changed notations and chose a particular vacuum as system's ground state to expand the Lagrangian around. This is why we may speak of a symmetry spontaneously broken: no external agent provoked it and the symmetry of the system is still preserved, although 'hidden' by a particular choice for the ground state.

Taking the limit l → 0 and plugging L f into the EuleroLagrange equations, we obtain the motion equation for fields s and q for a one-dimensional diatomic chain along spatial coordinate x as

$$\left ( \frac { \partial ^ { 2 } } { \partial t ^ { 2 } } - c _ { s } ^ { 2 } \frac { \partial ^ { 2 } } { \partial x ^ { 2 } } \right ) \theta ( t, x ) = 0$$

$$\left ( \frac { \partial ^ { 2 } } { \partial t ^ { 2 } } - c _ { s } ^ { 2 } \frac { \partial ^ { 2 } } { \partial x ^ { 2 } } + \omega _ { 0 } ^ { 2 } \right ) \sigma ( t, x ) = 0 \text{ \quad \ \ } ( 3 6 )$$

whose oscillatory solutions GLYPH&lt;181&gt; e i ( wq s , ( k t ) -kx ) yield the dispersion relations

$$\omega _ { \theta } ^ { 2 } = c _ { s } ^ { 2 } k ^ { 2 }$$

$$\omega _ { \sigma } ^ { 2 } = \omega _ { 0 } ^ { 2 } + c _ { s } ^ { 2 } k ^ { 2 }.$$

They describe respectively a gapless and a gapped mode of the displacement field f (and consequently a gapless and a gapped mode of the ensuing cell's dipole moment, if the solid is ionic), resembling respectively the acoustic and the optical phonon's dispersion law.

From a more formal point of view, we can state that the SSB has generated the real and positively defined mass-like term w s 2 0 2 for the amplitude or Higgs mode s (the optical phonon). From the present approach it is also clear that the optical phonon is not a gauge field, since the oscillation costs at least the energy ¯ h w 0, and it is not possible to change ground state changing s without spending energy.

In addition, still owing to the spontaneous breaking of a continuous Lagrangian symmetry, a scalar field (here q ) appears and results massless , yielding the gapless dispersion relation w q 2 = c s 2 k 2 . In this limit, it is possible to change the

ground state of q without spending energy, just operating a gauge transformation along the valley of the Mexican hat (see Figure 3(c)). It follows that the field q is a pure gauge field that we identify with the acoustic phonon , arising as the Goldstone mode associated to the breaking of a continuous symmetry.

However, it must be stressed that the present approach also describes in an unitary way several nonlinearities, in particular the coupling between optical and acoustic phonons, an argument addressed in section V.

## V. NONLINEARITIES AND PHONONS-COUPLING

An important outcome of the present SSB approach is the natural emergence of phonon-phonon interactions and selfinteractions, ensuing from the terms in l in the expressions of Eq. (31). The dynamics of q alone, described by L q in Eq. (32), is the same we obtained gauging the spatial symmetry T ( 4 ) in the framework of non-Abelian gauge theories (see Eq. (27)), revealing a deep link between the two approaches. The term ( l / 4 ) q 4 describes four-phonon processes, like a scattering between two acoustic phonons or a decay of one into three acoustic phonons (all them are impossible in an Abelian theory, like e.g. the electromagnetism). The same argument also applies to optical phonons, whose scattering is ruled by the corresponding term ( l / 4 ) s 4 in Eq. (33).

A further important point is the functional form of the Lagrangian L s -q , Eq. (34), neglected during the preliminary analysis. It consists of two terms, and it confirms a possible coupling between acoustic and optical phonons, as shown in Refs. [61 and 62] for bulk GaAs, GaN, ZnO, MoS2, and for BN monolayers. In Eq. (34), the term proportional to sq 2 describes the decay of an optical phonon s into two acoustic phonons q (or the reverse process), as outlined e.g. in Ref. [63] fitting the theory to silicon. In the same equation, the term proportional to q 2 s 2 describes a scattering between two acoustic and two optical phonons, a four-phonon scattering term whose role has been recently recognized 64 in boron arsenide to be responsible for a substantial reduction of its expected thermal conductivity.

At the lowest order in l , Eq. (34) provides additional terms in the coupled equations of motion for q and s , that become

$$\left ( \frac { \partial ^ { 2 } } { \partial t ^ { 2 } } - c _ { s } ^ { 2 } \frac { \partial ^ { 2 } } { \partial x ^ { 2 } } + \sqrt { 2 \lambda } c _ { s } \omega _ { 0 } \sigma ( t, x ) \right ) \theta ( t, x ) = 0 \quad ( 3 9 ) \quad \text{W}$$

$$\left ( \frac { \partial ^ { 2 } } { \partial t ^ { 2 } } - c _ { s } ^ { 2 } \frac { \partial ^ { 2 } } { \partial x ^ { 2 } } + \omega _ { 0 } ^ { 2 } \right ) \sigma ( t, x ) = 0. \quad \ \ ( 4 0 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

The equation for s ( t , x ) is the same Eq. (36): without loosing generality, we can write its oscillatory solution as s ( t , x ) GLYPH&lt;181&gt; cos ( ws t -kx ) and plug it in Eq. (39). Approximating ws ≈ w 0, separating the variables as q ( t , x ) = q t ( ) t q x ( x ) , and supposing q x ( x ) GLYPH&lt;181&gt; exp ( ikx ) , the equation of motion for q t ( ) t is

$$\frac { d ^ { 2 } \theta _ { t } ( t ) } { d t ^ { 2 } } + \left ( c _ { s } ^ { 2 } k _ { n } ^ { 2 } + q \, \omega _ { 0 } \cos \left ( \omega _ { 0 } t \right ) \right ) \theta _ { t } ( t ) = 0 \quad ( 4 1 ) \quad \text{$\mathfrak { r i c a l }$}$$

where q = √ 2 l . For simplicity, in the calculations it is also possible to set cs = 1. The phonon wavevector takes the discrete values kn = n p / ( Na ) , where n = 1 , ... N , and N is the number of atoms in the chain. Eq. (41) can be treated as a Mathieu's equation 65-68 , which describes time-dependent harmonic oscillators perturbed by a periodic load (the term in the cosine) representing the interaction with the optical phonon. For l = 0 the eigenfrequencies are wq = cskn , recovering the standard acoustic phonon dispersion law in the long wavelength limit.

glyph[negationslash]

In the general case q = 0, we treated the perturbation induced on the cell by the optical oscillation as a boundary value problem according to Ref. [69]. In short, Eq. (41) was written as q ′′ t +( bn + q w 0 cos ( w 0 t )) q t , where bn is an unknown parameter, to be regarded as a discrete eigenvalue. Defining the domain as t ∈ [ 0 , T ] and imposing as boundary conditions q t ( 0 ) = 1, and q ′ t ( 0 ) = q ′ t ( T ) = 0 (cosine-like solutions), the eigenvalues bn were found iteratively, starting from a guess solution q t = cos ( n p t / T ) and a guess value bn . The convergence is fast, and the obtained parameters bn can be regarded as the discrete eigenvalues w q 2 associated to eigenfunctions q t .

glyph[negationslash]

In order to see how the optical-acoustic interaction qualitatively works, since the energy of optical phonons is usually higher than acoustics', it makes sense to set w 0 = cskN . It is interesting that for q = 0 the present formulation predicts the existence of mini-gaps opening inside the Brillouin zone, for wq ≈ 0 5 . w 0, as visible in Figure 4. This is a feature typical of Mathieu's equation: a periodic modulation induces a Bragglike condition that prevents the phonon propagation. Similar behavior was predicted and experimentally observed by Raman spectroscopy in superlattices 70-73 : in both cases, there is a periodic modulation (caused by the optical mode according to the present work, or caused by the superlattice in the cited cases) that induces a Bragg condition which prevents the phonon propagation when its frequency is around half of the modulation frequency. Furthermore, piezoelectrical coupling of optical and acoustic phonon modes was found also in zincblende GaAs quantum-well slabs 61,74 and theoretically calculated for wurzites in the Lagrangian formalism in Ref. [75], confirming by different approaches the existence of the phenomenology described in the present work.

## VI. CONCLUSIONS

We investigated in depth the emergence of acoustic phonons in crystalline solids, both as Goldstone and as gauge bosons. The acoustic phonon was previously acknowledged as the gauge boson appearing when the electron's Lagrangian in a crystalline solid is requested to be locally invariant with respect to the group of spatial translations, whose generators appear broken due to the lattice itself.

However, the gauging of the translations group is not trivial: in order to better explain the context, first we reviewed the differences between the gauging of spatial (external) and internal symmetries, as in Yang-Mills' theories. Then, exploiting the mathematical similarities between the sonic (or acoustic) metric and the Lorentzian metric of the ordinary space-time, we

1

FIG. 4. Dispersion relation for the acoustic phonon: the pure acoustic, unperturbed case q = 0, and the effect of the interaction with optical phonons, for q = 2 5 10, where a mini-gap opens for , , wq ≈ 0 5 . w 0 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00918v3-figure-9f18bea6-da17-4349-a48e-59e566fa1df4.jpg)

gauged the translations group in solid crystals following the same tetrad (or vierbein) formalism employed to express the General Relativity as a gauge theory with space-time translations as structure group. This allowed to compare the dynamical role of the crystal elasticity tensor in solid state physics and the Ricci's tensor in General Relativity at a more formal and rigorous level than before. The subsequent analysis showed that, beyond the linear limit, the translations group's non-Abelianity may generate a gap in the frequency dispersion relation of the acoustic phonon, as if a mass-like term were present in its free field Lagrangian, although this is not the case.

In order to better investigate the latter point, we considered a different scenario, driving attention to the problem of understanding the nature of the optical phonons in standard solids. We showed that both the acoustic and optical phonon emerge respectively as the gapless Goldstone (phase) and the gapped Higgs (amplitude) fluctuation mode of an order parameter arising from the spontaneous breaking of a global symmetry, without invoking the gauge principle. In greater detail, following the linear-sigma model approach the Higgs' mechanism is shown to generate a massive amplitude (Higgs) mode that we identify with the optical phonon, with gapped frequency dispersion relation due to the mass-like term. At

the same time, the Higgs' mechanism generates a massless phase (Goldstone) mode, the acoustic phonon. A frequencygap only appears in the strong nonlinear regime, and it is due to an anharmonic term, the same that arises from the gauging of the spatial translations group, an approach which did not provide any description of the optical phonon, though.

Although the SSB is a well known phenomenon, it turns out that it is able to describe not only novel aspects of material science (like e.g. Cooper pair plasmon modes in superconductors layered cuprates and thin films 8,76 ), but also an old acquaintance of the solid state physics, as the acoustic and optical phonons in crystals. In particular, the Higgs mechanism describes all the phonon-phonon interactions, including a possible perturbation on the acoustic phonon's frequency dispersion relation induced by the eventual optical phonon, a peculiar behavior not described so far, at the best of author's knowledge, able to produce unexpected mini-gaps inside the Brillouin zone.

- 1 Y. Nambu, Phys. Rev. Lett. 1960 , 4 , 380.
- 2 J. Goldstone, Il Nuovo Cimento 1961 , 19 , 154.
- 3 P. W. Higgs, Phys. Rev. Lett. 1964 , 13 , 508.
- 4 F. Mandl and B. Shaw, Quantum Field Theory , 2 nd edn., Wiley, UK 2010 .
- 5 P. W. Anderson, Phys. Rev. 1963 , 130 , 439.
- 6 P. B. Littlewood and C. M. Varma, Phys. Rev. Lett. 1981 , 47 , 811.
- 7 P. B. Littlewood and C. M. Varma, Phys. Rev. B 1982 , 26 , 4883.
- 8 D. Podolsky, A. Auerbach, and D. P. Arovas, Phys. Rev. B 2011 , 84 , 174522.
- 9 M. Blaber, M. Arnold, and M. Ford, J. Phys. Condens. Matter 2010 , 22 , 143201.
- 10 A. Politano, L. Viti, and M. Vitiello, Appl. Phys. Materials 2017 , 5 , 035504.
- 11 M.Z. Hasan and C. L. Kane, Rev. Modern Phys. 2010 , 82 , 3045.
- 12 Y. Ando, J. Phys. Soc. Japan 2013 , 82 , 102001.
- 13 X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B 2011 , 83 , 205101.
- 14 A. A. Burkov and L. Balents, Phys. Rev. Lett. 2011 , 107 , 127205.
- 15 B. Singh, A. Sharma, H. Lin, M. Z. Hasan, R. Prasad, and A. Bansil, Phys. Rev. B 2012 , 86 , 115208.
- 16 Z. M. Raines, V. G. Stanev, and V. M. Galitski, Phys. Rev. B 2015 , 92 , 184511.
- 17 C. A. Kendziora, I. A. Sergienko, R. Jin, J. He, V . Keppens, B. C. Sales, and D. Mandrus, Phys. Rev. Lett. 2005 , 95 , 125503.
- 18 U. S. Pracht, Electrodynamics of quantum-critical conductors and superconductors , Springer 2017 .
- 19 U. S. Pracht, T. Cea, N. Bachar, G. Deutscher, E. Farber, M. Dressel, M. Scheffler, C. Castellani, A. M. Garcia-Garcia, and L. Benfatto, Phys. Rev. B 2017 , 96 , 094514.
- 20 N. W. Ashcroft and N. D. Mermin, Solid State Physics , Saunders College Publishing, Philadelphia, US 1976 .
- 21 H. Leutwyler, Phys. Rev. D 1994 , 49 , 3033.
- 22 H. Leutwyler, Helvetica Physica Acta 1997 , 70 , 275.
- 23 H. Weyl, Z. Phys. 1929 , 56 , 330.
- 24 C. Dartora and G. Cabrera, J. Phys. A 2014 , 47 , 035004.
- 25 D. Grensing and G. Grensing, Phys. Rev. D 1983 , 28 , 286.
- 26 F. Gronwald, Proc. Gauge theories of gravitation , Jadwisin, Poland 1998 , 1121.
- 27 I. Low and A. V. Manohar, Phys. Rev. Lett. 2002 , 88 , 101602.
- 28 H. Watanabe and T. Brauner, Phys. Rev. D 2012 , 85 , 085010.
- 29 M. Blagojevi´ c, F. Hehl, and T. Kibble, Gauge Theories of Gravitation , Imperial College Press, UK 2013 .
- 30 T. Brauner and H. Watanabe, Phys. Rev. D 2014 , 89 , 085004.
- 31 P. Sotirelis, P. von Allmen, and K. Hess, Phys. Rev. B 1993 , 47 , 12744 .
- 32 P. Sotirelis and K. Hess, Phys. Rev. B 1994 , 49 , 7543.
- 33 M. Vallone, F. Bertazzi, M. Goano, and G. Ghione, Proc. Numerical Simulation of Optoelectronic Devices (NUSOD) , Palma de Mallorca, Spain 2014 , 99.
- 34 M. Vallone, F. Bertazzi, M. Goano, and G. Ghione, Phys. Status Solidi B 2015 , 252 , 971.
- 35 M. Vallone, F. Bertazzi, M. Goano, and G. Ghione, J. Appl. Phys. 2017 , 121 , 123107.
- 36 S. Sternberg, Group theory and physics , Cambridge University Press, UK 1994 .
- 37 R. Gilmore, Lie groups, physics, and geometry: An introduction for physicists, engineers and chemists , Cambridge University Press, UK 2008 .
- 38 M.E. Peskin and D. V. Schroeder, An introduction to quantum field theory , Westview Press, Perseus Books Group, Reading, MA, US 1995 .
- 39 T. Brauner, Symmetry 2010 , 2 , 609 .
- 40 W.G. Unruh, Phys. Rev. Lett. 1981 , 46 , 1351.
- 41 C. Barceló, S. Liberati, and M. Visser, Living Reviews in Relativity 2011 , 14 , 3.
- 42 J. R. Schrieffer, Theory of superconductivity , Addison-Wesley, Reading, MA, US 1988 .
- 43 N. Nagaosa, Quantum field theory in condensed matter physics , Springer, Berlin, Germany 1999 .
- 44 G. 't Hooft, 50 Years of Yang-Mills theory , World Scientific Publishing Company 2005 .
- 45 R. Jackiw, Rev. Modern Phys. 1980 , 52 , 661.
- 46 M. Daniel and C. M. Viallet, Rev. Modern Phys. 1980 , 52 , 175.
- 47 P. Fré, in Gravity, a geometrical course , Springer, Dordrecht, Germany 2013 , Ch. 2.
- 48 M. Hohmann, L. Järv, M. Krššák, and C. Pfeifer, Phys. Rev. D 2018 , 97 , 104042.
- 49 R. Utiyama, Phys. Rev. 1956 , 101 , 1597.
- 50 T. Kibble, J. Math. Phys. 1961 , 2 , 212.
- 51 D. Sciama, Rev. Modern Phys. 1964 , 36 , 463.
- 52
- Y. Cho, Phys. Rev. D 1976 , 14 , 2521.
- 53 R. Tresguerres and E. W. Mielke, Phys. Rev. D 2000 , 62 , 044004.
- 54 J. Martín-Martín and A. Tiemblo, Int. J. Geom. Meth. Mod. Phys. 2010 , 7 , 323.
- 55 S. Weinberg, Gravitation and Cosmology: Principles and applications of the general theory of relativity , Wiley, New York, US 1972 .
- 56 J. Gegenberg, S. Rahmati, and S. S. Seahra, Phys. Rev. D 2016 , 93 , 064025.
- 57 M. Frasca, J. Nonlinear Math. Phys. 2011 , 18 , 291.
- 58 M. Frasca, Eur. Phys. J. Plus 2017 , 132 , 38.
- 59 M. Gell-Mann and M. Lévy, Il Nuovo Cimento 1960 , 16 , 705.
- 60 A. Nicolis and F. Piazza, Phys. Rev. Lett. 2013 , 110 , 011602.
- 61 M. Willatzen and Z. L. Wang, Phys. Rev. B 2015 , 92 , 224101.
- 62 Y. Lu, T. Sun, and D. B. Zhang, Phys. Rev. B 2018 , 97 , 174304.
- 63 P. Klemens, Phys. Rev. 1966 , 148 , 845.
- 64 T. Feng, L. Lindsay, and X. Ruan, Phys. Rev. B 2017 , 96 , 161201(R).
- 65 M. Abramowitz and I. Stegun, Handbook of mathematical functions , Dover Publ., New York, US 1972 .
- 66 P.M.Morse and H. Feshbach, Methods of Theoretical Physics , McGraw-Hill 1953 .
- 67 R. Coïsson, G. Vernizzi, and X. Yang, Proc. 2009 IEEE International Workshop on Open-source Software for Scientific Computation (OSSC) , US 2009 , 3.
- 68 A. Parra-Hinojosa and J. Gutiérrez-Vega, J. Phys. Conf. Ser. 2016 , 698 , 012005.
- 69 U. Ascher, R. Mattheij, and R. Russel, Numerical solution of boundary value problems for ordinary differential equations , Society for Industrial and Applied Mathematics (SIAM), 1953 .
- 70 S. Tamura and J. P. Wolfe, Phys. Rev. B 1987 , 35 , 2528.
- 71 F. Calle, M. Cardona, E. Richter, and D. Strauch, Solid State Commun. 1989 , 72 , 1153.
- 72 K. Shinokita, K. Reimann, M. Woerner, T. Elsaesser, R. Hey, and C. Flytzanis, Phys. Rev. Lett. 2016 , 116 , 075504.
- 73
- A. Bruchhausen, N. Lanzillotti-Kimura, A. Fainstein, A. Soukiassian, D. Tenne, D. Schlom, X. Xi, and A. Cantarero, J. Phys. Conf. Ser. 2007 , 92 , 012160.
- 74 L. Duggen and M. Willatzen, Phys. Rev. B 2017 , 95 , 035310.
- 75 D. Datta, M. Dutta, and M. Stroscio, IEEE 17th International Conference on Nanotechnology (IEEE-NANO) , Pittsburgh, PA, US 2017 , 189.
- 76 R. Côté and A. Griffin, Phys. Rev. B 1993 , 48 , 10404.
Marco Vallone
Higgs and Goldstone modes in crystalline solids
In crystalline solids the acoustic phonon is known to be the
frequency-gapless Goldstone boson emerging from the spontaneous breaking of the
continuous Galilean symmetry induced by the crystal lattice. It has also been
described as the gauge boson that appears when the free electrons' Lagrangian
in the crystal is requested to be locally gauge invariant with respect to T(3),
the group of the infinitesimal spatial translations. However, the
non-Abelianity of T(3) makes the acoustic phonon a frequency-gapped mode, in
contradiction with its description as Goldstone boson. A different perspective
overcomes this contradiction. In fact, we show that both the acoustic and
optical phonon - the latter never appearing following the other two approaches
- emerge respectively as the gapless Goldstone (phase) and the gapped Higgs
(amplitude) fluctuation mode of an order parameter arising from the spontaneous
breaking of a global symmetry, without invoking the gauge principle. The
optical phonon's frequency-gap is present in all regimes, and it arises from a
mass-like term in the Lagrangian due to the Higgs mechanism itself. Instead, an
eventual acoustic phonon's frequency-gap appears only in the strong nonlinear
regime, and it is due to an anharmonic term, the same term arising from the
gauging of T(3), an approach which did not provide any description of the
optical phonon, though. In addition, the Higgs mechanism describes all the
phonon-phonon interactions, including a possible perturbation on the acoustic
phonon's frequency dispersion relation induced by the eventual optical phonon,
a peculiar behavior not described so far in these terms.