In the current era of precision quantum many-body physics, one of the most scrutinized systems is the unitary limit of the nonrelativistic spin-1/2 Fermi gas, due to its simplicity and relevance for atomic, condensed matter, and nuclear physics. The thermodynamics of this strongly correlated system is determined by universal functions which, at high temperatures, are governed by universal virial coefficients b_{n} that capture the effects of the n-body system on the many-body dynamics. Currently, b_{2} and b_{3} are well understood, but the situation is less clear for b_{4}, and no predictions have been made for b_{5}. To answer these open questions, we implement a nonperturbative analytic approach based on the Trotter-Suzuki factorization of the imaginary-time evolution operator, using progressively finer temporal lattice spacings. By means of these factorizations and automated algebra codes, we obtain the interaction-induced change Δb_{n} from weak coupling to unitarity. At unitarity, we find that Δb_{3}=-0.356(4) in agreement with previous results, Δb_{4}=0.062(2), which is in agreement with all previous theoretical estimates but at odds with experimental determinations, and Δb_{5}=0.078(6), which is a prediction. We show the impact of those answers on the density equation of state and Tan contact, and trace their origin back to their polarized and unpolarized components.
In the current era of precision quantum many-body physics, one of the most scrutinized systems is the unitary limit of the nonrelativistic spin-1/2 Fermi gas, due to its simplicity and relevance for atomic, condensed matter, and nuclear physics. The thermodynamics of this strongly correlated system is determined by universal functions which, at high temperatures, are governed by universal virial coefficients b_{n} that capture the effects of the n-body system on the many-body dynamics. Currently, b_{2} and b_{3} are well understood, but the situation is less clear for b_{4}, and no predictions have been made for b_{5}. To answer these open questions, we implement a nonperturbative analytic approach based on the Trotter-Suzuki factorization of the imaginary-time evolution operator, using progressively finer temporal lattice spacings. By means of these factorizations and automated algebra codes, we obtain the interaction-induced change Δb_{n} from weak coupling to unitarity. At unitarity, we find that Δb_{3}=-0.356(4) in agreement with previous results, Δb_{4}=0.062(2), which is in agreement with all previous theoretical estimates but at odds with experimental determinations, and Δb_{5}=0.078(6), which is a prediction. We show the impact of those answers on the density equation of state and Tan contact, and trace their origin back to their polarized and unpolarized components.
The virial expansion characterizes the high-temperature approach to the quantum-classical crossover in any quantum many-body system. Here, we calculate the virial coefficients up to the fifth order of Fermi gases in one, two, and three dimensions, with attractive contact interactions, as relevant for a variety of applications in atomic and nuclear physics. To that end, we discretize the imaginary-time direction and calculate the relevant canonical partition functions. In coarse discretizations, we obtain analytic results featuring relationships between the interaction-induced changes $\mathrm{\ensuremath{\Delta}}{b}_{3}, \mathrm{\ensuremath{\Delta}}{b}_{4}$, and $\mathrm{\ensuremath{\Delta}}{b}_{5}$ as functions of $\mathrm{\ensuremath{\Delta}}{b}_{2}$, the latter being exactly known in many cases by virtue of the Beth-Uhlenbeck formula. Using automated-algebra methods, we push our calculations to progressively finer discretizations and extrapolate to the continuous-time limit. We find excellent agreement for $\mathrm{\ensuremath{\Delta}}{b}_{3}$ with previous calculations in all dimensions and we formulate predictions for $\mathrm{\ensuremath{\Delta}}{b}_{4}$ and $\mathrm{\ensuremath{\Delta}}{b}_{5}$ in one and two dimensions. We also provide, for a range of couplings, the subspace contributions $\mathrm{\ensuremath{\Delta}}{b}_{31}, \mathrm{\ensuremath{\Delta}}{b}_{22}, \mathrm{\ensuremath{\Delta}}{b}_{41}$, and $\mathrm{\ensuremath{\Delta}}{b}_{32}$, which determine the equation of state and static response of polarized systems at high temperature. As a performance check, we compare the density equation of state and Tan contact with quantum Monte Carlo calculations, diagrammatic approaches, and experimental data where available. Finally, we apply Pad\'e and Pad\'e-Borel resummation methods to extend the usefulness of the virial coefficients to approach and in some cases go beyond the unit-fugacity point.
By generalizing our automated algebra approach from homogeneous space to harmonically trapped systems, we have calculated the fourth- and fifth-order virial coefficients of universal spin-21 fermions in the unitary limit, confined in an isotropic harmonic potential. We present results for said coefficients as a function of trapping frequency (or, equivalently, temperature), which compare favorably with previous Monte Carlo calculations (available only at fourth order) as well as with our previous estimates in the untrapped limit (high temperature, low frequency). In addition to the conventional coefficients, we provide estimates for the contributions from subspaces with varying polarization. We use our results for the virial coefficients, together with resummation techniques, to calculate the compressibility and spin susceptibility.
Universality of strongly interacting fermions is a topic of great interest in diverse fields. Here we investigate theoretically the universal dynamic density response of resonantly interacting ultracold Fermi atoms in the limit of either high temperature or large frequency. At high temperature, we use quantum virial expansion to derive universal, nonperturbative expansion functions of the dynamic structure factor; at large momentum, we conjecture that the second-order expansion function gives the Wilson coefficient used in the operator product expansion method. The dynamic structure factor is therefore determined by its second-order expansion function with an overall normalization factor given by Tan's [S. Tan, Ann. Phys. (NY) 323, 2952 (2008); S. Tan, Ann. Phys. (NY) 323, 2971 (2008); S. Tan, Ann. Phys. (NY) 323, 2987 (2008)] contact parameter. We show that the spin-parallel and -antiparallel dynamic structure factors have, respectively, a tail of the form $\ensuremath{\sim}\ifmmode\pm\else\textpm\fi{}{\ensuremath{\omega}}^{\ensuremath{-}5/2}$ for $\ensuremath{\omega}\ensuremath{\rightarrow}\ensuremath{\infty}$, decaying slower than the total dynamic structure factor found previously ($\ensuremath{\sim}{\ensuremath{\omega}}^{\ensuremath{-}7/2}$). Our predictions for the dynamic structure factor at high temperature or large frequency are testable using Bragg spectroscopy for ultracold atomic Fermi gases.
By generalizing our automated algebra approach from homogeneous space to harmonically trapped systems, we have calculated the fourth- and fifth-order virial coefficients of universal spin- fermions in the unitary limit, confined in an isotropic harmonic potential. We present results for said coefficients as a function of trapping frequency (or, equivalently, temperature), which compare favorably with previous Monte Carlo calculations (available only at fourth order) as well as with our previous estimates in the untrapped limit (high temperature, low frequency). In addition to the conventional coefficients, we provide estimates for the contributions from subspaces with varying polarization. We use our results for the virial coefficients, together with resummation techniques, to calculate the compressibility and spin susceptibility.
The Unitary Fermi Gas (UFG) is one of the most strongly interacting systems known to date, as it saturates the unitarity bound on the quantum mechanical scattering cross section. The UFG corresponds to a two-component Fermi gas in the limit of short interaction range and large scattering length, and is currently realized in ultracold-atom experiments via Feshbach resonances. While easy to define, the UFG poses a challenging quantum many-body problem, as it lacks any characteristic scale other than the density. As a consequence, accurate quantitative predictions of the thermodynamic properties of the UFG require Monte Carlo calculations. However, significant progress has also been made with purely analytical methods. Notably, in 2005 Tan derived a set of exact thermodynamic relations in which a universal quantity known as the "contact" C plays a crucial role. Recently, C has also been found to determine the prefactor of the high- frequency power-law decay of correlators as well as the right-hand-sides of shear- and bulk viscosity sum rules. The contact is therefore a central piece of information on the UFG in equilibrium as well as away from equilibrium. In this talk we describe some of the known aspects of Fermi gases at and around unitarity, show our latest Monte Carlo results for the contact at finite temperature, and summarize the open questions in the field, some of which we are starting to answer using large-scale Monte Carlo calculations by adapting methods from Lattice QCD.
The Unitary Fermi Gas (UFG) is one of the most strongly interacting systems known to date, as it saturates the unitarity bound on the quantum mechanical scattering cross section.The UFG corresponds to a two-component Fermi gas in the limit of short interaction range and large scattering length, and is currently realized in ultracold-atom experiments via Feshbach resonances.While easy to define, the UFG poses a challenging quantum many-body problem, as it lacks any characteristic scale other than the density.As a consequence, accurate quantitative predictions of the thermodynamic properties of the UFG require Monte Carlo calculations.However, significant progress has also been made with purely analytical methods.Notably, in 2005 Tan derived a set of exact thermodynamic relations in which a universal quantity known as the "contact" C plays a crucial role.Recently, C has also been found to determine the prefactor of the high-frequency power-law decay of correlators as well as the right-hand-sides of shear-and bulk viscosity sum rules.The contact is therefore a central piece of information on the UFG in equilibrium as well as away from equilibrium.In this talk we describe some of the known aspects of Fermi gases at and around unitarity, show our latest Monte Carlo results for the contact at finite temperature, and summarize the open questions in the field, some of which we are starting to answer using large-scale Monte Carlo calculations by adapting methods from Lattice QCD.