Strong coupling effects on specific heat in the BCS-BEC crossover

TL;DR

This paper theoretically examines how strong-coupling effects influence the specific heat in a superfluid Fermi gas across the BCS-BEC crossover, highlighting the role of Goldstone modes in low-temperature behavior.

Contribution

It introduces a calculation of specific heat considering superfluid fluctuations at unitarity, revealing the dominance of Goldstone modes at low temperatures.

Findings

Low-temperature specific heat exhibits T^3 behavior.

Goldstone mode fluctuations dominate low-energy excitations.

Results enhance understanding of pairing fluctuations in superfluid Fermi gases.

Abstract

We theoretically investigate strong-coupling effects on specific heat at constant volume $C_{\rm V}$ in a superfluid Fermi gas with a tunable interaction associated with Feshbach resonance. Including fluctuations of the superfluid order parameter within the strong-coupling theory developed by Nozi\`eres and Schmitt-Rink, we calculate the temperature dependence of $C_{\rm V}$ at the unitarity limit in the superfluid phase. We show that, in the low temperature region, $T^3$-behavior is shown in the temperature dependence of $C_{\rm V}$. This result indicates that the low-lying excitations are dominated by the gapless Goldstone mode, associated with the phase fluctuations of the superfluid order parameter. Since the Goldstone mode is one of the most fundamental phenomena in the Fermionic superfluidity, our results are useful for further understanding how the pairing fluctuations affects physical properties in the BCS-BEC crossover physics below the superfluid transition temperature.

Equations

$$H=\sum _{\mathbf{p},\sigma }{\xi }_{\mathbf{p}}{c}_{\mathbf{p},\sigma }^{†}{c}_{\mathbf{p},\sigma }-U\sum _{\mathbf{p},{\mathbf{p}}^{\prime },\mathbf{q}}{c}_{\mathbf{p}+\mathbf{q}/2,\uparrow }^{†}{c}_{-\mathbf{p}+\mathbf{q}/2,\downarrow }^{†}{c}_{-{\mathbf{p}}^{\prime }+\mathbf{q}/2,\downarrow }{c}_{{\mathbf{p}}^{\prime }+\mathbf{q}/2,\uparrow }.$$

$$\frac{4\pi a_s}{m}=\frac{-U}{1-U{\sum }_{\mathbf{p}}^{p_{\mathrm{c}}}\frac{1}{2{\epsilon }_{\mathbf{p}}}}.$$

$$H=\sum _{\mathbf{p}}{\psi }_{\mathbf{p}}^{†}\left({\xi }_{\mathbf{p}}{\tau }_3-\Delta {\tau }_1\right){\psi }_{\mathbf{p}}-U\sum _{\mathbf{q}}{\rho }_{+}\left(\mathbf{q}\right){\rho }_{-}\left(-\mathbf{q}\right).$$

$${\rho }_{\pm }\left(\mathbf{q}\right)=\sum _{\mathbf{p}}{\psi }_{\mathbf{p}+\frac{\mathbf{q}}{2}}^{†}{\tau }_{\pm }{\psi }_{\mathbf{p}-\frac{\mathbf{q}}{2}}$$

$${\Omega }_{\mathrm{MF}}=-\frac{m{\Delta }^2}{4\pi a_s}+\sum _{\mathbf{p}}\left({\xi }_{\mathbf{p}}-E_{\mathbf{p}}+\frac{{\Delta }^2}{2{\epsilon }_{\mathbf{p}}}\right)-2T\sum _{\mathbf{p}}\ln \left(1+e^{-\frac{E_{\mathbf{p}}}{T}}\right),$$

$$\delta \Omega =-\frac{1}{2\beta }\sum _{\mathbf{q},i{\nu }_n}\mathrm{Tr}\left[\ln \hat{\Gamma }(\mathbf{q},i{\nu }_n)\right].$$

$$\hat{\Pi }(\mathbf{q},i{\nu }_n)$$

$${\Pi }_{s{s}^{\prime }}(\mathbf{q},i{\nu }_n)$$

$$G_0\left(\mathbf{p},i{\omega }_n\right)=\frac{1}{i{\omega }_n-{\xi }_{\mathbf{p}}{\tau }_3+\Delta {\tau }_1}$$

$$C_{\mathrm{V}}={\left(\frac{\partial E}{\partial T}\right)}_N.$$

$$E=\Omega +TS+\mu N=\Omega -T{\left(\frac{\partial \Omega }{\partial T}\right)}_{\mu }-\mu {\left(\frac{\partial \Omega }{\partial \mu }\right)}_T.$$

$$1=\frac{4\pi a_s}{m}\sum _{\mathbf{p}}\left(\frac{1}{2E_{\mathbf{p}}}\tanh \frac{2E_{\mathbf{p}}}{2T}-\frac{1}{2{\epsilon }_{\mathbf{p}}}\right),$$

$$N$$

$$v_{\varphi }$$

$$N_{\mathrm{G}}=-{\left(\frac{\partial \Omega }{\partial \mu }\right)}_{T,\Delta }$$

Full text

1 Introduction

Since the BCS-BEC crossover, where the superfluid properties continuously change from the weak-coupling BCS (Bardeen-Cooper-Schrieffer)-type to the BEC (Bose-Einstein condensation) of tightly binding molecular Bose gas as increasing the interaction strength, has been realized in 40KRegal and 6LiZwierlein2004 ; Kinast2004 ; Bartenstein2004 Fermi gases, strong-coupling effects on various physical quantities in this system, such as single-particle propertiesTsuchiya2009 ; Tsuchiya2011 , spin susceptibilityTajima2014 , and thermodynamic propertiesHaussmann ; Hu2006pra ; Pieter2016 , have extensively investigated. Recently, in ultracold Fermi gas, the great progress of experimental technique makes us enable to measure various thermodynamic quantitiesKinast2005 ; Luo2009 ; Horikochi2010 ; Nascimbene2010 ; Ku .

In the BCS-BEC crossover phenomena below the superfluid transition temperature, the low-energy excitations give us a useful information for further understanding this phenomenon. While in the weak-coupling BCS regime there are two kinds of fundamental excitations, i.e. the Fermionic single-particle excitations associated with pair breaking and the gapless Goldstone mode associated with the phase fluctuations of the superfluid order parameter $\Delta$, in the strong-coupling BEC regime the former is strongly suppressed due to the large binding energy of the pairs. Although in the low temperature region where $T\ll\Delta$ the low-lying excitations are always dominated by the Goldstone mode in the entire interaction regime, since the sound velocity of this mode strongly depends on the interaction strengthOhashi2003 , a characteristic temperature where the system starts to be dominated by the Goldstone mode can be used to determine the region where the phase fluctuations become remarkable.

For this purpose, we focus on the specific heat at constant volume $C_{\rm V}$, which has been known to be sensitive to pairing fluctuations above the superfluid transition temperature $T_{\rm c}$Pieter2016 . As well known, while within the weak-coupling mean-field theory in which the gapless Goldstone mode is ignored, $C_{\rm V}$ is dominated by the pair breaking leading $e^{-\Delta/T}$-behavior in the low temperature region, in the ideal Bose gas, $C_{\rm V}$ shows $T^{3/2}$-behavior. On the other hand, when the thermodynamic properties are dominated by the gapless Goldstone mode, $C_{\rm V}$ should be proportional to $T^{3}$Haussmann . Thus, the temperature dependence of $C_{\rm V}$ is sensitive to what kind of excitations is remarkable. In this sense, from the detailed temperature dependence of $C_{\rm V}$, we might determine the region where the gapless Goldstone mode associated with the phase fluctuations dominates the thermodynamic properties.

In this paper, we theoretically calculate $C_{\rm V}$ at the unitarity limit below the superfluid transition temperature $T_{\rm c}$ within a strong-coupling NSR theory developed by Nozières and Schmitt-RinkHu2006pra ; Ohashi2003 ; NSR ; Fukushima2007 ; Hu2006epl , in which the contributions from the gapless Goldstone mode are taken into account. We find that in the low temperature region where the thermal transfer from the gapless Goldstone mode to the single-particle excitations is sufficiently suppressed, $C_{\rm V}$ shows $T^{3}$-behavior. Throughout this paper, we take $\hbar=k_{\rm B}=1$, and the system volume $V$ is taken to be unity, for simplicity.

2 Formulation

We consider a two-component Fermi gas with a contact-type attractive interaction associated with Feshbach resonance, described by the Hamiltonian

[TABLE]

Here $c_{\bm{p},\sigma}$ denotes an annihilation operator of a Fermi atom with the momentum ${\bm{p}}$ and the pseudospin $\sigma=\uparrow,\downarrow$, and $\xi_{\bm{p}}=p^{2}/(2m)-\mu$ is the kinetic energy measured from the chemical potential $\mu$ (where $m$ is the atomic mass). The second term in the Hamiltonian Eq. (1) describes the contact-type attractive interaction with a tunable coupling constant $U>0$, which is conveniently measured in terms of the observable $s$-wave scattering length as

[TABLE]

Here, $p_{\rm c}$ is a cutoff momentum. Using this scale, the weak- and strong-coupling region are characterized by the region where $(k_{\rm F}a_{s})^{-1}<-1$ and $(k_{\rm F}a_{s})^{-1}>1$, respectively. The region between them ($-1<(k_{\rm F}a_{s})^{-1}<1$) is usually referred as the crossover region, in which pairing fluctuations might be remarkable.

To investigate the superfluid properties of this system, as usual, it is convenient to write the Hamiltonian Eq. (1) in the Nambu representationOhashi2003 by introducing the superfluid order parameter $\Delta=-U\sum_{\bm{p}}\left\langle c_{-{\bm{p}},\downarrow}c_{{\bm{p}},\uparrow}\right\rangle$, as

[TABLE]

Here $\psi_{\bm{p}}=\left(c_{{\bm{p}},\uparrow}c_{-{\bm{p}},\downarrow}^{\dagger}\right)^{T}$ is the two-component Nambu spinor operator, $\tau_{i}$ ($i=1,2,3$) is Pauli matrices acting on the particle-hole space, and

[TABLE]

is the generalized density operator describing fluctuations of $\Delta$, where $\tau_{\pm}=\frac{1}{2}\left(\tau_{1}\pm i\tau_{2}\right)$. In Eq. (3), $\Delta$ is taken to be real without loss of generality.

The strong-coupling corrections from the second term in Eq. (3) to the thermodynamic properties are taken into account within NSR theoryHu2006pra ; Ohashi2003 ; NSR ; Fukushima2007 ; Hu2006epl , in which the thermodynamic potential is given as the sum $\Omega=\Omega_{\rm{MF}}+\delta\Omega$ of the mean-field part $\Omega_{\rm MF}$ and the strong-coupling corrections $\delta\Omega$. The mean field part is given by

[TABLE]

where $E_{\bm{p}}=\sqrt{\xi_{\bm{p}}^{2}+|\Delta|^{2}}$ is the Bogoliubov single-particle excitation spectrum. The strong-coupling corrections to the thermodynamic potential $\delta\Omega$ are described by the Feynman diagrams shown in Fig. 1, which give

[TABLE]

Here, $\hat{\Gamma}({\bm{q}},i\nu_{n})=-U/[1+U\hat{\Pi}({\bm{q}},i\nu_{n})]$ is the 2$\times$2 particle-particle scattering matrix, where

[TABLE]

is the lowest-order 2$\times$2 matrix pair correlation function describing the fluctuation of the superfluid order parameter, and

[TABLE]

is the 2$\times$2 matrix single-particle Green’s function in the mean-field level.

The specific heat at constant volume $C_{\rm V}$ is obtained from a thermodynamic relation

[TABLE]

Here $E$ is the internal energy, which is given by the Legendre transformation from $\Omega$ as

[TABLE]

In this paper, by numerically carrying out the derivative of $E$ with respect to the temperature $T$ in Eq. (12), $C_{\rm V}$ is calculated. In this procedure, we first determine the chemical potential $\mu$ and the superfluid order parameter $\Delta$ as a function of $T$ by self-consistently solving the gap equation, which is obtained from the Thouless criterion ${\rm det}\left[\Gamma({\bm{q}}=0,i\nu_{n}=0)^{-1}\right]=0$ as

[TABLE]

together with the particle number equation

[TABLE]

We note that the Thouless criterion guarantees the existence of the gapless Goldstone mode, of which the dispersion relation is obtained by solving the equation ${\rm{det}}\left[\Gamma\left({\bm{q}},i\nu_{n}\to\omega+i\delta\right)^{-1}\right]=0$. In the low temperature limit ($T\ll\left\|\mu\right\|,\Delta$), indeed, we obtain the well-known Goldstone mode $\omega_{q}=v_{\phi}q$ in the low energy regionOhashi2003 , with the sound velocity

[TABLE]

which is expected to give a $T^{3}$ contribution to the specific heat. We also briefly mention that the superfluid transition temperature $T_{\rm c}$ is obtained by solving Eqs. (14) and (15) with $\Delta=0$, and $\mu$ above $T_{\rm c}$ is simply given by solving only Eq. (15) with $\Delta=0$. In the next section, in addition to the results in the superfluid phase, we also show the results in the normal phasePieter2016 , for comparison.

3 Result

Figure 2 shows the temperature dependence the chemical potential $\mu$, as well as the superfluid order parameter $\Delta$ in the unitarity limit, determined from Eqs. (14) and (15). Here, we note that $\Delta$ in Fig. 2 has finite value even at $T=T_{\rm c}$. However, it is known that this behavior indicating a first-order phase transition is an artifact in our approximation. Same problem has been reported in other diagrammatic strong-coupling theories, such as non-self-consistent $T$-matrix theoryTMA , as well as self-consistent $T$-matrix theorySCTMA . Since this unphysical behavior of $\Delta$ might affect our results just below $T_{\rm c}$, we will focus on the low-temperature behavior of $C_{\rm V}$, and we leave this problem as a future work. We also note that the obtained $\Delta(T)$ is always larger than $T$ below $T_{\rm c}$. Thus, in the superfluid phase at the unitarity limit, the single-particle excitations associated with the pair-breaking, which is characterized by an energy scale $2\Delta$, are thermally suppressed, and the thermodynamics of this system is expected to be dominated by the gapless Goldstone mode, as well as the thermal transfer from the collective mode to the single-particle excitations.

Figure 3 shows calculated temperature dependence of $C_{\rm V}$ at the unitarity limit. Starting from the high temperature region, in normal phase $C_{V}$ is found to increase as approaching $T_{\rm c}$, that is in contrast to one of the ideal Fermi gas where $C_{\rm V}$ monotonically decreases as decreasing $T$. As discussed in our previous paperPieter2016 , this enhancement of $C_{\rm V}$ near $T_{\rm c}$ originates from the appearance of the preformed Cooper pairs associated with a strong-attractive interaction. We also find that a jump of $C_{\rm V}$ across $T_{\rm c}=0.22T_{\rm F}$, as shown in the ordinary BCS superconductors. However, as mentioned above, since our results just below $T_{\rm c}$ are not reliable due to the artificial first-order phase transition, we do not further discuss this point.

In the superfluid phase below $T_{\rm c}$, $C_{\rm V}$ monotonically decreases as decreasing $T$, and eventually vanishes at $T=0$, as expected. As shown in Fig 3 (b), in the low temperature region ($T\lesssim 0.03T_{\rm F}$), we find that $C_{\rm V}$ is proportional to $T^{3}$, that indicates that $C_{\rm V}$ is dominated by the gapless Goldstone mode.

To more clearly see this, we conveniently write the particle number Eq. (15) as the sum $N=N_{0}+N_{\rm G}$ of the contribution from the Goldstone mode

[TABLE]

and the others $N_{0}$. Here, we mention that, although $N_{\rm G}$ includes the contribution from both the phase and the amplitude fluctuations of the superfluid order parameter, because the amplitude fluctuations are rapidly suppressed as developing $\Delta$, at least, in the low-temperature region $T\ll\Delta$, $N_{\rm G}$ can be regarded as the contribution from the Goldstone mode. Noting that the single-particle excitations associated with pair breaking can also be ignored in this low temperature region, the low-temperature behavior of $C_{\rm V}$ originates from the gapless Goldstone mode, as well as the thermal transfer from the Goldstone mode to the single-particle contributions. As shown in Fig. 4, as decreasing $T$, the latter effects are gradually suppressed, and when $T\lesssim 0.05T_{\rm F}$, $N_{\rm G}$ becomes almost constant. Then, $T^{3}$-dependence coming from the Goldstone mode becomes dominant in $C_{\rm V}$.

4 Conclusion

In this paper, we have theoretically investigated the effects of fluctuations of the superfluid order parameter on the specific heat at constant volume $C_{\rm V}$ at the unitarity limit within a strong-coupling NSR theory. We found that in the low-temperature region, where the thermal transfer from the gapless Goldstone mode to the single-particle excitations is sufficiently suppressed, $C_{\rm V}$ exhibits a $T^{3}$-dependence. Since the $T^{3}$-dependence comes from the gapless Goldstone mode associated with the phase fluctuations of the superfluid order parameter, our results indicate that the temperature dependence of $C_{\rm V}$ might be useful to determine the region where the phase fluctuations dominate the thermodynamic properties. Furthermore, since the gapless Goldstone mode always exists in the superfluid phase, but its properties remarkably depend on the interaction strength, it is our future problem how the temperature dependence of $C_{\rm V}$ changes as varying the interaction strength.

Acknowledgements.

This work was supported by KiPAS project in Keio University. DI was supported by Grant-in-aid for Scientific Research from JSPS in Japan (No.JP16K17773). YO was supported by Grant-in-aid for Scientific Research from MEXT and JSPS in Japan (No.JP18K11345, No.JP18H05406, No.JP16K05503).