Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0803.3238 · Physical Review A

Quantum Monte Carlo simulations of the BCS-BEC crossover at finite temperature

Aurel Bulgac, Joaquín E. Drut, Piotr Magierski

Published Aug 21, 2008• 173 citations
arXivHTMLPDF

Abstract

The quantum Monte Carlo method for spin-$\frac{1}{2}$ fermions at finite temperature is formulated for dilute systems with an $s$-wave interaction. The motivation and the formalism are discussed along with descriptions of the algorithm and various numerical issues. We report on results for the energy, entropy, and chemical potential as a function of temperature. We give upper bounds on the critical temperature ${T}_{c}$ for the onset of superfluidity, obtained by studying the finite-size scaling of the condensate fraction. All of these quantities were computed for couplings around the unitary regime in the range $\ensuremath{-}0.5\ensuremath{\leqslant}{({k}_{F}a)}^{\ensuremath{-}1}\ensuremath{\leqslant}0.2$, where $a$ is the $s$-wave scattering length and ${k}_{F}$ is the Fermi momentum of a noninteracting gas at the same density. In all cases our data are consistent with normal Fermi gas behavior above a characteristic temperature ${T}_{0}>{T}_{c}$, which depends on the coupling and is obtained by studying the deviation of the caloric curve from that of a free Fermi gas. For ${T}_{c}<T<{T}_{0}$ we find deviations from normal Fermi gas behavior that can be attributed to pairing effects. Low-temperature results for the energy and the pairing gap are shown and compared with Green-function Monte Carlo results by other groups.

DOI
10.1103/physreva.78.023625

Related papers

Pairing and Superfluid Properties of Dilute Fermion Gases at Unitarity

Nandini Trivedi

We study the system of a dilute gas of fermions in 3-dimensions, with attractive interactions tuned to the unitarity point, using the non-perturbative Restricted Path Integral Monte Carlo (R-PIMC) method. The pairing and superfluid properties of this system are calculated at finite temperature. The total energy at very low temperature from our results agrees closely with that of previous ground-state Quantum Monte Carlo calculations. We identify the temperature T � � 0.70ǫF below which pairing correlations develop, and estimate the critical temperature for the superfluid transition Tc � 0.25ǫF from a finite size scaling analysis of the superfluid density. PACS numbers: The recent experiments, starting with generating a degenerate gas of cold atoms[3], the use of Feshbach resonance to tune the effective interaction between the fermions[4, 5], the measurements of the gap to identify a pairing scale[6] and the measurements of vortices[7] to identify the superfluid state, have all given an impetus to the study of superfluidity in the BCS-BEC crossover regime. The nature of superconductivity or superfluidity in a many-particle system with an increasing pairing attraction was shown [1, 2] to interpolate smoothly between the BCS regime for weak attraction between the fermions, and the Bose-Einstein condensation (BEC) regime for strong coupling. In the weak coupling regime, the Cooper pair size is much larger than the interparticle spacing, and the simultaneous pairing and condensation of fermions is well described by the BCS theory. In the strong coupling regime, fermions form strongly bound bosonic molecules at a pairing temperature scale T ∗ that is considerably higher than the BEC temperature Tc. In a two species system of fermions with an attractive interaction between them, the unitary point is defined by the divergence of the zero energy s-channel scattering length, as, for two particles in free space. The unitary point is interesting since it describes a strongly correlated system with universal properties[8]. In general, there are two length scales that define the system: the interparticle spacing ∝ n −1/3 , and the scattering length as which contains information about the interaction potential. At the unitary point however, since as diverges, all the properties of the system are described by a single length scale k −1

Bulletin of the American Physical Society•Mar 6, 2007•0 citations

Quantum Monte Carlo study of dilute neutron matter at finite temperatures

Gabriel Wlazłowski, Piotr Magierski

