text:
# The low-energy phase-only action in a superconductor: a comparison with the 𝑋𝑌 model

## I Introduction

The description of the low-energy dynamics of the phase $`\theta `$ of the superconducting order parameter $`\mathrm{\Delta }=|\mathrm{\Delta }|e^{i\theta }`$ has recently attracted a renewed interest in connection to the phenomenology of high-$`T_c`$ superconducting cuprates. These materials are characterized by a strongly anisotropic (quasi-two-dimensional) crystal structure and by anomalous normal-state properties, which can be interpreted in terms of preformed Cooper pairs (non-zero amplitude $`|\mathrm{\Delta }|`$) without phase coherence (vanishing superfluid density $`\rho _s`$).review  Although in ordinary (weak-coupling) superconductors both $`|\mathrm{\Delta }|`$ and $`\rho _s`$ vanish at the critical temperature $`T_c`$, due to thermally excited quasiparticles, it is possible that $`\rho _s`$ vanishes when $`|\mathrm{\Delta }|`$ is finite, due to phase fluctuations. If such is the case, phase fluctuations should play a role more relevant than expected in ordinary superconductors, both near $`T_c`$ and deeply in the superconducting state.

A crucial point in the analyses of phase-fluctuation effects is clearly the choice of the effective model, the most common being the classical $`XY`$ model, as the paradigm for the universality class which is relevant to superconductivity.pf  The $`XY`$ model, which has also been extensively used in contexts such as systems of resistively shunted Josephson junctions or granular superconductors,chacra  assumes a phase field $`\theta _i`$ defined on a coarse-grained lattice with a lattice constant determined by the coherence length $`\xi _0`$, which sets the length scale above which the fluctuations of the amplitude $`|\mathrm{\Delta }|`$ become uncorrelated, and $`|\mathrm{\Delta }|`$ may be assumed to be fixed. Thus, at distances larger than $`\xi _0`$, phase fluctuations only are relevant, and can be described in terms of a Josephson-like interaction between neighboring sites $`<i,j>`$,

$$H_{cl}=\frac{D_0}{4}\underset{<i,j>}{}\left(1\mathrm{cos}\theta _{ij}\right)\frac{D_0}{8}𝑑𝐱\left[(\stackrel{}{}\theta )^2\frac{\xi _0^2}{12}\underset{\alpha =x,y}{}\left(\frac{\theta }{\alpha }\right)^4+\mathrm{}\right].$$
(1)
Here $`\theta _{ij}\theta _i\theta _j|\stackrel{}{}\theta |\xi _0`$ is the phase difference between nearest-neighboring sites, and $`D_0`$ is the bare coupling constant. Eq. (1) is written for a two-dimensional system (which is the case relevant for cuprate superconductors). In $`d=2`$, $`D_0`$ is an energy which measures the superfluid stiffness, i.e., the energetic cost to produce phase variations in the system. In the continuum (Galilean invariant) model $`D_0=\rho _s(T=0)/m=\rho /m`$, where $`m`$ is the electron mass, $`\rho `$ is the electron density, and we set $`\mathrm{}=1`$. In the case of a generic $`d`$-dimensional system, we still refer to $`D_0`$ as to the stiffness, even though the corresponding energy scale is $`D_0\xi _0^{d2}`$.

Within the $`XY`$ model, one can determine the depletion of the superfluid stiffness due to the collective-mode excitations, i.e., to the anharmonic terms in the gradient expansion of Eq. (1). At low temperature, the leading correction to the bare stiffness $`D_0`$ comes from the fourth-order term, and within perturbation theory with respect to the harmonic term one gets, e.g. in $`d=2`$,

$$D(T)=D_0\left[1\frac{\xi _0^2}{4}(\stackrel{}{}\theta )^2\right],$$
(2)
where $`D_0=D(T=0)`$, since $`(\stackrel{}{}\theta )^2=4T/\xi _0^2D_0`$ when evaluated at Gaussian level. This result led to the proposal that phase fluctuations, rather than the $`d`$ wave quasiparticles excitations characteristic of the superconductivity in the cuprates, are responsible for the linear thermal depletion of the superfluid densitypf  which is experimentally observed in these materials, down to very low temperatures (1—5 K).exp  However, one would expect that quantum effects are relevant at such low temperatures, and, even though the $`XY`$ model may be reasonably adopted in the classical regime, the investigation of the quantum regime is much more involved.

It has been proposedarun00 ; noi2  that the quantum effects can be partially included in the so-called quantum $`XY`$ model, by deriving the Gaussian phase-only effective action from a microscopic BCS model, while obtaining the anharmonic terms through the expansion the $`\mathrm{cos}\theta _{ij}`$ term of Eq. (1). As a consequence, the interaction terms in the phase, given by powers of $`\theta _{ij}\xi _0|\stackrel{}{}\theta |`$, are purely classical, whereas the Gaussian propagator is evaluated in the quantum regime. However, this approach turns out to be unsatisfactory within many respects.

Indeed, since in the quantum case $`(\stackrel{}{}\theta )^20`$ at $`T=0`$ (see below), Eq. (2) leads to a finite correction to $`D_0`$ even at $`T=0`$, which, in turn implies $`\rho _s(0)\rho `$ in the continuum model, thus explicitly violating Galilean invariance. Therefore, the description of the phase fluctuations by means of the quantum $`XY`$ model misses some important effect within the continuum model. A second issue, which is particularly relevant in connection to the description of high-$`T_c`$ superconductors, concerns the limit of strong pairing interaction. Indeed, while at weak coupling the amplitude fluctuations can be safely neglected in deriving the phase-only theory, as the interaction increases, the fluctuations of the modulus of the order parameter become intimately connected with the density fluctuations.depalo ; noi3  Since the phase and the density are conjugate variables, the description of phase fluctuations requires a proper treatment of the amplitude fluctuations in the strong-coupling regime.

In this paper we present a systematic comparison between the phase-only model derived microscopically and the quantum and classical $`XY`$ models. To this extent, after describing the formal steps to derive the effective action in various cases (continuum vs lattice model, neutral vs charged system), we evaluate the one-loop corrections to $`D_0`$ coming from the non-Gaussian (anharmonic) terms.

We find that quantum effects generate dynamical interaction terms, absent in the quantum $`XY`$ model, whose role is crucial, e.g., in preserving the Galilean invariance in the continuum model, and in reproducing the expected mapping of the lattice model onto the $`XY`$ model in the strong-coupling regime. Some of the results on quantum corrections in the weak-coupling limit where already presented in a short paper,noi3  and we provide here a detailed derivation, since the issue, which involves many subtleties, has been only partially addressed in the literature.arun00 ; noi2 ; depalo ; noi3 ; aitchison95 ; dattu98 ; aitchison00 ; kwon ; sharapov ; sharapov2 ; kim  We mainly focus on the quantum ($`T=0`$) corrections to the superfluid density, and we comment on the possible extension to the classical regime. Indeed, besides the restoration of the Galilean invariance in the continuum system, the microscopic derivation of a phase-only theory, developed in this paper, allows to demonstrate clearly that, in general, the $`XY`$ model leads to an overestimation of the depletion of the superfluid $`\delta D`$ stiffness, even when the realistic effects such as the Coulomb screening are taken into account, with the only exception of the extreme strong-coupling regime where the $`XY`$ and the microscopic estimate of $`\delta D`$ tend eventually to coincide. For sake of clarity we summarize below our final results, with reference to the corresponding Sections where they are derived, as a guideline for the reader. We report in tab. I and II (respectively for the continuum and the lattice models) the estimations for the ratio between the corrections to the superfluid density, derived within the microscopic models ($`\delta D`$) and within the $`XY`$ model ($`\delta D_{XY}`$) in all the cases considered in this work.

The specific issue of the temperature dependence of the superfluid density in high-$`T_c`$ superconductors has been analyzed in Ref. \[arun00, ; noi2, \]. In general, according to Eq. (1), one would expect that Coulomb effects, lifting the phase mode to the plasma frequency (see Sec. II B), completely suppress the contribution of phase fluctuations to $`D(T)`$. However, in Ref. \[noi2, \] it has been shown that the $`d`$-wave symmetry of the order parameter plays a crucial role in determining the behavior of the phase mode in the presence of both long-range Coulomb interactions and dissipation. Even though here we ignore dissipative effects, for the sake of completeness we still compare $`s`$-wave and $`d`$-wave superconductors. Once again, the formalism of the effective action makes this comparison, to some extent, particularly simple, and of transparent physical interpretation.

The plan of the paper is the following. In Sec. II.1 we derive the Gaussian phase-only action starting from the continuum model, in the absence of long-range Coulomb forces. The case of the charged system is dealt with in Sec. II.2. Although most of the results presented in Sec. II are standard, our scope here is to introduce the formalism and a classification scheme for the interaction vertices, which turns out to be particularly useful in dealing with non-Gaussian terms. In Sec. III we apply the same procedure to a lattice model, which does not appear as a trivial extension of the continuum model. In Sec. IV we calculate the anharmonic terms, providing a classification of the diagrams (Sec. IV.1), and determine the one-loop correction to the stiffness within the continuum model, both at $`T=0`$ (Sec. IV.2) and in the classical limit (Sec. IV.3). The results are then compared with the results of the $`XY`$ model. In Sec. V the same analysis is carried out in the lattice case. The strong-coupling regime of the lattice model is discussed in Sec. VI, and the results are summarized in Sec. VII. The details on the derivation of the phase-only action are reported in App. A, whereas in App. B we discuss the connection between the coefficients of the Gaussian action and the gauge-invariant electromagnetic response functions.

## II The effective action for the continuum model

### II.1 The neutral system

To introduce the general formalism adopted in this paper, and for the sake of completeness, we shortly discuss the standard formal steps to derive the phase-only action within the continuum BCS model $`H=H_0+H_I`$, with

$`H_0`$ $`=`$ $`{\displaystyle \underset{\sigma }{}}{\displaystyle 𝑑𝐫\frac{1}{2m}c_\sigma ^+(𝐫)^2c_\sigma (𝐫)},`$ (3)
$`H_I`$ $`=`$ $`{\displaystyle \frac{U}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤,𝐤^{^{}},𝐪}{}}w(𝐤)w(𝐤^{})c_{𝐤+\frac{𝐪}{\mathrm{𝟐}}}^+c_{𝐤+\frac{𝐪}{\mathrm{𝟐}}}^+c_{𝐤^{^{}}+\frac{𝐪}{\mathrm{𝟐}}}c_{𝐤^{}+\frac{𝐪}{\mathrm{𝟐}}}.`$ (4)
Most of our results are obtained for generic spatial dimension $`d`$ of the system, unless explicitly indicated. In Eqs. (3)-(4), $`c_\sigma ^{(+)}`$ is the annihilation (creation) operator of an electron of spin $`\sigma `$, $`\mathrm{\Omega }`$ is the volume, and $`U>0`$ is the pairing interaction strength. Hereafter we assume $`\mathrm{}=k_B=1`$. The factor $`w(𝐤)`$ controls the symmetry of the Cooper-pair wave function, e.g., $`w(𝐤)=1`$ for $`s`$-wave superconductors, and $`w(𝐤)=\mathrm{cos}2\varphi `$, with $`\varphi =\text{atan}(k_y/k_x)`$, for $`d`$-wave superconductors, in the continuum case and for $`d=2`$. Note that, in the case of a $`d`$-wave symmetry, four nodes are present in the quasiparticle excitation gap, and the existence of gapless excitations is relevant in determining the low-temperature thermodynamics of the system.

We first discuss the case of a neutral system, while the Coulomb interaction will be introduced in Sec. II.2, to deal with the charged case. As customary, the microscopic effective model for the collective modes is derived by considering the action corresponding to the Hamiltonian (3)-(4), within the finite-temperature Matsubara formalism,

$$S_{micro}=S_0+S_I=_0^\beta 𝑑\tau \left\{\underset{𝐤\sigma }{}c_{𝐤\sigma }^+(\tau )[_\tau +\xi _𝐤]c_{𝐤\sigma }(\tau )d\tau +H_I(\tau )\right\},$$
(5)
where $`\tau `$ is the imaginary time, $`\beta =1/T`$, and $`\xi _𝐤=ϵ_𝐤\mu `$ is the band dispersion with respect to the chemical potential $`\mu `$, with $`ϵ_𝐤=𝐤^2/2m`$ in the continuum case. To obtain the effective action in terms of the order-parameter collective degrees of freedom, the interaction is decoupled in the particle-particle channel by means of the Hubbard-Stratonovich transformation, introducing the auxiliary complex field $`\mathrm{\Delta }(x,\tau )`$, as explicitly reported in App. A. It is then possible to make the dependence on the phase of $`\mathrm{\Delta }(𝐪,\tau )`$ explicit in the action, by performing the gauge transformation Eq. (94) on the fermionic fields. Indeed the dependence on the phase is eliminated from $`S_I`$, and made explicit in the free part of the action $`S_0`$, which now reads

$$\stackrel{~}{S}_0=S_0+𝑑x\varphi ^+(x)\widehat{\mathrm{\Sigma }}\varphi (x),$$
(6)
with

$$\widehat{\mathrm{\Sigma }}=\{\frac{i}{2}_\tau \theta (x)+\frac{1}{8m}\left[\stackrel{}{}\theta (x)\right]^2\}\widehat{\tau }_3+\{\frac{i}{4m}\stackrel{}{}\theta (x)\stackrel{}{}\}\widehat{\tau }_0,$$
(7)
where $`\widehat{\tau }_i`$ are the Pauli matrices and the operator $`\stackrel{}{}(\stackrel{}{}\stackrel{}{})`$ acts on the fermionic Nambu spinor $`\varphi (x)`$ defined as the column vector $`(c_{}(x),c_{}^+(x))`$. It is worth noting that the field $`\theta `$ appears in the action only through its time and spatial derivatives, as it is expected for a Goldstone (massless) field. Moreover we point out that, while in the $`s`$-wave case the transformation (94) completely eliminates the phase $`\theta `$ from $`S_I`$, in the $`d`$-wave case a residual dependence survives. Nonetheless, in App. A, where we report some more detail, we show that this residual dependence can be safely neglected in the hydrodynamic limit, which is the relevant regime when discussing the low-energy properties of the collective mode. Thus one can now integrate out the fermions, which appear quadratically in $`\stackrel{~}{S}_0+S_I`$, leading to the effective action (99), written in terms of the collective variables only. This standard procedure allows one to analyze the starting problem (3)-(4) from a different point of view. Instead of perturbatively expanding the Hamiltonian in powers of the coupling $`U`$, the effective action (99) is expanded by assuming that the fluctuations of the fields around their saddle-point values are small, and that the variations of the phase in space and time are slow. Moreover, if one is interested in the dynamics of the phase at low temperature, the fluctuations of the modulus $`|\mathrm{\Delta }|`$ can be neglected. As far as the specific problem of high-$`T_c`$ superconductors is concerned, this approximation seems particularly appropriate for underdoped cuprates, which exhibit an anomalously large binding energy when compared to the critical temperature.review  On the other hand, one should expect that modulus fluctuations can still be relevant at distances smaller than the typical length scale over which $`|\mathrm{\Delta }|`$ fluctuates. It is then natural to introduce a spatial cut-off for phase fluctuations, of the order of the coherence length $`\xi _0`$. At distances greater than $`\xi _0`$ the phase fluctuations are the only relevant degrees of freedom, and a description in terms of the phase-only action becomes meaningful. Accordingly, the cut-off $`\xi _0`$ for the phase-only action is defined as the characteristic length scale of the spatial decay of the correlation function for $`|\mathrm{\Delta }|`$. This point is delicate, since the naive identification of the coherence length $`\xi _0`$ with the Cooper-pair size $`\xi _{pair}`$ is meaningless in $`d`$-wave superconductor, where $`\xi _{pair}`$ diverges at $`T=0`$, due to the presence of the gapless quasiparticle excitations mentioned above, as it is discussed in Ref. \[noi4, \].

Starting from Eq. (99) of App. A, with $`|\mathrm{\Delta }|=|\mathrm{\Delta }|(0)+\delta |\mathrm{\Delta }|(𝐪)`$, and neglecting the fluctuations of the modulus $`\delta |\mathrm{\Delta }|(𝐪)`$, the sum $`U^1_q|\mathrm{\Delta }|(q)^2`$ reduces to its $`q=0`$ value. As a consequence, the self-energy (100) does no longer depend on $`|\mathrm{\Delta }|`$. The action (99) is then decomposed as $`S_{neutral}=S_{MF}+S_{eff}(\theta )`$, where $`S_{MF}=|\mathrm{\Delta }|^2(0)/U\text{Tr}\left[\mathrm{ln}𝒢_0^1\right]`$ is the mean-field BCS action, and $`S_{eff}`$ is the phase-only

$$S_{eff}(\theta )=\text{Tr}\underset{N=1}{\overset{\mathrm{}}{}}\frac{1}{N}(\widehat{\mathrm{\Sigma }}\widehat{𝒢}_0)^N=\underset{N=1}{\overset{\mathrm{}}{}}S_{eff}^{(N)},$$
(8)
which takes into account phase fluctuations around the BCS saddle point. Here $`\widehat{𝒢}_0`$ is the BCS matrix defined in Eq. (101), and the trace Tr is explicitly given in $`k`$-space by

