Goldstone and Higgs Hydrodynamics in the BCS-BEC Crossover

We discuss the derivation of a low-energy effective field theory of phase (Goldstone) and amplitude (Higgs) modes of the pairing field from a microscopic theory of attractive fermions. The coupled equations for Goldstone and Higgs fields are critically analyzed in the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensate (BEC) crossover both in three spatial dimensions and in two spatial dimensions. The crucial role of pair fluctuations is investigated, and the beyond-mean-field Gaussian theory of the BCS-BEC crossover is compared with available experimental data of the two-dimensional ultracold Fermi superfluid.


Introduction
An important achievement in the physics of ultracold atoms was the experimental realization of the crossover from the Bardeen-Cooper-Schrieffer (BCS) superfluid phase of loosely-bound pairs of fermions to the Bose-Einstein condensate (BEC) of tightly-bound composite bosons [1,2]. Superfluid fermions in the BCS-BEC crossover are an ideal platform to investigate macroscopic quantum phenomena. For instance, macroscopic quantum tunneling and quantum self-trapping [3][4][5] have recently been observed across the BCS-BEC crossover [6].
Very recently, the BCS-BEC crossover has also been realized in quasi two-dimensional (2D) configurations [7][8][9][10]. Beyond-mean-field theoretical investigations of 2D Fermi gases in the full BCS-BEC crossover have been carried out both at zero and finite temperature [11][12][13][14][15][16][17]. These 2D results have clearly shown that contrary to the 3D case, mean-field theories are completely unreliable for the study of strongly-interacting superfluid fermions in two dimensions because of the huge increase of quantum fluctuations. The 2D BCS-BEC crossover is also interesting for high-T c superconductivity where the phase diagram of cuprate superconductors can be interpreted in terms of a BCS-BEC crossover as doping is varied. The critical temperature T c has a wide fluctuation region with pseudo-gap effects not yet fully understood [18,19]. Moreover, it has been suggested that iron-based superconductors have composite superconductivity, consisting of strong-coupling BEC in the electron band and weak-coupling BCS-like superconductivity in the hole band [20].
This year, we have used renormalization-group techniques to calculate the Berezinskii-Kosterlitz-Thouless (BKT) critical temperature [21,22], taking explicitly into account the formation of quantized vortices [23]. The presence of quantized vortices and antivortices renormalizes the superfluid density of the system as the temperature increases, and the renormalized superfluid density jumps from a finite value to zero as the temperature reaches the BKT critical temperature [24]. In this way, we have also obtained the BKT critical temperature in the full 2D BCS-BEC crossover for the uniform Fermi superfluid [23]. Quantized vortices in superfluids are a peculiar consequence of the existence of an underlying compact real field, whose spatial gradient determines the local superfuid velocity of the system [25][26][27]. This compact real field-the so-called Goldstone field-is the phase angle of the complex bosonic field which in the case of attractive fermions describes strongly-correlated Cooper pairs of fermions with opposite spins. The amplitude of this complex pairing field-sometimes also called a Higgs field [28]-is decoupled from the angular Goldstone field only in the deep BCS regime where the particle-hole symmetry is almost preserved.
In this paper, we review the low-energy effective field theory of Goldstone and Higgs modes derived from the microscopic theory of paired fermions. At zero temperature, we show that the velocity of the Goldstone mode obtained with and without amplitude fluctuations (in 2D but also in 3D) has a quite different behavior [11]. We find that amplitude fluctuations are necessary to identify the velocity of the Goldstone mode with the first sound velocity one gets from the mean-field equation of state by using familiar thermodynamics relationships [11]. Finally, we show that in the 2D case, the first sound velocity obtained from the beyond-mean-field equation of state is quite different with respect to the mean-field one-in particular in the BEC side of the BCS-BEC crossover [15]. For the sake of completeness, we also report our beyond-mean-field calculation [13][14][15][16] of the zero-temperature pressure of the Fermi superfluid in the 2D BCS-BEC crossover: these theoretical results are in very good agreement with the available experimental data [7].