We report results of fully nonperturbative, path integral Monte Carlo calculations for dilute neutron matter. The neutron-neutron interaction in the $s$ channel is parameterized by the scattering length and the effective range. We calculate the energy and the chemical potential as a function of temperature at density $\ensuremath{\rho}=0.003 {\text{fm}}^{\ensuremath{-}3}$. The critical temperature ${T}_{c}$ for the superfluid-normal phase transition is estimated from the finite size scaling of the condensate fraction. At low temperatures we extract the spectral weight function $A(\mathbf{p},\ensuremath{\omega})$ from the imaginary time propagator using the methods of maximum entropy and singular value decomposition. We determine the quasiparticle spectrum, which can be accurately parameterized by three parameters: an effective mass ${m}^{*}$, a mean-field potential $U$, and a gap $\ensuremath{\Delta}$. Large values of $\ensuremath{\Delta}/{T}_{c}$ indicate that the system is not a BCS-type superfluid at low temperatures.

arXiv:0912.0373•Physical Review C•Jan 31, 2011•13 citations

Quantum Monte Carlo studies of superfluid Fermi gases

Siyuan Chang, Vijay R. Pandharipande, Joseph Carlson +1 more

We report results of quantum Monte Carlo calculations of the ground state of dilute Fermi gases with attractive short-range two-body interactions. The strength of the interaction is varied to study different pairing regimes which are characterized by the product of the $s$-wave scattering length and the Fermi wave vector, $a{k}_{F}$. We report results for the ground-state energy, the pairing gap $\ensuremath{\Delta}$, and the quasiparticle spectrum. In the weak-coupling regime, $1∕a{k}_{F}<\ensuremath{-}1$, we obtain Bardeen-Cooper-Schrieffer (BCS) superfluid and the energy gap $\ensuremath{\Delta}$ is much smaller than the Fermi gas energy ${E}_{\mathrm{FG}}$. When $a>0$, the interaction is strong enough to form bound molecules with energy ${E}_{\mathrm{mol}}$. For $1∕a{k}_{F}\ensuremath{\gtrsim}0.5$, we find that weakly interacting composite bosons are formed in the superfluid gas with $\ensuremath{\Delta}$ and gas energy per particle approaching $\ensuremath{\mid}{E}_{\mathrm{mol}}\ensuremath{\mid}∕2$. In this region, we seem to have Bose-Einstein condensation (BEC) of molecules. The behavior of the energy and the gap in the BCS-to-BEC transition region, $\ensuremath{-}0.5<1∕a{k}_{F}<0.5$, is discussed.

arXiv:physics/0404115•Physical Review A•Oct 5, 2004•221 citations

Theory of ultracold atomic Fermi gases

S. Giorgini, Лев П. Питаевский, S. Stringari

The physics of quantum degenerate atomic Fermi gases in uniform as well as in harmonically trapped configurations is reviewed from a theoretical perspective. Emphasis is given to the effect of interactions that play a crucial role, bringing the gas into a superfluid phase at low temperature. In these dilute systems, interactions are characterized by a single parameter, the $s$-wave scattering length, whose value can be tuned using an external magnetic field near a broad Feshbach resonance. The BCS limit of ordinary Fermi superfluidity, the Bose-Einstein condensation (BEC) of dimers, and the unitary limit of large scattering length are important regimes exhibited by interacting Fermi gases. In particular, the BEC and the unitary regimes are characterized by a high value of the superfluid critical temperature, on the order of the Fermi temperature. Different physical properties are discussed, including the density profiles and the energy of the ground-state configurations, the momentum distribution, the fraction of condensed pairs, collective oscillations and pair-breaking effects, the expansion of the gas, the main thermodynamic properties, the behavior in the presence of optical lattices, and the signatures of superfluidity, such as the existence of quantized vortices, the quenching of the moment of inertia, and the consequences of spin polarization. Various theoretical approaches are considered, ranging from the mean-field description of the BCS-BEC crossover to nonperturbative methods based on quantum Monte Carlo techniques. A major goal of the review is to compare theoretical predictions with available experimental results.

