Superfluidity in the cold atomic two-species Fermi gas system in the unitary limit of infinite scattering length remains incompletely understood. In particular, a pseudogap phase has been proposed to exist above the superfluid critical temperature. Here we apply the auxiliary-field quantum Monte Carlo method to perform ab initio calculations of the temperature dependence of three quantities---the energy-staggering pairing gap, the condensate fraction, and the heat capacity---in a trapped finite-size cold atom system in the unitary limit of infinite scattering length. As the calculations of the energy-staggering pairing gap require the use of the canonical ensemble, we employ an efficient algorithm for the stabilization of particle-number projection that is essential for reaching convergence in the size of the model space. We observe clear signatures of the superfluid phase transition in all three quantities, including a signature of the recently measured $\ensuremath{\lambda}$ peak in the heat capacity. We also find that the energy-staggering gap does not lead the condensate fraction as the temperature decreases.
Superfluidity in the cold atomic two-species Fermi gas system in the unitary limit of infinite scattering length remains incompletely understood. In particular, a pseudogap phase has been proposed to exist above the superfluid critical temperature. Here we apply the auxiliary-field quantum Monte Carlo method to perform ab initio calculations of the temperature dependence of three quantities---the energy-staggering pairing gap, the condensate fraction, and the heat capacity---in a trapped finite-size cold atom system in the unitary limit of infinite scattering length. As the calculations of the energy-staggering pairing gap require the use of the canonical ensemble, we employ an efficient algorithm for the stabilization of particle-number projection that is essential for reaching convergence in the size of the model space. We observe clear signatures of the superfluid phase transition in all three quantities, including a signature of the recently measured \ensuremathλ peak in the heat capacity. We also find that the energy-staggering gap does not lead the condensate fraction as the temperature decreases.
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.
In the two-component Fermi gas with contact interactions a pseudogap regime, in which pairing correlations are present without superfluidity, can exist at temperatures between the superfluid critical temperature Tc and a temperature T∗>Tc. However, the existence of a pseudogap in the unitary limit of infinite scattering length is debated. To help address this issue, we have used finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods to study the thermodynamics of the spin-balanced homogeneous unitary Fermi gas on a lattice. We present results for the thermal energy, heat capacity, condensate fraction, a model-independent pairing gap, and spin susceptibility, and compare them to experimental data when available. 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. We find that the energy-staggering pairing gap vanishes above Tc and that T∗ at unitarity, as determined from the spin susceptibility, is lower than previously reported in AFMC simulations.
The unitary Fermi gas (UFG) is a strongly correlated system of two-species (spin-1/2) fermions with a short-range attractive interaction modeled by a contact interaction and has attracted much interest across different disciplines. The UFG is considered a paradigm for strongly correlated superfluids and has been investigated extensively, with generally good agreement found between theory and experiment. However, the extent of a pseudogap regime above the critical temperature Tc for superfluidity is still debated both theoretically and experimentally. Here we study thermodynamic properties of the UFG across the superfluid phase transition using finite-temperature lattice auxiliary-field quantum Monte Carlo (AFMC) methods in the canonical ensemble of fixed particle numbers. We extrapolate our lattice AFMC results to the continuous time and continuum limits, thus removing the systematic error associated with the finite filling factor of previous AFMC studies. We determine the critical temperature to be Tc=0.16(1)TF. For the largest particle number studied N=114, the energy-staggering pairing gap is suppressed above a pairing scale temperature of T∗≈0.2TF. The spin susceptibility displays moderate suppression above Tc with a spin gap temperature of Ts≈0.2TF. We also calculate a free energy-staggering pairing gap, which shows substantially reduced statistical errors when compared with the energy-staggering gap, allowing for a clear signature of pairing correlations in the finite-size system. All results indicate that the pseudogap regime is narrow, with pseudogap signatures emerging at temperatures below T∗≈0.2TF. The reduced statistical errors of the free energy gap enable an extrapolation at low temperatures, allowing an estimate of the zero-temperature pairing gap ΔE=0.576(24)εF.
The unitary Fermi gas is a many-body system of two-component fermions with zero-range interactions tuned to infinite scattering length. Despite much activity and interest in unitary Fermi gas and its universal properties, there have been great difficulties in performing accurate calculations of the superfluid condensate fraction and pairing wave function. In this paper, we present auxiliary-field lattice Monte Carlo simulations using a lattice interaction which accelerates the approach to the continuum limit, thereby allowing for robust calculations of these difficult observables. As a benchmark test, we compute the ground-state energy of 33 spin-up and 33 spin-down particles. As a fraction of the free Fermi gas energy EFG, we find E0/EFG=0.369(2),0.372(2), using two different definitions of the finite-system energy ratio, in agreement with the latest theoretical and experimental results. We then determine the condensate fraction by measuring off-diagonal long-range order in the two-body density matrix. We find that the fraction of condensed pairs is $\ensuremath{\alpha}=0.43(2)$. We also extract the pairing wave function and find the pair correlation length to be ${\ensuremath{\zeta}}_{p}{k}_{F}=1.8(3)\ensuremath{\hbar}$, where kF is the Fermi momentum. Provided that the simulations can be performed without severe sign oscillations, the methods we present here can be applied to superfluid neutron matter as well as more exotic P-wave and D-wave superfluids.
The two-species Fermi gas with attractive short-range interactions in two spatial dimensions provides a paradigmatic system for the understanding of strongly correlated Fermi superfluids in two dimensions. It is known to exhibit a BEC to BCS crossover as a function of ln(k_{F}a), where a is the scattering length, and to undergo a Berezinskii-Kosterlitz-Thouless superfluid transition below a critical temperature T_{c}. However, the extent of a pseudogap regime in the strongly correlated regime of ln(k_{F}a)∼1, in which pairing correlations persist above T_{c}, remains largely unexplored with controlled theoretical methods. Here, we use finite-temperature auxiliary-field quantum Monte Carlo methods on discrete lattices in the canonical ensemble formalism to calculate thermodynamical observables in the strongly correlated regime. We extrapolate to continuous time and the continuum limit to eliminate systematic errors and present results for particle numbers ranging from N=42 to N=162. We estimate T_{c} by a finite-size scaling analysis, and observe clear pseudogap signatures above T_{c} and below a temperature T^{*} in both the spin susceptibility and free-energy gap. We also present results for the contact, a fundamental thermodynamic property of quantum many-body systems with short-range interactions.
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