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:1507.08230 · Physical Review A

Numerical study of the unitary Fermi gas across the superfluid transition

Olga Goulko, Matthew Wingate A

Published May 4, 2016• 20 citations
arXivHTMLPDF

Abstract

We present results from Monte Carlo calculations investigating the properties of the homogeneous, spin-balanced unitary Fermi gas in three dimensions. The temperature is varied across the superfluid transition allowing us to determine the temperature dependence of the chemical potential, the energy per particle, and the contact density. Numerical artifacts due to finite volume and discretization are systematically studied, estimated, and reduced.

DOI
10.1103/physreva.93.053604

Related papers

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

S. Jensen, C. N. Gilbreth, Y. Alhassid

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. We also present results for the energy as a function of temperature in the continuum limit.

arXiv:1906.10117•Physical Review Letters•Jul 22, 2020•34 citations

Equation of state of the unitary Fermi gas: An update on lattice calculations

Joaquín E. Drut, Timo A. Lähde, Gabriel Wlazłowski +1 more

The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured to unprecedented accuracy at the MIT. In particular, these measurements provide an improved understanding of the regime below $T/{\ensuremath{\epsilon}}_{F}^{}\ensuremath{\simeq}0.20$, where a transition into a superfluid phase occurs. In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte Carlo (AFQMC) results for the UFG at finite temperature and compare them with the MIT data for the energy, chemical potential, and density. These AFQMC results have been obtained using methods based on the hybrid Monte Carlo (HMC) algorithm, which was first introduced within the context of lattice QCD.

arXiv:1111.5079•Physical Review A•May 1, 2012•41 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

Thermodynamic and hydrodynamic behaviour of interacting Fermi gases

Olga Goulko

Fermionic matter is ubiquitous in nature, from the electrons in metals and semiconductors or the neutrons in the inner crust of neutron stars, to gases of fermionic atoms, like 40K or 6Li that can be created and studied under laboratory conditions. It is especially interesting to study these systems at very low temperatures, where we enter the world of quantum mechanical phenomena. Due to the Fermi-Dirac statistics, a dilute system of spin-polarised fermions exhibits no interactions and can be viewed as an ideal Fermi gas. However, interactions play a crucial role for fermions of several spin species. This thesis addresses several questions concerning interacting Fermi gases, in particular the transition between the normal and the superfluid phase and dynamical properties at higher temperatures. First we will look at the unitary Fermi gas: a two-component system of fermions interacting with divergent scattering length. This system is particularly interesting as it exhibits universal behaviour. Due to the strong interactions perturbation theory is inapplicable and no exact theoretical description is available. I will describe the Determinant Diagrammatic Monte Carlo algorithm with which the unitary Fermi gas can be studied from first principles. This algorithm fails in the presence of a spin imbalance (unequal number of particles in the two components) due to a sign problem. I will show how to apply reweighting techniques to generalise the algorithm to the imbalanced case, and present results for the critical temperature and other thermodynamic observables at the critical point, namely the chemical potential, the energy per particle and the contact density. These are the first numerical results for the imbalanced unitary Fermi gas at finite temperature. I will also show how temperatures beyond the critical point can be accessed and present results for the equation of state and the temperature dependence of the contact density. At sufficiently high temperatures a semiclassical description captures all relevant physical features of the system. The dynamics of an interacting Fermi gas can then be studied via a numerical simulation of the Boltzmann equation. I will describe such a numerical setup and apply it to study the collision of two spin-polarised fermionic clouds. When the two components are separated in an elongated harmonic trap and then released, they collide and for sufficiently strong interactions can bounce off each other several times. I will discuss the different types of the qualitative behaviour, show how they can be interpreted in terms of the equilibrium properties of the system, and explain how they relate to the coupling between different excitation modes. I will also demonstrate how transport coefficients, for instance the spin drag, can be extracted from the numerical data.

•Jan 10, 2012•2 citations

Pairing and superfluid properties of dilute fermion gases at unitarity

Vamsi K. Akkineni, David M. Ceperley, Nandini Trivedi

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.

arXiv:cond-mat/0608154•Physical Review B•Oct 16, 2007•68 citations

Momentum Distribution and Contact of the Unitary Fermi Gas

Joaquín E. Drut, Timo A. Lähde, Timour B. Ten