arXiv:0706.3360•Reviews of Modern Physics•Oct 2, 2008•1,997 citations

Pairing Correlations across the Superfluid Phase Transition in the Unitary Fermi Gas

Scott M. Jensen, C. N. Gilbreth, Y. Alhassid

In the two-component Fermi gas with a contact interaction, a pseudogap regime can exist at temperatures between the superfluid critical temperature T_{c} and a temperature T^{*}>T_{c}. This regime is characterized by pairing correlations without superfluidity. However, in the unitary limit of infinite scattering length, the existence of this regime is still debated. To help address this, we have applied finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods to study the thermodynamics of the superfluid phase transition and signatures of the pseudogap in the spin-balanced homogeneous unitary Fermi gas. We present results at finite filling factor ν≃0.06 for the condensate fraction, an energy-staggering pairing gap, the spin susceptibility, and the heat capacity, and compare them to experimental data when available. In contrast to previous AFMC simulations, our model space consists of the complete first Brillouin zone of the lattice, and our calculations are performed in the canonical ensemble of fixed particle number. The canonical ensemble AFMC framework enables the calculation of a model-independent gap, providing direct information on pairing correlations without the need for numerical analytic continuation. We use finite-size scaling to estimate T_{c} at the corresponding filling factor. We find that the energy-staggering pairing gap vanishes above T_{c}, showing no pseudogap effects, and that the spin susceptibility shows a substantially reduced signature of a spin gap compared to previously reported AFMC simulations.

arXiv:1801.06163•Physical Review Letters•Mar 2, 2020•32 citations

Emergence of a Pseudogap in the BCS-BEC Crossover

Adam Richie-Halford, Joaquín E. Drut, Aurel Bulgac

Strongly correlated Fermi systems with pairing interactions become superfluid below a critical temperature T_{c}. The extent to which such pairing correlations alter the behavior of the liquid at temperatures T>T_{c} is a subtle issue that remains an area of debate, in particular regarding the appearance of the so-called pseudogap in the BCS-BEC crossover of unpolarized spin-1/2 nonrelativistic matter. To shed light on this, we extract several quantities of crucial importance at and around the unitary limit, namely, the odd-even staggering of the total energy, the spin susceptibility, the pairing correlation function, the condensate fraction, and the critical temperature T_{c}, using a nonperturbative, constrained-ensemble quantum Monte Carlo algorithm.