$$S_{eff}^{(N)}=\frac{1}{N}\underset{k_1,\mathrm{},k_N}{}\text{tr}[\widehat{𝒢_0}(k_1)\widehat{\mathrm{\Sigma }}(k_1k_2)\times \widehat{𝒢_0}(k_2)\widehat{\mathrm{\Sigma }}(k_2k_3)\times \widehat{𝒢_0}(k_N)\widehat{\mathrm{\Sigma }}(k_Nk_1)],$$

where $`\mathrm{tr}`$ is the trace in the Nambu space only and $`\widehat{\mathrm{\Sigma }}(k)`$ is the Fourier transform of the expression (7). Since the self-energy (7) contains powers of the field $`\theta `$ up to the second, the expansion (8) is not directly an expansion in powers of $`\theta `$ and the terms with the same power of $`\theta `$ must be recollected, leading to the expansion

$$S_{eff}(\theta )=\underset{k=2}{\overset{\mathrm{}}{}}S_k=\underset{k=2}{\overset{\mathrm{}}{}}\underset{q_1,\mathrm{},q_{k1}}{}A_k(q_1,\mathrm{},q_{k1})\theta (q_1)\mathrm{}\theta (q_{k1})\theta (q_1\mathrm{}q_{k1}),$$
(9)
which is completely equivalent to Eq. (8), and includes both the Gaussian ($`k=2`$) and the anharmonic ($`k>2`$) terms.

Each term in the effective action (9) is given by a closed fermionic loop $`A_k`$ with $`k`$ incoming lines corresponding to fluctuating $`\theta `$ fields. In forming the fermionic loops, according to Eq. (7), we have three kinds of fermion-$`\theta `$ vertices, depicted in Fig. 1:

* The vertex with a single incoming dotted line corresponds to the insertion of the time derivative of the phase, i.e., the term $`\widehat{\mathrm{\Sigma }}_B^1=(i_\tau \theta /2)\widehat{\tau }_3`$ of Eq. (7).
* The vertex with two incoming dashed lines corresponds to the insertion of the gradient squared of the phase, i.e., the term $`\widehat{\mathrm{\Sigma }}_B^2=[(\stackrel{}{}\theta )^2/8m]\widehat{\tau }_3`$ of Eq. (7).
* The vertex with a single incoming dashed line corresponds to the insertion of the spatial derivative of the phase associated to the spatial derivative of the incoming and outgoing Green functions, i.e., the term $`\widehat{\mathrm{\Sigma }}_F^1=i\{[\stackrel{}{}\theta (x)\stackrel{}{}]/4m\}\widehat{\tau }_0`$ of Eq. (7).

The terms generated in the effective action by the self-energy $`\widehat{\mathrm{\Sigma }}_B=\widehat{\mathrm{\Sigma }}_B^1+\widehat{\mathrm{\Sigma }}_B^2`$, associated to the $`\widehat{\tau }_3`$ matrix, are essentially the same which are met when deriving the phase-only action in a bosonic superfluid system.popov  We therefore refer to these terms as “bosonic” contributions. The $`\widehat{\mathrm{\Sigma }}_F^1`$ term generates instead contributions which are strictly related to the existence of quasiparticle excitations. Therefore, we refer to them as “fermionic” contributions. Clearly, the expansion (8) also generates “mixed” contributions, obtained by combinations of $`\widehat{\tau }_0`$ and $`\widehat{\tau }_3`$ terms of the self-energy $`\widehat{\mathrm{\Sigma }}`$. Such a distinction will be relevant in the following, while discussing the derivation of the anharmonic terms of the phase-only action.

At Gaussian level $`(k=2)`$ the effective action describes the collective mode of the system. Beyond the Gaussian level $`(k>2)`$, interaction terms in the field $`\theta `$ appear, and we obtain the quantum generalization of the anharmonic terms of the $`XY`$ model, Eq. (1). As in the classical case, the anharmonic terms are responsible for the reduction of the stiffness due to phase fluctuations. Before discussing this issue, we analyze the structure of the Gaussian phase propagator in the quantum regime. The Gaussian effective action is derived with long but straightforward calculations, and reads

$$S_{neutral}^G=\frac{1}{8}\underset{q}{}\left[\mathrm{\Omega }_m^2\mathrm{\Lambda }_{\rho \rho }(q)𝐪_a𝐪_b\mathrm{\Lambda }_{JJ}^{ab}(q)+2i\mathrm{\Omega }_m𝐪_a\mathrm{\Lambda }_{\rho J}^a(q)\right]|\theta (q)|^2.$$

The coefficients $`\mathrm{\Lambda }_{\rho \rho },\mathrm{\Lambda }_{\rho J}`$ and $`\mathrm{\Lambda }_{JJ}`$, which are depicted in Fig. 2 and whose explicit expressions are reported in App. A, represent respectively the density-density, density-current, and current-current bubbles, evaluated with the BCS Green functions. Since the BCS approximation explicitly breaks the gauge invariance,schrieffer  the $`\mathrm{\Lambda }`$ functions do not coincide with the physical correlation functions $`K`$ sharapov2 , which must obey the gauge invariance. However, as it is shown explicitly in App. B (where the general relationship between the $`\mathrm{\Lambda }`$ and the $`K`$ bubbles is obtained), it is possible to relate the coefficients of Gaussian effective action with measurable quantities such as the bare compressibility $`\kappa (T)`$ and the stiffness $`D(T)=\rho _s(T)/m`$ of the system, in the hydrodynamic limit.notacompr  Indeed, in this regime, we obtain the expression

$$S_{hydr}^G=\frac{1}{8}\underset{q}{}\left[\kappa _0(T)\mathrm{\Omega }_m^2+D(T)𝐪^2\right]\theta (q)\theta (q).$$
(10)
The Bogoljubov sound mode of the neutral system, with velocity $`c_s=\sqrt{D(T)/\kappa _0(T)}`$, is thus recovered. It is worth noting that, since the phase and the density are conjugate variables, the collective mode which appears in the phase-field propagator is the same as the density mode (see, e.g., Ref. \[anderson, \]), as it can be explicitly checked by deriving the physical density-density correlation function $`K_{\rho \rho }`$ according to Eq. (110).

By comparing Eq. (10) with the classical $`XY`$ model, one can recognize the same term $`\rho _s(\theta )^2`$ which appears in Eq. (1). However, the quantum Gaussian propagator for phase fluctuations contains also the $`\mathrm{\Omega }_m^2`$ term of Eq. (10) which takes into account the dynamic effects. Moreover, according to the results of App. B and to Eqs. (105)-(106), the temperature dependence of the bare stiffness $`D(T)`$ which appears in Eq. 10 is entirely due to the quasiparticle excitations. Indeed, the coefficient $`\mathrm{\Lambda }_{jj}^{ab}`$ in Eqs. (105)-(106) describes the normal-fluid component, whose contribution increases as the temperature increases

$$D(T)=\frac{\rho _{MF}}{m}\frac{2}{\mathrm{\Omega }}\underset{𝐤}{}\frac{𝐤^2}{m^2}\left(\frac{f}{E_𝐤}\right)=\frac{\rho _s(T)}{m}.$$
(11)
In the absence of dissipation the coefficients $`\mathrm{\Lambda }`$ are not analytic at $`q=0`$ and finite temperature, since in this case the static ($`\mathrm{\Omega }_m=0,𝐪0`$) and the dynamic limit ($`𝐪=0,\mathrm{\Omega }_m0`$) are different. This issue has been recently addressed in Ref. \[sharapov, \]. Since we are interested in the $`T=0`$ correction to the superfluid density we shall not discuss further the consequences of such non analyticity, which would disappear in any case in the presence of dissipation.

### II.2 The charged system

When dealing with a charged system, the long-range Coulomb interaction between the fermions,

$$H_{Coul}=\frac{1}{2\mathrm{\Omega }}\underset{𝐤,𝐤^{},𝐪}{}\underset{\sigma ,\sigma ^{^{}}}{}V(𝐪)c_{𝐤+𝐪,\sigma }^+c_{𝐤^{}𝐪,\sigma ^{}}^+c_{𝐤^{},\sigma ^{}}c_{𝐤,\sigma },$$
(12)
must be included in the Hamiltonian (3)-(4). $`V(𝐪)`$ is the Fourier transform of the three-dimensional Coulomb repulsion $`V(r)=e^2/\epsilon _Br`$, projected onto a $`d`$ dimensional system (e.g., $`d=2`$ for a single layer). Thus, for generic $`d`$, $`V(𝐪)=\lambda e^2/|𝐪|^{d1}`$, where $`\lambda `$ is a constant which depends on the dimension $`d`$, e.g., $`\lambda =4\pi /\epsilon _B`$ for $`d=3`$ (isotropic three-dimensional system) and $`\lambda =2\pi /\epsilon _B`$ for $`d=2`$ (single layer). Here $`\epsilon _B`$ is the dielectric constant of the ionic background. The background ensures overall charge neutrality and cancels out the apparent divergence of the Coulomb interaction $`V(𝐪)`$ at $`𝐪=0`$. We point out that Eq. (12) generically describes a density-density interaction mediated by the potential $`V(𝐪)`$. The particular choice of $`V(𝐪)`$ distinguishes the short-range from the long-range case.

After including the term (12) in Eq. (5), we follow the standard procedure and introduce a Hubbard-Stratonovich field $`\rho _{HS}`$ associated to the electron density, thus decoupling the Coulomb interaction in the particle-hole channel. Since the short-range pairing in Eq. (4) is important for small center-of-mass momentum, while the Coulomb effects are important for small momentum transfer, one can reasonably assume that the breakup of the actual interaction in this manner is physically sensible and does not lead to any “overcounting”. Note that in principle also the pairing term should be decoupled in the particle-hole channel, but these terms are safely negligible, at least for weak and intermediate $`U`$. Nevertheless, in Sec. VI we shall discuss the subtleties of such a procedure, in the case of strong pairing interaction.

As described in App. A, by separating $`\rho _{HS}(q)=\rho _{HS}(0)+\delta \rho _{HS}(q)`$, and by neglecting the fluctuations of the modulus $`\delta |\mathrm{\Delta }|`$, the expansion of the phase-only action around the BCS saddle point leads to $`S_{charged}=S_{MF}+S_{eff}(\theta ,\rho _{HS})`$, where now

$$S_{eff}(\theta ,\delta \rho _{HS})=\underset{q}{}\frac{\delta \rho _{HS}^2(𝐪)}{2V(𝐪)}+\text{Tr}\underset{N=1}{\overset{\mathrm{}}{}}\frac{1}{N}(\widehat{\mathrm{\Sigma }}_c\widehat{𝒢}_0)^N,$$
(13)
and

$$\widehat{\mathrm{\Sigma }}_c=[\frac{i_\tau \theta (x)}{2}+\frac{(\stackrel{}{}\theta (x))^2}{8m}+i\delta \rho _{HS}(x)]\widehat{\tau }_3+[\frac{i}{4m}\stackrel{}{}\theta (x)\stackrel{}{}]\widehat{\tau }_0.$$
(14)
The Hubbard-Stratonovich term $`\delta \rho _{HS}^2(0)/V(0)`$ does not appear in the mean-field action $`S_{MF}`$ as it is supposed to be canceled out, in order to preserve the charge neutrality of the system, by an equal (and opposite in sign) contribution coming from the interaction with the ionic background. This term survives instead in the case of short-range interaction. Since the self-energy $`\widehat{\mathrm{\Sigma }}_c`$ depends now also on the second Hubbard-Stratonovich field $`\delta \rho _{HS}`$, to obtain the phase-only action, one has to proceed in two steps, as is explicitly shown in App. A: first, the Gaussian action in both $`\theta `$ and $`\delta \rho _{HS}`$ is derived and, then, after integrating out the field $`\delta \rho _{HS}`$, one obtains the Gaussian phase-only action

$$S_{charged}^G(\theta )=\frac{1}{8}\underset{q}{}\left[\mathrm{\Omega }_m^2_{\rho \rho }(q)𝐪_a𝐪_b_{JJ}^{ab}(q)+2i𝐪_a\mathrm{\Omega }_m_{\rho J}^a(q)\right]\theta (q)\theta (q),$$
(15)
which generalizes Eq. (II.1) for a charged system. The coefficients $``$, whose definitions are reported in App.A, are simply the coefficients $`\mathrm{\Lambda }`$ of the unscreened phase-only action (II.1) dressed by the density fluctuations within the random-phase approximation (RPA). This represents the only difference between the Gaussian actions (II.1) and (15), which have therefore the same formal structure.

The inclusion of Coulomb interaction obviously reflects on the collective mode of the system. Indeed, the expression of the Gaussian action (15) in the hydrodynamic limit differs from Eq. (10) only in the coefficient of the $`\mathrm{\Omega }_m^2`$ term, since both $`\mathrm{\Lambda }_{\rho j}^a`$ and $`_{\rho j}^a`$ vanish in the static limit, and $`_{JJ}^{ab}`$ behaves as $`\mathrm{\Lambda }_{JJ}^{ab}`$, whereas $`_{\rho \rho }`$ gives the RPA compressibility of the charged system

$$_{\rho \rho }(𝐪,\mathrm{\Omega }_m=0)\begin{array}{c}𝐪0\hfill \end{array}\frac{\kappa _0(T)}{1+V(𝐪)\kappa _0(T)}\begin{array}{c}𝐪0\hfill \end{array}\frac{1}{V(𝐪)},$$
(16)
which vanishes for $`𝐪0`$. As a consequence, the collective mode described by the action

$$S_{charged}^G=\frac{1}{8}\underset{q}{}\left[\frac{\mathrm{\Omega }_m^2}{V(𝐪)}+D(T)𝐪^2\right]\theta (q)\theta (q),$$
(17)
is the plasma mode in $`d`$ dimensions. For example, in the three-dimensional case Eq. (17) reads

$$S_{charged}^G=\frac{1}{8}\underset{q}{}𝐪^2\left[\frac{1}{4\pi e^2}\mathrm{\Omega }_m^2+D(T)\right]\theta (q)\theta (q)$$
(18)
and the plasma mode is characterized by a finite energy at zero momentum, $`\mathrm{\Omega }_P=\sqrt{4\pi e^2D(T)}=\sqrt{4\pi e^2\rho _s(T)/m}`$. At $`T=0`$, where $`\rho _s`$ coincides with the particle density, $`\mathrm{\Omega }_P`$ is exactly the plasma frequency of a three-dimensional charged system.

The vanishing of the compressibility (16) in the presence of density-density interactions is peculiar of the Coulomb case. If one considers a short-range interaction $`V(𝐪)=\stackrel{~}{V}`$, Eq. (16) simply gives the RPA dressing of the bare compressibility

$$_{\rho \rho }(𝐪,\mathrm{\Omega }_m=0)\begin{array}{c}𝐪0\hfill \end{array}\kappa (T)=\frac{\kappa _0(T)}{1+\stackrel{~}{V}\kappa _0(T)},$$
(19)
which inserted into Eq. (15) leads to the Gaussian action

$$S_{short}^G=\frac{1}{8}\underset{q}{}\left[\frac{\kappa _0(T)}{1+\stackrel{~}{V}\kappa _0(T)}\mathrm{\Omega }_m^2+\frac{\rho _s(T)}{m}𝐪^2\right]\theta (q)\theta (q).$$

In the presence of short-range density-density interactions the collective mode is still a sound mode (known as Bogoljubov-Anderson modeanderson ), with velocity $`c_s^{}=\sqrt{\rho _s(T)/m\kappa (T)}=c_s\sqrt{1+\stackrel{~}{V}\kappa _0(T)}`$.

## III The effective action for the lattice model

We extend the standard procedures discussed in Sec. II, to deal with the more realistic case of a lattice model. This extension does not reduce to a trivial modifications of the continuum results, the most important differences being the appearance of vertices with more than two incoming dashed lines, and the impossibility to absorb the dressing of the fermionic lines in a chemical potential shift.

We first consider the neutral case, the generalization to the charged case being straightforward. The main difference between the continuum and the lattice model is in the kinetic term of the Hamiltonian. Here, we rewrite the free action (see Eqs. (3), (5)) introducing a hopping $`t`$ between nearest-neighboring sites on a lattice of spacing $`a=1`$,

$$S_{micro}=_0^\beta d\tau \{\underset{𝐱,\sigma }{}c_{𝐱\sigma }^+(\tau )(_\tau \mu )c_{𝐱\sigma }(\tau )t\underset{𝐱,𝐱^{},\sigma }{}[c_{𝐱\sigma }^+(\tau )c_{𝐱^{}\sigma }(\tau )+h.c.]+H_I(\tau )\},$$
(20)
so that the free-electron dispersion now reads $`\xi _𝐤=2t_j\mathrm{cos}k_j\mu `$, with $`j=x_1,\mathrm{},x_d`$. The different form of the kinetic term modifies the effect of the gauge transformation on the free action $`S_0`$. As we showed in Sec. II, after the transformation (94) the phase disappears from $`S_I`$ and modifies the free action $`S_0`$ according to Eq. (6). As a consequence, to obtain the self-energy $`\widehat{\mathrm{\Sigma }}`$ within the lattice model we analyze the effect of the gauge transformation on the free action (20), the subsequent steps being unchanged. We find two differences with respect to the continuum case:

* The terms $`\widehat{\mathrm{\Sigma }}_F^1`$ and $`\widehat{\mathrm{\Sigma }}_B^2`$ of the self-energy (7) arise, after the gauge transformation, from the kinetic term $`\stackrel{}{}^2/2m`$ of the free action. Since the analogous of the velocity $`𝐤/m`$ of the continuum case is $`\xi _𝐤/𝐤`$ in the lattice case, it follows that for a lattice model each insertion of a spatial derivative of $`\theta `$ in the fermionic loops is associated to a same-order $`𝐤`$-derivative of the band dispersion $`\xi _𝐤`$.
* On the lattice the coupling of the fermions to $`\theta `$ induced by the gauge transformation of the kinetic term generates higher than second-order derivatives of $`\theta `$ in the self-energy $`\widehat{\mathrm{\Sigma }}`$. In the language of Feynman diagrams this means that besides the vertices with one or two incoming $`\theta `$-lines depicted in Fig. 1, we should add all the vertices with $`n3`$ incoming $`\theta `$-lines, each associated with a $`n`$-order $`𝐤`$-derivative of $`\xi _𝐤`$, see Fig. 3. Observe that no mixing between different components $`(a=x,y,\mathrm{})`$ of the $`𝐤`$-derivative of $`\xi _𝐤`$ is produced, at least in the case of nearest-neighbor hopping.

As a consequence of (i) and (ii), the self-energy $`\widehat{\mathrm{\Sigma }}(k)`$ appearing in Eq. (8) for $`S_{eff}(\theta )`$ reads now

$`\widehat{\mathrm{\Sigma }}(k^{}k)=\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}{\displaystyle \frac{\omega _{(k^{}k)}}{2}}\theta (k^{}k)\widehat{\tau }_3+{\displaystyle \underset{n=2}{\overset{\mathrm{}}{}}}{\displaystyle \frac{i^n}{n!}}\left({\displaystyle \frac{T}{\mathrm{\Omega }}}\right)^{n1}{\displaystyle \underset{\genfrac{}{}{0pt}{}{q_1\mathrm{}q_n}{a=x_1\mathrm{}x_d}}{}}\mathrm{sin}\left({\displaystyle \frac{𝐪_{1}^{}{}_{a}{}^{}}{2}}\right)\mathrm{}\mathrm{sin}\left({\displaystyle \frac{𝐪_{n}^{}{}_{a}{}^{}}{2}}\right)\times `$
$`\theta (q_1)\mathrm{}\theta (q_n)\delta \left({\displaystyle \underset{i=1}{\overset{n}{}}}q_ik^{}+k\right){\displaystyle \frac{^n\xi _𝐤}{𝐤_a^n}}\widehat{\tau }(n)`$ (21)
where $`\delta (q)`$ is a delta-function and $`\widehat{\tau }(n)`$ corresponds to $`\widehat{\tau }_0`$ for $`n`$ odd (fermionic-like contributions) and to $`\widehat{\tau }_3`$ for $`n`$ even (bosonic-like contributions). The functions $`\mathrm{sin}(𝐪_i/2)`$ in Eq. (21) reduce to the factors $`𝐪_i/2`$ of the continuum case at small momenta, while preserving the lattice periodicity. Since we are interested in deriving the hydrodynamic action, in the following they will be approximated by $`𝐪_i/2`$.

In this section we focus on the generalization of the Gaussian action derived for the continuum model to the lattice case. To this purpose, it is sufficient to take into account the modifications described in the item (i). Indeed, at Gaussian level, only fermionic loops with two $`\theta `$ lines contribute, i.e., the same diagrams depicted in Fig. 2, which acquire however a different meaning. In particular, the first fermionic loop of Fig. 2, which corresponds to the diamagnetic contribution to the current, does not correspond to the particle density divided by the mass, resulting instead (in the specific case of a square lattice with nearest-neighbor hopping) proportional to the average kinetic energy $`E_{kin}`$ of the particles

$$𝒯=\frac{1}{d\mathrm{\Omega }}\underset{𝐤,a}{}\left(\frac{^2\xi _𝐤}{𝐤_a^2}\right)\left(1\frac{\xi _𝐤}{E_𝐤}\right)\mathrm{tanh}\left(\frac{\beta E_𝐤}{2}\right)=a^2E_{kin}$$
(22)
where we indicate with $`E_𝐤`$ the superconducting quasi-particle energy defined in App. A, and, for the sake of clarity, we made the dependence on the lattice spacing $`a`$ explicit. As a consequence, the current-current correlation function $`\mathrm{\Lambda }_{JJ}`$ becomes now $`\mathrm{\Lambda }_{JJ}^{ab}=𝒯\delta _{ab}+\mathrm{\Lambda }_{jj}^{ab}(q)`$, where it is assumed that in the definitions (104)-(106) each factor $`𝐤_a/m`$ is substituted with $`\xi _𝐤/𝐤_a`$ and the lattice dispersion $`\xi _𝐤`$ appears now in the quasiparticle energy $`E_𝐤`$. The stiffness, defined as the static limit of the transverse part of $`\mathrm{\Lambda }_{JJ}`$, is the lattice analogous of Eq. (11), i.e.

$$D(T)=\frac{1}{d\mathrm{\Omega }}\underset{𝐤,a}{}\left(\frac{^2\xi _𝐤}{𝐤_a^2}\right)\left(1\frac{\xi _𝐤}{E_𝐤}\right)\mathrm{tanh}\left(\frac{\beta E_𝐤}{2}\right)\frac{2}{d\mathrm{\Omega }}\underset{𝐤,a}{}\left(\frac{\xi _𝐤}{𝐤_a}\right)^2\left(\frac{f}{E_𝐤}\right).$$
(23)
At Gaussian level the only difference in the phase-only action between the continuum and the lattice case is in the definitions (11) or (23) of the superfluid stiffness, the remaining steps being the same. In the following, for the sake of simplicity, we write down the generic Gaussian actions as

$$S^G=\frac{1}{8}\underset{q}{}\left[\chi \mathrm{\Omega }_m^2+D(T)𝐪^2\right]\theta (q)\theta (q),$$
(24)
where $`\chi `$ indicates the static limit of the density-density correlation function $`\chi (q)\mathrm{\Lambda }_{\rho \rho }(q)`$, i.e., the compressibility, for neutral systems