Functional Integration for the BCS-BEC Crossover
We consider a D-dimensional system of two-spin-component fermions interacting through an attractive s-wave contact potential, contained in a volume V, at fixed chemical potential µ and temperature T. Within the path integral formalism, the partition function of the system can be written as [25][26][27] where ψ σ (r, τ),ψ σ (r, τ) are complex Grassmann fields (σ =↑, ↓), β = 1/k B T, with k B is Boltzmann constant, and the Euclidean Lagrangian density reads Here g is the attractive interaction strength (g < 0) of the s-wave coupling between fermions with opposite spin σ. Looking for analytical solutions, the quartic interaction cannot be treated exactly. The Hubbard-Stratonovich transformation [25][26][27] introduces an additional auxiliary pairing field ∆(r, τ) corresponding to a Cooper pair, decoupling the quartic interaction. The transformation is based on the following identity The partition function can then be written as with the following Lagrangian density The integration over the fermionic fields ψ σ (r, τ) andψ σ (r, τ) can now be carried out exactly, obtaining with G −1 the inverse Green's function, given by To investigate effects of quantum and thermal fluctuations of the gap field ∆(r, t) around its mean-field value ∆ 0 , we set where η(r, τ) is the complex field of pairing fluctuations. In this way, the inverse Green function G −1 is decomposed in a mean-field component −G −1 0 , where the pairing field ∆ is replaced by its uniform and constant saddle point value, plus a fluctuation part F:

Loop Expansion and Gaussian Approximation
The logarithm appearing in Equation (6) can be written as [29]: The Gaussian (one-loop) approximation consists of the following expansion for the second term in the right-hand-side of Equation (10): Within this Gaussian approximation, the partition function reads where is the mean-field partition function and is the Gaussian partition function characterized by the following Gaussian action: In this formula, we have introduced the Fourier transform of the fluctuation fields and the bosonic Matsubara frequencies Ω m = 2mπ/β. The matrix elements of the inverse pair fluctuation propagator M are given by [29] , (17) where is the spectrum of fermionic single-particle excitations.

Beyond-Mean-Field Grand Potential
The thermodynamic grand potential Ω of the fermionic superfluid is given by At the Gaussian one-loop level, one gets where the mean-field grand potential reads [25][26][27]] The Gaussian grand potential is instead The sum over Matsubara frequencies is quite complicated, and it does not give a simple expression. An approximate formula which is valid in the BEC regime of the crossover [30] is the following: where is the spectrum of bosonic collective excitations with ω(q) derived from In the Gaussian pair fluctuation (GPF) approach [31], given the grand potential the energy gap ∆ 0 is obtained from the mean-field gap equation The number density n is instead obtained from the beyond-mean-field number equation taking into account the gap equation (i.e., that ∆ 0 depends on µ and T: ∆ 0 (µ, T)). Notice that the Nozieres and Schmitt-Rink approach es [32] are quite similar, but in the number equation one forgets that ∆ 0 depends on µ.