We calculate the momentum distribution n(k) of the unitary Fermi gas by using quantum Monte Carlo calculations at finite temperature T/ϵ(F) as well as in the ground state. At large momenta k/k(F), we find that n(k) falls off as C/k⁴, in agreement with the Tan relations. From the asymptotics of n(k), we determine the contact C as a function of T/ϵ(F) and present a comparison with theory. At low T/ϵ(F), we find that C increases with temperature, and we tentatively identify a maximum around T/ϵ(F) ≃ 0.4. Our calculations are performed on lattices of spatial extent up to N(x) = 14 with a particle number per unit volume of ≃ 0.03-0.07.

arXiv:1012.5474•Physical Review Letters•May 20, 2011•79 citations
We present results from Monte Carlo calculations investigating the properties of the homogeneous, spin-balanced unitary Fermi gas in three dimensions. The temperature is varied across the superfluid transition allowing us to determine the temperature dependence of the chemical potential, the energy per particle, and the contact density. Numerical artifacts due to finite volume and discretization are systematically studied, estimated, and reduced.
A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. We also present results for the energy as a function of temperature in the continuum limit.
The thermodynamic properties of the unitary Fermi gas (UFG) have recently been measured to unprecedented accuracy at the MIT. In particular, these measurements provide an improved understanding of the regime below $T/{\ensuremath{\epsilon}}_{F}^{}\ensuremath{\simeq}0.20$, where a transition into a superfluid phase occurs. In light of this development, we present an overview of state-of-the-art auxiliary field quantum Monte Carlo (AFQMC) results for the UFG at finite temperature and compare them with the MIT data for the energy, chemical potential, and density. These AFQMC results have been obtained using methods based on the hybrid Monte Carlo (HMC) algorithm, which was first introduced within the context of lattice QCD.
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.
Fermionic matter is ubiquitous in nature, from the electrons in metals and semiconductors or the neutrons in the inner crust of neutron stars, to gases of fermionic atoms, like 40K or 6Li that can be created and studied under laboratory conditions. It is especially interesting to study these systems at very low temperatures, where we enter the world of quantum mechanical phenomena. Due to the Fermi-Dirac statistics, a dilute system of spin-polarised fermions exhibits no interactions and can be viewed as an ideal Fermi gas. However, interactions play a crucial role for fermions of several spin species. This thesis addresses several questions concerning interacting Fermi gases, in particular the transition between the normal and the superfluid phase and dynamical properties at higher temperatures. First we will look at the unitary Fermi gas: a two-component system of fermions interacting with divergent scattering length. This system is particularly interesting as it exhibits universal behaviour. Due to the strong interactions perturbation theory is inapplicable and no exact theoretical description is available. I will describe the Determinant Diagrammatic Monte Carlo algorithm with which the unitary Fermi gas can be studied from first principles. This algorithm fails in the presence of a spin imbalance (unequal number of particles in the two components) due to a sign problem. I will show how to apply reweighting techniques to generalise the algorithm to the imbalanced case, and present results for the critical temperature and other thermodynamic observables at the critical point, namely the chemical potential, the energy per particle and the contact density. These are the first numerical results for the imbalanced unitary Fermi gas at finite temperature. I will also show how temperatures beyond the critical point can be accessed and present results for the equation of state and the temperature dependence of the contact density. At sufficiently high temperatures a semiclassical description captures all relevant physical features of the system. The dynamics of an interacting Fermi gas can then be studied via a numerical simulation of the Boltzmann equation. I will describe such a numerical setup and apply it to study the collision of two spin-polarised fermionic clouds. When the two components are separated in an elongated harmonic trap and then released, they collide and for sufficiently strong interactions can bounce off each other several times. I will discuss the different types of the qualitative behaviour, show how they can be interpreted in terms of the equilibrium properties of the system, and explain how they relate to the coupling between different excitation modes. I will also demonstrate how transport coefficients, for instance the spin drag, can be extracted from the numerical data.
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 calculate the momentum distribution n(k) of the unitary Fermi gas by using quantum Monte Carlo calculations at finite temperature T/ϵ(F) as well as in the ground state. At large momenta k/k(F), we find that n(k) falls off as C/k⁴, in agreement with the Tan relations. From the asymptotics of n(k), we determine the contact C as a function of T/ϵ(F) and present a comparison with theory. At low T/ϵ(F), we find that C increases with temperature, and we tentatively identify a maximum around T/ϵ(F) ≃ 0.4. Our calculations are performed on lattices of spatial extent up to N(x) = 14 with a particle number per unit volume of ≃ 0.03-0.07.