$`\chi (q)`$ $`=`$ $`{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{k}{}}\mathrm{tr}[\widehat{𝒢}_0(k)\widehat{\tau }_3\widehat{𝒢}_0(k+q)\widehat{\tau }_3],`$ (25)
$`\chi `$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}\left[{\displaystyle \frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}}\mathrm{tanh}{\displaystyle \frac{\beta E_𝐤}{2}}2f^{}(E_𝐤){\displaystyle \frac{\xi _𝐤^2}{E_𝐤^2}}\right]\begin{array}{c}T0\hfill \end{array}{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}{\displaystyle \frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}}.`$ (28)
which has to be replaced by $`\chi _{LR}(𝐪)1/V(𝐪)`$ in the presence of Coulomb interaction. As a consequence, in the hydrodynamic regime, the Gaussian propagator for phase fluctuations will be generically given by

$$P(q)=\theta (q)\theta (q)=4\left[\chi (q)\mathrm{\Omega }_m^2+D𝐪^2\right]^1.$$
(29)
The stiffness which appears in the Gaussian propagator is the bare one, which we indicated with a subscript $`D_0`$ in Sec. I. However, to simplify the notation, in the following we drop the subscript, while keeping in mind that the Gaussian value $`D`$ is corrected by anharmonic phase fluctuation according to $`DD+\delta D`$, where $`\delta D`$ will be calculated within the one-loop approximation. When quantum corrections are finite, i.e. $`\delta D(T=0)0`$, the zero-temperature stiffness does not coincide with the Gaussian value, defined by the $`T=0`$ limit of the BCS formulas (11) and (23). The consequences of such an occurrence will be discussed in Sec. IV.

## IV One-loop corrections to $`D`$ within the continuum model

Before deriving the anharmonic terms in the phase-only action, we discuss the phase-fluctuation correction to the stiffness at $`T=0`$ in the quantum $`XY`$ model, which was briefly anticipated in Sec. I. The quantum $`XY`$ model includes quantum effects for the phase through a suitable modification of the Gaussian phase propagator only, while the interaction terms are powers of $`(\theta _i\theta _j)\xi _0|\stackrel{}{}\theta |`$, as in the expansion (1). Consequently $`\delta D=(D\xi _0^2/2d)(\stackrel{}{}\theta )^2`$, which generalizes Eq. (2) in $`d`$ dimensions. Since $`(\stackrel{}{}\theta )^2=𝐪^2P(q)`$, with the quantum propagator $`P(q)`$ defined in Eq. (29), we find

$`\delta D_{XY}`$ $`=`$ $`{\displaystyle \frac{D\xi _0^2}{2}}{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}𝐪^2P(q)=2D\xi _0^2{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}{\displaystyle \frac{𝐪^2}{\chi \mathrm{\Omega }_m^2+D𝐪^2}}=`$ (32)
$`=`$ $`D\xi _0^2{\displaystyle \frac{1}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪}{}}{\displaystyle \frac{𝐪^2}{\sqrt{\chi D}|𝐪|}}[1+2b(\beta c_s|𝐪|)]\begin{array}{c}T0\hfill \end{array}D\xi _0^2{\displaystyle \frac{1}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪}{}}{\displaystyle \frac{|𝐪|}{\sqrt{\chi D}}}={\displaystyle \frac{g}{\sqrt{\chi D}}}\xi _0^2D,`$
where $`b(x)=[e^x1]^1`$ is the Bose function, $`c_s=\sqrt{D/\chi }`$ is the sound velocity, $`g=𝒜_d\zeta ^{d+1}/d(d+1)`$, $`𝒜_d`$ is the solid angle in $`d`$ dimensions and $`\zeta 1/\xi _0`$ is the momentum cutoff.notacut  We considered the neutral case as an example, the generalization to the charged case being straightforward.

To compare this result with that expected for the microscopically derived phase-only action, a preliminary distinction arises between continuum and lattice models. The Gaussian phase propagator is given by Eq. (24) in both cases, and the temperature dependence of $`D(T)`$ is due to quasiparticles only. In spite of this analogy, the stiffness $`D`$ has a different meaning in the continuum or lattice case. In the former $`D=\rho _s/m`$, with the constraint $`\rho =\rho _s+\rho _n`$ at any $`T`$, whereas in the latter $`D`$ is related to the mean kinetic energy, see Eqs. (11) and (23), respectively. In particular, at $`T=0`$,

$`\mathrm{continuum}\mathrm{model}`$ $``$ $`D(0)={\displaystyle \frac{\rho _s(0)}{m}}={\displaystyle \frac{\rho }{m}},`$ (33)
$`\mathrm{lattice}\mathrm{model}`$ $``$ $`D(0)=𝒯(0)={\displaystyle \frac{1}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐤,a}{}}\left({\displaystyle \frac{^2\xi _𝐤}{𝐤_a^2}}\right)\left(1{\displaystyle \frac{\xi _𝐤}{E_𝐤}}\right).`$ (34)
The equality $`\rho _s=\rho `$ is expected to hold at $`T=0`$ in the continuum model, since Galilean invariance ensures that at $`T=0`$, where quasiparticle excitations are absent, there is no preferred reference frame in the system, and all the electrons contribute to the superfluid component.abrikosov  Thus, Eq. (33) states that the Gaussian approximations preserves that Galilean invariance of the microscopic model. Furthermore, despite the fact that, beyond the Gaussian level, anharmonic terms induce corrections to the bare Gaussian stiffness, such corrections are expected to vanish at $`T=0`$ in the continuum model, to preserve the Galilean-invariant relation $`\rho _s(0)=\rho `$. In the lattice case, instead, the presence of finite corrections to the Gaussian value $`D(0)`$ in Eq. (34) is not forbidden. This is exactly what happens in the quantum $`XY`$ model, which is a lattice model, according to Eq. (32). As a consequence, the quantum $`XY`$ model cannot be a good approximation for the phase-only model of the continuum system, since it misses the cancellation of the $`T=0`$ correction to $`D`$ due to anharmonic phase fluctuations. This requirement is instead fulfilled by expanding the microscopically derived $`S_{eff}(\theta )`$ beyond the Gaussian level, and deriving one-loop corrections to $`D`$, as we shall show in the following.

The failure of the quantum $`XY`$ model for continuum systems is not a purely formal question, as it rises the issue of its quantitative validity both in the lattice case at $`T=0`$, and in lattice and continuum systems in the classical regime, above some crossover temperature $`T_{cl}`$ (which depends on the parameter of the microscopic model). This issue is of particular relevance when high-$`T_c`$ superconductors are considered. Indeed, although in charged systems $`T_{cl}`$ can be a quite large energy scale, of the order of the plasma energy, dissipation in $`d`$-wave superconductors strongly reduces $`T_{cl}`$, as discussed in Ref. \[noi2, \], making the classical regime experimentally accessible at much lower temperature. This situation then requires a critical discussion of the corresponding classical model for phase fluctuations, with respect to the standard assumption of considering the $`XY`$ model. Here we do not discuss neither dissipative effects, nor the precise definition of $`T_{cl}`$, but we directly identify the classical regime by considering only the $`\mathrm{\Omega }_m=0`$ contribution in the Matsubara sums which define the phase-fluctuation corrections $`\delta D`$ and $`\delta D_{XY}`$. In the case of the $`XY`$ model, by considering Eq. (32), we find that in the classical regime

$$\delta D_{XY}2\frac{T}{d\xi _0^{d2}},$$
(35)
i.e., $`\delta D_{XY}(T)T`$, as already observed in Sec. I. We shall find below that also $`\delta D`$ is linear in $`T`$, but with a much smaller coefficient in the weak-to-intermediate-coupling regime (see also Ref. \[kwon, \]).

### IV.1 Classification of the diagrams

We start our analysis with the continuum neutral system, while the long-range Coulomb forces will be included later. To determine the one-loop correction $`\delta D`$ to the stiffness induced by anharmonic terms in $`S_{eff}(\theta )`$, we derive the contributions up to $`S_4`$ in Eq. (9). In general, if we separate the effective action in the Gaussian ($`S^G`$) and anharmonic ($`S_{anh}`$) part, and make a Taylor expansion of the latter, by means of the Wick theorem we can evaluate the two-field propagator as the resummation of all the connected diagrams obtained by fixing two $`\theta `$ legs and contracting the others with the Gaussian propagator $`P(q)`$,

$`\theta (q)\theta (q)`$ $`=`$ $`{\displaystyle \frac{\theta (q)\theta (q)e^{S^G(\theta )}e^{S_{anh}(\theta )}𝒟\theta }{Z}}`$
$`=`$ $`{\displaystyle \underset{n}{}}{\displaystyle \frac{(1)^n}{n!}}{\displaystyle \frac{\theta ^2e^{S^G(\theta )}[S_{anh}(\theta )]^n𝒟\theta }{Z}}={\displaystyle \underset{n}{}}{\displaystyle \frac{(1)^n}{n!}}\theta (q)\theta (q)[S_{anh}(\theta )]^n_{CON},`$
where $`\mathrm{}_{CON}`$ indicates connected diagrams only. All one-loop corrections to $`P(q)`$, i.e., diagrams with only one integration on bosonic momenta, are obtained by approximating $`_n(1)^n/n![S_{anh}]^nS_4+[S_3]^2/2`$. Moreover, since we are specifically interested in the correction to $`D`$, which is the coefficient of the gradient term in the effective Gaussian action, we select the diagrams which have at least two $`\stackrel{}{}\theta `$ lines, as depicted in Figs. 4 and 5.

To obtain one-loop corrections to $`D`$ we consider all the possible contractions of two $`\theta `$ lines in one of the fourth-order vertices in Fig. 5, or in the combination of two third-order vertices in Fig. 4. The above classification in bosonic, fermionic and mixed contributions turns out to be particularly convenient while evaluating all the resulting one-loop corrections to $`D`$. We can indeed distinguish:

* Fermionic corrections, which are given by the self-energy- and vertex-like diagrams obtained from the $`4a`$ loop in Fig. 5.
* Bosonic corrections, which arise from two loops of the $`3b`$ type and from the loop $`4e`$ of Fig. 5, see Fig. 6.
* Mixed corrections, which are the remaining combinations of third- and fourth-order diagrams.

The distinction between fermionic, bosonic, and mixed contributions to $`\delta D`$ is motivated a posteriori by the fact that all the contributions to $`\delta D`$ coming from the fermionic and mixed diagrams cancel out each other at $`T=0`$, whatever is the symmetry of the gap and both in continuum and lattice models. This result is connected to the fact that the temperature dependence of the fermionic and mixed bubbles is controlled by the thermal excitation of the quasiparticles, which are suppressed at low temperature. However, while for $`s`$-wave superconductors this suppression is exponential in temperature due to the presence of a finite gap in the excitation spectrum, for a $`d`$-wave superconductor the presence of nodal quasiparticles obliges a careful evaluation of the fermionic contribution to $`\delta D`$. For example, one can see that the fermionic contributions generated by the vertex (4a) in Fig. 5 can be written as

$$\delta D_𝒬=\frac{1}{8d}\frac{T}{\mathrm{\Omega }}\underset{𝐪,\mathrm{\Omega }_m}{}P(q)𝐪^2𝒬(𝐪,\mathrm{\Omega }_m),$$
(36)
where $`v_𝐤=\xi _𝐤/𝐤`$ and the hydrodynamic limit of the $`𝒬`$ bubble is

$$𝒬(𝐪0,\mathrm{\Omega }_m=0)=\frac{1}{\mathrm{\Omega }}\underset{𝐤}{}v_𝐤^4f^{\prime \prime \prime }(E_𝐤),$$
(37)
where $`f^{\prime \prime \prime }(x)`$ is the third-order derivative of the Fermi function. In the $`s`$-wave case the integral in Eq. (37) vanishes exponentially at $`T=0`$, while in the $`d`$-wave case, by evaluating the contribution of nodal quasiparticles to Eq. (37), which is the relevant contribution at low temperature, one finds that it diverges as $`1/T`$. This can be easily seen by considering that the density of states of nodal quasiparticles is linear in energy $`N(E)E`$, qp  so that $`𝒬(q=0)𝑑EN(E)f^{\prime \prime \prime }(E)1/T`$ as $`T0`$. It should be stressed that the divergence in the $`d`$-wave case of all fermionic loops with $`N4`$ Green functions, when evaluated at zero external momenta, is a characteristic feature of systems having Dirac quasiparticles. Nevertheless the one-loop expansion of $`D`$ is well defined: indeed by retaining all the $`𝐪,\mathrm{\Omega }_m`$ dependence on the fermionic bubble $`𝒬`$ one has no longer to deal with any divergence of fermionic corrections to $`D`$, which result even vanishing in the zero-temperature limit. samokhin

Special attention must be devoted also to the bosonic ($`d1`$ and $`d2`$) and fermionic ($`d3`$) diagrams depicted in Fig. 7. Indeed, these diagrams appear as corrections to the bare Green function in the diagrams of Fig. 2 which define the the mean-field current-current correlation function $`\mathrm{\Lambda }_{JJ}`$ in Eq. (105). In particular, $`d1`$ and $`d2`$ in Fig. 7 are phase-fluctuation corrections to the diamagnetic term $`\rho _{MF}/m`$, and $`d3`$ in Fig. 7 is a phase-fluctuation correction to the bubble $`\mathrm{\Lambda }_{jj}`$, which gives the quasiparticle depletion of $`\rho _s(T)`$ at finite temperature. One can easily see that these terms are due to the shift of $`\rho `$ with respect to the mean-field value $`\rho _{MF}`$, as induced by anharmonic phase fluctuations. On the other hand, working at fixed density, one assumes $`\rho _{MF}\rho `$ and takes into account these contributions by shifting the chemical potential $`\mu \mu +\delta \mu `$ with respect to the mean-field value, so that

$$\delta \rho =\rho \rho _{MF}=\frac{\mathrm{\Delta }F}{\mu }+\frac{\rho }{\mu }\delta \mu =0\delta \mu =\frac{1}{\chi }\frac{\mathrm{\Delta }F}{\mu },$$
(38)
where we use the definition $`\chi =\rho /\mu `$ for the compressibility, and the free-energy contribution $`\mathrm{\Delta }F`$ due to phase fluctuations is evaluated within the Gaussian approximation,

$$\mathrm{\Delta }F=\mathrm{ln}𝒟\theta e^{S^G(\theta )}=\frac{1}{2}\frac{T}{\mathrm{\Omega }}\underset{q}{}\mathrm{ln}[\chi \mathrm{\Omega }_m^2+D𝐪^2].$$
(39)
While evaluating one-loop corrections to $`D`$, one takes into account the chemical potential shift (38) by adding the contribution

$$\delta D_\mu =\frac{D}{\mu }\delta \mu ,$$
(40)
which cancels out exactly the contribution of the diagrams in Fig. (7). This can be easily understood, e.g., at $`T=0`$. Because at $`T=0`$ all the mixed diagrams vanish, only $`d1`$ and $`d2`$ contribute to $`\delta D`$. Let us consider the diagram $`d2`$. Its contribution to $`\delta D`$ is

$$\delta D_{d2}=\frac{1}{8m^2}\frac{T}{\mathrm{\Omega }}\underset{q}{}𝐪^2P(q)\left\{\frac{T}{\mathrm{\Omega }}\underset{k}{}\mathrm{tr}[\widehat{𝒢}_0(k)\widehat{\tau }_3\widehat{𝒢}_0(k)\widehat{\tau }_3]\right\}=\frac{\chi }{8m^2}\frac{T}{\mathrm{\Omega }}\underset{q}{}𝐪^2P(q).$$

Indeed, the fermionic loop of $`d2`$, given by the trace in the previous equation, corresponds to the hydrodynamic limit of the same density-density bubble $`\chi (q)`$ which appears as the coefficient of the $`\mathrm{\Omega }_m^2`$ term in the Gaussian propagator $`P(q)`$ (29), see Eqs.(25)-(28).notamat

On the other hand, while evaluating $`\mathrm{\Delta }F/\mu `$ according to Eq. (39), one finds that the $`D𝐪^2`$ term of Eq. (39) contributes with

$$\frac{\mathrm{\Delta }F}{\mu }=\mathrm{}+\frac{1}{2}\frac{T}{\mathrm{\Omega }}\underset{q}{}\frac{1}{\chi \mathrm{\Omega }_m^2+D𝐪^2}𝐪^2\frac{D}{\mu }=\mathrm{}+\frac{\chi }{8m}\frac{T}{\mathrm{\Omega }}\underset{q}{}𝐪^2P(q),$$

where we used Eq. (29) and $`D=\rho _s/m=[\rho \rho _n(T)]/m`$ so that, at $`T=0`$, $`D=\rho /m`$ and $`D/\mu =\chi /m`$. As a consequence, by means of Eqs. (38) and (40) we find that

$$\delta D_\mu =\frac{\chi }{m}\frac{1}{\chi }\frac{\mathrm{\Delta }F}{\mu }=\mathrm{}+\frac{\chi }{8m^2}\frac{T}{\mathrm{\Omega }}\underset{q}{}𝐪^2P(q),$$

which is exactly $`\delta D_{d2}`$. The same holds for $`\delta D_{d1}`$, $`\delta D_{d3}`$ and the contribution arising by differentiating $`\chi \mathrm{\Omega }_m^2`$ and $`\rho _n(T)`$ (at finite temperature) in Eq. (39), leading to a cancellation of $`\delta D_\mu `$ with the diagrams of Fig. 7.

### IV.2 One-loop corrections at $`T=0`$

The cancellation of fermionic and mixed contributions as $`T0`$ can be demonstrated with lengthly but straightforward calculations, so we will not discuss it here, where instead we analyze the behavior of the remaining bosonic contributions of Fig. 6, and in particular the role of dynamics in restoring Galilean invariance in the continuum model. We note in passing that the cancellation of the fermionic corrections has been already observed in Ref. \[kwon, \], where however the authors did not include both mixed terms and corrections to $`\delta D`$ coming from the third-order loops in $`S_3(\theta )`$. While neglecting mixed term at $`T=0`$ is not problematic, since they cancel out at $`T=0`$, neglecting the dynamic third-order vertex $`3b`$ of Fig. 4 leads to neglecting the bosonic diagrams $`2`$ and $`3`$ of Fig. 6, whose role at $`T=0`$ is crucial.

Roughly speaking, the diagrams depicted in Fig. 6, arising from the bosonic self-energy $`\widehat{\mathrm{\Sigma }}_B`$ only, correspond to

$$1\frac{\chi }{m^2}(\stackrel{}{}\theta )^2(\stackrel{}{}\theta )^2,2,3\frac{\chi ^2}{m^2}(\stackrel{}{}\theta )^2(\stackrel{}{}\theta )^2(_\tau \theta )^2,$$

where the angular brackets indicate the average with respect to the Gaussian action. As a consequence, the first diagram $`1`$, arising from a vertex $`(\stackrel{}{}\theta )^4`$ in the effective action $`S_{eff}`$, is the analogous of the term $`_{\alpha =1}^d(_\alpha \theta )^4`$ in the quantum $`XY`$ model (1). Indeed, the term $`1`$ generates the correction

$$\delta D_1=\frac{1}{4m^2}\frac{T}{d\mathrm{\Omega }}\underset{q}{}𝐪^2P(q)\chi (q),$$
(41)
where $`\chi (q)`$ is given by Eq. (25). By evaluating Eq. (41) in the hydrodynamic limit, in analogy with Eq. (32), we find that

$$\delta D_1=\frac{\chi }{m^2}\frac{T}{d\mathrm{\Omega }}\underset{𝐪,\mathrm{\Omega }_m}{}\frac{𝐪^2}{\chi \mathrm{\Omega }_m^2+D𝐪^2}=\frac{\chi }{2m^2}\frac{1}{d\mathrm{\Omega }}\underset{𝐪}{}\frac{𝐪^2}{\sqrt{\chi D}|𝐪|}[b(\beta c_s|𝐪|)b(\beta c_s|𝐪|)]\begin{array}{c}T0\hfill \end{array}\frac{g}{\sqrt{\chi D}}\frac{\chi }{2m^2},$$
(42)
which is again a finite correction, as in the quantum $`XY`$ model, see Eq. (32). By comparing Eq. (32) and (42), one finds that the strength of the term $`(\stackrel{}{}\theta )^4`$ in the continuum model is controlled by $`\chi /m^2`$, and not $`\xi _0^2D`$, as in the $`XY`$ model. Nevertheless, at $`T=0`$, the result (42) cannot be conclusive, since it would lead to $`\rho _s\rho `$ at zero temperature. Indeed, by adding the contribution of the remaining diagrams $`2`$ and $`3`$

$$\delta D_2+\delta D_3=\frac{1}{8m^2}\frac{T}{d\mathrm{\Omega }}\underset{q}{}𝐪^2\mathrm{\Omega }_m^2P^2(q)\chi ^2(q),$$

we find that the bosonic one-loop corrections to $`D`$ reduce to

$$\delta D=\frac{1}{8m^2d}\frac{T}{\mathrm{\Omega }}\underset{q}{}𝐪^2P(q)\chi (q)\left[2\mathrm{\Omega }_m^2P(q)\chi (q)\right].$$
(43)
In writing Eq. (43) we are relying on the fact that both the third- and fourth-order vertices needed to calculate $`\delta D`$ are expressed in terms of the fermionic bubble $`\chi (q)`$, Eq. (25). Indeed, both the static $`\widehat{\mathrm{\Sigma }}_B^2=((\stackrel{}{}\theta )^2/8m)\widehat{\tau }_3`$ and the dynamic $`\widehat{\mathrm{\Sigma }}_B^1=(i_\tau \theta /3)\widehat{t}_3`$ part of the “bosonic” self-energy $`\widehat{\mathrm{\Sigma }}_B`$ are tied to the Pauli matrix $`\widehat{\tau }_3`$, leading to the same fermionic bubble (25) as a coefficient of the phase-only action. As observed in Refs. \[aitchison95, ; depalo, \], this is the relevant consequence of the Galilean-invariant form of the bosonic $`\tau _3`$ term appearing in (7). Evaluating Eq. (43) in the hydrodynamic limit, we find

$`\delta D`$ $`=`$ $`{\displaystyle \frac{\chi }{2m^2}}{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}{\displaystyle \frac{𝐪^2}{\chi \mathrm{\Omega }_m^2+D𝐪^2}}\left[2{\displaystyle \frac{4\chi \mathrm{\Omega }_m^2}{\chi \mathrm{\Omega }_m^2+D𝐪^2}}\right]={\displaystyle \frac{1}{2m^2}}{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}𝐪^2{\displaystyle \frac{\mathrm{\Omega }_m^2+\omega _𝐪^2}{(\mathrm{\Omega }_m^2+\omega _𝐪^2)^2}}=`$ (44)
$`=`$ $`{\displaystyle \frac{1}{2m^2}}{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}𝐪^2{\displaystyle \frac{(i\mathrm{\Omega }_m+\omega _𝐪)^2}{[(i\mathrm{\Omega }_m)^2+\omega _𝐪^2]^2}}={\displaystyle \frac{1}{2m^2}}{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}𝐪^2{\displaystyle \frac{1}{[i\mathrm{\Omega }_m+\omega _𝐪]^2}}={\displaystyle \frac{1}{2m^2}}{\displaystyle \frac{1}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪}{}}𝐪^2b^{}(\omega _𝐪),`$
where $`\omega _𝐪=\sqrt{D/\chi }|𝐪|`$ is the sound mode, and $`b^{}(x)=\beta \mathrm{e}^{\beta x}/(\mathrm{e}^{\beta x}1)^2`$. As a consequence of the addition of the $`2`$ and $`3`$ contributions, the pole at $`i\mathrm{\Omega }_m=\omega _𝐪<0`$, responsible for a finite contribution at $`T=0`$ in Eq. (42), is canceled in favor of a double pole at $`i\mathrm{\Omega }_m=\omega _𝐪>0`$, leading to Eq. (44), where the standard Bogoljubov reduction of $`\rho _s/m`$ in a superfluid bosonic system is recognized popov . Thus, for $`T0`$ the integral in Eq. (44) vanishes, i.e., by fully including the dynamic structure of the interaction for the phase, we obtain that in a Galilean-invariant system $`\rho _s=\rho `$ at $`T=0`$. Notice that the result (44) arises from an exact compensation between static and dynamic bosonic vertices for phase fluctuations, both at Gaussian level, in determining the structure of the propagator $`P(q)`$, and beyond the Gaussian level, in weighting the relative contributions of static and dynamic parts.

### IV.3 One-loop corrections in the classical limit

According to Eq. (42), a further result of the microscopic derivation of the effective action is a reduction of the strength of the static interaction term $`(\stackrel{}{}\theta )^4`$ with respect to the $`XY`$ model, see Eq. (32). This reduction is relevant in the classical regime, where only the correction $`\delta D_{}1`$ to $`D`$ coming from the static interaction term $`1(\chi /m^2)(\stackrel{}{}\theta )^4`$ survives in Eq. (43), leading to

$$\delta D\frac{\chi }{m^2D\xi _0^2}\frac{T}{d\xi _0^{d2}},$$

which is qualitatively similar to the result of the classical $`XY`$ model, Eq. (35). To make a quantitative comparison, one can estimate the compressibility $`\chi `$ in the weak- or intermediate-coupling regime as the value of the density of states at the Fermi level,notachi  which, in the continuum model, is

$$\chi =N(\epsilon _F)=2\frac{d𝐤}{(2\pi )^d}\delta (\epsilon _F\epsilon _𝐤)=2𝒜_dmk_F^{d2},$$
(45)
where $`k_F`$ is the Fermi wave-vector. Nevertheless, also $`\rho _s`$, i.e., $`\rho `$ (if we neglect quasiparticle effects notabos ), can be defined through $`k_F`$, at least in the weak-to-intermediate-coupling regime (more precisely as long as $`T_c<E_F`$, i.e., $`U<t`$)

$$\rho _s=\rho =2_{|𝐤|<k_F}\frac{d𝐤}{(2\pi )^d}=2𝒜_d\frac{k_F^d}{d}.$$

It then follows that the ratio between the coefficient of the static interaction term $`(\stackrel{}{}\theta )^4`$ in the continuum or in the $`XY`$ model is

$$\frac{\chi }{m^2D\xi _0^2}=\frac{\chi }{m\rho _s\xi _0^2}=\frac{d}{(\xi _0k_F)^2},$$
(46)
which means that, apart from a factor $`d/2`$, we have

$$\mathrm{Classical}\mathrm{Regime}\frac{\delta D}{\delta D_{XY}}\frac{1}{\left(k_F\xi _0\right)^2}.$$
(47)
As a consequence, in the classical limit, the one-loop correction to $`D`$ derived within the effective action is smaller than the corresponding first-order correction in the $`XY`$ model as far as $`k_F^1<\xi _0`$, i.e., in the weak- and intermediate-coupling regime for the pairing interaction.

Finally, let us consider the effect of the Coulomb interaction between the electrons. As we showed in Sec. II.2, the presence of a density-density interaction in the microscopic Hamiltonian implies the dressing of the coefficients of the Gaussian action by the RPA series of the Coulomb potential $`V(𝐪)`$. As a consequence, in the hydrodynamic limit, the bare compressibility $`\chi `$ is replaced by the dressed bubble $`\chi _{LR}(𝐪)1/V(𝐪)`$, and the sound mode $`\omega _𝐪`$ of Eq. (44) is converted into the plasma mode $`\omega _𝐪^P`$ of the $`d`$-dimensional system (see Eq. (18) for the three-dimensional case). In the same way, when deriving the anharmonic terms in $`𝒮_{eff}`$, we must now include the RPA density fluctuations in all the vertices for $`S_3`$ and $`S_4`$. The one-loop corrections to $`D`$ are formally identical to Eqs. (43)-(44), with $`\chi \chi _{LR}`$. Thus at $`T=0`$ we recover again the cancellations (44) of the bosonic diagrams, with $`\omega _𝐪\omega _𝐪^P`$. At the same time, since $`\chi _{LR}(𝐪)`$ vanishes as $`𝐪0`$, the classical ($`\mathrm{\Omega }_n=0`$) term in Eq. (43), i.e., $`\delta D_1`$ given by Eq. (42), reads

$$\delta D_{LR}\frac{T}{m^2dD}\frac{1}{\mathrm{\Omega }}\underset{𝐪}{}\chi (𝐪)=\frac{T}{m^2dD\lambda e^2}\frac{𝒜_d\zeta ^{2d1}}{2d1},$$

whereas the quantum $`XY`$ model leads to the same result (35) of the neutral case. Thus we can roughly estimate

$$\mathrm{Classical}\mathrm{Regime}\mathrm{LR}\frac{\delta D_{LR}}{\delta D_{XY}}\frac{\epsilon _F}{\epsilon _C}\left(\frac{1}{k_F\xi _0}\right)^{d+1},$$

where $`\epsilon _F=k_F^2/2m`$ is the Fermi energy and $`\epsilon _C=\lambda e^2k_F`$ is a characteristic Coulomb energy scale. As a consequence, while within the $`XY`$ model the Coulomb interaction modifies only the low-temperature behavior of $`\rho _s`$, within the continuum model it affects also the high-temperature classical regime.

## V One-loop corrections to $`D`$ within a lattice model

To extend the previous results to a lattice model, we need to take into account the discussion of Sec. III, in which we derived the phase-only action starting from a nearest-neighbor tight-binding model for electrons in a lattice.

At Gaussian level, the item (i) of Sec. III allows us to define the superfluid stiffness of a lattice model according to Eq. (23). In particular, a factor $`\frac{1}{8}\mathrm{\Lambda }_a`$, with

$$\mathrm{\Lambda }_a=\frac{^2\xi _𝐤}{𝐤_a^2}.$$

appears in the fermionic bubble which defines the diamagnetic term, i.e., the coefficient of the $`(\stackrel{}{}\theta )^2`$ term in the Gaussian action for phase fluctuations. As a consequence, $`D`$ is controlled at $`T=0`$ by the mean kinetic energy $`𝒯`$, Eq. (22), instead of $`\rho /m`$ of the continuum case. In deriving anharmonic terms in $`S_{eff}(\theta )`$, we must take into account item (i) in the proper definition of the fermionic bubbles of Figs. 4 and 5, which give the vertices for the interacting phase-only model. Moreover, we need to include also the fourth-order vertex depicted in Fig. 8, which is peculiar of the lattice case, and arises from the $`n=4`$ term of Eq. (21) (see also Fig. 3).

At this point the calculation of one-loop correction to $`D`$ follows the same steps described in Sec. IV. The cancellation of fermionic and mixed one-loop corrections to $`D`$ at $`T=0`$ still holds, thus we can again focus on bosonic corrections only. As far as the corrections in Fig. 7 are concerned, they are no more canceled out by the contribution $`\delta D_\mu `$ due to the chemical-potential shift beyond the Gaussian level. In particular, while discussing in detail in Sec. IV the case of the term $`d2`$, we stressed that in the continuum model such a cancellation holds because the fermionic bubble of diagram $`d2`$ is the same bubble which one obtains by differentiating the diamagnetic part of $`D`$ with respect to $`\mu `$, so that $`D/\mu =\chi /m`$. Indeed, since in the continuum case the factor $`\mathrm{\Lambda }_\alpha =1/m`$ associated to the $`(\stackrel{}{}\theta )^2`$ insertion is a constant, one obtains the same two-line fermionic loop both by expanding $`S_{eff}`$ up to $`S_4`$, leading to $`d2`$, and by differentiating the single-line diagram which defines $`\rho `$. However, in the lattice case the factor $`\mathrm{\Lambda }_\alpha `$ is not a constant. As a consequence the diagram obtained by differentiating the diamagnetic term $`𝒯D(0)`$ with respect to $`\mu `$ has one insertion of $`\mathrm{\Lambda }_\alpha `$, while the the diagram $`d2`$, obtained by expanding $`S_{eff}`$ according to item (i), has a $`\mathrm{\Lambda }_\alpha `$ factor for each $`\theta `$-vertex (diagram $`_2`$ in Fig. 9). The same holds for the $`d1`$ diagram.

As a consequence, in the lattice case, the contributions to $`D`$ come from the bosonic diagrams depicted in Fig. 9, and from the correction $`\delta D_\mu `$ induced by the chemical-potential shift, which do not cancel anymore the diagrams $`_1`$ and $`_2`$ in Fig. 9. While evaluating $`\delta D_2`$ we must consider that we have $`2(d1)`$ bubbles with $`\mathrm{\Lambda }_\alpha \mathrm{\Lambda }_\beta `$ in the expansion of $`S_{anh}`$ in $`d`$ dimensions, so that

$`\delta D_1`$ $`=`$ $`{\displaystyle \frac{1}{8}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{q}{}}2P(q)\mathrm{\Omega }_m^2\eta _{\rho E}(q),`$ (48)
$`\delta D_2`$ $`=`$ $`{\displaystyle \frac{1}{8d}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{q}{}}P(q)𝐪^2[\chi _{EE}+(d1)\stackrel{~}{\chi }_{EE}],`$ (49)
$`\delta D_3`$ $`=`$ $`{\displaystyle \frac{1}{8d}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{q}{}}P(q)𝐪^2\stackrel{~}{𝒯},`$ (50)
$`\delta D_4`$ $`=`$ $`{\displaystyle \frac{1}{8d}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{q}{}}2P(q)𝐪^2\chi _{EE}(q),`$ (51)
$`\delta D_5+\delta D_6`$ $`=`$ $`{\displaystyle \frac{1}{8d}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{q}{}}P^2(q)𝐪^2\mathrm{\Omega }_m^2\chi _{\rho E}(q),`$ (52)
where $`\eta _{\rho E}(q)`$ is the fermionic bubble with three $`\widehat{𝒢}_0`$ lines and one insertion of $`\mathrm{\Lambda }_\alpha `$, and the $`\chi _{ab}(q)`$ bubbles correspond to the insertion of one ($`\chi _{\rho E}`$) or two equal/different ($`\chi _{EE}`$/$`\stackrel{~}{\chi }_{EE}`$) factors $`\mathrm{\Lambda }_\alpha `$ in the two-line bubble. As usual, when the dependence on $`q`$ is not explicitly indicated, we are considering the $`𝐪0,\mathrm{\Omega }_m=0`$ (static) limit. Notice that, since $`(^4\xi _𝐤/𝐤_\alpha ^4)=a^2(^2\xi _𝐤/𝐤_\alpha ^2)`$, the $`\stackrel{~}{𝒯}`$ which appears in $`\delta D_3`$ is proportional to the mean kinetic energy defined by Eq. (22), $`\stackrel{~}{𝒯}=a^2𝒯=a^4E_{kin}`$ (the dependence on the lattice parameter $`a`$ is made explicit, for the sake of clarity in the forthcoming discussion).

According to Eq. (38), the chemical-potential shift is

$$\delta \mu =\frac{1}{\chi }\frac{\mathrm{\Delta }F}{\mu }=\frac{1}{8\chi }\frac{T}{\mathrm{\Omega }}\underset{q}{}P(q)[\mathrm{\Omega }_m^2\eta (q)+𝐪^2\chi _{\rho E}],$$

where $`\eta (q)`$ is the three-line bubble without factor $`\mathrm{\Lambda }_\alpha `$, $`\chi `$ is the compressibility, which carries no $`\mathrm{\Lambda }_\alpha `$ insertion, and we used $`\chi _{\rho E}=D/\mu `$. By means of Eq. (40), we obtain

$$\delta D_\mu =\frac{D}{\mu }\delta \mu =\frac{\chi _{\rho E}}{\chi }\frac{1}{8}\frac{T}{\mathrm{\Omega }}\underset{q}{}P(q)[\mathrm{\Omega }_m^2\eta (q)+𝐪^2\chi _{\rho E}].$$
(53)
Finally, the overall one-loop correction to $`D`$ coming from Eqs. (48)-(52) and (53) is

$`\delta D=`$ $``$ $`{\displaystyle \frac{1}{8d}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{q}{}}P(q)𝐪^2\left[2\chi _{EE}(q)+\chi _{EE}+(d1)\stackrel{~}{\chi }_{EE}+\stackrel{~}{𝒯}d{\displaystyle \frac{\chi _{\rho E}^2}{\chi }}\right]`$ (54)
$``$ $`P(q)\mathrm{\Omega }_n^2\left[\chi _{\rho E}^2(q)𝐪^2P(q)2d\eta _{\rho E}(q)+d{\displaystyle \frac{\chi _{\rho E}}{\chi }}\eta (q)\right].`$
To make a comparison with the continuum case, Eq. (44), we evaluate the above equation in the hydrodynamic limit. This corresponds to calculate the $`\chi _{ab}`$ bubbles at zero incoming momentum and by carefully evaluating the $`𝐪0`$ limit for the terms which involve $`\eta _{\rho E}(q)`$ and $`\eta (q)`$. Indeed, by simply calculating $`_{\mathrm{\Omega }_m}P(q)\mathrm{\Omega }_m^2\eta (q)\eta _{\mathrm{\Omega }_m}P(q)\mathrm{\Omega }_m^2`$ one finds that $`_{\mathrm{\Omega }_m}P(q)\mathrm{\Omega }_m^2\mathrm{}`$. Instead, by retaining the $`\mathrm{\Omega }_m`$ dependence in the fermionic bubble one ensures the convergence of the Matsubara sum, calculating the $`𝐪0`$ limit afterward. Thus at $`T=0`$ limit we obtain that the overall correction to $`D`$ is finite in the lattice case,

$$\delta D=\frac{g}{4\sqrt{\chi D}}\left[3\chi _{EE}+(d1)\stackrel{~}{\chi }_{EE}+\stackrel{~}{𝒯}(d+2)\frac{\chi _{\rho E}^2}{\chi }\right]\frac{g(d+1)}{\zeta \chi }\left[\stackrel{~}{\eta }_{\rho E}\frac{\chi _{\rho E}}{\chi }\stackrel{~}{\eta }\right],$$
(55)
according to the following limiting values

$`{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}𝐪^2P(q)\begin{array}{c}T0\hfill \end{array}{\displaystyle \frac{2g}{\sqrt{\chi D}}},`$ (58)
$`{\displaystyle \frac{T}{d\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}𝐪^2\mathrm{\Omega }_m^2P^2(q)\begin{array}{c}T0\hfill \end{array}{\displaystyle \frac{4g}{\chi \sqrt{\chi D}}},`$ (61)
$`{\displaystyle \frac{T}{4\mathrm{\Omega }}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}\mathrm{\Omega }_m^2P(q)\eta _{\rho E}(q)\begin{array}{c}T0\hfill \end{array}{\displaystyle \frac{\stackrel{~}{\eta }_{\rho E}}{\chi }}{\displaystyle \frac{g(d+1)}{\zeta }}.`$ (64)
The $`T=0`$ value of $`\chi `$ is given by Eq. (28). The others bubbles are defined as

$`\chi _{\rho E}`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}\left({\displaystyle \frac{^2\xi _𝐤}{𝐤_\alpha ^2}}\right){\displaystyle \frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}},`$ (65)
$`\chi _{EE}`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}\left({\displaystyle \frac{^2\xi _𝐤}{𝐤_\alpha ^2}}\right)^2{\displaystyle \frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}},`$ (66)
$`\stackrel{~}{\chi }_{EE}`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}\left({\displaystyle \frac{^2\xi _𝐤}{𝐤_\alpha ^2}}\right)\left({\displaystyle \frac{^2\xi _𝐤}{𝐤_\beta ^2}}\right){\displaystyle \frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}},(\alpha \beta ),`$ (67)
and