Low-Energy Gaussian Action
We have seen that the analytical form of the inverse pair fluctuation propagator M(q, ω) (with ω = iΩ m ) is quite complicated, and one can find its matrix elements numerically. Here we use a series expansion of M(q, ω) up to the second order in q and ω [33][34][35]. Moreover, we decompose the fluctuation field as follows: where σ(r, τ) and θ(r, τ) are real and can be identified at the lowest order with amplitude and phase fluctuations, respectively [33][34][35]. In other words, σ(r, τ) is the Higgs field and θ(r, τ) is the Goldstone field [28].
In this way, after some calculations [11,34,35], one finds the low-energy real-time Gaussian action derived from the Euclidean Gaussian action (15) where t = −iτ is the real time. At zero temperature, the coefficients J, K θθ , K σσ , and k θσ are related to the partial derivatives of the zero-temperature mean-field grand potential Ω m f (µ, V, T = 0, ∆ 0 ) [11,35]. In particular, is the phase stiffness, is the phase-phase susceptibility, is the amplitude-amplitude susceptibility, and is the phase-amplitude susceptibility, that is the Goldstone-Higgs coupling constant. Equation (30) is practically the same action functional derived in [36,37] by using a two-channel model model for the BCS-BEC crossover. The Euler-Lagrange equations for the Goldstone field θ(r, t) and Higgs field σ(r, t) obtained from (30) are Calculating the time derivative of σ(r, t) from Equation (36), one finds Inserting this result in Equation (35), we finally get the d'Alambert equation of waves which admits the generic solution with the dispersion relation ω k = c s k and the velocity of propagation of the Goldstone mode, where is the effective susceptibility. Notice that only if the amplitude fluctuations are negligible (i.e., σ ≃ 0); from Equation (36), it follows K θσ ≃ 0 and consequently K ≃ K θθ .
In the 3D BCS-BEC crossover, the interaction strength g is usually written in terms of the to 3D s-wave scattering length a [27] 1 Instead, in the 2D BCS-BEC crossover, the interaction strength g is usually related to the binding energy of Cooper pairs by [16] 1 In fact, contrary to the 3D case, 2D realistic interatomic attractive potentials always have a bound state. Both Equations (43) and (44) are ultraviolet divergent, but they exactly compensate the divergence of the mean-field grand potential Ω mg , which depends of the bare interaction strength g [16]. Figure 1 shows that taking into account only phase fluctuations (i.e., setting K = K θθ ) leads to a quite different behaviour of the velocity c s from that obtained by considering both phase and amplitude fluctuations [11]. In the upper panel, there is the behaviour of c s in the 2D BCS-BEC crossover, while in the lower panel there are the results for the 3D BCS-BEC crossover. Only in the deep BCS regime (left side of the two panels of Figure 1) do the two approaches give the same results, while the phase-only sound velocity diverges in the BEC regime (right side of panels). Thus, the effect of Goldstone-Higgs coupling K θσ is crucial in the study of the BCS-BEC crossover. It is important to stress that the zero-temperature results of Figure 1 are obtained adopting the mean-field number equation; i.e., We used Equation (45) instead of Equation (28) in the calculation of c s in order to have, from Equation (31), n (46) which is the expected phase stiffness at zero temperature, where the total density n should be equal to the superfluid density n s . In this way, one also satisfies the compressibility sum rule [38]. Here ǫ F =h 2 k 2 F /(2m) is the Fermi energy and v F = √ 2ǫ F /2 the Fermi velocity. Adapted from [11]. BCS: Bardeen-Cooper-Schrieffer; BEC: Bose-Einstein condensate.

Connection with the Popov's Hydrodynamic Action Functional
Given the Goldstone-Higgs action functional (30) and performing functional integration over the Higgs field σ(r, t), one obtains the Goldstone action whose Euler-Lagrange equation is exactly Equation (38) with K given by Equation (42). Remarkably, this Goldstone action functional can be immediately derived from the Popov's hydrodynamic action [39] functional integrating over the field ρ(r, t) of density fluctuations. ρ(r, t) represents a small space-time-dependent perturbation with respect to the constant and uniform density n. In the zero-temperature action functional (48), both J and K depend on n through the relationship between µ and n. From the Euler-Lagrange equations of (48) with respect to θ(r, t) and ρ(r, t), and introducing the velocity field v(r, t) =h 2m ∇θ(r, t) , one finds the familiar linearized hydrodynamic equations of Euler from which one immediately finds the d'Alambert equation for density fluctuations Thus, one can identify the velocity c s of propagation of the Goldstone mode with the velocity c 1 of the first sound of the fermionic superfluid [25][26][27]. In fact, according to the two-fluid theory of Laudau [40], a superfluid is characterized by the presence of the first sound (where superfluid and normal components oscillate in phase), but also the second sound, where superfluid and normal components oscillate with opposite phases. For the sake of completeness, we observe that the Euler Equations (50) and (51) can also be re-written in terms of a nonlinear Schrödinger equation for the complex field Ψ(r, t) = ρ(r, t) 1/2 e iθ(r,t) [36,37,41,42].

