{'fraction_non_alphanumeric': 0.03235775176561923, 'fraction_numerical': 0.007085088902791073, 'mean_word_length': 4.490791706327856, 'pattern_counts': {'":': 1, '<': 8, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 25, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Spontaneous Symmetry Breaking and Nambu-Goldstone Bosons in Quantum Many-Body Systems

Spontaneous symmetry breaking is a general principle, that constitutes the underlying concept of a vast number of physical phenomena ranging from ferromagnetism and superconductivity in condensed matter physics to the Higgs mechanism in the standard model of elementary particles. I focus on manifestations of spontaneously broken symmetries in systems that are not Lorentz invariant, which include both, nonrelativistic systems as well as relativistic systems at nonzero density, providing a self-contained review of the properties of spontaneously broken symmetries specific to such theories. Topics covered include: (i) Introduction to the mathematics of spontaneous symmetry breaking and the Goldstone theorem. (ii) Minimization of Higgs-type potentials for higher-dimensional representations. (iii) Counting rules for Nambu-Goldstone bosons and their dispersion relations. (iv) Construction of effective Lagrangians. Specific examples in both relativistic and nonrelativistic physics are worked out in detail.


Introduction
Symmetry considerations play a key role in our understanding of Nature. Ever since the birth of science, the aesthetic beauty of symmetry has appealed to the minds of natural philosophers, and later physicists, seeking for the origin of the laws of Nature. With the advent of quantum mechanics in the twentieth century, symmetry techniques based on group theory became an indispensable tool of a theoretical physicist. While in classical physics symmetries are straightforwardly incorporated and directly connected to physical observables by means of the Noether theorem, the situation in the quantum theory is more subtle. A fundamental theorem due to Wigner (see Sec. 2.2 in [1]) states that invariance of observables under certain transformation implies the existence of a unitary operator on the Hilbert space of states. If the symmetry transformation is moreover compatible with the dynamics of the system, this unitary operator commutes with the Hamiltonian and gives rise to a characteristic multiplet structure in its spectrum. However, there are systems whose dynamics is invariant under a symmetry transformation, yet this symmetry is not manifest in the spectrum or physical observables. We speak of spontaneous symmetry breaking (SSB). This phenomenon is ubiquitous in condensed matter physics where it is responsible, for example, for the peculiar behavior of superconductors, superfluids, and ferromagnets. In high energy physics, it underlies much of our current understanding of the fundamental interactions and the origin of masses of elementary particles.
This paper provides a pedagogical review of the properties of spontaneously broken continuous global internal symmetries with special emphasis on systems lacking Lorentz invariance that is, either intrinsically nonrelativistic systems or relativistic systems at nonzero density. (In the following, all these systems will for simplicity be referred to as nonrelativistic, unless explicitly indicated otherwise.) It aims to be a self-contained introduction into the subject, requiring just moderate knowledge of quantum field theory. Nevertheless, it will hopefully also offer some new material to the experts.
The plan of the paper is as follows. In Section 2 the notion of SSB is introduced and the mathematical subtleties associated with its implementation on the Hilbert space of a quantum system are discussed to some depth. A specific example, a free nonrelativistic field theory, is analyzed in detail in Section 3 illustrating explicitly the general material of the previous section, and at the same time showing some nontrivial features of SSB in nonrelativistic systems to be discussed in what follows. The technically easiest way to achieve SSB in an interacting field theory is to introduce an effective scalar field and adjust its phenomenological potential so that it has a symmetry-breaking minimum at the tree level. While in the case of the Higgs mechanism in the standard model of elementary particles or the Ginzburg-Landau theory for ferromagnets this is rather trivial, it becomes a difficult problem with the increasing dimension of the representation in which the scalar field transforms. Section 4 describes a general method to minimize Higgs potentials; the specific example of a spin-one color superconductor is worked out in detail. Section 5 presents, from the theoretical point of view, the most important part of the paper. I first formulate and prove the Goldstone theorem and point out its limitations in case of nonrelativistic systems, essentially following the classic review [2]. Then the issue of the counting rules for Nambu-Goldstone (NG) bosons is introduced, starting with the seminal contribution of Nielsen and Chadha [3], and then discussing the recent developments. In Section 6 the general results are elucidated using several explicit examples. Finally, Section 7 is devoted to the model-independent description of SSB based on low-energy effective field theory, worked out in the general nonrelativistic case by Leutwyler [4]. Section 8 provides a summary and offers some comments on spontaneous breaking of spacetime symmetries, not covered by this paper.
What I have regretfully omitted from this paper is a thorough discussion of applications of SSB and detailed examples of its appearance in realistic interacting quantum field theories. Even though this would certainly help the reader close the gap between the general principles and their concrete applications in current literature, I had several reasons for this omission. First, the list of problems to discuss would strongly depend on personal taste, and even with a very limited selection would increase the volume of the paper inadequately. Second, I prefer to stay at an as general level as possible so that the paper may be found useful by a wider audience. I thus discuss the principles and the model-independent approach of effective field theory, and only use specific model examples to illustrate the general results. Finally, dealing with realistic field theories would almost inevitably mean using suitable approximation schemes. One then has to make extra effort to show that the obtained results are physical, and not just artifacts of the approximation. On the contrary, most of the examples presented in the text are, sometimes in a restricted sense, exactly solvable.
Since this is a review, most of its content is essentially compiled from existing literature. However, some of the presented material is new. This concerns in particular Section 6: to the best of my knowledge this is the first time to give an explicit, and exactly solvable, example of a system exhibiting a NG boson whose energy is proportional to a noninteger power of momentum. The nonanalytic long-distance behavior of such a system clearly precludes any description based on conventional effective field theory. Furthermore, I have improved upon the proof of Theorem 4 on the NG boson counting in absence of charge density. Finally, as far as I know, even the identification of a free nonrelativistic particle as a type-II NG boson presented in Section 3 is novel. Throughout the text I use, for the sake of simplicity, the natural units in which = c = 1. Also, I use the Einstein summation convention and the Lorentz indices are contracted with the timelike Minkowski metric, g µν = diag(1, −1, −1, −1).
The literature dealing with various applications of the general principle of SSB is vast, and I therefore apologize to all authors not mentioned explicitly in the list of references for having omitted their valuable contributions. I have, at least, tried to include papers discussing to some extent the general properties of spontaneously broken symmetries, in particular those focusing on nonrelativistic systems. Besides papers covering specific subtopics to be discussed in the following sections, I would also like to point out several good reviews of the physics of Nambu-Goldstone bosons and their effective field theory description [5][6][7][8][9] as well as some specialized in the relativistic chiral perturbation theory [10][11][12].

Basic properties of spontaneously broken symmetries
Consider a quantum field theory whose action is invariant under a continuous group of symmetry transformations. By Noether's theorem, this gives rise to the existence of a conserved current j µ (x), satisfying the continuity equation, ∂ µ j µ (x) = 0. In the standard, Wigner-Weyl realization of the symmetry, the associated integral charge Q serves as a generator of the symmetry transformations, which are implemented on the Hilbert space of states by unitary operators. The ground state is assumed to be a discrete, nondegenerate eigenstate of the Hamiltonian. Consequently, it bears a one-dimensional representation of the symmetry group, and therefore also is an eigenstate of the charge Q. The spectrum of all eigenstates splits into multiplets of the symmetry, corresponding to irreducible representations of the symmetry group.

Ground state and finite symmetry transformations
A spontaneously broken symmetry can be naively characterized as such symmetry that the ground state is not an eigenstate of its generator. In order to give a more formal definition, let us follow the review [2]. Given the charge density j 0 (x), one introduces for an arbitrary finite space domain Ω the operator In the following I will use the same letter Ω to denote the domain and its volume; the precise meaning of the symbol will always be clear from the context. The symmetry breaking condition can be restated as the existence of a (not necessarily local) operator Φ such that where |0 is a translationally invariant ground state. This expectation value is known as the order parameter. Clearly, this formal definition immediately implies the previous one: if the vacuum were an eigenstate of the charge operator, the expectation value of this commutator would have to be zero. It is customary to identify Q(t) = lim Ω→∞ Q Ω (t) formally with the integral charge operator. However, this operator strictly speaking does not exist [13]. Indeed, translational invariance of the vacuum would then imply translational invariance of the state Q(t)|0 , and consequently which diverges unless Q(t)|0 = 0. This would, however, be in contradiction with the symmetry breaking condition (2). The intuitive picture of spontaneous symmetry breaking, based on the observation that a symmetry transformation does not leave the ground state intact, suggests high degeneracy of equivalent ground states. Indeed, since the charge operator commutes with the Hamiltonian, so will a finite symmetry transformation generated by this operator. It will therefore transform the ground state into another state with the same energy. As long as the symmetry group is continuous, we will find infinitely many degenerate ground states. On account of the fact that they are all connected by symmetry transformations, they must be physically equivalent and any one of them can serve as a starting point for the construction of the spectrum of excited states. The issue of the choice of a ground state will be discussed in more detail in the next subsection.
However, mathematical implementation of these ideas is subtle. The finite volume charge operator Q Ω (t) induces a "finite symmetry transformation," U Ω (θ, t) = exp[iθQ Ω (t)], which in turn gives rise to a "rotated ground state," |θ, t Ω = U † Ω (θ, t)|0 . However, very much like the limit lim Ω→∞ Q Ω (t) does not exist, the operator exp[iθQ(t)] is not well defined either. In fact, it can be proved that lim as will be shown later on explicit examples. It means that in the infinite volume (thermodynamic) limit, any two ground states, formally connected by a symmetry transformation, are actually orthogonal. The same conclusion holds for excited states constructed above these vacua. All these states therefore cannot be accommodated in a single separable Hilbert space, forming rather two distinct Hilbert spaces of their own. Any of these Hilbert spaces can, nevertheless, be taken as a basis for an equivalent description of the system, and the choice has no observable physical consequences.
Unlike the transformations of physical states, finite symmetry transformations of observables can be defined consistently in a theory which is sufficiently causal. Using the Baker-Campbell-Hausdorff formula one obtains for any operator A that As long as the theory satisfies the microcausality condition, that is, the commutator of any two local operators separated by a spacelike interval vanishes, and as long as the operator A is localized in a finite domain of spacetime, there will be a region Ω 0 such that the charge density outside this region does not contribute to the commutator, The transformation U Ω (θ, t)AU † Ω (θ, t) therefore has a well-defined limit as Ω → ∞. The expectation value of the rotated operator A θ,t;Ω in the vacuum |0 can then be interpreted as the expectation value of A in the rotated vacuum |θ, t Ω . One should nevertheless keep in mind that in the limit of infinite volume the correct formal definition proceeds as above.
Let me recall at this point that in Lorentz invariant theories the microcausality condition is automatically guaranteed. This is, of course, not the case for nonrelativistic field theories where interaction is often modeled by a nonlocal instantaneous potential rather than by an exchange of a propagating mode. However, the above argument can still be applied provided the range of the interaction is short enough [14]. Typically, it is sufficient when the correlations decay exponentially at long distance. If they decrease just as some power of distance, the conditions for applicability of the general arguments concerning SSB have to be inspected case by case.

Explicit symmetry breaking and the choice of the ground state
In the preceding subsection it was demonstrated that SSB gives rise to multiple degenerate ground states in which the order parameter acquires different values. I also stressed that any of these degenerate vacua can be chosen as the physical ground state. A natural question then arises: why does the system actually choose one of the states with a definite value of the order parameter, and not their superposition in which the expectation value of the order parameter could be even made to vanish? The answer to this question is tightly connected to the presence of external perturbations and the thermodynamic limit. A thorough discussion of this issue is given in Sec. 19.1 of [15].
In finite volume the ground state of a quantum system is usually (though not necessarily always) nondegenerate. One still has a large number of states formally connected by symmetry transformations, but their exact degeneracy is lifted by boundary conditions. (An exception is the frequently used periodic boundary condition that preserves translational invariance and thus inherits many of the properties of the infinite volume system. Nevertheless, this boundary condition seems somewhat artificial from a physical point of view.) The unique ground state is then given by a symmetric superposition of all the states, very much like the ground state of the quantum mechanical particle in a potential double well. The distance of the individual energy levels is essentially determined by the tunneling amplitude for the transition from one state with a definite value of the order parameter to another. In the infinite volume limit, this amplitude is exponentially suppressed and the ground states become perfectly degenerate. In principle it is indeed possible, though technically not so advantageous, to construct a quantum mechanical description of a spontaneously broken symmetry using a symmetric ground state. SSB is then manifested by long range correlations rather than nonzero vacuum expectation values [16].
In practice, the symmetry is almost never exact but is usually disturbed by small perturbations such as external fields. A typical example is the ferromagnet, where the role of the order parameter is played by the spontaneous magnetization. In principle, the magnetization can take any direction, but an arbitrarily weak external magnetic field will cause the magnetization to align with it. This phenomenon is general. As soon as the system volume is large enough, the energy difference induced by the external perturbation will be much larger than the tiny intrinsic splitting of the energy levels. The ground state will then be determined solely by the perturbation; this is called vacuum alignment. By choosing an appropriate perturbation one can select the corresponding vacuum. The key step is now that there is a basis in the space of degenerate ground states in which all observables become diagonal operators in the limit of infinite volume (see Sec. 19.1 in [15]). Therefore, it is the same basis of states in which the order parameter takes definite values, and one of which is selected by the external perturbation, regardless of what operator is used as the perturbation. This concludes the argument that the physical vacua are those with a definite value of the order parameter; their superpositions do not survive the infinite volume limit. One practical means of selecting the ground state thus is to switch on an external perturbation and then go to infinite volume. After this is done, the perturbation can be adiabatically switched off without disturbing the vacuum. It is symptomatic of SSB that the order of these two limits (infinite volume versus zero perturbation) cannot be reversed.
Before concluding the section, let me summarize the main features of SSB. First and most importantly, the symmetry is not realized by unitary operators on the Hilbert space, so it does not give rise to multiplets in the spectrum. A hallmark of SSB is the existence of an order parameter, that is, nonzero vacuum expectation value of an operator which transforms nontrivially under the symmetry group. There is a continuum of degenerate ground states, labeled by different values of the order parameter. Each of these ground states constitutes a basis of a distinct Hilbert space, all of them bearing unitarily inequivalent representations of the broken symmetry. This is the socalled Nambu-Goldstone realization of symmetry. From the phenomenological point of view, the most important consequence of SSB is the existence of soft modes in the spectrum, the Nambu-Goldstone bosons. Section 5 will be devoted to a detailed discussion of this phenomenon.

Example: Free nonrelativistic particle
As an illustration of the general abstract material of the previous section, let us consider the simplest example imaginable, the free nonrelativistic field theory (see also Chap. 2 of [17]). In order to demonstrate the variety of technical approaches to SSB, I will follow a different path: instead of defining the charge operator by an integral of the charge density over a finite space domain, the whole system will be quantized in finite volume Ω.

6
The dynamics of the system is governed by the Lagrangian density Obviously, the action defined by this Lagrangian is invariant under the following three kinds of transformations, ψ with a real parameter θ. (Complex conjugated transformations for ψ † are implied.) Note that the transformation #3, corresponding to the usual conservation of particle number, leaves the Lagrangian (7) itself invariant, while the transformations #1, #2 change it by a total derivative. Taking this into account, one infers from the Noether theorem the following list of conserved charge and current densities, The three independent types of transformations listed in Eq. (8) generate the symmetry group ISO (2), that is, the two-dimensional Euclidean group. Transformations #1, #2 correspond to translations in the plane defined by the real and imaginary part of ψ, while #3 yields the rotation in this plane. Demanding as usual that the Schrödinger picture fields ψ(x), ψ † (x) satisfy the canonical commutation relation, one easily derives the algebra of the Noether charges in finite volume, Q a = Ω d 3 x a (x), While the first two relations are identical to their classical counterparts, the two translations at the classical level commute. The last commutator in (11) hence shows that after quantization, the algebra ISO(2) develops a central charge, and the representation of finite group transformations becomes projective.
As an aside let me remark that a projective representation of a symmetry group of a quantum system may arise either as a result of nontrivial global topology of the symmetry group, or as a consequence of a central charge in the Lie algebra as above. For semi-simple Lie algebras the central charges can always be removed by a proper redefinition of generators. In the case of the non-semi-simple Euclidean algebras ISO(N), the central charges vanish for all N but N = 2, which is exactly our case. Further details on the projective representations of symmetry groups may be found in [1] (Sec. 2.7) or [18] (Chap. 13).

Hilbert space and inequivalent ground states
In finite volume one can always define the (countable) basis of the physical Hilbert space using the Fock construction. In the following, we will implicitly assume periodic boundary conditions so that (discrete) translational invariance is preserved and the basis of one-particles states may be chosen as the eigenstates of the momentum operator. These one-particle states, labeled by the three-momentum k, can thus be obtained from the Fock vacuum |0 by the action of creation operators a † k , |k ≡ a † k |0 . Let us suppose for the sake of generality that the one-particle states are normalized as k |k = N k δ kk . Then the annihilation and creation operators satisfy the corresponding relation [a k , a † k ] = N k δ kk . The canonical commutator of the fields (10) subsequently fixes the normalization of their expansion in the annihilation and creation operators, Using this expansion of the fields, one finds an explicit representation of the symmetry generators, From here we immediately see that Q 1,2 |0 2 = Ω → ∞ in the infinite volume (thermodynamic) limit. This clearly indicates that the operators Q 1,2 are ill-defined on the Hilbert space built above the Fock vacuum in the limit that the space volume is sent to infinity. The Hamiltonian of the system is in the second quantization expressed as Obviously, the Fock vacuum |0 represents one possible ground state with zero energy. However, we already know that other, equivalent, ground states may be obtained from |0 by applying symmetry transformations. Let us therefore denote Substituting from (13) yields and a simple manipulation shows that all these states are in fact coherent states, corresponding to the eigenvalue z √ N 0 Ω of the annihilation operator a 0 . Naturally, they are annihilated by all a k with nonzero momentum k. They have the same energy as the Fock vacuum since adding an arbitrary number of zero-momentum quanta does not change the total energy of the system. The magnitude of the scalar product of two such states is equal to This demonstrates explicitly the general fact that in the infinite volume limit, any two ground states connected by a broken symmetry transformation become orthogonal.

Spontaneous symmetry breaking and the Nambu-Goldstone boson
Let us first analyze the pattern of symmetry breaking when the ground state is chosen as the Fock vacuum |0 . Apparently, out of the three generators Q a only the third one preserves the ground state, the other two are spontaneously broken. In the field space this corresponds to the fact that the full Euclidean group ISO(2) is spontaneously broken to its rotation subgroup, SO(2). The Goldstone theorem (to be discussed in detail in Section 5) now asserts that for each broken generator there should be a state in the spectrum which couples to the associated current. In this case there is obviously only one such state, namely the one-particle state |k , and the corresponding amplitudes read All these formulas are valid in the Schrödinger picture where the fields are time independent. However, transition to the Heisenberg picture is trivial in this noninteracting field theory. Due to the dispersion relation, E k = k 2 /2m, the state |k represents a NG boson with momentum k. Since its dispersion relation is quadratic, it is the simplest possible example of a type-II NG boson, to be discussed in depth later. Note that there is a single NG boson that couples to two broken generators. As will become clear in Section 5 this is only possible when the commutator of the two broken generators develops nonzero vacuum expectation value. A glance at Eq. (11) shows that this is indeed the case. Moreover, our simple noninteracting field theory provides a rather nontrivial realization of this condition. The vacuum expectation value of the commutator [Q 1 , Q 2 ] is guaranteed by the existence of a central charge in the quantized theory, rather than by a nonzero charge density as one might naively expect based on the Lie algebraic structure of the symmetry.
What now if we choose as the ground state one of the coherent states |z with nonzero z? In field space this means shifting the ground state off the origin, Since the new ground state |z is (in finite volume) connected to the Fock vacuum |0 by a symmetry transformation, they are physically equivalent. The corresponding symmetry breaking patterns thus also have to be the same, just with a different basis of broken and unbroken generators. The geometric picture is clear: any expectation value of the field always breaks the two translations in the field space; it only preserves the rotations about the point z|ψ(x)|z = z. From a physical point of view, the states |z are very interesting (see Sec. 2.3 in [17]). Note that the expectation of the particle number operator Q 3 is The operator Q 3 itself is apparently not well defined in the infinite volume limit. However, the density z| 3 |z remains finite in this limit. The states |z therefore describe a finite density free nonrelativistic many-body system with all particles occupying the lowest energy level, that is, a Bose-Einstein condensate (this issue is discussed with much greater mathematical rigor in [19]). Different choices of z correspond to condensates with a different density and/or phase. A cautious reader may now be confused by my previous statement, based on general considerations, that all ground states |z are equivalent since they are connected by symmetry transformations. Yet, systems with different density do not look physically equivalent. The resolution of this seeming paradox is tightly bound with measurement theory [2]. Any device that can measure the system density must be able to distinguish states with different density, and therefore explicitly break those symmetry transformations that connect these states. On the other hand, if we only allow for a measuring apparatus that preserves the symmetry of the system, then the different ground states will really be indistinguishable.
Let us finally comment on the general technique that allows one to pick a single ground state out of the continuum of degenerate ones. One introduces an explicit symmetry breaking term in the Hamiltonian which splits the energy of the ground states. Then the infinite volume limit is performed and only afterwards the symmetry breaking term is removed. In our case such an explicit symmetry breaking term can be proportional to the particle number operator Q 3 . One replaces the Hamiltonian (14) with H µ = H − µQ 3 . The parameter µ plays the role of a chemical potential. In view of Eq. (20) the energy of the state |z is then z|H µ |z = −µΩ |z| 2 . In order that the system has a ground state, the chemical potential must be non-positive. For any strictly negative value of µ the Fock vacuum |0 becomes the single nondegenerate ground state, while the energy of all other coherent states goes to infinity in the thermodynamic limit. In other words, in presence of a Bose-Einstein condensate the chemical potential must be zero, as is of course well known. Obviously, any of the coherent states |z can be selected as the preferred one by a proper choice of the explicit symmetry breaking term.
Even though this section was devoted to a free nonrelativistic particle, it is amusing to compare to the free relativistic massless complex scalar field theory. There, the action possesses the same symmetry group ISO(2), while the spectrum contains two NG bosons with energy proportional to momentum: the quantum of the field and its antiparticle. The spectrum is therefore dramatically different even though the symmetry breaking pattern is the same. The reason is that in the relativistic case the shift generators Q 1 and Q 2 commute even after quantization since the contributions to the commutator from the particle and antiparticle sectors cancel each other. The algebra ISO(2) thus does not develop a central charge and the implementation of broken symmetry as well as the spectrum of NG bosons is usual.

Achieving spontaneous symmetry breaking: Minimization of Higgs potentials
Demonstrating that a given physical system exhibits SSB requires two key steps. First, one has to identify a suitable order parameter. This is usually very difficult to do from first principles. The reason is that SSB is a nonperturbative phenomenon, so it cannot be achieved at any finite order of perturbation theory based on a ground state which preserves the symmetry. One then often relies on observations (such as in the case of ferromagnets where the order parameter, the spontaneous magnetization, is obvious) or physical insight (such as in the case of superconductors, where pair correlations are crucial).
The second step is the actual calculation of the order parameter. From the qualitative point of view, there are essentially two approaches. The first, physically more satisfactory but technically much more difficult, is to find a self-consistent (thus nonperturbative) symmetry-breaking solution to the equations of motion, starting from whatever degrees of freedom are present in the theory. The NG bosons then typically appear as collective modes of these elementary degrees of freedom. This is the case, for instance, of the Heisenberg model of a ferromagnet to be discussed in Section 6, the celebrated Nambu-Jona-Lasinio model [20,21], or the technicolor theories in high energy physics. The most serious technical drawback of this approach is that one usually has to introduce an Ansatz for the symmetry-breaking solution and show in turn that it is consistent with the equation of motion. Thus one only gets such solutions that are put in by hand from the beginning.
The second, more phenomenological approach is to introduce an effective field with suitable transformation properties under the symmetry in question, which creates the NG bosons and whose ground state expectation value at the same time serves as the order parameter. This is done, for example, in the Higgs mechanism in the standard model of elementary particles or in the Ginzburg-Landau theory of phase transitions. While the minimization of a quartic potential with one complex doublet scalar field as in the standard model is straightforward, it can become a complicated task provided the order parameter transforms in a higher-dimensional representation of the symmetry group such as in some models of grand unification. It is therefore worthwhile to ask what can be said about the minima of a potential just from symmetry considerations. The investigation of this issue was initiated by Michel [22,23], and in the next subsection I will review some basic results. A heuristic method to minimize general Higgs-type potentials, developed by Kim [24], will then be presented and illustrated on an explicit nontrivial example.

Group action on the order parameter space
Let me introduce, in a non-rigorous manner following Sec. 6.2 of [25], some basic mathematical terms needed for the discussion of invariant potentials and their minima. Let G be a compact Lie group, acting smoothly on an infinitely differentiable manifold M. The points on this manifold may be thought of as values of the order parameter, and the action of a group element, g ∈ G, on a point φ ∈ M will be denoted simply as gφ. For each point we also denote as H φ the set of all group elements that leave φ invariant, that is, H φ = {g ∈ G | gφ = φ}. This set forms a subgroup of G and is called the little (or isotropy) group of φ. In physical terms, it consists exactly of those transformations left unbroken when the order parameter takes the value φ.
A very important notion that we will further use is that of the orbit G(φ) of a given point φ. It consists of all points of M which can be reached from φ by a (naturally, broken) symmetry transformation, G(φ) = {gφ | g ∈ G}. The relation defined by the condition that two points be connected by a group transformation is an equivalence, and the group orbits then define a decomposition of the manifold M into equivalence classes. Any potential V on the manifold which is invariant under the group action, V(φ) = V(gφ) for all φ ∈ M and g ∈ G, may be thought of as a function on the orbits. The minimization problem for a given potential on M can therefore be reformulated as a minimization of a function on the space of orbits.
It is clear from the definition of the orbit that two points on the same orbit have isomorphic little groups. This is in fact an immediate consequence of the stronger statement that the two little groups are conjugate, H gφ = gH φ g −1 . In addition to the orbit G(φ) there may be other points on the manifold whose little group is conjugate to H φ . For example, when M is a linear space such as in the Higgs mechanism, multiplying φ by any (nonzero) number we obviously get a point with the same little group as φ. The set of all points with little groups conjugate to H φ is called a stratum, S (φ). Intuitively, a stratum consists of all points of the same symmetry "class": they have the same unbroken subgroup and the same symmetry-breaking pattern. In the phase diagram, a stratum would be associated with a particular phase. At a phase transition, the order parameter moves from one stratum to another, and the symmetry class changes.
We are now ready to formulate a fundamental theorem which constrains possible stationary points of group-invariant potentials [22]: (Michel). Let G be a compact Lie group acting smoothly on the real manifold M and let φ ∈ M. Then the orbit G(φ) is critical, that is, every smooth real G-invariant function on M is It is worth emphasizing that the theorem makes no particular assumption about the form of the invariant function, so it may be a Higgs-type quartic potential as well as the full quantum effective potential whose power expansion would contain terms of arbitrarily high orders in φ. In case M is a linear space Theorem 1 is actually not very useful. For every point φ ∈ M its stratum S (φ) also contains all points λφ with nonzero λ, and the orbit G(φ) is therefore never isolated in its stratum unless φ = 0. So in this case, the theorem simply asserts a rather obvious fact that the origin φ = 0 is a stationary point of every G-invariant function on the space M.
On the other hand when the manifold M is not a linear space, Theorem 1 may have remarkably strong consequences. Let us consider a manifold M made of the values of the order parameter of fixed "length" and arbitrary "direction", that is, constrain the linear space of the previous paragraph by the requirement φ = 1, where · is a suitably defined norm. This is a situation one deals with in the low energy effective descriptions of spontaneously broken symmetries (such as the so-called nonlinear sigma model). As I will show in the next subsection, also the minimization of a quartic Higgs potential can under certain circumstances be reduced to minimization over directions of the order parameter, its norm being fixed. When classifying possible shapes of the order parameter by their symmetry properties, one often encounters the situation that the form of the order parameter is fixed up to a symmetry transformation (and trivial rescaling). These states are called inert in condensed matter physics (see Sec. 6.2.3 in [25]). By definition, their stratum consists of a single (hence isolated) orbit. The inert states are therefore found among stationary states of any G-invariant function on the manifold of order parameters of fixed norm. Explicit examples will be shown later.
Before concluding the general discussion let me stress that an invariant function can, of course, have other stationary states than those guaranteed by Theorem 1. Indeed, any state can be realized as a minimum of a suitably defined potential. To see this, note that a given group orbit can be uniquely specified by the values c α of a set of group invariants, say P α (φ). The desired potential then reads f (φ) = α [P α (φ) − c α ] 2 . Such potentials will, however, typically involve terms of high order. If we restrict ourselves to Higgs-type polynomials of fourth order, the list of allowed stationary states may become more narrow. Another point to emphasize is that Michel's theorem 1 makes a claim about stationary states of a given function, while in physical applications one would usually like to know the (absolute) minimum. To that end, Michel conjectured [23] that if the representation of the symmetry group G of the Higgs-type fourth order potential is irreducible on the real, its absolute minima are realized by orbits with maximal little groups. Even though this rule turned out not be hold in general (see [26,27] and references therein), it provides a useful guide for locating the minimum of Higgs potentials. A more rigorous and complete analysis of this problem may be found in [28].

Minimization of Higgs potentials
An elegant general method to minimize Higgs potentials was developed by Kim [24]. I will restrict to the simplest case where the order parameter transforms in an irreducible representation of the symmetry group. Then there is a single invariant quadratic term in the potential, proportional (We can without lack of generality assume that the components φ i of the order parameter φ are real.) Assuming further reflection symmetry or any other constraint which rules out a cubic term, the most general potential up to fourth order in the field reads where P (4) α is a set of algebraically independent fourth-order invariants in φ. (The mass term m 2 was written suggestively with a minus sign to ensure that a nontrivial minimum exists.) One of these invariants is always P (4) Factoring it out, one gets where A α (φ) = P (4) α (φ)/P (4) 1 (φ). These quantities are dimensionless and represent the "angles", that is, the orientation of the condensate in field space. The main observation now is that minimization of the potential with respect to φ is equivalent to successive minimization with respect to the norm φ and the angles A α . For fixed values of the angles, the norm of the field and the value of the potential in the minimum are given by In order to find the absolute minimum of the potential, one has to maximize φ min , and hence minimize the expression Ξ(A α ) ≡ α 1 λ α A α for given couplings λ α . This may be done using an appealing geometric picture. Consider for simplicity the case of three independent quartic invariants P (4) α (φ), that is, two angles, A 2 and A 3 . These cannot acquire arbitrary values, but rather span some domain in a two-dimensional plane which is indicated by the gray shaded area in Fig. 1; it will be referred to as the target space. The shape of the target space is independent of the couplings λ α , and is a sole characteristics of the algebraic structure of the symmetry group and its particular representation. The set of constant Ξ(A α ) is in the (A 2 , A 3 ) plane represented by a straight line with the normal vector (λ 2 , λ 3 ). If we choose too low a value of Ξ, the line will not intersect the target space, that is, there is no φ that would yield the desired value of Ξ. This corresponds to the dashed red line. As we increase Ξ, the line shifts parallel until it touches for the first time the target space (solid line, point K). This point determines the absolute minimum of Ξ(A α ) and thus also the absolute minimum of the potential V(φ). It is then a matter of mere algebra do find out which shapes of the order parameter map onto the point K in the target space.
This geometric picture can be conveniently used to scan the whole phase diagram as the couplings λ α are varied. First, it is obvious that the shape of the order parameter in the minimum does not depend on λ 1 ; this merely affects its overall magnitude, and at the same time must be large enough to ensure the boundedness of the potential from below. Coming back to the example with three quartic invariants, the phase diagram can then be plotted in the (λ 2 , λ 3 ) plane. For any values of the couplings, the ground state will always be represented by a point on the boundary of the target space. If the boundary is concave, as is the case of its edges adjacent to the point K, then by tuning the couplings continuously the ground state will reside at K, and then change abruptly to L or O. In other words, the system will undergo a first-order phase transition. If, on the other hand, the boundary has a convex segment, such as the edge between the M, N vertices, then tuning the couplings smoothly will result in a continuous change of the ground state. In the phase diagram, one will then observe an M-phase, an MN-phase (where the ground state travels from M to N and typically has a different symmetry structure), and finally an N-phase, the three of them separated by a sequence of two second-order transitions. An explicit example will be given in the next subsection.
When the number of quartic invariants differs from three, the situation is very similar. First of all, in the case of two invariants only the solution is trivial. The target space is one-dimensional  and there are just two phases, corresponding to the minimum and maximum allowed values of A 2 (φ). This happens, for instance, in the Color-Flavor Locked phase of dense three-flavor quark matter [29]. For more quartic invariants the target space is multidimensional, and keeping track of all faces and edges of its boundary, which determine the phases appearing in the phase diagram, may become increasingly difficult. Yet the presented technique seems to be the most elegant way to accomplish this task. The same method can also be used, with necessary complications, when the representation of the order parameter is reducible, or when other (cubic or higher-order) terms are present in the potential.

Example: spin-one color superconductor
As an explicit nontrivial example, I will now discuss the phases of a spin-one color superconductor. The details of all derivations may be found in [30], see also [31,32] for a partial discussion of the problem. A similar analysis for the case of d-wave pairing, that is, a theory of a traceless complex symmetric matrix with SO(3) × U(1) symmetry, was performed in [33].
Spin-one color superconductivity is a viable candidate phase for the ground state of cold dense quark matter at moderate densities (see [34] for a review). The most favored pairing pattern involves quarks of a single flavor, pairing in the color SU(3) antitriplet and spin SO(3) triplet channel. All we need to know for the purposes of the present paper is that the order parameter ∆ is a complex 3 × 3 matrix, transforming as ∆ → U∆R where U ∈ U(3) L is unitary and R ∈ SO(3) R is orthogonal. (The overall phase transformations, complementing the color SU(3), stem from conservation of baryon number.) Given the fact that a complex 3 × 3 matrix involves eighteen degrees of freedom, twelve of which can be transformed away, the order parameter is fully specified by a set of six independent real numbers. It can always be cast in the form This expression is very convenient for the analysis of the possible inequivalent structures of the order parameter, that is, the group orbits and strata. Focusing for simplicity on continuous symmetry transformations only, there are altogether eight nontrivial strata, summarized in Tab. 1. Three of them correspond to inert states: A, polar, and CSL (abbreviated from Color-Spin Locking). It is worth emphasizing that the isomorphy of little groups of two different values of the order parameter does not necessarily imply that they lie in the same stratum. Indeed, the cylindrical and ε, as well as the polar and A, phases do have isomorphic little groups.
To decide which of these phases occupy a part of the phase diagram one needs to write down the most general U(3) L × SO(3) R invariant potential of the Higgs type (21), where 16 The three independent quartic invariants satisfy the following inequalities which establish the shape of the target space: and finally which holds only for those values of ∆ satisfying P (4) 3 ≤ P (4) 1 /9. The conditions for the saturation of these inequalities determine which strata appear on the boundary of the target space, that is, what are the candidate ground states.
The target space of the spin-one color superconductor is plotted in the left panel of Fig. 2. The three inert states are represented by single points, forming the corners of the target space. This is in agreement with the general discussion since their strata consist of a single orbit each. Accordingly, they occupy the majority of the phase diagram, shown in the right panel of Fig. 2. However, a part of the phase diagram belongs to the noninert ε phase, which corresponds to the AX segment of the boundary of the target space. This phase was missed in a previous analysis [31]. In contrast to the sketch in Fig. 1, the remaining segments of the target space boundary are not strictly convex, but straight. Consequently, right at the phase transition between two adjacent phases, many more states can actually coexist than just the two. Thus, the A-polar line is occupied by all matrices of rank one, the polar-CSL line by all real matrices (up to a symmetry transformation), and the CSL-X segment by all matrices of the oblate type with ∆ 2 = ∆ 2 1 + a 2 . In view of Michel's conjecture mentioned at the end of Section 4.1 it is useful to investigate the little groups of the eight strata listed in Tab. 1. Their hierarchy is displayed in Fig. 3. The three inert states all have maximal little groups, which is actually a general feature of the inert states. On the other hand, the little group of the ε phase is non-maximal, and its presence in the phase diagram thus provides a counterexample to Michel's conjecture. Note that this does not follow from the mere fact that the ε state is noninert, as even such states can in principle have maximal little groups (see Sec. 6.2.4 in [25]). The ε and A phases are connected by a second-order phase transition; the fact that the little group of the ε phase is a subgroup of that of the A phase is in accord with the general Landau theory of phase transitions. At the other end of the hierarchy of symmetry breaking stand the oblate and N 1 states which have no nontrivial (continuous) subgroups, hence representing the most general states preserving a part of the symmetry.

Goldstone theorem
One of the most striking consequences of SSB is the existence of soft modes in the spectrum whose energy vanishes in the long-wavelength limit. This is the celebrated Goldstone theorem [35,36] and the soft modes are usually referred to as the Nambu-Goldstone bosons. Since the aim of the paper is to discuss SSB in nonrelativistic systems, which involves a number of subtleties, Oblate N 1 Figure 3: Hierarchy of the little groups in a spin-one color superconductor. Arrows indicate successive breaking into smaller and smaller subgroups. The relation G 1 → G 2 means that the group G 2 is conjugate to a subgroup of G 1 , or in other words, the orientation of the order parameter in the two phases can be chosen so that G 2 ⊂ G 1 .
I will now briefly review the standard proof of the Goldstone theorem in a version applicable to these systems. From a physical point of view, the most important ingredient responsible for the presence of NG bosons is, apart from the symmetry breaking itself, the existence of a conserved charge. In Section 2.1 the integral charge Q Ω (t) was defined in finite volume and it was shown that it does not exist in the limit Ω → ∞. Since SSB is based on the existence of the order parameter (2) which, being a commutator of the broken charge, is well defined even in infinite volume, let us inspect its time dependence. To that end, one takes the commutator [∂ µ j µ (x), Φ] and integrates it over the domain Ω. Using the continuity equation and the Gauss theorem one arrives at The key technical assumption of the Goldstone theorem, apart from translational invariance of the vacuum, is the vanishing of the surface integral in the infinite volume limit. This, upon taking the vacuum expectation value, guarantees that the order parameter is time independent. The same remark as at the end of Section 2.1 applies here: in causal theories the vanishing of the surface term is guaranteed as long as the operator Φ is localized to a finite domain of spacetime. (In practice, it is often even strictly local.) In acausal theories such as some nonrelativistic models with instantaneous interaction, the surface integral tends to zero in the infinite volume limit when the interaction is of finite range or decreases exponentially with distance. In case of long-range interactions, however, the disappearance of the surface term must be checked case by case.
With these remarks in mind, I will from now on simply assume that the surface term in Eq. (29) vanishes. This makes sure that the order parameter is time independent as Ω → ∞. It is given by the spatial integral of 0|[ j 0 (x, t), Φ]|0 . Inserting partition of unity in terms of eigenstates of the Hamiltonian and using translational invariance of the vacuum, this becomes The summation runs over discrete labels distinguishing different excitation branches in the spectrum as well as continuous internal variables of the multiparticle states. The states are normalized by n k |m q = (2π) 3 δ mn δ 3 (k − q). Integral over total three-momentum of the modes is indicated explicitly. Integrating now over the domain Ω yields In a large volume, this function is strongly peaked around k = 0 and thus only states with low momentum will contribute to the right hand side. At the same time we know that as Ω → ∞ the left hand side becomes time independent. This is only possible when the energy of the contributing states vanishes in the limit of zero momentum. Finally, we use the broken symmetry assumption which states that the left hand side actually is nonzero, so that there must be at least one state that couples to both the broken current j 0 (x) and the interpolating field Φ. The most general formulation of the Goldstone theorem therefore is: Theorem 2 (Goldstone). Spontaneous breaking of a continuous global internal symmetry implies the existence of a mode in the spectrum such that lim Several remarks to the theorem, its technical assumptions and derivation are in order. First of all, the theorem guarantees existence of a NG mode in the spectrum. In the most general formulation it does not tell us how many NG bosons there are. In Lorentz invariant theories there turns out to be exactly one NG boson for each broken generator. (Also, due to the strong constraints on the form of the dispersion relation, a NG boson is then simply a massless particle.) Perhaps this is the reason why, somewhat unfortunately, many textbooks focused on relativistic quantum field theory and its applications to particle physics overlook the issue of NG boson counting. Here, it will be discussed in detail in the following subsection for the case of internal symmetries. It is also quite nontrivial when spacetime symmetries are spontaneously broken. For more details, the reader is referred to [37]. Second, the Goldstone theorem provides information about the low energy behavior of the NG boson. Since it constitutes a whole excitation branch in the spectrum, it will (but does not necessarily need to) presumably also exist at high momentum. However, in this range its behavior is determined by details of short distance physics, hence it is nonuniversal and cannot be predicted solely from the broken symmetry. Examples are the acoustic phonon in crystalline solids whose dispersion relation is linear at low momentum, but gradually flattens until the group velocity becomes zero at the edge of the Brillouin zone, or the Bogolyubov-Anderson mode in Bose-Einstein condensates whose dispersion at high momentum becomes that of a free boson, essentially insensitive to the presence of the condensate.
Third, as was repeatedly stressed above, a sufficient technical condition for the integral charge Q to be time independent, and therefore for the Goldstone theorem to hold, is the causality which is inherent in all Lorentz invariant theories. In such a case the consequences of the Goldstone theorem cannot be escaped. Yet it does not mean that the predicted gapless NG boson can actually be observed, for it may appear in the unphysical sector of the Hilbert space. A distinguished class of systems where this may happen are gauge theories. In continuum, they are ill-defined unless the gauge is fixed appropriately. The remaining global symmetry can then be broken spontaneously. In covariant gauges that preserve Lorentz invariance and hence causality, a NG boson appears, but it is in the unphysical part of the Hilbert space, as mentioned above [38][39][40]. On the other hand, in noncovariant (such as Coulomb) gauges, the presence of a long-range interaction invalidates the Goldstone theorem. This is well known to happen, for instance, in electric superconductors where the soft mode acquires nonzero energy even in the long-wavelength limit. In the lattice formulation where there is no need for gauge fixing, spontaneous breaking of the gauge symmetry is ruled out by Elitzur's theorem [41].
Fourth, the assumption of unbroken translational invariance of the vacuum is natural. Without translational invariance, there is no conserved momentum to label the excitations so the whole quasiparticle picture loses its meaning. Note that for the sake of the low-energy physics that the Goldstone theorem is concerned with, only discrete translational invariance like in crystalline solids is in principle sufficient. The values of momentum are then restricted to a finite range, but this still allows for well-defined low-momentum excitations.
Finally, it should be made clear that the Goldstone theorem speaks of a limit of the energy as momentum goes to zero. This is ensured by the derivation above where integration is first performed over a finite space domain Ω. The function ϕ Ω (k) is thus smeared around the origin and only in the limit Ω → ∞ goes to (2π) 3 δ 3 (k). This distinction is very important since in principle there can be isolated states with zero momentum and energy which do not represent NG bosons, but still give a spurious contribution to the right hand side of Eq. (31). Thanks to the fact that the integrand is smeared by the factor ϕ Ω (k), these states fortunately do not affect the value of the momentum integral. In finite volume, such states may be, for example, the other, degenerate ground states. One can then show explicitly [14] that their contribution is suppressed when the limit of infinite volume is taken.
Concerning the observable consequences of SSB, I have focused exclusively on the presence of massless particles in the spectrum so far. However, the Goldstone theorem has in fact much broader consequences. The reason is that we can consider not just the (quasi-)local field Φ, but in general an n-point Green's function constructed from such fields as the interpolating field that gives rise to an order parameter. An immediate corollary of the Goldstone theorem then is that if two Green's functions, connected by a symmetry transformation, are different then there must be a NG mode in the spectrum. For example, if the masses of two particles, created by fields that lie in the same multiplet of the symmetry, are different then the spectrum exhibits a NG state [42]. This makes it clear that a spontaneously broken symmetry by no means resembles symmetry that is broken explicitly. Indeed, the former still gives rise to exact constraints on the Green's functions of the theory, conveniently summarized in terms of a set of Ward identities.

Goldstone boson counting: Dispersion relations
As already suggested in the previous subsection, it is in general a complicated task to determine the precise number of NG bosons and its connection to the number of broken symmetry generators. In order to appreciate how frequently one deals with nontrivial realizations of broken symmetry, let me list several examples from nonrelativistic as well as relativistic physics. The most profound example is that of ferromagnets versus antiferromagnets. In both the SU(2) symmetry of spin rotations is spontaneously broken by mutual spin alignment to the U(1) subgroup of rotations about the direction of total magnetic moment (staggered, in the case of antiferromagnets). This means that two generators are spontaneously broken. While an antiferromagnet possesses two NG modes (spin waves or magnons), both with a linear dispersion relation at low momentum, there is just one NG mode in a ferromagnet, its dispersion relation being quadratic. Similarly, in the so-called canted phase of ferromagnets, the SU(2) symmetry is completely broken, and there is one NG mode of either of the above mentioned types [43].
The same phenomenon was predicted to occur in multicomponent Bose-Einstein condensates of alkalic atomic gases [44,45]. In particular with three spin polarizations the global symmetry in question is SO(3) × U(1), corresponding to rotational invariance and conservation of particle number. There are two phases, roughly analogous to the polar and A phases of superfluid 3 He or the spin-one color superconductor analyzed in Section 4.3, leaving unbroken SO(2) and U(1) subgroups which are isomorphic, but realized differently due to a different structure of the ground state. (The respective values of the order parameter lie in different strata.) The three broken generators give rise to three NG bosons with linear dispersions in the polar phase, and to one linear and one quadratic NG mode in the A phase. Similar conclusions have recently been reported for spin-two Bose-Einstein condensates [46].
In a three-component Fermi gas [47][48][49], the global SU(3) × U(1) symmetry is spontaneously broken in presence of Cooper pairs built from fundamental fermions. Due to Pauli principle, these pairs transform as an antisymmetric tensor (antitriplet) of SU(3), leaving an SU(2) × U(1) symmetry unbroken. The five broken generators lead to one linear and two quadratic NG bosons. One should perhaps add that while in (anti)ferromagnets the dispersion relation of magnons can be probed directly by means of neutron scattering, similar direct measurements for superfluid atomic gases are beyond the reach of current experiments.
In relativistic physics, the issue of NG boson counting has been discussed primarily in the context of dense nuclear and quark matter. For instance, a condensate of neutral kaons can accompany the Color-Flavor Locked (CFL) phase of three-flavor quark matter at very high densities. Such a system exhibits a global SU(2) × U(1) symmetry stemming from isospin and strangeness conservation. The condensate leaves just its U(1) subgroup intact, thus breaking spontaneously three generators. These give rise to one linear and one quadratic NG boson [50][51][52]. The same behavior as in condensed matter ferromagnets can be observed in nuclear matter [53]. Yet other examples can be found among the various color superconducting phases of dense quark matter. In a two-flavor color superconductor, known as the 2SC phase, the symmetry as well as its breaking pattern are the same as in the above mentioned three-component Fermi gas [54]. However, in this case the symmetry is gauged, and as will become clear later, this result may thus be just an artifact of the model treatment based on a global color invariance. On the other hand, in a spin-one color superconductor, one can find spin waves and the "abnormal" NG bosons with quadratic dispersion relation are then physical [55].
All the preceding examples suggest that there is a deeper connection between the number of NG bosons and their dispersion relations. Quite generally, when there are some modes with quadratic dispersion relation, the total number of NG bosons seems less than the number of broken generators, in contrast to the naive expectation based on experience with Lorentz invariant systems. Incidentally, a great deal of understanding of this problem was achieved long before many of the above listed examples were even known. More than thirty years ago, Nielsen and Chadha [3] formulated the following theorem: Theorem 3 (Nielsen and Chadha). Assume that translational invariance is not completely broken spontaneously and that there are no long-range interactions. Then the energy of a NG boson is analytic in momentum. Denoting NG modes whose energy is proportional to an odd power of momentum as type-I, and those whose energy is proportional to an even power of momentum as type-II, the number of type-I NG bosons plus twice the number of type-II NG bosons is greater than or equal to the number of broken symmetry generators.
The requirement on the absence of long-range forces was technically formulated as a commutativity condition: for any two local operators A(x) and B(0) there is a real positive number τ such that It is this assumption which is responsible for the analyticity of the Fourier transform of the commutator (30), and in turn for the analyticity of the dispersion relation. Obviously it is a stronger assumption than is necessary for the Goldstone theorem itself. Indeed, in Section 6 I will show examples of systems with long-range forces where the Goldstone theorem is still valid, but the NG dispersion relation is not analytic. The actual proof of the claim about the number of NG bosons is algebraic and somewhat involved, so I will not repeat details and instead refer the reader to the original paper [3]. It is notable that the statement of the Nielsen-Chadha theorem is rather general. First, it does not specify the power of momentum to which energy is proportional. Second, it only gives an inequality for the number of NG bosons. On the other hand, I am not aware of any example of a NG boson whose energy would be proportional to a higher power of momentum than two. Also, the question of the possible general saturation of the inequality for the number of NG bosons is yet to be understood. Nevertheless, note that the counting may be obscured by the presence of gapless non-NG modes, as recognized recently in the context of spin-two Bose-Einstein condensates [46].

Goldstone boson counting: Charge densities
What I have not mentioned before in the list of examples exhibiting type-II NG bosons was that they are all accompanied by nonvanishing density of some of the conserved charges. Indeed, this is exactly the sought property that distinguishes ferromagnets and antiferromagnets: the ferromagnetic ground state features nonzero net spin density. The connection of the NG boson counting and the presence of charge density in the ground state is rather general and a number of (partial) results exists in this respect. I will start with a theorem due to Schäfer et al. [50]: Note that in [50] the theorem was formulated as a strict equality. However, this does not follow from the proof presented there. In general it is difficult to place an upper bound on the number of NG modes; see also the remark at the end of the previous subsection in this respect. The original version of the proof was based on the identification of vectors Q a |0 with the NG states. This may be troublesome because as explained in Section 2 these vectors are not well defined in the infinite volume limit. Also, it is not clear a priori that the number of linearly independent vectors Q a |0 is equal to the number of NG states in the spectrum. Below I present a slightly more rigorous modification of the proof which avoids this step.
Take j 0 a (x) as the broken charge density in Eq. (30) and j 0 b (0) as the interpolating field for the NG mode. Integrating over the space in the infinite volume limit then yields Denoting M an = 0| j 0 a (0)|n 0 , this equals (MM † ) ab − (MM † ) ba . The assumption of the theorem is equivalent to the requirement that MM † be symmetric, and hence real (it is automatically Hermitian). It may therefore be diagonalized by a real orthogonal transformation, MM † → RMM † R T , which is equivalent to M → RM, that is, to a change of basis of the broken generators. (It is essential that the transformation matrix R is real, for otherwise R ab Q b would not necessarily be a generator of the symmetry.) Assume now that the number of NG states, n NG , is smaller than the number of broken currents, n BC . The rank of the n BC × n NG matrix M then cannot be larger than n NG , and so cannot the rank of the n BC × n BC diagonal matrix MM † . The basis of broken generators may be ordered so that the first n NG diagonal elements of MM † are nonzero while the remaining n BC − n NG ones are zero. Since MM † is just the matrix of scalar products of rows of M, this means that the last n BC − n NG rows are zero. Consequently, the generators Q n NG +1 , . . . , Q n BC have zero matrix elements with all NG states |n 0 , which is in contradiction with the assumption that they are broken. The theorem is thus proved.
A short precaution is appropriate regarding the interpretation of Theorem 4. Since the symmetry generators furnish a Lie algebra, it is tempting to conclude that nonzero density of some of the conserved charges is a necessary condition for the number of NG bosons to be smaller than the number of broken generators [50,56]. This is indeed the case in all systems listed in the previous subsection. However, the otherwise trivial example of a free nonrelativistic particle (see Section 3) teaches us that it is possible to achieve nonzero vacuum expectation value of a commutator of two generators even if the expectations of all generators themselves vanish. The resolution is that the Lie algebra of symmetry generators picks a central charge upon quantization of the theory. One therefore cannot generally conclude that nonzero charge density is necessary in order to have fewer NG bosons than broken generators, unless additional assumptions are made such as that the Lie algebra of symmetry generators is semi-simple, in which case there are no nontrivial central charges (see Sec. 2.7 in [1]).
Looking back at Eq. (33), one observes that in case that the left hand side is actually nonzero, there must be a NG mode that couples to both broken generators Q a , Q b [56]. This indicates dispersion relations number of NGBs charge densities [3] [4] [57] [50] at a very elementary level that the opposite to Theorem 4 should also hold: nonzero vacuum expectation value of a commutator of two broken generators implies that the number of NG bosons is smaller than the number of broken generators. However, it would be desirable to prove this statement at a similar level of rigor as Theorem 4 itself.
While the Nielsen-Chadha theorem 3 relates the number of NG bosons to their dispersion relations, Theorem 4 of Schäfer et al. provides a (partial) connection of the number of NG bosons to charge density. (Keeping in mind the example of Section 3 I use the term "charge density" here as a synonym for "expectation value of a commutator of two broken generators.") However, there is also a rather straightforward connection between the latter two, that is, NG dispersion relations and charge density. This was shown by Leutwyler [4] in the framework of low-energy effective field theory. A general argument based just on basic symmetry assumptions and analyticity is presented below. For a summary of all these relations, see Fig. 4.
Consider now two broken generators, Q a and Q b , such that their commutator has a nonzero vacuum expectation value. I will demonstrate that, quite generally, this implies existence of one NG boson with nonlinear (typically quadratic) dispersion relation at low momentum. For the sake of simplicity I will assume rotational invariance. This may seem rather restrictive, note however that even the low-energy behavior of many solid state systems such as ferromagnets or superconductors is, to leading order, described by a rotationally invariant effective field theory. Keeping this in mind, the transition amplitude for the annihilation of the NG boson by the current operator can be parameterized as 0| j µ a (0)|n k = ik µ on F an (|k|) + iδ µ0 G an (|k|) (34) where k µ on is the on-shell wave vector of the one-particle state, that is, k 0 on is given by the dispersion relation, k 0 on = E n,k . Using current conservation, this leads to the following equation for the NG dispersion relation, (E 2 n,k − k 2 )F an + E n,k G an = 0 A few special cases are worth mentioning. When G an = 0 the standard Lorentz invariant dispersion relation is recovered. When the ratio G an /F an is small of order O(|k|) in the limit k → 0, the dispersion relation comes out linear with phase velocity differing from one. Finally, when G an /F an has a nonzero limit as k → 0, the dispersion relation takes the low-momentum form This is the case of most interest since it turns out to be implied by the presence of nonzero charge density.
In order to decide which of the possibilities outlined above is actually realized, one needs more information about the dynamics of the system. To that end, consider the time-ordered currentcurrent correlation function, D µν Let us also assume the following form of the equal-time commutator of two charge density operators, [ j 0 a (x, t), j 0 b (y, t)] = iδ 3 (x − y)C ab (x, t). By this general notation the possibility of a nontrivial central charge is taken into account. Taking the divergence of the correlation function and employing the Källén-Lehmann spectral representation, one arrives at the following Ward identity in momentum space, where k µ = (k 0 , −k). The index ν was set to zero since only then the left-hand side can have nonzero vacuum expectation value due to the assumed rotational invariance. In order that the particle and antiparticle poles on the right-hand side are canceled, they must appear at energies satisfying Eq. (35), which is just another derivation of the NG dispersion relation. Upon cancelation of the poles, the right-hand side of Eq. (37) takes the form 2i Im E 2 n,k F an F * bn + E n,k (F an G * bn + F * bn G an ) + G an G * bn (38) Taking now the limit of zero momentum, the broken symmetry condition expressed in terms of the nonvanishing expectation value of C ab (0) implies a "density rule," In this expression, G an stands for the zero-momentum limit of the amplitude function defined by Eq. (34). This is the desired result. It means that when the commutator of Q a and Q b has nonzero vacuum expectation value, the amplitudes G a , G b have finite limits as k → 0. According to Eq. (36) the dispersion relation of the NG boson that couples to both broken generators is then nonlinear. In fact, in all cases I am aware of it is quadratic. However, this does not strictly speaking follow from the general argument presented here.

Linear sigma model
The above arguments connecting the number of NG bosons, their dispersion relations and charge densities are general, yet not as conclusive as one might desire. This last part of the section will therefore be devoted to a discussion that is less general, but allows much stronger statements. For many purposes it is often sufficient to introduce an effective scalar field whose vacuum expectation value plays the role of an order parameter. This was actually the underlying motivation behind Section 4 where the general problem of minimization of potential for such a field was considered. Let me first recall one of the standard proofs of the Goldstone theorem that applies to Lorentz invariant theories [36].
Consider a set of real scalar fields φ i , transforming linearly in a real representation of a symmetry group with generators (T a ) i j , supposed to be pure imaginary and Hermitian. This means in practice that the symmetry group is realized by orthogonal transformations, which can always be ensured as long as the group is compact (see Chap. 7 of [18]). The restriction to real fields is not essential for every complex field can be decomposed into two real ones, and I do not make any assumptions about (ir)reducibility of the group representation.
For linearly realized symmetry, the quantum effective potential V eff (φ) of the theory is invariant under the same transformations as the classical Lagrangian (see Sec. 16.4 in [15]). The condition of invariance therefore reads for all generators T a . Differentiating with respect to φ and setting it to its vacuum expectation value ϕ, satisfying ∂V eff (ϕ)/∂φ i = 0, one obtains The matrix of second partial derivatives is the mass matrix of the theory and as long as Lorentz invariance is preserved, it completely determines the spectrum of the system. This is the Goldstone theorem at work: spontaneous breaking of the generator T a by definition means that it does not annihilate the vacuum, that is, (T a ) i j ϕ j 0. But then the real vector χ ai = i(T a ) i j ϕ j is a zero mode of the mass matrix, hence corresponding to a NG boson. How many NG bosons are there? It seems almost obvious that there is one for each broken generator, but let us be a bit more precise at this point. The vectors χ a give rise to a matrix of scalar products, h ab = χ T a χ b = ϕ T T a T b ϕ. This matrix is real, symmetric and positive semi-definite, and may be diagonalized by an orthogonal transformation, h → RhR T . This in turn corresponds to the change of basis of the Lie algebra of generators, T a → R ab T b . After the transformation, the nonzero vectors χ a , equal in number to the number of broken generators, form an orthogonal system. Since they are all zero modes of the mass matrix, it is now obvious that there is exactly one NG mode for each broken generator.
In order to see what difference lack of Lorentz invariance brings in, I will now switch to a class of theories with "minimal" breaking of Lorentz symmetry, induced by chemical potential. These describe many-body systems with relativistic dynamics. Lorentz invariance is broken explicitly by medium effects. The main message will be that the number of zero modes of the mass matrix is still equal to the number of broken generators, but in general no longer coincides with the number of NG modes. Consider a general model for a complex field φ, defined by the Lagrangian [56] where V(φ) is a renormalizable potential with terms up to fourth order in the field φ. The covariant derivative reads D µ φ = (∂ µ −iA µ )φ and makes the Lagrangian formally gauge invariant (see Sec. 2.4 in [58]). In the end the background gauge field is set to A µ = ( a µ a T a , 0), where µ a are chemical potentials associated with a set of mutually commuting conserved charges represented by matrices T a . Upon expanding the covariant derivatives, the Lagrangian becomes with the modified potentialṼ The chemical potential thus gives a negative contribution to the mass matrix, and when it exceeds certain critical value, the perturbative vacuum ϕ = 0 will no longer be stable and the field will condense. As soon as the field develops nonzero expectation value, some of the symmetry generators become spontaneously broken. Keeping in mind that the zero modes of the mass matrix are iT a ϕ, one reparameterizes the field as Here Π(x) is a linear combination of broken generators and represents the NG fields, while H(x) stands for the massive, "Higgs" modes. After some manipulation, the bilinear part of the Lagrangian which governs the spectrum at tree level, becomes whereṼ bilin is the bilinear part of the modified potential, which only depends on the Higgs field H. This is the manifestation of the NG nature of the field Π: it drops from the mass part of the Lagrangian. The Noether currents in the model are j µ a = −i(D µ φ † T a φ − φ † T a D µ φ), and the charge densities in the ground state read j 0 a = 2ϕ † A 0 T a ϕ. Consequently, the last term in the Lagrangian (45) is proportional to the density of a commutator of two generators. Not surprisingly, it plays a crucial role in the structure of the NG spectrum, in accordance with Theorem 4. In order to understand the spectrum of the system governed by the Lagrangian (45), it is desirable to simplify the problem as much as possible; it looks like one has to deal with a very complicated mixing of all the modes. Fortunately, it turns out that when the fields H, Π are properly decomposed into irreducible representations of the unbroken subgroup, the whole Lagrangian splits into sectors in which just two modes can mix with each other at a time (see [56] for details). A key role in the argument is played by the statement, whose validity is not restricted to the framework of the linear sigma model, formulated here as a simple theorem: Theorem 5. It is possible to choose a basis of the symmetry generators in such a way that only mutually commuting charges have nonzero vacuum expectation value.
Unfortunately, I am not aware of a completely general proof of this claim. However, for the wide class of systems whose symmetry group is the special unitary group and its products, the proof is simple. Denote the expectations of the charge operators as q a = 0|Q a |0 . Forming a linear combination with the generator matrices T a one gets a Hermitian matrix q a T a . Under the unitary group transformation U on the Hilbert space, the charge operators transform in the adjoint representation, Q a → UQ a U † . This in turn translates into the transformation rule for the charge expectations q a . However, under the adjoint action of the unitary group, the matrix q a T a can always be diagonalized. In such a diagonal representation, only mutually commuting charges have nonzero expectation values, as was to be proved.
The search for charges with nonzero expectation value is greatly facilitated by the fact that they must transform as singlets under the unbroken symmetry. The candidates can therefore be found by a decomposition of the Lie algebra of generators into irreducible representations of the unbroken subgroup. Once all charges that have nonzero vacuum expectation value have been identified, they can be completed to form the Cartan subalgebra of the Lie algebra. The standard root decomposition of Lie algebras (see Chaps. 6 and 8 of [59]) then guarantees that the remaining generators split up into pairs whose commutators lie in the Cartan subalgebra. Only such pairs can be mixed by the last term in the Lagrangian (45). The task to determine the NG spectrum thus boils down to the problem of mixing of two fields, captured by the Lagrangian As the notation suggests, π is the NG field while H is a Higgs field as long as its mass function f (µ) is nonzero. The coefficient g(µ) accounts for the last two terms in Eq. (45), which have just a single time derivative. The spectrum of this Lagrangian consists of two modes with squared energies at low momentum given by the following formulas, If f (µ) 0, the Lagrangian indeed describes a massive state and a type-I NG boson with linear dispersion relation at low momentum. On the other hand, if f (µ) = 0, that is, if both π and H have no mass term, the spectrum contains a massive mode and only one NG boson of type-II with quadratic dispersion relation at low momentum, E 2,k = k 2 / |g(µ)|. This is a sheer consequence of the term in the Lagrangian with a single time derivative, and in turn of the nonzero expectation value of a commutator of two broken generators. One can thus conclude with the strong statement that there is one type-II NG boson with quadratic dispersion relation for each pair of generators whose commutator has nonzero vacuum expectation value. Explicit examples, further elucidating this result, will be analyzed in the next section.

Further examples
In order to shed some light on the general arguments given in the preceding section, I will now work out in detail some examples of systems featuring type-II NG bosons. Interestingly enough, these examples turn out to be nontrivial and thus hopefully contribute to deeper understanding of the underlying physics. I will start with a nonrelativistic, completely solvable toy model.

Nonrelativistic Boulware-Gilbert model
Following [2], consider a nonrelativistic version of the model introduced in [38], defined in finite volume Ω by the Hamiltonian The canonical coordinate and momentum satisfy the usual commutation relation, [φ(x), π(y)] = iδ 3 (x − y). Since the Hamiltonian is bilinear, this is essentially a free field theory. The theory can be solved by switching to Fourier transformed variables, φ k = Ω d 3 x φ(x)e −ik·x , and analogously for π k and V k . The Hamiltonian thus becomes From here it is already straightforward to derive the Lagrangian, Inverting the kernel of this bilinear Lagrangian one obtains the time-ordered Green's function (propagator) of the theory, where E 2 k = k 2 (1 + V k ) is the dispersion relation of the one-particle excitation of the model. Apparently, when the potential V(x) falls off sufficiently fast at large distance so that a finite limit lim k→0 V k exists, the theory has a gapless excitation. This is a NG boson associated with the invariance of the Lagrangian (50) under the constant shift, φ(x) → φ(x) + θ. The associated current is j = −∇φ and charge density = π. The integral charge then is Q Ω = Ω d 3 x (x) = π 0 , which indeed commutes with the Hamiltonian.
Let us now have a closer look at the NG dispersion relation in the case of a long-range interaction. Assume that the potential is radially symmetric, regular at the origin and drops as 1/r α at large distance r, considering the following approximation, for r r 0 γ/r α for r r 0 (52) The momentum integration in where k = |k|, can then be split into three regions, 0 < x < kr 0 , kr 0 < x < 1, and 1 < x < ∞, which are well separated provided the momentum is small enough. In the integration over the first region we approximate the potential by V 0 and the sine by its argument, getting 4 3 πV 0 r 3 0 . In the 29 second region we approximate the potential by γ/r α and the sine by its argument, which yields 4πγ k α−3 − r 3−α 0 /(3 − α). In the last region we approximate just the potential, leading to the integral 4πγk α−3 ∞ 1 dx sin x/x α−1 . We do not need to evaluate this integral explicitly, it is enough to know that it converges as long as α > 1. Putting all the pieces together, we arrive at the asymptotic scaling of the energy at low momentum, For α ≤ 1 the Fourier transform of the potential does not exist and the theory is ill-defined in three spatial dimensions. The conclusion is that for α > 3 the range of the interaction is short enough to ensure applicability of the Nielsen-Chadha theorem 3; the energy is analytic in momentum. On the other hand, for 1 < α ≤ 3 the interaction is long-ranged, but there is still a NG boson. It is, however, a very nontrivial one. The nonanalytic structure of the dispersion relation clearly precludes a low-energy description in terms of an effective Lagrangian. Expansion in powers of momentum is impossible in such a case.
Another interesting special case is the screened Coulomb interaction, V(r) = γe −µr /4πr [2]. Recalling that V k = γ/(k 2 + µ 2 ), one observes that the system can be stable even with an attractive interaction, provided it is not too strong so that V 0 ≥ −1. In the limiting case, γ = −µ 2 , we find Energy is quadratic in momentum and the NG boson is type-II! Does it mean that we have found an explicit example of a system exhibiting a strict inequality in the Nielsen-Chadha theorem 3? Not really, because there are in fact two broken conserved charges. The reason is that for V 0 = −1, the operator φ 0 also commutes with the Hamiltonian (49). There are thus two conserved charges, π 0 and φ 0 , both being spontaneously broken. Moreover, their commutator is proportional to the unit operator, and hence obviously has nonzero vacuum expectation value in accord with the general arguments of the previous section. Yet another nontrivial fact is that the new conserved charge is not of the Noether type; it is not associated with a symmetry of the Lagrangian. This emphasizes the fact that the Goldstone theorem assumes just the existence of a conserved charge, and it is not important where this charge comes from.

Heisenberg ferromagnet
The next example of SSB is already quite realistic and represents a model for isotropic ferromagnets [14]. The degrees of freedom are spin-1 2 particles fixed to nodes of a crystal lattice, represented by operators s i . (The Roman subscript labels the lattice sites.) The Hamiltonian of the model reads It is invariant under simultaneous rotations of all the spins, forming the group SU (2). As long as all the couplings J i j are positive, the ground state of the Hamiltonian is obviously formed by all spins aligned in the same direction. Formally, one rewrites s i · s j = 1 2 (s i + s j ) 2 − 3 4 . With all spins aligned, each pair resides in a state with total spin one so that the ground state energy is − 1 Any state of the ith particle may be labeled by a direction n as |i, n , meaning that it is an eigenstate of the projection n · s i with the eigenvalue 1/2. The degenerate ground states of the Hamiltonian (55) are then |0 n = N i=1 |i, n , where N is the total number of lattice sites. Using an explicit expression for the eigenvectors of spin-1 2 operators, one finds the scalar product of two states associated with different directions, | n 1 |n 2 | = cos(θ n 1 ,n 2 /2), where θ n 1 ,n 2 is the angle between the two directions. The scalar product of two ground states of the Heisenberg model therefore is and it apparently goes to zero as N → ∞ unless n 1 and n 2 are identical. This is in agreement with the discussion in Section 2 where it was pointed out as one of the characteristic features of SSB.
To construct the whole Hilbert space above the ground state |0 n , one can conveniently use the formalism of creation and annihilation operators. To that end, recall that the two-dimensional space of spin 1 2 may be viewed as the Fock space of the fermionic oscillator. One defines annihilation operators a i (n) and creation operators a † i (n) so that a i (n)|i, n = 0 and {a(n), a † (n)} = 1. These are nothing else than the lowering and raising operators familiar from the theory of angular momentum. In addition to their defining relations above, they satisfy [a i (n), a † i (n)] = 2n · s i . Note that, in this setting, annihilation and creation operators at different lattice sites commute rather than anticommute as usual. The change of sign induced by the interchange of two distinguishable fermions is, however, merely a convention.
The Hilbert space of the Heisenberg ferromagnet is set up as a Fock space above the vacuum |0 n . In the ground state all spins point in the direction n, while the excited states are obtained by the action of the creation operators a † i (n) that flip the spin at the ith lattice site to the opposite direction. The countable basis of the space Hilbert space contains all vectors of the form a † i 1 (n)a † i 2 (n) · · · |0 n where a finite number of spins are flipped. It is now obvious that in the infinite volume limit, all basis vectors from the Hilbert space corresponding to a direction n 1 are orthogonal to all basis vectors from the Hilbert space characterized by a different direction n 2 . This means that these two spaces are completely orthogonal.
The symmetry group SU(2) is generated by the operator of the total spin, S = i s i . Only rotations about the direction n of the ground state, that is, the spontaneous magnetization, are unbroken. The unbroken group U(1) is generated by the spin projection n · S. The two remaining generators of the symmetry group are spontaneously broken. They are obviously ill-defined operators, as are the induced finite transformations, for they take any state out of the above constructed separable Hilbert space.
In order to identify the NG boson state(s), one rewrites the Hamiltonian (55) in terms of the annihilation and creation operators. Since the orientation of the ground state is now fixed, the argument (n) of these operators will be for simplicity omitted. First observe that The Hamiltonian preserves the "particle number," that is, the number of flipped spins generated by the operator i a † i a i . This is of course, up to irrelevant constants, nothing but the n-projected component of the total spin, which is not spontaneously broken and thus can be used to label physical states. Let us restrict our attention to the "one-particle" space, spanned on the basis |i = a † i |0 n . The physical reason behind this restriction is that the sought NG boson turns out to be the spin wave, a traveling perturbation induced by flipping a single spin.
On the one-particle space, the second term on the right hand side of Eq. (57) gives zero unless i = j. The one-particle Hamiltonian thus reads, up to an irrelevant constant, and acts on the basis states as It is worth emphasizing that so far we have nowhere used the discrete translational invariance implied by the symmetries of the crystal lattice. In fact, the indices i, j have been just labels distinguishing different degrees of freedom. However, as we already know, translational invariance is needed in order to have well defined quasiparticle excitations carrying conserved momentum. The stationary states can then be sought as the common eigenstates of the discrete translation operators, that is, the plane waves, |k = i e ik·x i |i . Assuming finally that the couplings depend just on the distance of the lattice sites, J i j = J( x i − x j ), one finds by direct substitution that the plane wave |k indeed is an eigenstate of the one-particle Hamiltonian (59), with the energy given by where J k is the discrete Fourier transform, J k = i J(x i )e −ik·x i . As a concrete example, the nearest neighbor interaction of strength J on a square lattice of spacing would give the dispersion E k = J(3 − cos k x − cos k y − cos k z ) = 2J sin 2 k x 2 + sin 2 k y 2 + sin 2 k z 2 . The dispersion relation of the NG boson is quadratic at low momentum, hence it is of type-II. In fact, this follows directly from the Nielsen-Chadha theorem 3 since there are two broken generators and only one NG state. In this case one can easily see that acting with the two broken generators on the ground state |0 n formally produces states that differ just by a phase factor.
At the end, let us again look more closely at the NG dispersion relation (60) and assume a long-range potential as in Eq. (52). Since we are only interested in the low-momentum behavior, the discrete Fourier transform of the interaction can be replaced with the continuous radial integral (53) with an additional prefactor 1/ 3 taking into account the volume of the elementary lattice cell. In order that the spectrum be well defined at all, the Fourier transform at zero momentum must exist, which is true only for α > 3. By an analysis similar to that for the Boulware-Gilbert model one finds that In conclusion, for α ≥ 5 the interaction drops sufficiently rapidly so that the Nielsen-Chadha theorem applies and the dispersion relation is analytic, concretely quadratic since the two broken generators demand a type-II NG boson. For 3 < α ≤ 5 the NG boson still exists, but its energy is once again nonanalytic in momentum. Remember that even though we have not solved the Heisenberg model (55) completely, Eq. (60) still represents the exact energy of exact eigenstates of the Hamiltonian, and therefore this conclusion about the NG energy is not just an artifact of some approximation. Finally, note that for α ≥ 5 the Fourier transform at low momentum may be evaluated by Taylor expansion up to order k 2 , resulting in the dispersion relation 6.3. Linear sigma model For a final example let me get back to the linear sigma model. Unlike the previous two examples where the NG state and its dispersion relation were constructed exactly, here just the classical approximation will be used. As a trade-off one is able to gain deeper insight in the nature of type-II NG bosons.
For illustration purpose, I will restrict to the model with SU(2) × U(1) symmetry, defined by the Lagrangian [50,51] L Here φ is a complex doublet of scalar fields and the covariant derivative includes chemical potential µ associated with the U(1) factor of the symmetry group, D ν φ = (∂ ν − iδ ν0 µ)φ. For µ > M the perturbative Fock vacuum becomes unstable with respect to fluctuations of the field φ, and the field develops nonzero vacuum expectation value, ϕ. At tree level, The four Noether currents of the theory have the form where τ a are the Pauli matrices. Choosing the orientation of the ground state so that v is real and positive and resides solely in the lower component of ϕ, the only unbroken generator will be τ + = 1 2 (1 1 + τ 3 ). The broken generators can be conveniently chosen as τ − = 1 2 (1 1 − τ 3 ) and τ 1,2 . Then, only the generator τ − has nonzero expectation value in the ground state. Therefore, the commutator [τ 1 , τ 2 ] also has nonzero expectation value and we expect one type-II NG boson which couples to both j 1 ν and j 2 ν . In addition, there should be one type-I NG boson coupled to τ − . Remember from Section 5.4 that the NG fields can be identified as iT a ϕ, using the broken generators T a . This means that the upper component of φ, φ 1 , represents the type-II NG boson while the imaginary part of the lower component, φ 2 , represents the type-I NG boson. The real part of φ 2 corresponds to fluctuations of the magnitude of the condensate, and is expected to be massive. Moreover, φ 1 carries unit charge of the unbroken U(1) symmetry generated by τ + while φ 2 does not. Therefore, these two components do not mix and their bilinear Lagrangians and propagators can be evaluated separately.
One thus finds that φ 1 excites modes with dispersion relations E ∓,k = k 2 + µ 2 ∓ µ. The upper sign corresponds to the type-II NG boson; indeed, at low momentum the energy is approximately E −,k ≈ k 2 /2µ. The lower sign describes a mode with mass 2µ. As predicted in Section 5.4 the presence of nonzero vacuum expectation value of the commutator of two broken generators turns the two associated modes into one massive state and one NG boson of type-II. Similarly, the dispersion relations of the two modes in φ 2 are found to be Again, the upper sign refers to the NG boson, whose dispersion relation at low momentum this time is E −,k ≈ |k| (µ 2 − M 2 )/(3µ 2 − M 2 ). Apart from the dispersion relations, the knowledge of the propagators allows us to determine the couplings of the modes to the fields from the residua of the respective poles using the Källén-Lehmann spectral representation [60]. Denoting the type-I and type-II boson states as |π(k) and |G(k) , one obtains and and finally 0|φ † 1 (0)|G(k) = 0 thanks to the unbroken U(1) charge. These expressions allow one to evaluate, with the help of Eq. (64), the couplings of the NG bosons to the broken currents [57], With all these formulas at hand one readily checks all general formulas derived in Section 5.3 such as the density rule (39). In order to gain further insight into the nature of the type-II NG boson, let us investigate the solution of the classical equation of motion. The bilinear part of the Lagrangian (63) containing just the field φ 1 is up to a total derivative. This has plane wave solutions of the type φ 1 (x) = ae −ik·x with k 0 = k 2 + µ 2 − µ, which is exactly the dispersion relation of the type-II NG boson in the model. Substituting into Eq. (64) one arrives at the expressions for the currents carried by this plane wave, The second line shows that the plane wave is associated with a uniform current of the unbroken charge, τ + . This corresponds to the observation made above that the field φ 1 , and hence its quantum |G(k) , carries unit charge of the unbroken group U(1) . Also, there is uniform background density of the τ − charge, induced by the condensate v. It is independent of the plane wave since the associated current does not couple to |G(k) . The first line of Eq. (70) reveals that the isospin vector of charge densities rotates in the (1, 2) plane perpendicular to the direction of the condensate. In other words, the type-II NG boson is a circularly polarized isospin wave. This is, of course, also in accord with the fact that it carries the unbroken charge. One may then naturally ask where the plane wave with opposite polarization is. In this model, this "antiparticle" of the type-II NG boson is heavy; its mass is 2µ.
Interestingly, the answer to this question is quite different in a nonrelativistic ferromagnet, analyzed in the previous subsection. There, too, the magnon is a spin wave which carries (minus) one unit of the unbroken charge, that is, projection of spin into the direction of total magnetization. However, in this case, the spin wave with opposite circular polarization does not even exist. The explanation is twofold. First, quantum mechanically, the ground state is the state with the highest possible spin, and the magnon corresponds to a propagating perturbation caused by flipping of one of the spins. The spin wave polarized in the opposite way would have to have total spin one unit higher than the ground state, which is not possible. Classically, imagine the magnon as a single test spin propagating on a background of other spins. The Heisenberg Hamiltonian (55) shows that this background acts on the test spin as a homogeneous magnetic field. The evolution of a spin in a magnetic field is well known as the Larmor precession. Its sense is fixed by the magnetic moment of the spin; rotation of the spin in the opposite direction is not possible.

Low-energy effective field theory for NG bosons
Effective field theory (EFT) is a powerful approach to problems involving several energy scales that was developed in the full generality in particle physics [61], but has grown into an indispensable tool in essentially all branches of physics that use the methods of quantum field theory. The basic assumption is that the physical degrees of freedom split into groups well separated in energy. In case one is interested only in low-energy, or long-distance, observables, EFT may be used to describe the system in a way which involves just the low-energy degrees of freedom. The effects of the microscopic, short-distance dynamics are incorporated in the values of coupling constants of the effective theory. For more details see the review papers [5][6][7][8].
There are fortunate cases where the low-energy EFT can be derived directly from the underlying microscopic theory by "integrating out" the heavy, high-energy degrees of freedom. A profound example is the Euler-Heisenberg Lagrangian for the electromagnetic field, which can be inferred from quantum electrodynamics by eliminating electrons or any other massive charged particles present. In such cases EFT still provides a very convenient tool which dispenses with the degrees of freedom irrelevant for the description of low-energy physics. However, the true power of EFT lies in its successful application to systems where the high-energy degrees of freedom cannot be simply integrated out of the microscopic theory, either because of its nonperturbative nature, or because the low-energy modes are not even present in the microscopic theory (such as the hadrons that are not the fundamental degrees of freedom in quantum chromodynamics). The utility of EFT then relies on the observation that field theory is merely a convenient way to incorporate the general principles such as unitarity and cluster decomposition but contains no further dynamical assumptions, and it can therefore reproduce the predictions at low energy of any microscopic theory as long as it involves the appropriate low-energy degrees of freedom and symmetries, see e.g. [9] or Sec. 19.5 in [15].
A distinguished class of systems featuring separation of scales is represented by those with SSB. The NG bosons, guaranteed by the Goldstone theorem, are then the low-energy degrees of freedom. In the most frequent case that there are no other gapless states in the spectrum, they are in fact the only degrees of freedom in the low-energy EFT. A systematic approximation scheme for calculations of low-energy observables is then provided by expansion in powers of momentum, or derivatives. (This requires absence of long-range interactions, as shown in Sections 6.1 and 6.2 on explicit examples with non-analytic low-energy behavior that makes momentum expansion impossible.) For a detailed explanation how to set up a consistent power counting scheme, see [7,9]. The standard procedure to construct invariant Lagrangians for NG bosons [62,63] is briefly reviewed in the following subsection. After then I will point out some specific issues concerning its application to nonrelativistic systems.

Coset construction of effective Lagrangians
As was shown already in Section 5 NG states are created from the vacuum by the action of the broken symmetry generators. Formally, a global broken symmetry transformation, which merely moves the vacuum into an equivalent ground state, may be viewed as a zero-momentum NG boson. The finite-momentum NG state can then be excited by a local infinitesimal broken symmetry transformation. This corresponds to the intuitive picture of NG modes as small local fluctuations of the order parameter. Denoting the vacuum expectation value of the order parameter field φ(x) as ϕ, all possible values of the order parameter that can be reached from ϕ by a symmetry transformation span a manifold, known as the coset space of the broken symmetry, G/H. (Strictly speaking, the coset space may in general consist of several orbits of the symmetry group G. For the construction of the effective Lagrangian, only the orbit containing the vacuum ϕ is important.) The elements of the coset space may be thought of as the sets gH = {gh | h ∈ H} with arbitrary g ∈ G. Thanks to the group properties, two sets g 1 H and g 2 H are either identical or disjoint, and G/H is the space of these classes.
While the ground state is represented by a single point in the coset space, a general NG field configuration can be viewed as a map from the spacetime to the coset space, φ(x) : R 4 → G/H. The task to construct the most general G-invariant Lagrangian for the NG bosons is equivalent to the geometric problem of constructing a G-invariant function on G/H, given the action of the group on the order parameter.
Let the symmetry group G be compact and semi-simple. Then one can choose a basis T a of its Lie algebra such that Tr(T a T b ) = δ ab (see Chap. 1 of [18]). Splitting the generators T a into the unbroken ones, V α , generating the unbroken subgroup H, and the broken ones, A i , it follows that Tr(V α V β ) = δ αβ , Tr(A i A j ) = δ i j , and Tr(V α A i ) = 0. Furthermore, the structure constants f abc of the Lie algebra of G are fully antisymmetric and the broken generators span a representation of the unbroken group, hence the commutation relations Any element g ∈ G can be, at least in the neighborhood of unity, decomposed as g = e iπ i A i e iv α V α . Since e iv α V α leaves the ground state intact, the points of the coset space can be parameterized as φ(π) = U(π)ϕ, where U(π) = e iπ i A i . Thus, the NG fields π i serve as coordinates on G/H. The action of the symmetry group on G/H is defined by left multiplication, φ(π) g − → gφ(π) = φ(π). Decomposing gU(π) as U(π)h(g, π), where h(g, π) = e iv α (g,π)V α belongs to the unbroken subgroup H, one arrives at the transformation rule This defines a nonlinear realization of the group G on G/H. In order that the group structure be preserved, the function h : G × G/H → H must satisfy basic constraints such as h(1 1, π) = 1 1 and the associativity, h(g g, π) = h(g ,π(g, π))h(g, π). Also, for transformations from the unbroken subgroup it is independent of π, h(h, π) = h, so that the unbroken group acts on G/H linearly, The effective Lagrangian is most easily constructed using the Maurer-Cartan form, Employing the very useful formula, one easily realizes that the Maurer-Cartan form lies in the Lie algebra of G. Under the group action (72) it transforms as ω Since h dh † lies in the Lie algebra of H, it is useful to decompose the Maurer-Cartan form into projections on the spaces of unbroken and broken generators, ω = ω + ω ⊥ . The "longitudinal" part transforms as a gauge field, ω g − → hω h † − ih dh † , while the "transverse" part transforms covariantly ω ⊥ → hω ⊥ h † .
Should it be necessary to include other, non-NG modes in the effective theory, one proceeds as follows. Taking advantage of the fact that the unbroken subgroup H is realized linearly even on the NG bosons, consider a set of "matter" fields ψ that transform in the representation R of H, ψ h − → h R ψ. This is promoted to a nonlinear realization of the full group G, defined by Any Lagrangian invariant under H can be made invariant under the full group G using this procedure. Note that since the transformation matrix h R (g, π) depends on the NG fields π and thus implicitly on the spacetime coordinates, the transformation becomes local. In order that terms containing derivatives of ψ be still G-invariant, one uses the longitudinal component of the Maurer-Cartan form, ω = ω µ dx µ , to promote derivatives to covariant ones, ∂ µ ψ → (∂ µ + iω µ )ψ. The most general G-invariant Lagrangian can then be constructed [9] from the objects ω ⊥ , ψ and their covariant derivatives. For instance, the expression obviously includes the kinetic term for the NG bosons, as can be seen with the help of Eq. (74). The f 2 π is a Hermitian, positive definite matrix that must commute with all matrices from H in order that L bilin be G-invariant. If the broken generators transform irreducibly under H, it is proportional to the unit matrix by Schur's lemma, hence providing merely an overall normalization of the Lagrangian. In general it contains one arbitrary parameter per each irreducible representation, thus giving an independent normalization of the kinetic term for each multiplet of NG bosons.
One often meets the situation that the symmetry in question is local rather than global. This can happen both, when the system interacts with a dynamical gauge (such as electromagnetic) field and when external sources are coupled to the conserved currents of the theory [64]. Consider a gauge field A transforming under local G transformations as usual as A g − → gA g † + ig dg † . It is straightforward to see that the Maurer-Cartan form, defined now by iω = U † (π)DU(π) = U † (π)(d−iA )U(π), satisfies the same rule (75) as in the case of a global symmetry.
This concludes the brief review of the coset construction of effective Lagrangians. It is useful to note that while the choice of the coset parameterization as U(π) = e iπ i A i is very convenient in particular because it renders the action of the unbroken subgroup linear, the construction is not limited to this choice. In fact, the key elements are only the geometry of the coset space and its symmetry. The choice of coordinates π i has no observable effects; one speaks of reparameterization invariance of physical observables [62].

Geometric interpretation
The kinetic term (77) can be evaluated more or less explicitly. Using the orthogonality of the generators, the Maurer-Cartan form can be expanded in V α , A i . Denoting ω ⊥ = ω ⊥i A i , one finds The kinetic term acquires the form g i j dπ i dπ j with g i j = Σ ki Σ l j Tr( f 2 π A k A l ). (I deliberately use differentials instead of spacetime derivatives in order to emphasize the geometric nature of the quantities, and to avoid having to treat space and time derivatives separately.) The matrix g i j can be interpreted as a metric on the coset space. The construction of the most general term in the effective Lagrangian with two derivatives therefore reduces to the geometric problem of finding all G-invariant metrics on G/H [4,64].
In order to see what are the coordinate transformations that leave the metric g i j invariant, one must explicitly evaluate the change of the fields π i under a G transformation. Writing Eq. (72) as one sets g = e iθ a T a and expands both sides to first order in the infinitesimal transformation parameters θ a . (No assumption on the smallness of π i is made though.) Denoting further Comparing coefficients at A i , a formula for the shift of NG fields, δπ i , follows, The set of objects ξ i a that depend just on the coset geometry define the Killing vectors on G/H. They represent infinitesimal motions on G/H induced by the symmetry group transformations. A similar expression for the parameters v α of the compensating transformation h(g, π) is implied.
The above formulas were derived for infinitesimal transformation parameters θ a , but hold for arbitrary values of the NG fields π i . In practice, the vacuum is usually set at π i = 0 and for calculation of low-energy observables such as scattering amplitudes, one gets by with a power expansion in π i . All expressions of this subsection can then be evaluated explicitly in terms of the structure constants of the Lie algebra of symmetry generators. To first order in the NG fields one finds, for instance, This essentially confirms the picture that underlies all discussions of effective Lagrangians for NG bosons: unbroken symmetry is realized linearly on the NG fields whereas broken transformations are equivalent, to lowest order, to a mere shift. This makes sure that the NG fields correspond to gapless modes in the spectrum. It is worthwhile to inquire how the Maurer-Cartan form changes in presence of a gauge field. According to our discussion above, it picks an additional contribution, One then immediately finds that the transverse component of the Maurer-Cartan form that determines the kinetic term of the NG bosons, modifies according to ω ⊥i = Σ i j dπ j → Σ i j Dπ j , where Dπ i = dπ i − ξ i a A a . This shows that even in this very general setting where the symmetry group is realized nonlinearly, the infinitesimal shift of the fields under group motion serves to construct their gauge-covariant derivative.
In the above derivation, the metric g i j was derived from the Maurer-Cartan form and its Ginvariance followed from the G-invariance of the Lagrangian. It is instructive to reverse the argument and use an invariant metric on the coset space to construct Eq. (77) in a purely geometric way. To this end, recall that we deal with compact symmetry groups that can be represented by unitary matrices. The space of complex matrices possesses a natural distance function, invariant under both left and right multiplication by unitary matrices, defined by d(x, y) = x − y , where x = Tr(x † x). (I use the term distance instead of metric to distinguish this object from the metric tensor, which roughly speaking measures the distance on the tangent space at unity.) This induces a distance function on the coset G/H, manifestly invariant under the action (72) of G, It is easy to show that this prescription indeed yields a well defined distance. First, the unbroken subgroup H is compact so that the minimum in Eq. (85) exists. Furthermore, d G/H is independent of the choice of the representative elements x, y for the cosets xH, yH and it is zero if and only if x and y lie in the same coset. Finally, it is straightforward to prove the triangle inequality as a consequence of that for the distance d.
To determine the metric tensor following from Eq. (85), one takes the points x, y infinitesimally far apart, x = U(π) and y = U(π + dπ), whence The minimum will be realized for h infinitesimally close to unity, hence 1 1 − h will lie in the Lie algebra of H, and affect only the longitudinal part of the Maurer-Cartan form, ω (π). Clearly, the minimum is achieved when 1 1 − h exactly cancels ω so that Consequently, the bilinear Lagrangian (77) is directly related to a G-invariant metric on the coset space G/H. The only difference is the factor f 2 π present in Eq. (77). Mathematically, this comes from the fact that Eq. (87) represents an intrinsic metric of the coset space, invariant under both left and right G transformations, while Eq. (77) is by construction invariant under left G multiplication, but only under H acting from the right [64].

Nonrelativistic effective Lagrangians
In the preceding two subsections a general method how to build effective Lagrangians was presented. The construction is based on an implicit assumption that the Lagrangian is invariant under the symmetry group of the system. Global symmetry of the theory can be conveniently formulated in terms of invariance of the generating functional of Green's functions of the conserved currents, j µ a (x), with respect to local transformations of external sources, A a µ (x), coupled to these currents [65,66]. This implies a set of Ward identities that the theory must satisfy in order to preserve the symmetry. In Lorentz invariant theories without quantum anomalies, the most general solution of the Ward identities can indeed always be obtained from a Lagrangian invariant under the symmetry group [64]. On the other hand, in systems that lack Lorentz invariance ("nonrelativistic" systems) this conclusion is no longer valid. The most general solution of the Ward identities to the lowest order of the derivative expansion was given by Leutwyler [4] who showed that in general it is just the action, not the Lagrangian, what is invariant.
The Lagrangian of a low-energy EFT in principle contains an infinite number of terms. The predictive power of EFT lies in the fact that at any desired accuracy of an observable to be calculated, only a finite number of terms contribute. The Lagrangian is systematically expanded in powers of space (s) and time (t) derivatives, L = s,t L (s,t) . In order to achieve a consistent power counting scheme, it is convenient to assign the external source A a µ the same order as the derivative ∂ µ . However, the scaling of spatial and temporal derivatives may in general differ. For the sake of simplicity, I will assume that the system is invariant under the continuous rotation group SO(3). This is certainly plausible for fluids, and even for some crystalline (such the plain cubic) structures, where the anisotropy induced by the lattice only appears at higher orders in the derivative expansion. The first few terms in the effective Lagrangian then read, where the dot denotes a time derivative, and the covariant derivative, D µ π i = ∂ µ π i − ξ i a A a µ , was introduced below Eq. (84). Both g i j (π) andḡ i j (π) are G-invariant metrics on the coset space. It follows from the above discussion that if the NG fields transform irreducibly under the unbroken symmetry, they must be identical functions on G/H up to a trivial scale factor. Both L (0,2) and L (2,0) are obviously invariant under simultaneous gauge transformations of the NG fields and the external sources.
It would be tempting to conclude that the coefficient functions C i (π) and a (π) are related by a = −C i ξ i a so that L (0,1) = C i (π)D 0 π i . However, this does not follow from general symmetry considerations. Instead, the Ward identities encoded in the gauge invariance of the generating functional imply a set of constraints for all coefficients C i (π), a (π), g i j (π),ḡ i j (π), ξ i a (π) [4]. The functions ξ i a (π), defined by Eq. (82), are Killing vectors of both invariant metrics, g i j (π) andḡ i j (π), that is, they obey the Killing equation ∇ i ξ ja + ∇ j ξ ia = 0. The Riemannian covariant derivative ∇ i is defined with the help of the Christoffel symbol, ∇ i ξ ja = ∂ i ξ ja − Γ k i j ξ ka . (Following [4] I use the notation ∂ i ≡ ∂/∂π i , not to be mixed up with a derivative with respect to the space coordinate.) Note that using the two metrics to lower the index in ξ i a , one obtains two different objects, ξ ia and ξ ia , that satisfy Killing equations with respective Christoffel symbols.
Since the Killing vectors ξ i a generate the symmetry transformation (82) on the coset space, they must reproduce the structure of the Lie algebra of the symmetry generators, d a ξ Finally, the coefficients C i (π) and a (π) satisfy the identities While the former asserts that a (π) transforms in the adjoint representation of the symmetry group and fixes its value at any point of G/H in terms of that at the origin (π = 0), the latter determines the function C i (π) up to a gradient. This ambiguity is intrinsic to EFT and has no observable consequences. Indeed, changing C i (π) by a gradient merely modifies the Lagrangian by a total time derivative, hence leaving the action intact. Several remarks are in order here. First, L (0,1) is in general not gauge invariant: under the transformation (82) it changes by a total derivative, δL (0,1) = d dt θ a (C i ξ i a + a ) . Second, evaluating the action at its extremum in order to establish the leading contribution to the generating functional, one finds by differentiation with respect to A a 0 that a (π) at the origin determines the density of the conserved charge, 0| j 0 a (x)|0 = a (0). When all charge densities are zero, it follows from Eq. (89) that both a (π) and C i (π) vanish identically and the lowest order Lagrangian is essentially the same as in Lorentz invariant theories, up to a possible redefinition of the phase velocity of the NG bosons. The energy of the NG bosons is linear in momentum so that one assigns space and time derivatives the same degree in the power expansion, and the leading order Lagrangian reads L = L (0,2) + L (2,0) . Also, when the symmetry is Abelian Eq. (89) guarantees that a (π) = a (0) on the whole coset space, and in turn C i (π) = 0. Since constant a (π) does not contribute to the dynamics of NG bosons, merely generating background charge density, the equation of motion and dispersion relation of the NG bosons will be the same as above.
On the other hand, nonzero density of a non-Abelian charge leads to nonzero C i (π). The Lagrangian is gauge invariant only up to a total derivative, and the term with the single time derivative leads to a quadratic dispersion relation of the NG boson. One must then count a time derivative as two space derivatives so that the lowest order Lagrangian is L = L (0,1) + L (2,0) . This confirms the observation that density of a non-Abelian charge implies a type-II NG boson in the spectrum, made in Section 5 on a general ground.
It is instructive to illustrate the construction of effective Lagrangians on an explicit example. The most convenient one is provided by the (anti)ferromagnet. Both these systems possess an SU(2) symmetry group of spin rotations, which is spontaneously broken by an ordered ground state to its U(1) subgroup. As already mentioned in Section 5.2 in a ferromagnet all spins sitting on the crystal lattice are aligned in the same direction and there is a net total magnetization. Accordingly, the two broken generators give rise to one type-II NG boson with quadratic dispersion relation. On the contrary, in an antiferromagnet antiparallel alignment of neighboring spins is favored and the ground state, though ordered, does not develop nonzero spin density. There are two NG bosons with linear dispersion relations, corresponding to two linearly polarized spin waves.
The coset space is SU(2)/U(1) SO(3)/SO(2) in both cases and it is equivalent to a twosphere, S 2 . It can therefore be most conveniently represented by a unit vector, n. The symmetry group SU(2) acts on n through rotations, which can be parameterized by a vector g whose direction and magnitude represent the axis and angle of the rotation, respectively. An infinitesimal symmetry transformation therefore reads δn = n × g, and the covariant derivative in presence of an external source A µ consequently D µ n = ∂ µ n + A µ × n. Since all charge densities are zero in the antiferromagnet, the term L (0,1) is missing in the Lagrangian. The leading order Lagrangian thus is Obviously, for an arbitrary set of coordinates π 1 , π 2 on the coset, the invariant metrics will be g i j = F 2 s ∂ i n·∂ j n andḡ i j = F 2 t ∂ i n·∂ j n. From the Lagrangian one can also extract the Killing vectors, h i = −F 2 s n × ∂ i n. Indices are raised with the inverse metric, g i j = ε ik ε j ∂ k n · ∂ n/F 2 s ∂ 1 n × ∂ 2 n 2 . The case of the ferromagnet is more subtle. The ground state has nonzero charge density and accordingly the term C i (π) appears in the Lagrangian. First observe that the zero component of the external field, A a 0 , couples to the Noether charge density, which is proportional to the total magnetization. Hence A a 0 is interpreted as an external magnetic field. The direction of the magnetization in the ground state is aligned with it so that the vector of couplings a (π) is simply proportional to the order parameter n(π) up to a constant, = αn. By explicitly solving the constraint (89), the function C i (π) can be shown to have a topological nature, related to the Brouwer degree of the map n(π) [4]. Constructing a pathπ i (π, λ), 0 ≤ λ ≤ 1 such thatπ i (π, 0) = 0 andπ i (π, 1) = π i , and setting n = n π(π, λ) , one has C i (π) = α 1 0 dλ (∂ i n × ∂ λ n) · n (91) The integral is independent of the choice of the path, connecting the origin with the point π i . This 42 concludes the construction of the leading order effective Lagrangian for the ferromagnet. Details can be found in [4]. One comment concerning the difference between the antiferromagnet and the ferromagnet is worthwhile here. Due to the presence of a term with a single time derivative, the effective Lagrangian of the ferromagnet is sometimes claimed to break time reversal invariance. However, this is somewhat misleading. What breaks the time reversal is the ferromagnetic ground state, thanks to the nonzero charge density, that is, a zero component of a conserved current. On the other hand, any microscopic description of a ferromagnet will be time reversal invariant. Indeed, for example the Hamiltonian (55) can describe both the ferromagnetic and the antiferromagnetic state, depending on the sign of the coupling constants J i j . Both systems therefore must have the same symmetry. Since one of the basic requirements on the low-energy EFT is that it should reproduce the symmetry of the underlying theory, the effective Lagrangian of the ferromagnet must actually inherit this time reversal invariance. What really distinguishes ferromagnets from antiferromagnets is the way time reversal is implemented on the coset space, which is natural because this implementation does depend on the structure of the ground state. The transformation of the NG fields under time reversal in the antiferromagnet is such that it prohibits the term with a single time derivative, whereas in the ferromagnet it is allowed [67].

Applications of effective field theory
Effective field theory has become an invaluable tool of a theoretical physicist in all branches of physics where the methods of quantum field theory are used. It often facilitates calculations that would otherwise be very complicated, or hard to carry out at all. It would be hopeless to even try to give a representative list of all its applications. I will therefore just mentioned briefly a few explicit examples that are related to other material presented in this paper. The application of EFT to ferromagnets was worked out in [4,[67][68][69], and in a similar spirit also for antiferromagnets [70][71][72][73]. The utility of EFT of course goes beyond the mere derivation of the equation of motion and dispersion relation of spin waves. One can calculate other observables as well such as the cross section for scattering of magnons, or of slow neutrons off the magnons [9]. The use of EFT as well as other techniques for the nonrelativistic weakly interacting Bose gas and Bose-Einstein condensation is reviewed in [74]. See also [75] for the calculation of the damping rate of the corresponding NG boson and the construction of the effective Lagrangian based on Galilean invariance. Subtle features associated with a consistent definition of thermodynamics in various approximation schemes for Bose condensed systems are discussed in detail in [76]. Extension of Galilean invariance in nonrelativistic physics to general coordinate invariance was worked out in [77] and applied to strongly attractive Fermi gas at unitarity.
Among the varied applications of EFT, I would like to discuss in a little more detail the technique to construct the effective action for relativistic superfluids from the equation of state, developed by Son [78]. It is based on the observation that in a relativistic many-body system at zero temperature, the chemical potential is the sole source of Lorentz violation. Lorentz invariance can therefore be restored when the chemical potential is treated as a zero component of a background gauge field, very much like in Section 5.4 where the linear sigma model was investigated. Consider a relativistic system with Abelian U(1) symmetry that acts on the order parameter φ as a phase transformation: φ → φe iθ . Assigning a chemical potential µ to the Noether charge of this symmetry and denoting as in Section 5.4 A µ = (µ, 0), the quantum effective action of the system, Γ eff [A µ , φ], will be invariant under a simultaneous local transformation of φ and A µ , A µ → A µ + ∂ µ θ. As soon as the order parameter develops nonzero vacuum expectation value, ϕ, it can be conveniently parameterized as φ(x) = e iπ(x) [ϕ + H(x)]. The phase transformation of φ is equivalent to a shift of the NG field, π → π + θ.
At low energy the modes created by the "Higgs" field cannot be excited. In order to obtain an effective action for the NG boson only, one "integrates out" the Higgs mode by minimizing the effective action with respect to the modulus |φ|. The leading term in the derivative expansion of the resulting effective action, Γ eff [A µ , π], depends just on A µ and ∂ µ π. The underlying gauge invariance makes sure that they appear in the effective Lagrangian only in the combination D µ π = ∂ µ π − A µ . At the minimum of the effective action, the effective Lagrangian is equal to minus the energy density of the ground state in the background field A µ (see Sec. 16.3 in [15]). Recalling that the energy density is equal to minus the thermodynamic pressure, p(µ), and observing that µ = A µ A µ = D µ πD µ π for a constant field π, the full dependence on D µ π can be restored with the knowledge of the function p(µ), that is, the equation of state, Γ eff [µ, π] = d 4 x p( D µ πD µ π) = d 4 x p (∂ 0 π − µ) 2 − (∇π) 2 (92) One therefore arrives at the conclusion that to leading order in the derivative expansion, the full quantum effective action of the superfluid NG boson is completely determined by the equation of state. Since the equation of state is usually straightforward to calculate perturbatively, this is a rare example of a system where the low-energy EFT can be derived explicitly from the underlying microscopic theory. The effective Lagrangian can be used to derive transport coefficients and other observables associated with low-energy processes in which only the NG bosons can be excited [79].

Conclusions
In the present paper, I have provided a self-contained introduction to spontaneous symmetry breaking and the properties of the associated NG bosons. I discussed to some detail the mathematical subtleties, connected with the implementation of a broken symmetry on the Hilbert space of states. The focus of the paper was on nontrivial features of SSB specific to nonrelativistic systems. These include unconventional dispersion relations of NG bosons as well as improved rules for their counting. The general results were demonstrated on several, mostly exactly solvable, examples.
Due to limited space, several issues were touched only superficially, or not at all. The first of them is the possibility of explicit symmetry breaking. If weak enough, spontaneous breaking of such symmetry will give rise to a pseudo-NG boson. Its energy at zero momentum is not strictly zero, but it is still small compared to the characteristic energy scale of the system. However, the counting of NG bosons then becomes blurred, since it may not be possible to distinguish an "approximately linear" dispersion relation from an "approximately quadratic" one. Only when the explicit symmetry breaking is weak, one can still distinguish type-I and type-II NG bosons [57].
A special case of explicit symmetry breaking is represented by a wide class of theories where a part of the Lagrangian, responsible for the NG boson spectrum at the tree level or in a similar approximation, has a higher symmetry than the rest of it. This typically happens in relativistic scalar theories coupled to gauge fields [80,81], but it can also occur in some condensed matter systems such as superfluid 3 He [82,83]. In the appropriate approximation, the system then exhibits more gapless states than would correspond to the symmetry of the full theory, some of them stemming from the extended symmetry of the part of the Lagrangian. These spurious NG bosons receive a gap once quantum corrections are taken into account. This mechanism can be responsible for the presence of naturally light states in the spectrum.
Second, in this paper I have considered exclusively continuous internal symmetries. While spontaneous breaking of a discrete symmetry does not give rise to NG bosons, spontaneous breaking of spacetime symmetries is subtle. First of all, the NG field configurations can be imagined as small local fluctuations of the order parameter. Since the local versions of different spacetime transformations may coincide, there are typically fewer NG bosons than broken generators [37]. Examples of spontaneous breaking of spacetime symmetries include homogeneous but anisotropic states such as in superfluid 3 He (see [25] for an extensive review), spin-one color superconductors [32], or even imbalanced spin-zero color superconductors [84], and in Bose-Einstein condensates of relativistic vector fields [85,86], as well as inhomogeneous states such as crystalline solids [87] or superconductors with inhomogeneous pairing [88]. Accordingly, the behavior of NG bosons may be highly nontrivial. For example, in helical ferromagnets [89,90] the local magnetization field forms a spiral structure. Along the axis of the helix, the average magnetization is zero and the NG boson is type-I. In the transverse directions, the NG boson feels the uniform magnetization and is type-II like in ordinary ferromagnets.
My main motivation for writing this paper was to provide a review on spontaneous symmetry breaking that would be general enough to cover both, relativistic many-body systems as well as intrinsically nonrelativistic systems, and thus to bridge the gap between the communities. If a reader with expertise in high energy physics discovers that there is more to SSB and NG bosons than usually presented in the textbooks on particle physics, and if a reader with background in condensed matter or atomic physics finds here a general framework for the variety of fascinating phenomena he or she is familiar with, my goal will be achieved.