$`\stackrel{~}{\eta }`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}{\displaystyle \frac{\mathrm{\Delta }_𝐤^2\xi _𝐤}{E_𝐤^4}},`$ (68)
$`\stackrel{~}{\eta }_{\rho E}`$ $`=`$ $`{\displaystyle \frac{1}{\mathrm{\Omega }}}{\displaystyle \underset{𝐤}{}}\left({\displaystyle \frac{^2\xi _𝐤}{𝐤_\alpha ^2}}\right){\displaystyle \frac{\mathrm{\Delta }_𝐤^2\xi _𝐤}{E_𝐤^4}}.`$ (69)
Thus in the lattice case the $`T=0`$ mean-field value of the superfluid density is corrected by one-loop phase fluctuations according to Eq. (55). It is interesting to observe that the continuum case is recovered from Eq. (55) by performing the limit $`a0,t\mathrm{}`$, keeping $`2a^2t=1/m`$ finite. Indeed, since $`\xi _𝐤=2tϵ_𝐤\mu `$, with $`ϵ_𝐤=_a\mathrm{cos}(ak_a)`$, $`\mathrm{\Lambda }_\alpha =(^2\xi _𝐤/k_a^2)=2ta^2ϵ_𝐤1/m`$ in this limit, and according to Eqs. (65)-(69), one finds $`\chi _{EE},\stackrel{~}{\chi }_{EE}\chi /m^2`$, $`\chi _{\rho E}\chi /m`$, $`\stackrel{~}{\eta }_{\rho E}\stackrel{~}{\eta }/m`$, while $`\stackrel{~}{𝒯}=a^2D0`$. Thus in this limit, Eq. (55) correctly reduces to the continuum (Galilean-invariant) result $`\delta D=0`$ of Eq. (44).

We now compare $`\delta D`$ and $`\delta D_{XY}`$ in the lattice case. Differently from the continuum case, this comparison is now meaningful both in the classical limit and at $`T=0`$, since now we find finite corrections to $`D`$. By comparing Eq. (55) and Eq. (32) at $`T=0`$, we see that roughly speaking the strength $`\xi _0^2D`$ of the interaction in the quantum $`XY`$ model is given by the quantity in square brackets in the first line of Eq. (55), i.e., it is of order $`\chi _{EE},\stackrel{~}{\chi }_{EE},\chi _{\rho E}^2/\chi `$. Following the same arguments which led us to Eq. (46), we can estimate

$$\chi _{EE},\stackrel{~}{\chi }_{EE},\frac{\chi _{\rho E}^2}{\chi }\frac{\xi _0^2D}{(k_F\xi _0)^2},$$

where $`\chi _{EE},\stackrel{~}{\chi }_{EE},\chi _{\rho E}^2/\chi `$ have the same role as $`\chi /m^2`$ in Eq. (46), and $`k_F1/a`$ in the lattice case. However, in the lattice case the dynamic contribution controlled by $`\chi _{\rho E}/\chi ^2`$ does not cancel with the static contribution, and both contribute to $`\delta D`$ at $`T=0`$. Analogously, one finds that $`𝒯=a^2DD/k_F^2=\xi _0^2D/(k_F\xi _0)^2`$. Moreover,

$$\stackrel{~}{\eta }\frac{N(\epsilon _F)}{1+(\mu /\mathrm{\Delta })^2}\frac{N(\epsilon _F)}{(k_F\xi _0)^2}$$

where $`N(\epsilon _F)`$ is evaluated in Eq. (45), $`\mathrm{\Delta }`$ is the maximum gap value, and we used the fact that at weak coupling $`k_F\xi _0k_Fv_F/\mathrm{\Delta }=\mu /\mathrm{\Delta }1`$.noi4  Taking into account the prefactor $`g/\zeta \chi `$ in the second line of Eq. (55), we have

$$\frac{g}{\zeta \chi }\left[\stackrel{~}{\eta }_{\rho E}\frac{\chi _{\rho E}}{\chi }\stackrel{~}{\eta }\right]\delta D_{XY}\frac{1}{(k_F\xi _0)^3}$$
(70)
where $`\stackrel{~}{\eta }_{\rho E}\stackrel{~}{\eta }/m`$. Thus, the contributions coming from the three-line fermionic bubbles are subleading, and at $`T=0`$ Eq. (55) leads to an estimate

$$\frac{\delta D}{\delta D_{XY}}\frac{1}{\left(k_F\xi _0\right)^2}.$$
(71)
of the same order of Eq. (47), within a numerical factor. The presence of Coulomb forces does not change qualitatively this conclusion, and introduces only minor quantitative corrections. This is due to the fact that the RPA expressions for $`\chi _{EE},\stackrel{~}{\chi }_{EE}`$ have a finite limit for $`𝐪0`$, contrary to the continuum case,

$$\chi _{EE}^{LR}=\chi _{EE}\frac{\chi _{\rho E}^2V(𝐪)}{1+V(𝐪)\chi }\left(\chi _{EE}\frac{\chi _{\rho E}^2}{\chi }\right)+\frac{\chi _{\rho E}^2}{\chi }\frac{|𝐪|^{d1}}{\lambda e^2}.$$
(72)
When $`\delta D`$ is evaluated with the long-range phase propagator $`P(q)`$, one sees that the leading terms in powers of $`1/(\xi _0k_F)`$ come from the finite $`𝐪=0`$ limit of $`\chi _{EE},\stackrel{~}{\chi }_{EE}`$ in Eq. (72), and from the $`\stackrel{~}{𝒯}`$ term, all of order $`1/(\xi _0k_F)^2`$ with respect to the coefficient $`\xi _0^2D`$ of the $`XY`$ model, whereas

$$\chi _{\rho E}^{LR}=\frac{\chi _{\rho E}}{1+V(𝐪)\chi _{\rho E}}\begin{array}{c}𝐪0\hfill \end{array}\frac{1}{V(𝐪)},$$

so that dynamical terms are subleading at $`T=0`$ in the long-range case. This is clearly true as far as the constant in the parenthesis of Eq. (72) is a finite quantity. Indeed, when it vanishes by approaching the continuum limit, all (static and dynamics) terms are of the same order, canceling each other as discussed above.

Finally, the same result (71) holds in the classical limit, both for the neutral and for the charged system. As we discussed above, at high temperature only the static correction to $`D`$ contributes, with a coefficient controlled in the neutral case by $`\chi _{EE},\stackrel{~}{\chi }_{EE},\stackrel{~}{𝒯}`$, or by their $`𝐪=0`$ finite limit in the charged case, which play the same role (and have the same estimate) as $`\chi /m^2`$ in the continuum case.

As far as the effect of bosonic degrees of freedom is concerned, our result (71) suggests that the $`XY`$ model can eventually lead to an overestimate of phase-fluctuation effects with respect to the result obtained by means of the microscopically derived effective action. Indeed, according to Eq. (71) the depletion of $`D`$ due to phase fluctuations can be quite small at weak and intermediate coupling (particularly in the BCS limit) both in the quantum and in the classical regime.

On the other hand, as far as the temperature dependence of $`\delta D`$ for temperatures $`0<T<T_{cl}`$ is considered, one should take into account also the effects of fermionic and mixed contributions. As we discussed above, quasiparticles contribute to the depletion of $`D`$ at finite temperature already at Gaussian level, with a temperature dependence determined by the symmetry of the order parameter, which is linear in $`T`$ in $`d`$-wave superconductors. One would expect that at $`T>0`$ the one-loop corrections to $`D`$ coming from fermionic and mixed diagrams correct quantitatively but not qualitatively the $`D(T)`$ dependence due to quasiparticles at Gaussian level. However, a detailed analysis of the temperature dependence of $`D(T)`$ as arising from bosonic, mixed and fermionic contribution at low temperature is a much harder task. From one side, one has to face the difficulty of evaluating the $`q=0`$ limit of these expressions, as we discussed in Sec. IV A. On the other side, at finite temperature, non-local terms appear when deriving the phase-only action, making the hydrodynamic expansion not well-defined. These are the so-called Landau singularities, arising from the non analytic behavior at small $`q`$ of the $`\mathrm{\Lambda }_{JJ}`$.sharapov  In any case the analysis of these issues requires separate investigation and is not addressed in the present work.

## VI The strong-coupling regime

The results (71) for the lattice case was derived by assuming $`1/(k_F\xi _0)`$ as a small expansion parameter. In particular, this assumption allowed us to estimate the fermionic bubble according to Eq. (45), i.e., $`\chi N(\epsilon _F)`$, and to neglect the contribution to $`\delta D`$ of the $`\stackrel{~}{\eta }_{ab}`$ bubbles, see Eq. (70). One should expect that such approximations are reliable in the weak- or intermediate-coupling regime, where the correlation length $`\xi _0`$ is larger than the lattice parameter $`a`$, which sets the scale of the Fermi wave-vector $`k_F1/a`$. However, in the strong-coupling regime for the pairing interaction several modifications should be included. The fermionic bubbles can no longer be estimated as the density of states at the Fermi level; the correlation length is generically expected to be of the same order of magnitude as the lattice parameter; one should include both the RPA fluctuations, also induced by $`U`$, in the particle-hole channel, and the fluctuation of $`|\mathrm{\Delta }|(U/2)\sqrt{\rho (2\rho )}`$, which fluctuates because $`\rho `$ fluctuates.depalo  This last issue can be better understood by analogy with the superfluid bosonic systems, where the order parameter is directly the square-root of the density, and one expects that, in the strong-coupling regime, the fermionic system can be mapped onto an effective bosonic system, where Cooper pairs act as bosons with a weak residual repulsion. The above points are intimately connected, as we already discussed in Ref. noi4, . There, the behavior of the correlation length, properly defined as the characteristic length scale for the fluctuations of the modulus of the superconducting order parameter, was analyzed in the strong-coupling regime. The main finding is that $`\xi _0`$ attains a value of the order of the lattice parameter $`a`$ at finite densities, and diverges as $`1/\sqrt{\rho }`$ in the low-density regime $`\rho 0`$, recovering the well-known behavior of a bosonic system.bos

In the following we consider the $`s`$-wave case, which allows us for a transparent analytical treatment. We first address the neutral case, i.e. the model (3)-(4) with $`w(𝐤)=1`$ (negative-$`U`$ Hubbard model), and then we extend our results to the charged case. In the presence of $`d`$-wave pairing the analysis becomes more difficult, even though the physics is not expected to be different from the $`s`$-wave case, as we briefly discuss at the end of this section.

As we said at the beginning of Sec. II.1, the Hubbard-Stratonovich decoupling is a useful tool which allows us to substitute the microscopic model with an effective action written in terms of the relevant collective degrees of freedom. However, the choice of the relevant variables depends on the physical properties of interest and on the energetic scales involved in the problem. For example, the negative-$`U`$ Hubbard model can be decoupled both in the particle-particle channel, as we did, and in the particle-hole channel. The latter decoupling would lead to an effective action depending on the particle density as well, $`S_{eff}=S_{eff}(\theta ,\delta |\mathrm{\Delta }|,\delta \rho )`$.notaspin  The resulting action is formally the same that we obtained by decoupling the additional interaction Hamiltonian (12) in the particle-hole channel, but with $`V(𝐪)`$ replaced by the same $`U/2`$ which appears in the Hubbard term. The relevance of density fluctuations within the Hubbard model at strong coupling is addressed in Refs. \[depalo, ; noi4, \]. In deriving $`S_{eff}`$ up to Gaussian level a term proportional to the product of amplitude and density fluctuation, $`\delta |\mathrm{\Delta }|\delta \rho `$, appears. The coefficient of such a term in the Gaussian action is vanishingly small at weak coupling. This means that density and modulus fluctuations decouple at weak and intermediate coupling, and one can safely neglect modulus fluctuations without loosing information on the contribution of the density fluctuations. As a consequence, one can derive the phase-only action according to the procedure described in the previous sections. However, at strong coupling the situation is more involved: as it is shown in Refs. \[depalo, ; noi4, \], by increasing $`U`$ the density-modulus coupling become sizable, and one finds that the density of particles and the modulus of the order parameter experience the same fluctuations. Since the phase mode has the same behavior of the density mode, in the strong-coupling regime we should also take into account amplitude fluctuations to recover a consistent description of the phase fluctuations.

In Sec.II.2 we found that in deriving the effective action for the phase of the order parameter the inclusion of density fluctuations reflects in the replacement of the “bare” coefficients $`\mathrm{\Lambda }`$ of the Gaussian action with their counterparts $``$ evaluated within RPA resummation in the particle-hole channel. Beyond Gaussian level one expects the RPA dressing of the various fermionic loops which appear as coefficients of the effective action. As far as the $`\delta D`$ correction in Eq. (55) is concerned, this leads, e.g., to substitute the $`\chi _{ab}`$ fermionic bubbles with the RPA resummation with the potential $`U/2`$ of the corresponding density or current irreducible correlation functions $`\chi _{ab}^{irr}(q)`$. We are then left with the problem of including modulus fluctuations. In principle, to take them into account one should derive the Gaussian effective action $`S_G(\theta ,\delta |\mathrm{\Delta }|,\delta \rho )`$ by including all the fields and then integrate out modulus fluctuations as well, see, e.g., Refs. \[depalo, ; noi4, \]. However, the contribution of amplitude fluctuations to the static long-wavelength limit of the fermionic bubbles $`\chi _{ab}^{irr}(𝐪0,\omega _n=0)`$ can be easily derived in a different way. Let us consider, e.g., the density-density correlation function $`\chi ^{irr}`$. Since $`\chi ^{irr}=\rho /\mu `$ (see also Sec. IV), we can deduce it from the mean-field expression (107) of the density of particles, by taking into account the dependence on $`\mu `$ of the gap amplitude $`|\mathrm{\Delta }|(0)\mathrm{\Delta }_0`$ itself, which appears in the definition of $`E_𝐤=\sqrt{\xi _𝐤^2+\mathrm{\Delta }_0^2}`$ (see also Eqs. (102)). At weak coupling $`\mathrm{\Delta }_0/\mu `$ can be neglected, and one finds the definition (28) of $`\chi ^{irr}`$. At strong coupling, where modulus and density fluctuations are proportional, by working at fixed number of particles the chemical-potential variations needed to preserve the particle number reflect into amplitude variations, and $`\mathrm{\Delta }_0/\mu `$ gives a significant contribution to the irreducible bubbles. We then have (at $`T=0`$):

$$\chi ^{irr}=\frac{\rho _{MF}}{\mu }=\frac{1}{\mathrm{\Omega }}\underset{𝐤}{}\frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}+\frac{1}{\mathrm{\Omega }}\underset{𝐤}{}\frac{\xi _𝐤}{E_𝐤^3}\mathrm{\Delta }_0\frac{\mathrm{\Delta }_0}{\mu },$$

where $`\mathrm{\Delta }_0/\mu `$ can be obtained from the self-consistency saddle-point equation

$$\frac{2}{U}=\frac{1}{\mathrm{\Omega }}\underset{𝐤}{}\frac{1}{E_𝐤},$$
(73)
by differentiating both sides with respect to $`\mu `$,

$$\mathrm{\Delta }_0\frac{\mathrm{\Delta }_0}{\mu }=\frac{_𝐤(\xi _𝐤/E_𝐤^3)}{_𝐤(1/E_𝐤^3)},$$

so that

$$\chi ^{irr}=\frac{1}{\mathrm{\Omega }}\underset{𝐤}{}\frac{\mathrm{\Delta }_𝐤^2}{E_𝐤^3}+\frac{1}{\mathrm{\Omega }}\frac{\left[_𝐤(\xi _𝐤/E_𝐤^3)\right]^2}{_𝐤(1/E_𝐤^3)}.$$
(74)
In the limit of strong pairing interaction, we can evaluate the previous bubble for $`t/U1`$, by taking into account the strong-coupling solution of the coupled saddle-point equations for the amplitude $`\mathrm{\Delta }_0`$ Eq. (73) and for the particle number Eq. (107),

$`\mathrm{\Delta }_0`$ $`=`$ $`{\displaystyle \frac{U}{2}}(1d\alpha ^2)\sqrt{1\delta ^2}+𝒪\left({\displaystyle \frac{t^4}{U^3}}\right),`$
$`\mu `$ $`=`$ $`{\displaystyle \frac{U}{2}}\delta (1+2d\alpha ^2)+𝒪\left({\displaystyle \frac{t^4}{U^3}}\right),`$
where we introduced the short-hand notation $`\delta 1\rho `$ and the small parameter $`\alpha 2t/U`$. We then obtain

$$\chi ^{irr}=\frac{2}{U}\left(12d\alpha ^2\right).$$

Thus, including also the RPA resummation of $`\chi ^{irr}`$ in the particle-hole channel, we finally get

$$\chi =\frac{\chi ^{irr}}{1(U/2)\chi ^{irr}}=\frac{1}{2d\alpha t}.$$
(75)
The expression (75) for the density-density electronic correlation function at $`q=0`$ is the same which one would obtain by explicitly integrating out both the $`\delta \rho `$ and the $`\delta |\mathrm{\Delta }|`$ field in the action $`S(\theta ,\delta |\mathrm{\Delta }|,\delta \rho )`$. The latter procedure was indeed adopted in Ref. \[noi4, \] while deriving the propagator for the Hubbard-Stratonovich field associated to the density. The result (75) shows a large increase of $`\chi `$ at large $`U`$: this could be expected, since the electrons are strongly paired, naturally increasing the compressibility of the system.

As far as the other fermionic bubbles are concerned, we follow the same steps. In analogy with the previous discussion

$`\chi _{\rho E}^{irr}`$ $``$ $`{\displaystyle \frac{1}{d}}{\displaystyle \frac{D}{\mu }},`$
$`\chi _{EE}^{irr}+(d1)\stackrel{~}{\chi }_{EE}^{irr}`$ $``$ $`t{\displaystyle \frac{D_{EE}}{t}},`$
where the stiffness is defined in Eq. (23), and

$$D_{EE}=\frac{1}{2t\mathrm{\Omega }}\underset{𝐤}{}\left(\frac{^2\xi _𝐤}{𝐤_\alpha ^2}\right)\left(1\frac{\xi _𝐤}{E_𝐤}\right),$$

at $`T=0`$. By taking the derivative with respect to the hopping $`t`$ one includes also $`\mathrm{\Delta }_0/t`$, as obtained again from Eq. (73). In the strong-coupling limit

$`\chi _{\rho E}^{irr}`$ $`=`$ $`4\alpha ^2\delta ,`$ (76)
$`\chi _{EE}^{irr}`$ $`=`$ $`2t\alpha (1\delta ^2),`$ (77)
$`\stackrel{~}{\chi }_{EE}^{irr}`$ $`=`$ $`𝒪(t\alpha ^3).`$ (78)
By introducing the dressed potential $`\stackrel{~}{U}=(U/2)/[1(U/2)\chi ^{irr}]=U/4d\alpha ^2`$, in analogy with Eq. (75), we find

$`\chi _{\rho E}=\chi _{\rho E}^{irr}[1+\chi _{\rho \rho }^{irr}\stackrel{~}{U}]={\displaystyle \frac{2\delta }{d}},`$ (79)
$`\chi _{EE}=\chi _{EE}^{irr}+(\chi _{\rho E}^{irr})^2\stackrel{~}{U}=2t\alpha (1\delta ^2)+{\displaystyle \frac{8t\alpha }{d}}\delta ^2,`$ (80)
$`\stackrel{~}{\chi }_{EE}=\stackrel{~}{\chi }_{EE}^{irr}+(\chi _{\rho E}^{irr})^2\stackrel{~}{U}={\displaystyle \frac{8t\alpha }{d}}\delta ^2.`$ (81)
The inclusion of the modulus fluctuations always adds a second term to the weak-coupling definition of the irreducible bubbles, see, e.g., Eq. (74). This additional term is never of lower order in $`\alpha `$ with respect to the term obtained with the weak-coupling definitions (65)-(67). Since one expects that the same holds for the $`\eta `$ bubbles, let us first estimate the RPA resummation of the weak-coupling definitions (68) and (69), which do not include modulus-fluctuation contribution. The resulting corrections to $`D`$, arising from the second line of Eq. (55), are subleading with respect to the contribution coming from the $`\chi _{ab}`$ bubbles, first line of Eq. (55). Since modulus-fluctuation contributions do not change this result, we can neglect the correction $`\delta D`$ coming from the $`\eta `$ bubbles.

By using the expressions (79)-(81), the strong-coupling expression

$$D=a^{2d}2t\alpha (1\delta ^2),$$

and the fact that $`\stackrel{~}{𝒯}=a^2D`$, we find that Eq. (55), at large $`U/t`$, reduces to

$$\delta D=\frac{g}{\sqrt{\chi D}}8a^{4d}t\alpha (1\delta ^2)=\frac{4g}{\sqrt{\chi D}}a^2D,$$
(82)
where we reintroduced explicitly the dependence on the lattice parameter $`a`$, to allow one for a direct comparison with the result of the $`XY`$ model, Eq. (32). We then obtain that

$$\delta D=\delta D_{XY}\left(\frac{a}{\xi _0}\right)^2.$$
(83)
By taking into account the previous results for $`\xi _0`$,noi4  we find that, since at finite densities $`\xi _0a`$, the $`XY`$ model appears as the proper “strong-coupling” effective quantum model for phase fluctuations in a lattice. A noticeable exception is found in the extreme low-density regime, where $`\xi _0`$ diverges as $`1/\sqrt{\rho }`$. Moreover, since $`g1/\xi _0^{d+1}`$, Eq. (82) shows that in the low-density limit the relative superfluid-stiffness correction $`\delta D/D`$ vanishes as $`\rho ^{d/2}`$. Observe that instead the $`XY`$ model would lead to a relative correction $`\delta D_{XY}/D\rho ^{d/21}`$ which does not vanish at $`\rho =0`$ in $`d=2`$.

In the classical regime, as we already discussed in the weak-coupling case, only static anharmonic terms survive, leading to

$$\frac{\delta D}{\delta D_{XY}}=\left[1+\frac{2\delta ^2}{d(1\delta ^2)}\right]\left(\frac{a}{\xi _0}\right)^2.$$

Since the quantity in parenthesis behaves as the microscopically derived $`(\xi _0/a)^2`$ at all densities,noi4  up to a numerical constant of order one, we find that the in the classical regime the $`XY`$ model turns out to be appropriate, regardless of the density. Observe that the numerical difference comes out from the fact that the $`XY`$ model attributes a somewhat arbitrary coefficient to the anharmonic term $`(\stackrel{}{}\theta )^4`$. If one assumed in Eq. (1) higher order harmonics as

$$\frac{D_n}{4}\underset{<i,j>}{}(1\mathrm{cos}\theta _{ij})^n,(n2),$$

which do not modify the Gaussian term, by adjusting the $`D_n`$ coefficients one could reproduce the correct $`(\stackrel{}{}\theta )^{2n}`$ term as it arises in the microscopic expansion. In such a case the classical limit of this extended $`XY`$ model and of the microscopic model would lead to exactly the same result. As it has been noted recently in Ref. kim , the mapping of the classical phase-only action of a lattice systems on an extended $`XY`$ model can also lead in two dimensions to an enhanced fluctiation region near the Kostelitz-Thouless transition.

Let us consider now the case when also the Coulomb potential is present. The effect of modulus fluctuations on the irreducible fermionic bubbles is the same as already discussed above , so Eqs. (74) and (76)-(78) are still valid. As far as density fluctuations are concerned, we just have to take into account that the RPA resummation of the irreducible bubbles in the particle-hole channel must be performed with the full potential $`V(𝐪)U/2`$. As a consequence, even though we are considering a strong-coupling (i.e. large $`U`$) limit, as $`𝐪`$ goes to zero the Coulomb interaction $`V(𝐪)`$ is always predominant. Thus in the hydrodynamic limit the contributions to the correction $`\delta D`$ only come from the static interaction vertices, as already discussed in the weak-coupling limit, since the $`\chi _{\rho E}`$ bubble vanishes for small $`𝐪`$, while the current-current correlation functions attain a finite value. More explicitly, by performing the RPA resummations (79)-(81) with $`\stackrel{~}{U}`$ replaced by the full dressed potential $`\stackrel{~}{U}\stackrel{~}{V}(𝐪)=[V(𝐪)+U/2]/[1+(V(𝐪)U/2)\chi ^{irr}]`$, it can be easily seen that

$`\chi _{\rho E}`$  $`{\displaystyle \frac{1}{V(𝐪)}},`$ (86)
$`\chi _{EE}`$  $`\chi _{EE}^{irr}{\displaystyle \frac{(\chi _{\rho E}^{irr})^2}{\chi ^{irr}}}=2t\alpha (1\delta ^2)+𝒪(\alpha ^3t),`$ (89)
$`\stackrel{~}{\chi }_{EE}`$  $`\stackrel{~}{\chi }_{EE}^{irr}{\displaystyle \frac{(\chi _{\rho E}^{irr})^2}{\chi ^{irr}}}=𝒪(\alpha ^3t).`$ (92)
As a consequence, in Eq. (54) for $`\delta D`$ a coefficient $`(3\chi _{EE}^{irr}+\stackrel{~}{𝒯})/8=a^2D/2`$ survives. This coefficient plays the same role as $`\xi _0^2D/2`$ in the $`XY`$-model, see Eq. (32) (which however should be evaluated with the long-range propagator $`P(q)`$). Thus, Eq. (83) still holds, both at zero temperature and in the classical regime, where by definition only static anharmonic terms survive. Observe that this result is again not trivial, because only properly including modulus fluctuations one can find the equivalence between $`\chi _{EE}^{irr},\stackrel{~}{𝒯}`$ and $`a^2D`$ in the strong-coupling regime for the pairing potential $`U`$.

Finally, we briefly comment on the extension of the previous results to the case of $`d`$-wave symmetry of the order parameter. As far as modulus fluctuations in the particle-hole channel are concerned, in the $`d`$-wave case this would reflect in a Hartree-Fock-like correction to the band dispersion. Since this makes the analytical treatment not viable, both for the coherence length and for $`\delta D`$ itself we do not address this issue in the present work. However, at generic filling the physics is not expected to be different from the $`s`$-wave case, and the strong-coupling mapping of the effective action in the $`XY`$ model should be preserved.

## VII Conclusions

In this paper we analyzed in detail the issue of the collective-mode description of a superconductor in the quantum regime. In particular, we addressed the issue from the point of view of the phase-fluctuation correction to the superfluid density, and we extended the previously known results for continuum microscopic models to the more realistic lattice case.

Usually, the evaluation of the phase-fluctuation contribution to the superfluid-density depletion is addressed within the quantum generalizations of the classical $`XY`$ model, which however fails in providing a complete description of dynamical effects. Specifically, in the continuum system the finite correction to $`D`$ down to $`T=0`$ derived within the $`XY`$ model implies the explicit violation of Galilean invariance. As we showed, the inadequacy of the quantum $`XY`$ model can be understood as one derives the phase-only action starting from the microscopic model. This approach allows us to treat on the same footing both static and dynamic interaction terms in the phase field: besides the classical (static) anharmonic terms $`(\stackrel{}{}\theta )^4`$, one has third- and fourth-order quantum (dynamic) interaction terms which contain the time derivative of $`\theta `$. These quantum terms are absent in the quantum $`XY`$ model, where the dynamics only appears at the Gaussian level, and induce a correction to $`D`$ which cancels exactly, in the continuum case, the contribution due to the classical interaction, restoring the equality $`\rho _s=\rho `$ at $`T=0`$. We point out that this evident failure of the quantum $`XY`$ in the continuum case naturally poses the question of the quantitative validity of the quantum $`XY`$ model estimates even in the case of a lattice microscopic model.

The detailed description of the main steps required to derive the phase-only action was motivated by the fact that this formalism also provides a more general framework for the description of different microscopic systems. As a consequence, we were able to investigate the phase mode both in the presence and in the absence of long-range Coulomb interactions, and for both $`s`$-wave and $`d`$-wave symmetry of the order parameter. As far as the lattice case is concerned, we also showed that the phase-only action acquires a more complex structure due to the fact that the minimal substitution couples the fermions to the electromagnetic field at all the orders. Moreover, since the stiffness at $`T=0`$ is no longer related to the particle density but rather proportional to the average kinetic energy (for nearest-neighbor hopping), no conservation law is violated by finding that classical and quantum interaction terms lead now to a finite one-loop correction $`\delta D`$. When compared with the result of the quantum $`XY`$ model, one also sees that $`\delta D`$ is of order $`1/(k_F\xi _0)^2`$ with respect to $`\delta D_{XY}`$. A similar result also holds in the classical regime, where we found that the phase-fluctuation correction to $`D`$ is smaller (in the weak-to-intermediate-coupling regime) than within the classical $`XY`$ model by the same factor $`1/(k_F\xi _0)^2`$, in both the continuum and the lattice model. The reduction of the phase-fluctuation effects for $`k_F\xi _01`$ is made even more pronounced in the continuum case by the inclusion of long-range Coulomb forces.

We also devoted particular emphasis to the discussion of the strong-coupling regime. As the interaction strength increases the fermionic model evolves toward a bosonic one, where the Cooper pairs act as individual bosons. Since in the superfluid the order parameter is connected to density, one expects that the description of the strongly-interacting electronic system requires the inclusion of both density and modulus fluctuations. We explicitly considered the negative-$`U`$ Hubbard ($`s`$-wave) model. This model and its extensions are widely considered in the literature, as candidates to capture the main features of the intermediate-to-strong-coupling superconductors. This issue has attracted renewed interest due to its possible relevance to the description of the underdoped regime of the superconducting cuprates.

The inclusion of density fluctuations within the negative-$`U`$ Hubbard model is performed by RPA resummation in the particle-hole channel of the fermionic bubbles which appear as coefficients of the phase-only action. As far as the modulus fluctuations are concerned, their inclusion in the fermionic bubbles is directly derived by the definition of the susceptibilities in terms of the physical quantities, although the formal treatment of modulus fluctuation at Gaussian level is possible. The discussion of the result strongly involves the knowledge of the behavior of the coherence length $`\xi _0`$ at strong coupling as a function of the fermion density, since it sets the spatial cut-off for phase fluctuations. This issue was addressed elsewhere,noi4  and allows us to conclude that at strong coupling, both in the presence and in the absence of long-range Coulomb forces, the fermionic system maps onto the quantum $`XY`$ model, except in the very low-density quantum regime.

Even though our attention was devoted to the correction to the superfluid density, our results for the structure of the phase-only action also suggest the possible outcomes of the analysis of phase-fluctuation effects on other quantities. The connection between the coefficients of the phase-only action and the various physical response functions (in specific regimes) allows us to relate the correction to the Gaussian phase-only action to the phase-fluctuation induced correction to the corresponding physical quantity. This has been done, e.g., for the optical conductivity, as discussed in Ref. \[noicond, \] in connection to the physics of high-$`T_c`$ superconductors. A second possibility concerns the phase-fluctuation effect on the thermal conductivity, which is presently under investigation thermal .

###### Acknowledgements.

This work comes out from a fruitful collaboration with C. Castellani, to whom we are particularly indebted. We also thank S. De Palo, C. Di Castro, M. Grilli, S. Sharapov for many useful discussions and suggestions. L.B. acknowledges partial financial support by the Swiss National Science Fundation under Grant No. 620-62868.00. SC acknowledges financial support from the Italian MIUR, Cofin 2001, prot. 20010203848, and from INFM, PA-G0-4.

## Appendix A Deduction of the phase-only action by means of Fermionic variable integrations

In this appendix we report some details on the derivation of the phase-only action. As a first point we explicitly perform here the Hubbard-Stratonovich transformation on the microscopic action of Eq. (5), introducing the auxiliary field $`\mathrm{\Delta }(𝐪,\tau )`$. In this way we get the decoupling of the interacting part $`S_I`$ of the action (5) in the particle-particle channel, which now reads

$`S_I(c_\sigma ,c_\sigma ^+,\mathrm{\Delta },\mathrm{\Delta }^{})`$ $`=`$ $`{\displaystyle \frac{1}{U}}{\displaystyle \underset{𝐪,\mathrm{\Omega }_m}{}}|\mathrm{\Delta }(𝐪,\mathrm{\Omega }_m)|^2\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}{\displaystyle \underset{𝐤,𝐪,\omega _n,\mathrm{\Omega }_m}{}}\mathrm{\Delta }(𝐪,\mathrm{\Omega }_m)w(𝐤)c_{𝐤+\frac{𝐪}{2}}^+(\omega _n)c_{𝐤+\frac{𝐪}{2}}^+(\omega _n\mathrm{\Omega }_m)`$ (93)
$``$ $`\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}{\displaystyle \underset{𝐤,𝐪,\omega _n,\mathrm{\Omega }_m}{}}[\mathrm{\Delta }(𝐪,\mathrm{\Omega }_m)]^{}w(𝐤)c_{𝐤+\frac{𝐪}{\mathrm{𝟐}}}(\omega _n)c_{𝐤+\frac{𝐪}{2}}(\omega _n\mathrm{\Omega }_m).`$
Here $`\mathrm{\Delta }(𝐪,\mathrm{\Omega }_m)`$ is the Fourier transform of the auxiliary field, $`𝐪`$ is the momentum and $`\mathrm{\Omega }_m=2\pi m/\beta `$ is the Matsubara frequency.

A second issue to be discussed is the effect of the gauge transformation

$`c_\sigma (𝐫,\tau )c_\sigma ^{}(𝐫,\tau )`$ $`=`$ $`c_\sigma (𝐫,\tau )e^{i\theta (𝐫,\tau )/2},`$
$`c_\sigma ^+(𝐫,\tau )c_{}^{+}{}_{\sigma }{}^{}(𝐫,\tau )`$ $`=`$ $`c_\sigma ^+(𝐫,\tau )e^{i\theta (𝐫,\tau )/2},`$ (94)
on the fermionic field appearing $`S_0`$ and $`S_I`$. We start by rewriting Eq. (93) in real space

$`S_I(c_\sigma ,c_\sigma ^+,\mathrm{\Delta },\mathrm{\Delta }^{})={\displaystyle \frac{|\mathrm{\Delta }(𝐑,\tau )|^2}{U}𝑑𝐑𝑑\tau }`$ $``$ $`{\displaystyle \mathrm{\Delta }(𝐑,\tau )w(𝐫)c_{}^+(𝐑+\frac{𝐫}{2},\tau )c_{}^+(𝐑\frac{𝐫}{2},\tau )𝑑𝐑𝑑𝐫𝑑\tau }`$ (95)
$``$ $`{\displaystyle \mathrm{\Delta }^{}(𝐑,\tau )w(𝐫)c_{}(𝐑\frac{𝐫}{2},\tau )c_{}(𝐑+\frac{𝐫}{2},\tau )𝑑𝐑𝑑𝐫𝑑\tau },`$
where $`w(𝐫)=_𝐤w(𝐤)e^{i\mathrm{𝐤𝐫}}`$ is the Fourier transform of $`w(𝐤)`$. After the transformation (94), we obtain, e.g., for the first term linear in $`\mathrm{\Delta }`$ in Eq. (95)

$`S_I`$ $`=`$ $`\mathrm{}{\displaystyle _\mathrm{\Omega }}{\displaystyle _0^\beta }|\mathrm{\Delta }(𝐑,\tau )|e^{i\theta (𝐑,\tau )}w(𝐫)c_{}^+(𝐑+{\displaystyle \frac{𝐫}{2}},\tau )c_{}^+(𝐑{\displaystyle \frac{𝐫}{2}},\tau )𝑑𝐑𝑑𝐫𝑑\tau \mathrm{}`$
$``$ $`{\displaystyle _0^\beta }`$ $`\left[|\mathrm{\Delta }(𝐑,\tau )|e^{i[\theta (𝐑,\tau )\frac{\theta (𝐑\frac{𝐫}{2},\tau )}{2}\frac{\theta (𝐑+\frac{𝐫}{2},\tau )}{2}]}w(𝐫)c_{}^{+}{}_{}{}^{}(𝐑+{\displaystyle \frac{𝐫}{2}},\tau )c_{}^{+}{}_{}{}^{}(𝐑{\displaystyle \frac{𝐫}{2}},\tau )\right]d𝐑d𝐫d\tau .`$ (96)
If the interaction is local in real space, i.e., $`w(𝐫)=\delta (𝐫)`$, as in the (isotropic) $`s`$-wave case, the phase $`\theta `$ disappears from the exponential in Eq. (96). In the $`d`$-wave case instead, one can expand $`e^{i[\theta (𝐑,\tau )\frac{\theta (𝐑\frac{𝐫}{2},\tau )}{2}\frac{\theta (𝐑+\frac{𝐫}{2},\tau )}{2}]}1i𝐫_\alpha 𝐫_\beta \stackrel{}{}_\alpha \stackrel{}{}_\beta \theta (𝐑)+\mathrm{}`$. The residual dependence on the phase affects the form of the phase-fluctuation propagator by introducing terms which are irrelevant in the hydrodynamical limit. As a consequence, in the following we completely discard the dependence on the phase field in $`S_I`$, assuming that after the gauge transformation (94) $`S_I`$ has the form (95) with both $`\mathrm{\Delta }(x)`$ and $`\mathrm{\Delta }^{}(x)`$ substituted by $`|\mathrm{\Delta }(x)|`$.

As far as $`S_0`$ in Eq. (5) is concerned, the effect of the gauge transformation can be better expressed by first rewriting

$$S_0=𝑑x\varphi ^+(x)\left(\begin{array}{cc}_\tau \frac{\stackrel{}{}^2}{2m}\mu & 0\\ 0& _\tau +\frac{\stackrel{}{}^2}{2m}+\mu \end{array}\right)\varphi (x),$$

where we introduced the standard Nambu spinor $`\varphi ^+(𝐫,\tau )=\left(c_{}^+(𝐫,\tau )c_{}(𝐫,\tau )\right)`$. After the transformation (94), $`S_0`$ gets modified into

$$\stackrel{~}{S_0}=S_0+dx\varphi ^+(x)\{[\frac{i_\tau \theta (x)}{2}+\frac{(\stackrel{}{}\theta (x))^2}{8m}]\widehat{\tau }_3+[\frac{i}{4m}\stackrel{}{}\theta (x)\stackrel{}{}]\widehat{\tau }_0\}\varphi (x),$$
(97)
where the operator $`\stackrel{}{}(\stackrel{}{}\stackrel{}{})`$ acts on the fermionic variables.

Introducing the Nambu notation also for the interacting part of the action, after the gauge transformation we find

$$S=\underset{k^{}k}{}\underset{ij}{}\varphi _i^+(k^{})A_{k^{}k}^{ij}\varphi _j(k)+\underset{q}{}\frac{|\mathrm{\Delta }(q)|^2}{U},$$
(98)
where we put $`k=(𝐤,\omega _n)`$ and $`q=(𝐪,\mathrm{\Omega }_m)`$, and the matrix $`A_{k^{}k}^{ij}`$ can be deduced from Eqs. (93) and (97). As we observed in Sec. II.1, the fermionic variables appear in the action (98) only in the Gaussian term, whose coefficient $`A_{k^{}k}^{ij}`$ contains the collective variables $`|\mathrm{\Delta }|`$ and $`\theta `$. By integrating out the fermions,negele

$$𝒟\varphi 𝒟\varphi ^+e^{_{k^{^{}}k}_{ij}\varphi _i^+(k^{^{}})A_{k^{^{}}k}^{ij}\varphi _j(k)}=\text{Det}A_{k^{^{}}k}^{ij},$$

and using the well-known identity: $`\mathrm{ln}\text{Det}A_{k^{^{}}k}^{ij}=\text{Tr}\mathrm{ln}A_{k^{^{}}k}^{ij}`$, where the trace Tr acts in the Nambu space and on the four-momenta $`k,k^{}`$, we finally get the effective action for the bosonic degrees of freedom

$$S_{eff}(|\mathrm{\Delta }|,\theta )=\underset{q}{}\frac{|\mathrm{\Delta }(q)|^2}{U}\text{Tr}\mathrm{ln}A_{k^{^{}}k}^{ij}.$$
(99)
The second term of Eq. (99) acquires a more readable form by separating the part of the matrix $`A_{k,k^{}}`$ with the explicit structure of $`\delta _{kk^{}}`$ from the rest. This leads to

$$\text{Tr}\mathrm{ln}A_{k^{^{}}k}^{ij}=\text{Tr}\mathrm{ln}\left[\widehat{𝒢}_0^1\widehat{\mathrm{\Sigma }}\right]=\text{Tr}\mathrm{ln}\widehat{𝒢}_0^1+\text{Tr}\mathrm{ln}\left[\widehat{1}\widehat{𝒢}_0\widehat{\mathrm{\Sigma }}\right]=\text{Tr}\mathrm{ln}\widehat{𝒢_0^1}\underset{N}{}S_{eff}^N,$$

where

$$S_{eff}^N=\frac{1}{N}\text{Tr}\stackrel{N\text{ times}}{\stackrel{}{\left[\widehat{𝒢_0}\widehat{\mathrm{\Sigma }}\times \mathrm{}\times \widehat{𝒢_0}\widehat{\mathrm{\Sigma }}\right]}},$$

and

$`\widehat{\mathrm{\Sigma }}`$ $`=`$ $`\left[\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}\left({\displaystyle \frac{\omega _n^{^{}}\omega _n}{2}}\right)\theta (k^{^{}}k){\displaystyle \frac{1}{8m}}{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{s}{}}(𝐤^{}𝐤𝐬)𝐬\theta (k^{^{}}ks)\theta (s)\right]\widehat{\tau _3}`$ (100)
$`+`$ $`\left[{\displaystyle \frac{i}{4m}}\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}(|𝐤^{}||𝐤|)(|𝐤^{}|+|𝐤|)\theta (k^{^{}}k)\right]\widehat{\tau }_0\left[\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}w\left({\displaystyle \frac{𝐤^{}𝐤}{2}}\right)|\mathrm{\Delta }|(k^{^{}}k)\right]\widehat{\tau _1}.`$
The BCS matrix $`𝒢_0`$ is

$$\widehat{𝒢}_0=\left(\begin{array}{cc}𝒢_0(𝐤,\omega _n)& _0(𝐤,\omega _n)\\ _0(𝐤,\omega _n)& 𝒢_0(𝐤,\omega _n)\end{array}\right),$$
(101)
where

$$𝒢_0(𝐤,\omega _n)=\frac{i\omega _n\xi _𝐤}{\omega _n^2+E_𝐤^2},_0(𝐤,\omega _n)=\frac{w(𝐤)\sqrt{\frac{T}{\mathrm{\Omega }}}|\mathrm{\Delta }|(0)}{\omega _n^2+E_𝐤^2},$$
(102)
and $`E_𝐤=\sqrt{\xi _𝐤^2+w^2(𝐤)\mathrm{\Delta }_0^2}`$. The value of the gap amplitude $`\mathrm{\Delta }_0=\sqrt{T/\mathrm{\Omega }}|\mathrm{\Delta }|(0)`$ in Eqs. (102) is determined by the solution of the saddle-point equation for the effective action (99), $`(S/|\mathrm{\Delta }|(q=0))=0`$, and coincides with the mean-field value of the superconducting order parameter. The saddle point equation imposes a constraint only on the modulus of the order parameter, since a constant value of the Goldstone field $`\theta `$ does not affect the action. As a consequence, $`𝒢_0`$ and $`_0`$ are exactly the BCS normal and anomalous Green functions for the Hamiltonian (3)-(4).schrieffer  By separating the field $`|\mathrm{\Delta }|(q)=|\mathrm{\Delta }|(0)+\delta |\mathrm{\Delta }|(q)`$ in its constant and fluctuating parts, the effective action (99) can be expanded in terms of phase and modulus fluctuations with respect to the BCS state. Observe that since one separates modulus and phase $`|\mathrm{\Delta }|(q)`$ is the Fourier transform of the real field $`|\mathrm{\Delta }|(x)`$, and consequently $`[|\mathrm{\Delta }|(q)]^{}=|\mathrm{\Delta }|(q)`$ and $`_q|\mathrm{\Delta }(q)|^2=_q||\mathrm{\Delta }|(q)|^2=_q|\mathrm{\Delta }|(q)|\mathrm{\Delta }|(q)`$, where $`\mathrm{\Delta }(q)`$ is the Fourier transform of the complex field $`\mathrm{\Delta }(x)`$.

When, at the end, the gaussian expansion in the phase fluctuations around the saddle-point solution is performed, the Gaussian effective action reported in Eq. (II.1) is obtained,

$$S_{neutral}^G=\frac{1}{8}\underset{q}{}\left[\mathrm{\Omega }_m^2\mathrm{\Lambda }_{\rho \rho }(q)𝐪_a𝐪_b\mathrm{\Lambda }_{JJ}^{ab}(q)+2i\mathrm{\Omega }_m𝐪_a\mathrm{\Lambda }_{\rho J}^a(q)\right]|\theta (q)|^2.$$

where the explicit expressions of coefficient $`\mathrm{\Lambda }`$ are the following

$`\mathrm{\Lambda }_{\rho \rho }(q)`$ $`=`$ $`{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{k}{}}\mathrm{tr}[\widehat{𝒢}_0(k)\widehat{\tau }_3\widehat{𝒢}_0(k+q)\widehat{\tau }_3]`$ (103)
$`\mathrm{\Lambda }_{\rho J}^a(q)`$ $`=`$ $`{\displaystyle \frac{T}{\mathrm{\Omega }}}{\displaystyle \underset{k}{}}{\displaystyle \frac{𝐤^a}{m}}\mathrm{tr}\left[\widehat{𝒢}_0(𝐤+{\displaystyle \frac{𝐪}{2}},\omega _n)\widehat{\tau }_0\widehat{𝒢}_0(𝐤{\displaystyle \frac{𝐪}{2}},\omega _n\mathrm{\Omega }_m)\widehat{\tau }_3\right]`$ (104)
$`\mathrm{\Lambda }_{JJ}^{ab}`$ $`=`$ $`{\displaystyle \frac{\rho _{MF}}{m}}\delta _{ab}+\mathrm{\Lambda }_{jj}^{ab}(q),`$ (105)
with

$$\mathrm{\Lambda }_{jj}^{ab}(q)=\frac{T}{\mathrm{\Omega }}\underset{k}{}\frac{𝐤^a}{m}\frac{𝐤^b}{m}\mathrm{tr}\left[\widehat{𝒢}_0(𝐤+\frac{𝐪}{2},\omega _n)\widehat{\tau }_0\widehat{𝒢}_0(𝐤\frac{𝐪}{2},\omega _n\mathrm{\Omega }_m)\widehat{\tau }_0\right]$$
(106)
and

$$\rho _{MF}=\frac{T}{\mathrm{\Omega }}\underset{k}{}\mathrm{tr}\left[\widehat{𝒢}_0(k)\widehat{\tau }_3\right]=1\frac{1}{\mathrm{\Omega }}\underset{𝐤}{}\frac{\xi _𝐤}{E_𝐤}\mathrm{tanh}\left(\frac{\beta E_𝐤}{2}\right),$$
(107)
is the saddle-point fermion density.

In the presence of long-range Coulomb interaction, the Hubbard-Stratonovich transformation of the microscopic action (5) reads

$`S_I`$ $`(c_\sigma ,c_\sigma ^+,\mathrm{\Delta },\mathrm{\Delta }^{},\rho _{HS})={\displaystyle \underset{q}{}}\{{\displaystyle \frac{|\mathrm{\Delta }(q)|^2}{U}}+{\displaystyle \frac{|\rho _{HS}(q)|^2}{2V(𝐪)}}\mathrm{\Delta }(q)\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}{\displaystyle \underset{𝐤,\omega _n}{}}w(𝐤)c_{𝐤+\frac{𝐪}{2}}^+(\omega _n)c_{𝐤+\frac{𝐪}{2}}^+(\omega _n\mathrm{\Omega }_m)`$
$``$ $`\mathrm{\Delta }^{}(q)\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}{\displaystyle \underset{𝐤,\omega _n}{}}w(𝐤)c_{𝐤+\frac{𝐪}{\mathrm{𝟐}}}(\omega _n)c_{𝐤+\frac{𝐪}{2}}(\omega _n\mathrm{\Omega }_m)+i\rho _{HS}(q)\sqrt{{\displaystyle \frac{T}{\mathrm{\Omega }}}}{\displaystyle \underset{𝐤,\omega _n,\sigma }{}}c_{𝐤𝐪,\sigma }^+(\omega _n\mathrm{\Omega }_m)c_{𝐤,\sigma }(\omega _n)\},`$
where the first, the third and the fourth terms were already present in Eq. (93). The last term describes the interaction of the fermions with the $`\rho _{HS}`$. This term is not affected by the gauge transformation (94), since it is of the form $`c^+c`$. Integrating out the fermions the effective action for the collective variables(99) is now

$$S_{eff}(|\mathrm{\Delta }|,\theta ,\rho _{HS})=\underset{q}{}\left\{\frac{|\mathrm{\Delta }(q)|^2}{U}+\frac{|\rho _{HS}(q)|^2}{2V(𝐪)}\right\}\text{Tr}\mathrm{ln}A_{k^{^{}}k}^{ij},$$
(108)
where the matrix elements $`A_{k^{}k}^{ij}`$ include now the density field $`\rho _{HS}`$, $`A_{k^{^{}}k}A_{k^{^{}}k}+i\sqrt{T/\mathrm{\Omega }}\rho _{HS}(kk^{^{}})\widehat{\tau }_3`$. As a consequence, the density contributes to both the diagonal part of $`A_{k^{}k}^{ij}`$, whence

$$\widehat{𝒢}_0^1\widehat{𝒢}_0^1i\sqrt{\frac{T}{\mathrm{\Omega }}}\rho _{HS}(0)\widehat{\tau }_3,$$
(109)
and to the self-energy matrix $`\widehat{\mathrm{\Sigma }}\widehat{\mathrm{\Sigma }}+i\rho _{HS}(x)\widehat{\tau }_3.`$ Again, we separate $`\rho _{HS}(q)=\rho _{HS}(0)+\delta \rho _{HS}(q)`$. The value of $`\rho _{HS}(0)`$ is determined by the solution of the saddle-point equation for the action (108), i.e.

$$\frac{i\rho _{HS}(0)}{V(0)}=\frac{T}{\mathrm{\Omega }}\underset{k}{}\left[𝒢_0(k)+𝒢_0(k)\right]=\rho _{MF}$$

where $`\rho _{MF}`$ is the particle density evaluated at mean-field level according to Eq. (107). Therefore, the correction (109) to the Green function is the standard Hartree shift of the chemical potential $`\mu \mu +V(0)\rho _{MF}`$.notamu .

In the charged case, as anticipated in Sec. II B, in order to get the phase-only effective action at the Gaussian level, it is necessary at first to evaluate the Gaussian expansion of $`S_{eff}`$ both in the phase and the density fluctuations

$$S_{eff}^G(\theta ,\delta \rho _{HS})=\frac{1}{8}\underset{q}{}\left(\begin{array}{cc}\theta (q)& \delta \rho _{HS}(q)\end{array}\right)\widehat{B}(q)\left(\begin{array}{c}\theta (q)\\ \delta \rho _{HS}(q)\end{array}\right),$$

where

$$\widehat{B}(q)=\left(\begin{array}{cc}\mathrm{\Omega }_m^2\mathrm{\Lambda }_{\rho \rho }(q)𝐪_a𝐪_b\mathrm{\Lambda }_{JJ}^{ab}(q)+2i\mathrm{\Omega }_m𝐪_a\mathrm{\Lambda }_{\rho J}^a& 2𝐪_a\mathrm{\Lambda }_{\rho J}^a(q)+2i\mathrm{\Omega }_m\mathrm{\Lambda }_{\rho \rho }(q)\\ 2𝐪_a\mathrm{\Lambda }_{\rho J}^a(q)2i\mathrm{\Omega }_m\mathrm{\Lambda }_{\rho \rho }(q)& 4\left(\frac{1}{V(𝐪)}+\mathrm{\Lambda }_{\rho \rho }(q)\right)\end{array}\right),$$

with the same definitions (103)-(105) of the coefficients $`\mathrm{\Lambda }`$, and, in a second time, to perform the Gaussian integral over the field $`\delta \rho _{HS}`$.

As a final result of this procedure, the following expression (15) for the effective action is obtained

$$S_{charged}^G(\theta )=\frac{1}{8}\underset{q}{}\left[\mathrm{\Omega }_m^2_{\rho \rho }(q)𝐪_a𝐪_b_{JJ}^{ab}(q)+2i𝐪_a\mathrm{\Omega }_m_{\rho J}^a(q)\right]\theta (q)\theta (q),$$

where the coefficients $``$ are given by

$`_{\rho \rho }(q)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Lambda }_{\rho \rho }(q)}{1+V(𝐪)\mathrm{\Lambda }_{\rho \rho }(q)}},`$
$`_{\rho J}^a(q)`$ $`=`$ $`{\displaystyle \frac{\mathrm{\Lambda }_{\rho J}^a(q)}{1+V(𝐪)\mathrm{\Lambda }_{\rho \rho }(q)}},`$
$`_{JJ}^{ab}(q)`$ $`=`$ $`\mathrm{\Lambda }_{JJ}^{ab}(q)V(𝐪){\displaystyle \frac{\mathrm{\Lambda }_{\rho J}^a(q)\mathrm{\Lambda }_{\rho J}^b(q)}{1+V(𝐪)\mathrm{\Lambda }_{\rho \rho }(q)}}.`$
## Appendix B Effective action and gauge invariance

In Sec. II.1 we observed that the BCS bubbles $`\mathrm{\Lambda }(q)`$ which appear as coefficients of the Gaussian effective action (II.1) break, as it is well known schrieffer , the gauge invariance of the theory. In this appendix we derive the general relationship between the coefficents $`\mathrm{\Lambda }(q)`$ and the corresponding gauge-invariant electromagnetic response functions $`K(q)`$, which control the physically accessible quantities (see also Ref. \[arun00, ; sharapov, \]).

We first introduce the electromagnetic potential $`A_\mu =(\varphi ,𝐀)`$ into the Gaussian effective action (II.1) via the minimal substitution $`(\theta /t)(\theta /t)+2e\phi `$, $`\stackrel{}{}\theta \stackrel{}{}\theta (2e/c)𝐀`$, which gives

$$S_G(\theta )S_G(\theta ,A)=S_G(\theta )+\frac{e^2}{2}\underset{q}{}A_\mu (q)\mathrm{\Lambda }^{\mu \nu }(q)A_\nu (q)+\underset{q}{}ie^2[q_\mu \mathrm{\Lambda }^{\mu \nu }(q)A_\nu (q)\theta (q)q_\mu \mathrm{\Lambda }^{\mu \nu }(q)A_\nu (q)\theta (q)],$$

where the four-dimensional notation for $`\mathrm{\Lambda }^{\mu \nu }`$ was introduced, so that $`\mathrm{\Lambda }_{\rho \rho }=\mathrm{\Lambda }^{00}`$, $`\mathrm{\Lambda }_{\rho J}^a=\mathrm{\Lambda }^{0a}`$ and $`\mathrm{\Lambda }_{JJ}^{ab}=\mathrm{\Lambda }^{ab}`$. We can now integrate out the $`\theta `$ field to obtain the partition function $`Z[A_\mu ]`$ which allows one to define the electromagnetic response functions of the system as

$$K^{\mu \nu }(q)\frac{\delta ^2\mathrm{ln}Z[A]}{\delta A_\mu (q)\delta A_\nu (q)}|_{A_\mu (q)=0,A_\nu (q)=0}=\mathrm{\Lambda }^{\mu \nu }(q)+\frac{\mathrm{\Lambda }^{\mu \alpha }(q)q_\alpha q_\beta \mathrm{\Lambda }^{\beta \nu }(q)}{q_\gamma \mathrm{\Lambda }^{\gamma \delta }(q)q_\delta }.$$
(110)
This equation provides the relation between the coefficients $`\mathrm{\Lambda }^{\mu \nu }`$ of Eq. (II.1) and the physical response functions $`K^{\mu \nu }`$. Now it can be easily seen that the difference between the $`\mathrm{\Lambda }`$ and the $`K`$ functions (i.e., the second term in the right-hand side of Eq. (110)) has a purely longitudinal structure. In other words phase fluctuations affect only longitudinal correlation functions schrieffer . As a consequence, we can safely identify the transverse part of the Gaussian bubbles ($`\mathrm{\Lambda }^T`$) with the transverse part of the physical response functions ($`K^T`$).

This fact turns out to be crucial in interpreting the Gaussian action (II.1) in the hydrodynamic limit, which is the physically interesting case while evaluating the low-temperature behavior of the collective mode. Indeed, by calculating the static limit $`\mathrm{\Omega }_m=0,𝐪0`$ of the coefficients (103)-(105) we observe that:

* In the Gaussian action (II.1) the longitudinal part of $`\mathrm{\Lambda }_{JJ}^{ab}`$ appears, which, as we just showed, is not the physical longitudinal correlation function. However, in the $`q0`$ limit, the transverse $`\mathrm{\Lambda }_{JJ}^T`$ and longitudinal $`\mathrm{\Lambda }_{JJ}^L`$ part of $`\mathrm{\Lambda }_{JJ}^{ab}`$ do coincide, leading to the identification of the coefficient of $`𝐪^2`$ with the stiffness, defined through the transverse physical function (see, e.g., Ref. \[arun00, \]),

$$\underset{𝐪0}{lim}\mathrm{\Lambda }_{JJ}^L(𝐪,\mathrm{\Omega }_m=0)=\underset{𝐪0}{lim}\mathrm{\Lambda }_{JJ}^T(𝐪,\mathrm{\Omega }_m=0)=\underset{𝐪0}{lim}K_{JJ}^T(𝐪,\mathrm{\Omega }_m=0)=D(T).$$
* The coefficient $`\mathrm{\Lambda }_{\rho J}^a(q)`$ vanishes in the static limit, as it can be seen by considering its definition (104).
* From Eq. (110) it turns out that the coefficient $`\mathrm{\Lambda }_{\rho \rho }`$ coincides in the static limit with the physical density-density correlation function $`K_{\rho \rho }`$, whose static limit is by definition the (bare) compressibility $`\kappa _0(T)`$ of the system, $`lim_{𝐪0}\mathrm{\Lambda }_{\rho \rho }(𝐪,\mathrm{\Omega }_m=0)=lim_{𝐪0}K_{\rho \rho }(𝐪,\mathrm{\Omega }_m=0)=\kappa _0(T)`$.

As a consequence of the above results, the Gaussian phase-only action in the hydrodynamic regime as exactly the expression presented in Eq. (II.1).