First Sound Velocity from Thermodynamics
On the basis of the two-fluid theory [40,43], at zero temperature the first sound velocity c 1 of a superfluid is simply given by where µ is the chemical potential, n is the number density, and V is the volume. Refs. [11,44] have shown that the velocity c s of the Goldstone mode given by Equation (41) coincides with the first sound velocity c 1 , given by Equation (53), if one uses the mean-field equation of state (45). At zero temperature, the difference between Equations (28) and (45) is small in the 3D case, but it is instead very large in the 2D case, and in particular in the BEC regime of the crossover [15]. This effect is clearly shown in Figure 2, where we report the zero-temperature first sound velocity c 1 in the 2D BCS-BEC crossover, regulated by the binding energy ǫ B of Cooper pairs. In Figure 2, the dot-dashed line is obtained by using the mean-field number Equation (45), while the solid line is based on the beyond-mean-field number Equation (28). In the strong coupling BEC regime, the beyond-mean-field equation of state is needed to accurately describe the thermodynamic quantities, and the corresponding first found velocity c 1 correctly goes to the correct composite boson limit [12]. Bighin and Salasnich [15] have shown that the first sound velocity c 1 obtained with the beyond-mean-field Gaussian theory (solid line of Figure 2) is in good agreement with very preliminary experimental data of 6 Li atoms [45]. Unfortunately, there are not yet fully reliable published experimental data on the sound velocity c 1 . However, there are reliable experimental data of the pressure P at very low temperature for the dilute gas of 6 Li atoms [7]. In Figure 3, we plot these experimental data (filled circles with error bars) and compare them with our beyond-mean-field theory at zero temperature (solid line). The pressure is immediately obtained as with Ω given by Equation (26). Figure 3 clearly shows that, at zero temperature, the agreement between experimental data and our beyond-mean-field theory is very good in the full BCS-BEC crossover. In the deep BEC regime, the beyond-mean-field pressure becomes [12] This formula (dashed line of Figure 3 is nothing else than the Popov equation of state [46] of weakly-interacting repulsive bosons where m B = 2m is the mass of composite bosons and µ B = 2(µ + ǫ B /2) is the bosonic chemical potential. Note that very recently we have obtained non-universal corrections to the Popov equation of state, taking account of finite-range effects of the inter-atomic potential [47]. Fig. 3. Adimensional pressure P/P id in the 2D BCS-BEC crossover, calculated using Gaussian pair fluctuation (GPF) theory (solid line) and using the Popov bosonic theory (dashed line). Experimental data (filled circles with error bars) are taken from [7]. P id is the pressure of the ideal 2D Fermi gas.

Conclusions
We have analyzed the derivation of a low-energy effective field theory of Goldstone and Higgs fields from the beyond-mean-field BCS theory of attractive fermions. We have shown that across the BCS-BEC crossover, the inclusion of the Goldstone-Higgs coupling is crucial to identify the velocity of the Goldstone mode with the first sound velocity one gets from the mean-field equation of state. Moreover, we have explicitly shown that in the BEC side of the 2D BCS-BEC crossover, the first sound velocity obtained from the beyond-mean-field equation of state is quite different with respect to the mean-field one. Finally, comparing our theoretical results with 2D experimental data, we have found that at zero temperature the beyond-mean-field theory based on Gaussian pair fluctuations seems reliable in the full 2D BCS-BEC crossover, giving the Popov equation of state in the deep BEC regime. However, as shown in Refs. [48,49], at low temperatures and in the deep BCS regime, quantized vortices and dark solitons obtained with effective field approaches based on low-energy expansion may contradict with the Bogoliubov-de Gennes theory, which is well-founded in the BCS regime.

Conflicts of Interest:
The author declares no conflict of interest.