We study the system of a dilute gas of fermions in three dimensions, with attractive interactions tuned to the unitarity point, using the nonperturbative restricted path-integral Monte Carlo 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}^{*}\ensuremath{\approx}0.70{ϵ}_{F}$, below which pairing correlations develop, and estimate the critical temperature for the superfluid transition ${T}_{c}\ensuremath{\approx}0.245{ϵ}_{F}$ from a finite size scaling analysis of the superfluid density.
We study the system of a dilute gas of fermions in three dimensions, with attractive interactions tuned to the unitarity point, using the nonperturbative restricted path-integral Monte Carlo 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∗\ensuremath≈0.70ϵF, below which pairing correlations develop, and estimate the critical temperature for the superfluid transition Tc\ensuremath≈0.245ϵF from a finite size scaling analysis of the superfluid density.
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
By using exact path-integral Monte Carlo methods we calculate the equation of state of an interacting Bose gas as a function of temperature both below and above the superfluid transition. The universal character of the equation of state for dilute systems and low temperatures is investigated by modeling the interatomic interactions using different repulsive potentials corresponding to the same s-wave scattering length. The results obtained for the energy and the pressure are compared to the virial expansion for temperatures larger than the critical temperature. At very low temperatures we find agreement with the ground-state energy calculated using the diffusion Monte Carlo method.
The quantum Monte Carlo method for spin-21 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 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{\leqslant}{({k}_{F}a)}^{\ensuremath{-}1}\ensuremath{\leqslant}0.2$, where a is the s-wave scattering length and 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, 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 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.
Quantum degenerate Fermi gases can be created in the laboratories using alkali atoms. These gases can be in different regimes of density and interaction strength and provide an ideal test bed for the basic properties of the quantum and statistical mechanics. Also some astrophysical objects such as neutron stars can be idealized as interacting fermion systems. In this thesis work, the ground state properties of dilute Fermi gases with attractive short range two-body interactions are reported. Main results of this work are produced by application of quantum Monte Carlo methods. We introduce variational nodal structure to the trial wave function in order to deal with the fermion sign problem. We consider the possibility of Cooper pairing between fermions of different species. The strength of the interaction is varied to study different pairing regimes which are characterized by the product of the s-wave scattering length a and the Fermi wave vector kF. We also review some of the conventional variational methods such as the Lowest Order Constrained Variational which is suitable for Fermi fluids without pairing correlation. The ground state energy, pairing gap Δ, quasi-particle spectrum, two-body distribution function, and momentum distribution of the superfluid were studied with greater accuracy using ab initio quantum Monte Carlo techniques. In the case of the two component Fermi fluid, we find that in the weak coupling regime, 1/akF < -1, the system is a Bardeen-Cooper-Schrieffer type superfluid and the energy gap Δ is much smaller than the Fermi gas energy EFG. For 1/akF ≥ 0:5 we find that weakly interacting composite bosons or molecules are formed in the superfluid gas with Δ and gas energy per particle approaching half of the molecular energy jEmolj=2. In this region we seem to have Bose-Einstein condensation of molecules. The behavior of the energy, chemical potential and energy gap in the transition region, -0:5 < 1/akF < 0:5, is discussed. The study is extended to the three component Fermi fluid in relation to the spin polarized alkali gases and color superconductivity of the quark matter. Unlike the two component fluids, zero potential range limit cannot be reached without causing collapse of the system. We also study the scaling behavior of the three-body system. We calculate the threshold of the three-body bound state as a function of the potential range and strength. For the many-body systems, we find that ground state can be achieved by broken symmetry pairing. However, the restored symmetry pairing scenario in the strongly interacting regime is not conclusive. We also discuss the pairing gap as well as the structural properties such as two-body distribution functions.
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.
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 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∗, 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.