Arxiv:Hep-Ph/0607102V2 22 Sep 2006

Home , Goldstone boson

arXiv:hep-ph/0607102v2 22 Sep 2006 8 oBs–isencnesdaoi ae 9 0 11] 10, [9, gases ferromagnet atomic nuclear condensed relativistic the Bose–Einstein and to 7] [8] 6, 5, of [4, phases high-density QCD various however, from have pro- ranging decade, bosons studied, is Goldstone past been type-II happens the with systems In this other pro- several where ferromagnet. A the system be by generators. a vided may broken of bosons of example number Goldstone found the of than number smaller the boson, stone fboe ymtygnrtr ytefloiginequal- number following the the by to ity, generators related Goldstone symmetry are the broken types of of respective numbers the The type-II, of and respectively. bosons odd, even, is is power it are this if if bosons type-I, Goldstone as proportional The classified momentum. limit, then of long-wavelength power the some to in symmetry Gold- spontaneous is, the the breaking of from stemming energy under bosons the that, stone manner assumptions, [3]: Chadha nonperturbative technical and a certain Nielsen in by showed achieved result They was fundamental The respect this intricate. in quite be may Lorentz-noninvariantsystems in breaking symmetry spontaneous transport in viscosity. or role conductivity major heat a as play such modes phenomena soft par- the these In (and ticular, system. dynamics the low-energy of Goldstone thermodynamics) the the low-temperature for of crucial properties and are number bosons the- The Goldstone 2]. celebrated [1, gap- the orem of by guaranteed existence modes, the less is breaking symmetry spontaneous ∗ lcrncades [email protected] address: Electronic q 1 hw htweee hr satp-IGold- type-II a is there whenever that shows (1) Eq. yeIGs+2 + GBs type-I # of physics the that time long a for known been has It of consequences general and striking most the of One ewrs iersgamdl pnaeu ymtybreaki symmetry Spontaneous model, sigma Linear Keywords: i 11.30.Qc determine numbers: and PACS allowe momentum, is high — at phonon relation superfluid dispersion the its conden — of Bose–Einstein One-loo boson the Goldstone in type-I calculated matter. appears the are that quark propagator boson w its dense Goldstone as which II of well phase potential, as potential CFL chemical effective the finite in at condensation symmetry internal U(1) eivsiaesotnossmer raigwti h li the within breaking symmetry spontaneous investigate We .INTRODUCTION I. eateto hoeia hsc,NcerPyisInstit Physics Nuclear Physics, Theoretical of Department pnaeu ymtybekn ntelna im model sigma linear the in breaking symmetry Spontaneous × yeI GBs type-II # tfiieceia oeta:Oelo corrections One-loop potential: chemical finite at ≥ rkngenerators broken # Tomáˇs Brauner . (1) n h eea ffc frltvsi etrcondensation vector relativistic 13]. of [12, effect general the and is,i sol tteoelo ee htthe that level one-loop the at pa- necessary. reasons only corrections present several is loop are The it the There First, of inclusion gap. level. the this tree make fill that the to intended at is time per the all statements. worked these of we proof model, detailed sigma ‘acci- a linear give an the to it able Within type-I were making a one. thus of type-II vanish, dental’ velocity may phase boson the the like Goldstone cases where exceptional transitions, to up phase saturated, is (1) law with equality boson dispersion Goldstone quadratic type-II a one implies generators ken count- general rule: a conjectured ing den- and nonzero charges of Noether Gold- of possibility previous of sities the connection on with the counting Based clarified sigma boson we stone 15] linear potential. [5, Lorentz- results chemical relativistic of partial finite class the at particular — model a systems in noninvariant counting boson stone nodrt hc h auaino h Nielsen–Chadha the of necessary saturation the therefore [3]. check inequality are to corrections order quantum in correctly The qualitatively) (even sufficient[16]. spectrum not is the level determine tree the to that symme- Lagrangian, higher full known has the well Lagrangian than is the try classical of it part the static Third, the where when trusted. transition, be phase cannot in the analysis especially of role, the vicinity significant Second, a phase, merely the play potential. normal it may field corrections the scalar phase loop the in condensed stabilize spectrum Bose–Einstein to the serves the on in effect level, while tree no the at has Indeed, it role. nontrivial a plays teraction hs fdneQD[,5.I sas pca aeof case special a also is It 5]. [4, QCD Color-Flavor-Locked dense the of in phase global condensation the kaon with describe model par- sigma a linear analyze SU(2) the we — Instead, model [14]. ticular Ref. of pattern breaking u nlssws oee,prl lsia we — classical purely however, was, analysis Our Gold- the investigated we [14] paper recent our In o ipiiy ed o olwtegnrlsymmetry- general the follow not do we simplicity, For atclratnini adt h type- the to paid is attention Particular . ∗ g yeI odtn bosons. Goldstone Type-II ng, × e hs.Frhroe eso that show we Furthermore, phase. sed sdcywidth. decay ts ()smer.Ti oe a enue to used been has model This symmetry. U(1) o-eodniyo omttro w bro- two of commutator a of density Non-zero t,25068 ute, ersgamdlwt h SU(2) the with model sigma near ssgetda oe o kaon for model a as suggested as odcydet h nonlinearity the to due decay to d orcin otesaa field scalar the to corrections p ez zc Republic Reˇz, Czech ˇ sacneune h in- the consequence, a As . × λ ( φ † φ ) 2 in- 2 a class of models of relativistic Bose–Einstein condensa- are two Goldstone bosons, one type-I and one type-II. tion investigated by Andersen [17]. Their low-energy dispersion relations read The plan of the paper is as follows. Using the manifestly invariant formalism of generating functionals, we µ2 M 2 p2 determine in Sec. II the one-loop effective potential and E = − p and E = , (3) s3µ2 M 2 | | 2µ numerically find its minimum as a function of the chem- − ical potential. In Sec. III, we calculate the one-loop cor- respectively. In addition to the Goldstone bosons, rection to the scalar field propagator. We demonstrate there are two massive excitations with energy gaps analytically the presence of gapless (Goldstone) poles and 2(3µ2 M 2) and 2µ. then compute numerically the corrections to the tree- − level dispersion relations of the type-II Goldstone boson p and its massive counterpart. In Sec. IV, we calculate the decay rate of the superfluid phonon [18]. The paper A. One-loop correction is concluded with a summary and a discussion of the results. Some necessary technical details are deferred to Let us now proceed to the one-loop calculation. In the appendices. order to account for the breaking of the U(1) symmetry associated with the particle number (or strangeness, in the context of kaon condensation) we introduce the II. EFFECTIVE POTENTIAL formal Nambu doublet

φi We work with the model of Miransky and Shovkovy [4] Φi = † . and Schaefer et al. [5] which is defined by the Lagrangian φi [The index i refers to the doublet representation of the SU(2).] The one-loop effective action is then given by the = D φ†Dµφ M 2φ†φ λ(φ†φ)2. (2) textbook formula [20], L µ − − i Here φ is a complex doublet of the global SU(2) sym- Γ [φ, φ†]= S[φ, φ†]+ log det ∆−1 (4) metry and the chemical potential µ, associated with the 1L 2 global U(1) symmetry (particle number), is included in the covariant derivative, D φ = (∂ iδ µ)φ [19]. This plus counterterms. Here S[φ, φ†]= d4x (x) is the clas- µ µ µ0 L model has been analyzed, at the tree− level, in detail in sical action and ∆ is the tree-level matrix propagator R Refs. [4, 5], so we only summarize the main results for defined by later reference. 2 When µ>M the static potential develops a nontrivial −1 δ S ∆ij (x, y)= † . minimum and the scalar field condenses. As a result, δΦi (x)δΦj (y) the SU(2) U(1) symmetry of the Lagrangian (2) breaks down to its× U(1) subgroup [different from the original Displaying explicitly the matrix structure in the Nambu U(1)]. Associated with the three broken generators there space, it reads, in the momentum representation,

(p + µ)2 p2 M 2 2λφ†φ δ 2λφ†φ 2λφ φ ∆−1(p)= 0 − − − ij − j i − i j . (5) ij 2λφ†φ† (p µ)2 p2 M 2 2λφ†φ δ 2λφ†φ − i j 0 − − − − ij − i j !

Here the classical field φ is already assumed to be con- In this form, the effective potential is manifestly stant. There is no lack of generality in this requirement SU(2) U(1) invariant. For detailed calculations it is, as long as the vacuum is translationally invariant. however,× more convenient to fix the direction of φ in the With the constant classical field φ we may evaluate the SU(2) ‘flavor’ space. As usual, we set one-loop effective potential as 1 0 4 φ = . (7) i d k −1 √2 v V1L = Vcl log det ∆ (k) (6) − 2 (2π)4 Z With this choice, ∆−1 becomes diagonal in the flavor plus counterterms, where V = (M 2 µ2)φ†φ + λ(φ†φ)2. space and may easily be inverted, yielding cl − 3

−1 (p + µ)2 p2 M 2 λv2 0 ∆ (p)= 0 , 11 − − − 2 2 2 2 −1 0 (p0 µ) p M λv ! − − − − (8) 2 2 2 2 2 1 (p0 µ) p M 2λv λv ∆22(p)= 2 2 2 2 − − 2− − 2 2 2 2 , [p E (p)][p E (p)] λv (p0 + µ) p M 2λv 0 − + 0 − − − − − where allows to change the sign of the integration variable so that the contribution of the particles circulating in the 2 2 2 2 2 E±(p)= p + µ + M +2λv loop is the same as that of antiparticles. For the later calculation of the scalar field propagator it will also be 4µ2(p2 µ2 + M 2 + λv2)+(2µ2 + λv2)2. ± − useful to be able to shift the integration variable. In concrete computations we shall use the analytic integration By minimizingp the classical potential V we find that, at cl over energy (frequency) in combination with dimensional tree level, λv2 = µ2 M 2, which yields the dispersion regularization around d = 3 space dimensions [17]. relations − Eq. (10) may be rewritten in a particularly useful way. E2 = p2 +3µ2 M 2 4µ2 p2 + (3µ2 M 2)2, (9) With a decent use of Eq. (5) we find (for the sake of − ± − brevity, we temporarily drop out the arguments of the in accord with the result ofp Refs. [4, 5]. propagators as well as the integration measure)

∂ ∂ † † † † det ∆−1 = ∆−1φφ ∆−1φ φ ∆−1φφ∆−1φ φ B. Renormalization ∂v2 22 ∂v2 22 22 − 22 22 −1φφ† −1φ†φ −1φφ = 2λ ∆22 + ∆22 ∆22 The loop integral in Eq. (6) is divergent. To renormal- − − ize it, we add to the eﬀective potential the counterterm so that

2 † † 2 ∂ −1 Vcounterterm = δv + δM φ φ + δλ(φ φ) . 2 det ∆ † † ∂v 22 =2λ ∆φ φ + ∆φφ + ∆φφ − det ∆−1 22 22 22 This term generates, by means of its functional deriva- Z 22 Z † tives, also the counterterm for the scalar propagator. The =2λ 2∆φφ + ∆φφ . constant δv is irrelevant for our purposes as it just fixes 22 22 Z the offset of the effective potential at φ = 0. The counterterms in the effective potential are sufficient to absorb Eq. (10) thus acquires the form all infinities for, as is well known, there is no wave func- 4 µ2 (M 2 + δM 2) (λ + δλ)v2 tion renormalization in the scalar φ theory at one loop − − 4 [27]. d k † † =2iλ ∆φφ (k)+2∆φφ (k) + ∆φφ(k) (11) Actually, we shall not need to evaluate the effective (2π)4 11 22 22 potential at all. What we are interested in is just the Z h i vacuum expectation value of the scalar field, or v. We Even though this result has been derived by a not very shall be therefore solving the ‘gap equation’ ∂V1L/∂v = transparent manipulation, it could have been expected: 0. Now since the effective potential is a function of v2, It is equivalent, up to an overall factor, to the require- there is always the solution v = 0. In the Bose–Einstein ment that the one-particle-irreducible tadpole contribu- condensed phase this is, however, not the only solution tions to the vacuum expectation value of φ2 vanish. For and not the global minimum. The nontrivial solution is the diagrammatic representation of the individual terms 2 found by means of ∂V1L/∂v = 0. Formulas (6) and (8) see Appendix B. give For the concrete calculations, we adopted the following renormalization conditions. The scalar mass was renor- ∂V 1 1 malized by subtracting the whole one-loop correction at 0= 1L = (M 2 + δM 2 µ2)+ (λ + δλ)v2 ∂v2 2 − 2 µ = 0 and v = 0 [see also Eq. (15) later], i.e., ∂ −1 4 † i d k φφ ∂v2 det ∆22 (k) d + 2λ∆ (k) . (10) 2 3−d d k 1 2 (2π)4 11 − −1 δM = 3λΛ , (12) " det ∆22 (k) # − (2π)d √k2 + M 2 Z Z The momentum integral in Eq. (10) is only formal and where the integral is regularized by minimal subtraction implicitly involves a convenient regularization. In deriv- in d dimensions and Λ is the renormalization scale. This ing Eq. (10) we already assumed that the regularization procedure guarantees that the parameter M keeps its 4

interpretation as the physical (pole) mass of the scalar ¼º¼½ at zero chemical potential. The coupling constant λ was renormalized by modiﬁed minimal subtraction (MS) so that ¼º¼¼8

3λ2 2 ½ ¼º¼¼6 δλ = 2 γ + log 4π , 4π ǫ −

¼ ¼º¼¼4 where ǫ =3 d is the expansion parameter of dimensional Ú =Ú regularization− [28]. (The pole does not depend on the chemical potential and hence could be extracted from the analytically calculable loop integral at µ = 0.) The ¼º¼¼¾ renormalization scale Λ was set equal to the characteristic

scale of the system: we chose Λ = M in the normal phase ¼

and Λ = µ in the Bose–Einstein condensed phase. Since

½ ¾ ¿ 4 5 the phase transition occurs at µ = M (see the following discussion), the subtraction point is a connected function =Å of the chemical potential. For sake of numerical solution the gap equation (11) was rewritten, with the help of Eq. (8) and upon the FIG. 1: Relative increase of the vacuum expectation value µ2−M2 0 integration over frequencies, as v with respect to the tree-level value v = λ . The numerical data were obtained with the renormalizedq coupling µ2 M 2 d3k 1 λ set to 1. v2 = − λ − (2π)3 √k2 + M 2 + λv2 Z 2 k2 + M 2 + 3 λv2 µ2 + 1+ 2 − III. SCALAR PROPAGATOR E (k)+ E (k) E (k)E (k) + − + − δM 2 + δλ v2 The one-loop-corrected propagator is obtained as . (13) D − λ a second functional derivative of the one-loop effective action (4), This equation was solved iteratively with the initial 2 2 µ −M 2 ansatz given by the tree-level value, v0 = . To δ Γ λ −1(x, y)= 1L . q ij † implement the MS scheme within the numerical compu- D δΦi (x)δΦj (y) tation, we subtracted the pole part in the form After differentiating the logarithm of the determinant we 3λ2 Λ2 ddk 1 δλ = log +3λ2Λ3−d . arrive at −4π2 M 2 (2π)d (k2 + M 2)3/2 Z −1(x, y) = ∆−1(x, y) The integral here reproduces the pole and is, in fact, Dij ij equal to the full one-loop correction to λ at µ = 0 and v = i δ2∆−1 δ∆−1 δ∆−1 0. Upon the subtraction of δM 2 and δλ, the integral in + Tr ∆ ∆ ∆ 2 † − † δΦ (y) (13) is rendered finite and was evaluated approximately " δΦi (x)δΦj (y) δΦi (x) j # with a simple cutoff. The numerical results are displayed in Fig. 1. We em- plus counterterms. Here ‘Tr’ denotes the full functional phasize the fact that the phase transition to the Bose– trace, over both spacetime and internal degrees of free- Einstein condensed phase happens at the chemical poten- dom. −1 tial equal to the renormalized scalar mass. This is of The derivatives of the inverse bare propagator ∆ are course not surprising from the physical point of view and nothing but the cubic and the quartic interaction vertex. in fact is an important check of consistency of our calcu- They are local and may be represented by the coupling (i) (ij) lations. matrices T and Q , defined as This conclusion is also easily proved analytically. Just −1 set v = 0 and then immediately observe that the right δ∆ (s,t) (i) kl = T δ4(s t)δ4(s x) hand side of Eq. (11) is exactly canceled by the coun- † kl δΦi (x) − − terterm δM 2 i.e., v = 0 is a solution of Eq. (11) exactly δ2∆−1(s,t) for µ = M. kl (ij) 4 4 4 † = Qkl δ (s t)δ (s x)δ (t y). Finally, note that the convergence of our iterative so- δΦi (x)δΦj (y) − − − lution of the gap equation guarantees that the found solution is, as required, a (at least local) minimum of the The explicit formulas for the matrices T (i) and Q(ij) are effective potential. listed in the Appendix A. With this notation, the inverse 5 propagator acquires the form Note that the loop correction in Eq. (15) is momentum-independent, as usual in a φ4 theory. The −1 −1 renormalization is therefore trivial and corresponds to ij (p) = ∆ij (p) D a mere redefinition of the mass. i d4k (ij) (i) (j)† Since in dimensional regularization we may freely shift + 4 tr ∆(k)Q ∆(k)T ∆(p + k)T 2 (2π) − the integration variable, the loop integral in Eq. (15) Z h i (14) is µ-independent and thus is completely canceled by the mass counterterm (12). Therefore, in the normal phase, plus counterterms. The symbol ‘tr’ here refers to a trace the scalar field propagator gets no correction and the over internal degrees of freedom only, that is, over the mass spectrum simply consists of two doubly degenerate flavor and Nambu spaces. levels at M µ [4, 5]. ±

A. Normal phase B. Bose–Einstein condensed phase When the chemical potential is small enough, the system is in the normal phase: v = 0 and the symmetry When the scalar condenses i.e., v > 0, we have to is not spontaneously broken. In this phase the one-loop work with the full formulas (8). The matrix propagator has altogether 16 components (four flavor times four correction to the propagator is particularly simple. The D bare propagator ∆ is diagonal in both the flavor and the Nambu). Fortunately, several of them are actually zero Nambu space, the only nonzero components being due to the clever choice of the vacuum, Eq. (7). With this choice, the unbroken flavor U(1)Q symmetry is generated φφ† 2 2 2 −1 1 ∆ (p) = [(p + µ) p M ] , by the matrix Q = (11+ τ3). This means that the upper ij 0 − − 2 † flavor component of φ, φ1, carries the unbroken charge ∆φ φ(p) = [(p µ)2 p2 M 2]−1. ij 0 − − − while the lower component φ2 does not. The charge conservation then immediately implies that the propagator Since in the normal phase there is no cubic coupling is diagonal in the flavor space, i.e., 12(p)= 21(p) = 0. T (i), the final formula for the one-loop propagator is eas- D D Moreover, the propagator of φ1 is diagonal in the Nambu ily found to be φφ φ†φ† space, that is, 11 (p) = 11 (p) = 0. [All these relations could, of course,D alsoD be demonstrated explicitly by −1 (p)= (p + µ)2 p2 (M 2 + δM 2) δ Dijφφ† 0 − − ij a proper analysis of the formula (14).] d4k 1 The remaining nonzero components of the propagator 6iλδij 4 2 2 2 . (15) − (2π) (k0 + µ) k M are strongly constrained by Eq. (16) so that only three Z − − φφ† φφ† φφ of them are independent: 11 , 22 and 22 . The other nonzero component of the propagator, φ†φ, D D D is related to Eq. (15) by the first of the identities, D

1. Propagator of φ1 φ†φ(p)= φφ† ( p), D D∗ − ∗ (16) † † (p)= (p)= ( p). Dφ φ Dφφ Dφφ − φφ† The only independent component of 11 is 11 . which follow from the very deﬁnition of the propagator. A straightforward application of Eq. (14) yieldsD D

† −1φφ (p) = (p + µ)2 p2 (M 2 + δM 2) (λ + δλ)v2 D11 0 − − − 4 4 d k † † d k † † † † † 2iλ 2∆φφ (k) + ∆φφ (k) 2iλ2v2 ∆φφ (p + k) ∆φφ (k) + ∆φφ(k) + ∆φ φ (k) + ∆φ φ(k) . − (2π)4 11 22 − (2π)4 11 22 22 22 22 Z h i Z h i(17)

The ﬁrst integral represents the tadpoles and is the propagator is given in Appendix B. momentum-independent. A qualitative diﬀerence in † The propagator φφ (p) is expected to have a mass- comparison with the normal phase comes in the second D11 integral, which contains the contribution of the cubic ver- less pole corresponding to the Goldstone boson. Indeed, a short glance at Eq. (8) reveals that at tree level, when tices. The diagrammatic representation of all terms in † λv2 = µ2 M 2, there is a gapless pole in ∆φφ (p). The − 11 6

field φ1 annihilates a particle with dispersion relation remaining momentum integral. E = p2 + µ2 µ, i.e., a type-II Goldstone boson. The At the leading order of the momentum power expan- − † sion, the inverse propagator reads otherp excitation, annihilated by φ1, with a tree-level gap φ†φ 2µ is manifested as a pole in 11 (p). † D −1φφ (p)=2µp (1 + Z ) p2(1 + Z ). The analytic proof of the existence of a massless pole D11 0 1 − 2 φφ† in 11 (p) is provided in Appendix C 1. Here we just D −1φφ† The coefficients Z1,Z2 determine the renormalization of observe that once the relation 11 (0) = 0 is proved, both the wave function and the dispersion relation of −1φφ†D we may subtract it from 11 (p) to obtain a convenient the Goldstone boson. Their dependence on the chemical expression D potential is displayed in Fig. 2. Now we turn our attention to the massive partner of † −1φφ (p)= p2 +2µp p2 the type-II Goldstone boson. Analogously to Ref. [5], D11 0 0 − 4 we shall calculate just the leading order of its disper- d k † † 2iλ2v2 ∆φφ (p + k) ∆φφ (k) sion relation, i.e., its gap. The corresponding propa- 4 11 11 † − (2π) − gator, φ φ, is obtained from Eq. (18) via the rela- † Z h † † † i 11 φφ φφ φ φ φ φ D p 0 ∆22 (k) + ∆22 (k) + ∆22 (k) + ∆22 (k) . (18) tions (16). Next we set = . However, the explicit × form of the tree-level propagators, Eq. (8), reveals that h i φφ† φ†φ Note that the propagator no longer depends explicitly on ∆11 (k0 2µ, k) = ∆11 (k0, k) so that the difference of the renormalization counterterms. Also, the momentum- − φφ† propagators in the loop correction to ( p0, 0) ex- independent tadpole graphs disappeared and the remain- D11 − actly vanishes at p0 =2µ, see Eq. (18). As a result, the ing integral is finite and perfectly well defined even when renormalized gap is equal to its tree-level value 2µ and the regularization is removed. receives no radiative corrections. In order to establish the radiative corrections to the type-II Goldstone boson dispersion relation, we per- formed the integration over frequency and expanded the result in powers of the external momentum p, up to or- 2. Propagator of φ2 2 der p0 and p , respectively. The final formula is rather cumbersome, so we do not write it out and instead re- In this case, Eq. (14) yields, after a reasonable appli- port just the result of the numerical integration of the cation of the rules (16),

4 † d k † † −1φφ (p) = (p + µ)2 p2 (M 2 + δM 2) 2(λ + δλ)v2 2iλ ∆φφ (k)+2∆φφ (k) D22 0 − − − − (2π)4 11 22 Z 4 h i d k † † 2iλ2v2 ∆φφ (k)∆φφ (p + k) − (2π)4 11 11 Z 4 d k † † † † 1 † † † 8iλ2v2 ∆φφ (p + k) ∆φφ (k) + ∆φφ(k) + ∆φ φ (k)+ ∆φ φ(k) + ∆φφ(p + k)∆φ φ (k) , (19) − (2π)4 22 22 22 22 2 22 22 22 Z

4 4 d k d k † † −1φφ(p)= (λ + δλ)v2 2iλ ∆φφ(k) 2iλ2v2 ∆φφ (k)∆φφ (p + k) D22 − − (2π)4 22 − (2π)4 11 11 Z Z 4 d k † 1 † † † † † 8iλ2v2 ∆φφ(p + k) ∆φφ (k)+ ∆φφ(k) + ∆φ φ (k) + ∆φ φ(k) + ∆φφ (p + k)∆φφ (k) . (20) − (2π)4 22 22 2 22 22 22 22 22 Z

Here the propagator is not diagonal in the Nambu Fortunately, we need not calculate the determinant ex- space and the situation is thus more complicated. We plicitly, at least to prove the existence of a Goldstone cannot look for the zeros of the individual components of boson. Recalling the relations (16), the determinant at the inverse propagator. Instead, we have to solve the con- zero momentum simpliﬁes to −1 dition det 22 (p) = 0, which seems a rather formidable task in viewD of Eqs. (19) and (20).

† 2 2 det −1(0) = −1φφ (0) −1φφ(0) . D22 D22 − D22 h i

¼º¼¼5 strongly aﬀects the transport properties of the system Z ½ such as the shear viscosity [22]. Here we calculate the two-particle decay rate of the ¼º¼¼4 phonon at the leading order of perturbation theory. Such a quantity could in principle be extracted, via the opti-

¼º¼¼¿ cal theorem, from the imaginary part of the one-loop Z ¾ propagator calculated in Sec. III B 2. This would, how-

ever, be rather complicated for two reasons. First, the ¼º¼¼¾ phonon is, in the propagator of φ2, entangled with the massive ‘radial’ mode. Second, we would have to deal

with the full momentum dependence of the propagator; ¼º¼¼½ the low-momentum limit studied before would not be suf- ﬁcient for our purpose. Since the upcoming calculation

¼ of the phonon decay rate involves a few rather nonstan-

¾ ¿ 4 5 ½ dard pieces, we regard it instructive to perform it step

by step, using ordinary perturbation theory. =Å Before starting the calculation we brieﬂy comment why we have chosen to study such a particular process. First, FIG. 2: Renormalization constants of the type-II Goldstone as has already been mentioned above, the low-energy boson propagator. The numerical results indicate that their properties of the Goldstone bosons are essential for the dependence on the chemical potential is linear near the phase thermodynamics of the system at low temperature. Two transition and saturates at large µ. The data were obtained factors are determinant for the possibility of the decay of with λ = 1. the Goldstone bosons: The conservation of the unbroken U(1)Q charge and phase space (kinematic) restrictions. The type-II Goldstone boson, G, carries the unbroken The secular equation thus reduces to charge. Since at low energy it can only decay to gapless † particles, the only processes allowed by charge conser- −1φφ (0) = −1φφ(0) . D22 ± D22 vation are G Gπ, G Gππ, etc. (Here and in the following, π denotes→ the→ type-I Goldstone boson i.e., the

In view of Eq. (5) we see that the existence of a Goldstone phonon.) Taking into account the dispersion relations boson is proved once we verify the relation (3), these decay modes, however, fail to satisfy the en- † ergy and momentum conservation simultaneously. Thus, −1φφ (0) = −1φφ(0). (21) D22 D22 the type-II Goldstone boson is stable at low momentum. A detailed proof of this statement is provided in Ap- On the other hand, the phonon may decay, by the processes π ππ, π πππ etc. Only the two processes ex- pendix C 2. → → The analysis of Eqs. (19) and (20) would yield the ra- plicitly written here receive nonzero contribution at the diative correction to the type-I Goldstone boson disper- first order of perturbation theory. As will become clear sion relation, i.e., the speed of sound in the Bose–Einstein later, the three-particle decay is suppressed, at least in condensate. We do not perform the detailed calculation certain ranges of momentum and chemical potential. We here since our main concern was the type-II Goldstone shall therefore calculate just the decay rate for the process π ππ and comment on the relevance of our result boson studied above. Moreover, such a straightforward → approach is unnecessarily complicated because the one- afterwards. loop correction to the speed of sound has been determined by Andersen using a different method [17]. In the following, we shall instead concentrate on another spe- A. Vacuum transition amplitudes cific property of the superfluid phonon in our system — its decay. The decay rate cannot be straightforwardly deduced from the Lagrangian (2) because it is expressed in terms of the scalar field φ rather than the physical modes π, G IV. PHONON DECAY and their massive counterparts. The problem is that the kinetic term in the Lagrangian cannot be diagonalized at Upon inspection of Eq. (9), the full tree-level dis- the level of fields. Indeed, Eq. (8) shows that the di- persion relation of the type-I Goldstone boson (or, the agonalizing unitary transformation is necessarily energy- phonon) turns out to be convex. This fact has far- dependent. We therefore have to calculate the three- reaching consequences for it kinematically allows the point Green’s function of φ and use the LSZ reduction phonon to decay. (A similar, convex, dispersion law is formula to extract the physical scattering amplitude. found e.g. for the superfluid phonon in the Bose–Einstein Alternatively, we may directly calculate the S-matrix; condensed weakly interacting Bose gas [21].) This in turn the couplings of φ to the physical states follow from 8 the Källén–Lehmann representation of the propagator. where E(p) is the energy of the Hamiltonian eigenstate Let A(x) and B(x) be two local bosonic operators. n, p . The summation index n is discrete for one- The general form of the Källén–Lehmann representa- |particlei intermediate states, and continuous for multi- tion of the time-ordered Green’s function, AB(x, y) = particle states. Also, the states are assumed to be nor- D 3 i 0 T A(x)B(y) 0 , reads, in the momentum represen- malized as n, p m, k = δmnδ (p k). Note that we use tation,− h | { }| i the generalh Källén–Lehmann| i representation− (22) rather than the covariant form usual in literature [20]. 0 A(0) n, p n, p B(0) 0 AB(p,ω)=(2π)3 h | | ih | | i D ω E(p)+ iǫ n X − The poles in the propagator of φ correspond to the 0 B(0) n, p n, p A(0) 0 2 h | | − ih − | | i , (22) dispersion relations given in Eq. (9). The contribution − ω + E(p) iǫ of the massless pole at ω = E (p) iǫ is − − −

1 1 1 ∆pole(p,ω)= 22 2E (p) E2 (p) E2 (p) ω E (p)+ iǫ − − − + − − 2 2 2 2 2 [E−(p) µ] p M 2λv λv − − 2 − − 2 2 2 2 . × λv [E−(p)+ µ] p M 2λv − − − This expression immediately yields the transition amplitudes

1 1 1 1/2 p p 2 p2 2 2 0 φ2(0) π( ) = 3/2 [E−( ) µ] + + M +2λv , h | | i (2π) 2E−(p) 2 p 2 p − − E+( ) E−( )n o − (23) 1 p 1 q 1 1/2 † p p 2 p2 2 2 0 φ2(0) π( ) = 3/2 [E−( )+ µ] + + M +2λv . h | | i −(2π) 2E−(p) E2 (p) E2 (p) − + − − n o p q

The relative sign of the amplitudes follows from the non- where the scattering amplitude is given by diagonal part of the propagator. Note that the ampli- M tudes are particularly simple in the limit v 0 (i.e., at 1 √ 2 → 2 i 3 = 2 2iλv the phase boundary) where E±(p) = p + M µ = − M 3/2 − ± (2π) ǫp µ. It follows that, in this limit, ± p π(k) φ† 0 π(q) φ† 0 0 φ π(p) × h | 2| ih | 2| ih | 2| i 1 0 φ (0) π(p) = , 0 φ† (0) π(p) =0. +h π(k) φ† 0 π(q) φ 0 0 φ† π(p) 2 3/2 2 2 2 2 h | | i (2π) 2ǫp h | | i h | | ih | | ih | | i + π(k) φ 0 π(q) φ† 0 0 φ† π(p) (24) h | 2| ih | 2| ih | 2| i Thus the result typical forp Lorentz-invariant theories is † + π(k) φ2 0 π(q) φ2 0 0 φ2 π(p) recovered in the normal phase. h | | ih | | ih | | i + π(k) φ 0 π(q) φ† 0 0 φ π(p) At the first order of perturbation theory, the ampli- h | 2| ih | 2| ih | 2| i tude for the two-particle decay of π follows from the La- † + π(k) φ2 0 π(q) φ2 0 0 φ2 π(p) . (25) grangian (2) upon shifting φ2 by its vacuum expectation h | | ih | | ih | | i √ i value, v/ 2. The relevant cubic terms read For legibility we omitted the arguments ‘(0)’ of the field operators. Also, the initial state is labeled by the mo- √2λv(φ† φ φ† + φ† φ φ ). mentum p, the momenta of the products are k, q. − 2 2 2 2 2 2 These two terms generate three diagrams each, contribut- ing to the given decay process, that differ just by a per- B. Kinematics mutation of the lines. The result for the S-matrix ele- ment is Our analysis of the two-particle decay is further complicated by the fact that the Lorentz invariance is broken 4 4 1 by the chemical potential. We therefore have to evalu- S = (2π) δ (p k q)i 3 , − − − M (2π)3/2 ate the decay rate as a function of the momentum of the 9

decaying particle. After the common trick of squaring ¼º¼4

the S-matrix in ﬁnite volume and performing the trivial 4¼ three-momentum integration we arrive at the formula for ¼º¼¿

the decay rate, ¿¼ ¾¼ 1 1 ¼º¼¾ Γ= k2d3k 2δ E (p) E (k) E (q) .

2 4π2 |M| − − − − − ½¼ Z ¼º¼½ 1 The factor 2 arises from the fact that the two particles in the ﬁnal state are identical. The δ-function ¼

ﬁxes the moduli of the momenta k and q as a func- ¼

tion of the scattering angle θ, which measures the de- ¹¼º¼½ ﬂection of k from the direction of p. The q is given by 2 2 | |

q = p + k 2 p k cos θ and k is then ﬁxed by ¹¼º¼¾

| | − | || | | |

¼º¾ ¼º4 ¼º6 ¼º8 ½ energy conservation. Due to the form of the dispersion, ¼

Eq. (9),p the solution of the energy conservation condition

kj=Å is, however, rather involved, so we work it out analyti- j cally only in a special case (see Sec. IV C). In the general case, the decay rate is computed numerically. FIG. 3: Sample plot of the expression −E−(p)+ E−(k)+ Using the chain rule to differentiate E−(q) in the δ- E−(q) (in units of M) as a function of |k|, for µ = 2M and function, we arrive at the final formula for the decay rate, |p| = M. The curves are labeled by the value of the angle θ in degrees. The critical angle for the chosen set of parameters θm 1 is θ ≈ 26.3 deg. Γ= sin θ dθ k2(θ) m 4π 0 Z 2 |M| , (26) C. Decay rate × ′ ′ |k(θ)|−|p| cos θ E− k(θ) + E− q(θ) |q(θ)| The phonon decay rate may be determined analytically where θm is the maximum angle to which π can de- in the limit v 0. Strictly speaking, the following cal- cay at a given energy. Fig. 3 displays the expression culation is justified→ only for p2 λv2. However, later ≫ E−(p)+ E−(k)+ E−(q) as a function of k , for sev- the analytic result will be checked by a direct numerical − | | eral values of the angle θ. The critical angle θm is ap- integration of the full formula (26). For the time being, parently determined by the condition that the deriva- we simply set v = 0 everywhere except, of course, the tive of the argument of the energy-conserving δ-function overall factor v in the amplitude . [i.e., the expression in the denominator in Eq. (26)] Eq. (24) shows that only the firstM term survives in Eq. vanishes at k = 0. Since for k = 0 both energy and (25) so that momentum conservation require q = p, this reduces to ′ ′ 1 1 1 cos θm = E−(0)/E−(p). Using explicitly the dispersion =2√2λv . relation (9), we have M 2ǫp √2ǫk 2ǫq

Setting µ = M in Eq. (9)p gives E±(pp)= ǫp M, hence 2 2 E−(p) µ M ′ ± E−(p)= p /ǫp and Eq. (26) reduces to cos θm = 2− 2 | | p s3µ M | | − −1 λ2v2 π/2 2µ2 Γ= sin θ dθ 1 . (27) 4πǫp 0 × " − 4µ2p2 + (3µ2 M 2)2 # Z − k2(θ) . (28) p × k(θ) (ǫ + ǫ ) p ǫ cos θ Formula (27) admits several particularly simple lim- | | k(θ) q(θ) − | | k(θ) its. For p 0 we find θ = 0. This is easily un- | | → m With the simplified dispersion relation, the energy and derstandable: At low momentum the phonon dispersion momentum conservation may easily be solved to yield relation is linear and the decay thus kinematically for- 2 2 bidden. For p , on the other hand, we have 2M p E+(p)cos θ 2Mp cos θ | |2 →2 ∞ k = | | , E−(k)= . µ −M 2 p p2 2 2 p p2 2 θ = arccos 2 2 . Finally, for µ M+ at fixed | | E+( ) cos θ E+( ) cos θ m 3µ −M → − − momentum i.e.,q at the phase transition, θm = π/2. This The angular integration in Eq. (28) is then elementary shows that one has to be careful when performing these and we arrive at the final formula for the decay rate, limits. Near the phase transition there may be two small scales, p and µ M, and the result depends on their λ2v2 p | | − Γ= | | . (29) hierarchy. 4π ǫp(ǫp + M) 10

In Sec. IIIB1, we determined the radiative corrections

:¼¼½

¼ to the dispersion relations of the type-II Goldstone bo-

:¼¼¿

¼ son and its massive counterpart. The numerical results

¼º¼¾

:¼½ ¼ suggest that the dependence of the renormalization con- ¾

stants Z1 and Z2 on the chemical potential saturates at

:¼¿ ¼ large µ. Since in this region the two massive modes in

Å =Ú the model become very heavy, the asymptotic behavior

of the type-II Goldstone boson dispersion relation should

¼:½ ¼º¼½ be governed by a low-energy eﬀective theory for the Gold-

stone bosons, presumably analogous to that for the non-

:¿ ¼ relativistic ferromagnet [15, 23]. Furthermore, we showed that the massive counterpart ½ of the type-II Goldstone boson has gap 2µ, which re-

¼ ceives no radiative corrections. This assertion has a sim-

½ ¾ ¿ 4 5

¼ ple physical interpretation. Recall that µ is, originally,

Ôj=Å j a chemical potential associated with the U(1) symmetry of the Lagrangian (2). Once this symmetry is spontaneously broken, µ may be reinterpreted as a chemical FIG. 4: Decay rate for the process π → ππ as a function of potential of the unbroken U(1) . Since φ and φ† carry momentum (solid lines). The numerical data were obtained Q 1 1 the charges Q = 1, the corresponding excitations differ with λ = 1. The curves are labeled by the dimensionless pa- ± rameter µ/M −1. The dashed line corresponds to the analytic by 2µ in energy. result, Eq. (29). The decay rates for different values of the In a sense, this is a remnant of the Goldstone theo- chemical potential are scaled by the factor λv2 so that their rem: It guarantees nonrenormalization of the Goldstone convergence to the analytic result is comparable. boson mass to all orders in the loop expansion. Once the chemical potential is switched on, the excitation annihi- † lated by φ1 can no longer be massless, but the Goldstone Such a dependence of the decay rate on the momentum theorem still assures that its gap does not renormalize. of the initial phonon might have been expected: The de- In Sec. IV, we determined the two-particle decay rate cay rate tends to zero at both small and large momenta, of the superfluid phonon. Such a quantity as well as the where the dispersion relation becomes linear. A compar- very possibility of the phonon decay may be important ison with the exact numerical integration of Eq. (26) is for the transport phenomena in the Color-Flavor-Locked displayed in Fig. 4. quark matter with a kaon condensate [22]. As shown in Refs. [24, 25, 26], the CFL phase with a meson conden- V. CONCLUSIONS sate appears to be energetically preferred to the simple CFL phase, thus making our results of a possible relevance for the phenomenology of compact stars. In this paper we studied spontaneous symmetry break- To summarize, in this paper we investigated sponta- ing by Bose–Einstein condensation in scalar φ4 theory neous symmetry breaking within the linear sigma model with two charged scalar flavors. Our primary goal was at the one-loop level. We verified the saturation of the to refine our previous results [14] by including radiative Nielsen–Chadha inequality (1) and hence also the count- corrections. ing rule proposed (and proved at the tree level) in our pre- Within the one-loop calculation we verified that the vious work [14]. We thus fulfilled one part of the program Bose–Einstein condensation sets on at the chemical po- outlined in Ref. [14]. The other, more ambitious part is tential equal to the renormalized scalar mass (at µ = 0). to extend the results achieved so far to other relativistic This assertion, though perfectly physical, may seem sur- systems with finite chemical potential. Our recent work prising within the approximation used since, as is well shows that at least the type-II Goldstone bosons may be known [20], in massless φ4 theory the one-loop correction described, to a large extent, in a model-independent way. draws the scalar field expectation value off the origin to These results will be reported elsewhere. a nonzero value. The spontaneous symmetry breaking thus achieved, however, eventually turns out to be an ar- tifact of the one-loop effective potential and disappears upon the renormalization group improvement. Fortunately, such a fake vacuum apparently does not APPENDIX A: CUBIC AND QUARTIC COUPLING MATRICES occur in our case for the phase transition point is exactly where it should be. This gives us hope that the renormalization group improvement, however important for In this appendix the explicit form of the coupling ma- critical phenomena, should not change the conclusions of trices T (i) and Q(ij) is given. The indices φ and φ† refer this paper qualitatively. to the upper and lower components of the Nambu dou- 11 blet, respectively. We start with the interaction vertices. These all come from the λ(φ†φ)2 term in the Lagrangian. Upon shift- 2λ(δ φ + δ φ ) 0 − (i) il k kl i ing φ2 by v/√2, the interaction yields the following ver- Tklφ = − † † , 2λ(δilφ + δikφ ) 2λ(δikφl + δklφi) tices, − k l − † † (i) 2λ(δikφl + δklφi ) 2λ(δilφk + δikφl) T † = − − , klφ 0 2λ(δ φ† + δ φ†) − il k kl i = 4iλ, = 2iλ, (ij) 2λ(δjkδil + δklδij ) 0 − − Q † = − , klφφ 0 2λ(δikδjl + δklδij ) − 0 0 Q(ij) = , klφφ 2λ(δ δ + δ δ ) 0 − ik jl jk il = = √2iλv, − (ij) 0 2λ(δ δ + δ δ ) Q = ik jl jk il , klφ†φ† 0− 0 (ij) 2λ(δikδjl + δklδij ) 0 Q † = − . = = 2√2iλv, = 4iλ. klφ φ 0 2λ(δjkδil + δklδij ) − − −

APPENDIX B: DIAGRAMMATIC REPRESENTATION One-particle irreducible contributions to the vacuum expectation value of the scalar field φ2, see Eq. (11): In the following, we summarize the diagrammatic representation of the individual terms in the effective potential as well as the propagators. The flavor components = + φ1 and φ2 are denoted by oriented dashed and dotted lines, respectively. Since the 22 propagator is not di- 1 D + + agonal in the Nambu space, it has to carry two arrows, 2 specifying the operators it connects at its ends. The bare propagators (8) are represented by the following lines. 1 We have explicitly indicated the symmetry factor 2 at the † † ∆φφ = ∆φ φ = last graph. Note that in our formalism all such factors are 11 11 automatically generated by the functional differentiation φφ† φφ ∆22 = ∆22 = of the effective action. φ†φ† φ†φ Next is the diagrammatic representation of the propa- ∆ = ∆ = 22 22 gators. The loop graphs are in one-to-one correspondence with the terms in the formulas (17), (19) and (20).

† † −1φφ = ∆−1φφ + + + + + + D11 11

† † −1φφ = ∆−1φφ + + + D22 22

1 + + + + + 2

1 −1φφ = ∆−1φφ + + D22 22 2

1 + + + + + 2 12

APPENDIX C: PROOF OF THE EXISTENCE OF GAPLESS POLES

Here we provide a detail proof of the existence of massless poles in the propagators 11 and 22. The key ingredients are the formulas (17), (19) and (20), as well as the gap equation (11). D D

1. Pole in D11

† We are going to prove that −1φφ (0) = 0. Eq, (17) together with the expressions (8) for the bare propagators give D11 4 4 2 2 2 † d k 1 d k (k µ) ǫ 2λv −1φφ (0) = µ2 (M 2+δM 2) (λ+δλ)v2 4iλ 2iλ 0 − − k − D22 − − − (2π)4 (k + µ)2 ǫ2 λv2 − (2π)4 [k2 E2 (k)][k2 E2 (k)] Z 0 − k − Z 0 − + 0 − − d4k 1 [(k + µ)2 ǫ2 λv2] + [(k µ)2 ǫ2 λv2] 2iλ2v2 0 − k − 0 − − k − , (C1) − (2π)4 (k + µ)2 ǫ2 λv2 [k2 E2 (k)][k2 E2 (k)] Z 0 − k − 0 − + 0 − − 2 2 where we abbreviate ǫk = √k + M . Next we expand the gap equation (11), using the explicit form of the propagators (8), d4k 1 µ2 (M 2 + δM 2) (λ + δλ)v2 =2iλ − − (2π)4 (k + µ)2 ǫ2 λv2 Z 0 − k − d4k (k + µ)2 ǫ2 2λv2 d4k 1 +4iλ 0 − k − +2iλ2v2 . (2π)4 [k2 E2 (k)][k2 E2 (k)] (2π)4 [k2 E2 (k)][k2 E2 (k)] Z 0 − + 0 − − Z 0 − + 0 − − This expression is used to substitute for µ2 (M 2 + δM 2) (λ + δλ)v2 in Eq. (C1), which ﬁnally yields − − 4 2 2 2 2 † d k 1 (k µ) ǫ 2λv λv −1φφ (0) = 2iλ 0 − − k − 1 . D11 − (2π)4 (k + µ)2 ǫ2 λv2 − [k2 E2 (k)][k2 E2 (k)] − (k + µ)2 ǫ2 2λv2 Z 0 − k − 0 − + 0 − − 0 − k − Using the identity det ∆−1(k) = [k2 E2 (k)][k2 E2 (k)] = [(k + µ)2 ǫ2 2λv2][(k µ)2 ǫ2 2λv2] (λv2)2, 22 0 − + 0 − − 0 − k − 0 − − k − − which follows from Eq. (8), it is now straightforward to show that the integrand in the last expression vanishes.

2. Pole in D22

In this case we are going to prove Eq. (21). We shall use a more symbolic method, without referring to the explicit form of the bare propagators (8). Upon a comparison of Eqs. (19) and (20) we can see that several terms immediately cancel out. After a mild rearrangement we are to prove the identity

† † † † µ2 (M 2 + δM 2) (λ + δλ)v2 =2iλ ∆φφ + 2∆φφ ∆φφ +4iλ2v2 ∆φφ ∆φ φ ∆φφ∆φφ . (C2) − − 11 22 − 22 22 22 − 22 22 Z Z det∆22 For the sake of legibility we dropped out the integration measure and arguments of| the propagators.{z Also,} we used of φφ φ†φ† the fact that ∆22 = ∆22 . 2 2 −1 φφ On the right hand side of Eq. (C2), we rewrite λv det ∆22 as λv / det ∆22 = ∆22 and immediately arrive at an expression identical to the right hand side of Eq. (11). Eq. (C2) is thus proved and hence also the existence of a massless pole in the propagator . D22

ACKNOWLEDGMENTS manuscript. The present work was supported in part by the Institutional Research Plan AV0Z10480505, and by The author is indebted to J. Hoˇsek and J. Novotn´y the GACR grants No. 202/06/0734 and 202/05/H003. for helpful discussions and/or critical reading of the

[1] J. Goldstone, Nuovo Cim. 19, 154 (1961). [2] J. Goldstone, A. Salam, and S. Weinberg, Phys. Rev. 13

127, 965 (1962). [18] D. T. Son (2002), hep-ph/0204199. [3] H. B. Nielsen and S. Chadha, Nucl. Phys. B105, 445 [19] J. I. Kapusta, Finite-Temperature Field Theory, Cam- (1976). bridge Monographs on Mathematical Physics (Cam- [4] V. A. Miransky and I. A. Shovkovy, Phys. Rev. Lett. 88, bridge University Press, Cambridge, 1989). 111601 (2002), hep-ph/0108178. [20] M. E. Peskin and D. V. Schroeder, An Introduction [5] T. Schaefer, D. T. Son, M. A. Stephanov, D. Toublan, to Quantum Field Theory (Addison–Wesley, Reading, and J. J. M. Verbaarschot, Phys. Lett. B522, 67 (2001), 1995). hep-ph/0108210. [21] J. O. Andersen, Rev. Mod. Phys. 76, 599 (2004), cond- [6] M. Buballa, J. Hoˇsek, and M. Oertel, Phys. Rev. Lett. mat/0305138. 90, 182002 (2003), hep-ph/0204275. [22] C. Manuel, A. Dobado, and F. J. Llanes-Estrada, JHEP [7] D. Blaschke, D. Ebert, K. G. Klimenko, M. K. Volkov, 09, 076 (2005), hep-ph/0406058. and V. L. Yudichev, Phys. Rev. D70, 014006 (2004), [23] C. P. Hofmann, Phys. Rev. B60, 388 (1999), cond- hep-ph/0403151. mat/9805277. [8] A. Beraudo, A. De Pace, M. Martini, and A. Molinari, [24] P. F. Bedaque and T. Schaefer, Nucl. Phys. A697, 802 Ann. Phys. 317, 444 (2005), nucl-th/0409039. (2002), hep-ph/0105150. [9] T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). [25] D. B. Kaplan and S. Reddy, Phys. Rev. D65, 054042 [10] T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 (2002), hep-ph/0107265. (1998). [26] H. J. Warringa (2006), hep-ph/0606063. [11] L. He, M. Jin, and P. Zhuang (2006), cond-mat/0604580. [27] Strictly speaking, this is only true for the normal phase. [12] F. Sannino and W. Schäfer, Phys. Lett. B527, 142 In the Bose–Einstein condensed phase, there is a finite (2002), hep-ph/0111098. wave function renormalization induced by the second [13] F. Sannino, Phys. Rev. D67, 054006 (2003), hep- term in the integral in Eq. (14). This wave function renor- ph/0211367. malization is, in fact, responsible for the renormalization [14] T. Brauner, Phys. Rev. D72, 076002 (2005), hep- of the Goldstone boson dispersion relations. ph/0508011. [28] This is of course a standard textbook result — see, for [15] H. Leutwyler, Phys. Rev. D49, 3033 (1994), hep- instance, Eq. (11.76) in Ref. [20]. It should be noted, how- ph/9311264. ever, that we use a slightly different regularization than [16] S. R. Coleman and E. Weinberg, Phys. Rev. D7, 1888 usual, namely the dimensional regularization in three (1973). spatial dimensions. [17] J. O. Andersen (2005), hep-ph/0501094.