arXiv:2004.05014•Physical Review Letters•Aug 4, 2020•30 citations
The quantum Monte Carlo method for spin-12\frac{1}{2}21​ fermions at finite temperature is formulated for dilute systems with an sss-wave interaction. The motivation and the formalism are discussed along with descriptions of the algorithm and various numerical issues. We report on results for the energy, entropy, and chemical potential as a function of temperature. We give upper bounds on the critical temperature Tc{T}_{c}Tc​ for the onset of superfluidity, obtained by studying the finite-size scaling of the condensate fraction. All of these quantities were computed for couplings around the unitary regime in the range \ensuremath−0.5\ensuremath⩽(kFa)\ensuremath−1\ensuremath⩽0.2\ensuremath{-}0.5\ensuremath{\leqslant}{({k}_{F}a)}^{\ensuremath{-}1}\ensuremath{\leqslant}0.2\ensuremath−0.5\ensuremath⩽(kF​a)\ensuremath−1\ensuremath⩽0.2, where aaa is the sss-wave scattering length and kF{k}_{F}kF​ is the Fermi momentum of a noninteracting gas at the same density. In all cases our data are consistent with normal Fermi gas behavior above a characteristic temperature T0>Tc{T}_{0}>{T}_{c}T0​>Tc​, which depends on the coupling and is obtained by studying the deviation of the caloric curve from that of a free Fermi gas. For Tc<T<T0{T}_{c}<T<{T}_{0}Tc​<T<T0​ we find deviations from normal Fermi gas behavior that can be attributed to pairing effects. Low-temperature results for the energy and the pairing gap are shown and compared with Green-function Monte Carlo results by other groups.
We study the system of a dilute gas of fermions in 3-dimensions, with attractive interactions tuned to the unitarity point, using the non-perturbative Restricted Path Integral Monte Carlo (R-PIMC) method. The pairing and superfluid properties of this system are calculated at finite temperature. The total energy at very low temperature from our results agrees closely with that of previous ground-state Quantum Monte Carlo calculations. We identify the temperature T � � 0.70ǫF below which pairing correlations develop, and estimate the critical temperature for the superfluid transition Tc � 0.25ǫF from a finite size scaling analysis of the superfluid density. PACS numbers: The recent experiments, starting with generating a degenerate gas of cold atoms[3], the use of Feshbach resonance to tune the effective interaction between the fermions[4, 5], the measurements of the gap to identify a pairing scale[6] and the measurements of vortices[7] to identify the superfluid state, have all given an impetus to the study of superfluidity in the BCS-BEC crossover regime. The nature of superconductivity or superfluidity in a many-particle system with an increasing pairing attraction was shown [1, 2] to interpolate smoothly between the BCS regime for weak attraction between the fermions, and the Bose-Einstein condensation (BEC) regime for strong coupling. In the weak coupling regime, the Cooper pair size is much larger than the interparticle spacing, and the simultaneous pairing and condensation of fermions is well described by the BCS theory. In the strong coupling regime, fermions form strongly bound bosonic molecules at a pairing temperature scale T ∗ that is considerably higher than the BEC temperature Tc. In a two species system of fermions with an attractive interaction between them, the unitary point is defined by the divergence of the zero energy s-channel scattering length, as, for two particles in free space. The unitary point is interesting since it describes a strongly correlated system with universal properties[8]. In general, there are two length scales that define the system: the interparticle spacing ∝ n −1/3 , and the scattering length as which contains information about the interaction potential. At the unitary point however, since as diverges, all the properties of the system are described by a single length scale k −1
We report results of fully nonperturbative, path integral Monte Carlo calculations for dilute neutron matter. The neutron-neutron interaction in the sss channel is parameterized by the scattering length and the effective range. We calculate the energy and the chemical potential as a function of temperature at density $\ensuremath{\rho}=0.003 {\text{fm}}^{\ensuremath{-}3}$. The critical temperature Tc{T}_{c}Tc​ for the superfluid-normal phase transition is estimated from the finite size scaling of the condensate fraction. At low temperatures we extract the spectral weight function $A(\mathbf{p},\ensuremath{\omega})$ from the imaginary time propagator using the methods of maximum entropy and singular value decomposition. We determine the quasiparticle spectrum, which can be accurately parameterized by three parameters: an effective mass m∗{m}^{*}m∗, a mean-field potential UUU, and a gap $\ensuremath{\Delta}$. Large values of $\ensuremath{\Delta}/{T}_{c}$ indicate that the system is not a BCS-type superfluid at low temperatures.
We report results of quantum Monte Carlo calculations of the ground state of dilute Fermi gases with attractive short-range two-body interactions. The strength of the interaction is varied to study different pairing regimes which are characterized by the product of the sss-wave scattering length and the Fermi wave vector, akFa{k}_{F}akF​. We report results for the ground-state energy, the pairing gap $\ensuremath{\Delta}$, and the quasiparticle spectrum. In the weak-coupling regime, $1∕a{k}_{F}<\ensuremath{-}1$, we obtain Bardeen-Cooper-Schrieffer (BCS) superfluid and the energy gap $\ensuremath{\Delta}$ is much smaller than the Fermi gas energy EFG{E}_{\mathrm{FG}}EFG​. When a>0a>0a>0, the interaction is strong enough to form bound molecules with energy Emol{E}_{\mathrm{mol}}Emol​. For $1∕a{k}_{F}\ensuremath{\gtrsim}0.5$, we find that weakly interacting composite bosons are formed in the superfluid gas with $\ensuremath{\Delta}$ and gas energy per particle approaching $\ensuremath{\mid}{E}_{\mathrm{mol}}\ensuremath{\mid}∕2$. In this region, we seem to have Bose-Einstein condensation (BEC) of molecules. The behavior of the energy and the gap in the BCS-to-BEC transition region, $\ensuremath{-}0.5<1∕a{k}_{F}<0.5$, is discussed.
The physics of quantum degenerate atomic Fermi gases in uniform as well as in harmonically trapped configurations is reviewed from a theoretical perspective. Emphasis is given to the effect of interactions that play a crucial role, bringing the gas into a superfluid phase at low temperature. In these dilute systems, interactions are characterized by a single parameter, the sss-wave scattering length, whose value can be tuned using an external magnetic field near a broad Feshbach resonance. The BCS limit of ordinary Fermi superfluidity, the Bose-Einstein condensation (BEC) of dimers, and the unitary limit of large scattering length are important regimes exhibited by interacting Fermi gases. In particular, the BEC and the unitary regimes are characterized by a high value of the superfluid critical temperature, on the order of the Fermi temperature. Different physical properties are discussed, including the density profiles and the energy of the ground-state configurations, the momentum distribution, the fraction of condensed pairs, collective oscillations and pair-breaking effects, the expansion of the gas, the main thermodynamic properties, the behavior in the presence of optical lattices, and the signatures of superfluidity, such as the existence of quantized vortices, the quenching of the moment of inertia, and the consequences of spin polarization. Various theoretical approaches are considered, ranging from the mean-field description of the BCS-BEC crossover to nonperturbative methods based on quantum Monte Carlo techniques. A major goal of the review is to compare theoretical predictions with available experimental results.
In the two-component Fermi gas with a contact interaction, a pseudogap regime can exist at temperatures between the superfluid critical temperature T_{c} and a temperature T^{*}>T_{c}. This regime is characterized by pairing correlations without superfluidity. However, in the unitary limit of infinite scattering length, the existence of this regime is still debated. To help address this, we have applied finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods to study the thermodynamics of the superfluid phase transition and signatures of the pseudogap in the spin-balanced homogeneous unitary Fermi gas. We present results at finite filling factor ν≃0.06 for the condensate fraction, an energy-staggering pairing gap, the spin susceptibility, and the heat capacity, and compare them to experimental data when available. In contrast to previous AFMC simulations, our model space consists of the complete first Brillouin zone of the lattice, and our calculations are performed in the canonical ensemble of fixed particle number. The canonical ensemble AFMC framework enables the calculation of a model-independent gap, providing direct information on pairing correlations without the need for numerical analytic continuation. We use finite-size scaling to estimate T_{c} at the corresponding filling factor. We find that the energy-staggering pairing gap vanishes above T_{c}, showing no pseudogap effects, and that the spin susceptibility shows a substantially reduced signature of a spin gap compared to previously reported AFMC simulations.
Strongly correlated Fermi systems with pairing interactions become superfluid below a critical temperature T_{c}. The extent to which such pairing correlations alter the behavior of the liquid at temperatures T>T_{c} is a subtle issue that remains an area of debate, in particular regarding the appearance of the so-called pseudogap in the BCS-BEC crossover of unpolarized spin-1/2 nonrelativistic matter. To shed light on this, we extract several quantities of crucial importance at and around the unitary limit, namely, the odd-even staggering of the total energy, the spin susceptibility, the pairing correlation function, the condensate fraction, and the critical temperature T_{c}, using a nonperturbative, constrained-ensemble quantum Monte Carlo algorithm.