We map out the interaction effects on the first six virial coefficients of one-dimensional Fermi gases with zero-range attractive and repulsive interactions and the first four virial coefficients of the two-dimensional analog with attractive interactions. To that end, we use two nonperturbative stochastic methods: projection by complex stochastic quantization, which allows us to determine high-order coefficients at weak coupling and estimate the radius of convergence of the virial expansion, and a path-integral representation of the virial coefficients. To complement our numerical calculations, we present leading-order results in a semiclassical lattice approximation, which we find to be surprisingly close to the expected answers.
We map out the interaction effects on the first six virial coefficients of one-dimensional Fermi gases with zero-range attractive and repulsive interactions and the first four virial coefficients of the two-dimensional analog with attractive interactions. To that end, we use two nonperturbative stochastic methods: projection by complex stochastic quantization, which allows us to determine high-order coefficients at weak coupling and estimate the radius of convergence of the virial expansion, and a path-integral representation of the virial coefficients. To complement our numerical calculations, we present leading-order results in a semiclassical lattice approximation, which we find to be surprisingly close to the expected answers.
Many-body quantum systems provide an interesting playground for experimentalists and theorists alike. In particular, ultracold atomic gases provide not only a realm of study for which experiments have a high degree of control but also provides theorists with intriguing, yet challenging systems that may be accessed computationally. It is for these reasons that it is here where the majority of our work will focus. We discuss the necessary physics, thermodynamics, and computational methods for carrying out this work, and present several methods for computing coefficients of the quantum virial expansion across dimensions for Fermi gases with zero-range interactions. The techniques discussed are non-perturbative and involve stochastic methods as well as a semi-classical lattice approximation. Through these algorithms, we have computed the interaction dependence (both attractive and repulsive where applicable) of the virial coefficients to a relatively high order as compared to current literature in one, two, and three spatial dimensions. Through one of said methods, we are able to provide a prediction for the radius of convergence of the virial expansion in one spatial dimension.
Following up on recent calculations, we investigate the leading- and next-to-leading-order semiclassical approximation to the virial coefficients of a two-species fermion system with a contact interaction. Using the analytic result for the second-order virial coefficient as a renormalization condition, we derive expressions for up to the seventh-order virial coefficient $\mathrm{\ensuremath{\Delta}}{b}_{7}$. Our results at leading order, though approximate, furnish simple analytic formulas that relate $\mathrm{\ensuremath{\Delta}}{b}_{n}$ to $\mathrm{\ensuremath{\Delta}}{b}_{2}$ for arbitrary dimension, providing a glimpse into the behavior of the virial expansion across dimensions and coupling strengths. As an application, we calculate the pressure and Tan's contact of the two-dimensional attractive Fermi gas and examine the radius of convergence of the virial expansion as a function of the coupling strength.
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.
We use the virial expansion to investigate the behavior of the two-component, attractive Fermi gas in the high-temperature limit, where the system smoothly evolves from weakly attractive fermions to weakly repulsive bosonic dimers as the short-range attraction is increased. We present a formalism for computing the virial coefficients that employs a diagrammatic approach to the grand potential and allows one to easily include an effective range R∗ in the interaction. In the limit where the thermal wavelength $\ensuremath{\lambda}\ensuremath{\ll}{R}^{*}$, the calculation of the virial coefficients is perturbative even at unitarity and the system resembles a weakly interacting Bose-Fermi mixture for all scattering lengths a. By interpolating from the perturbative limits $\ensuremath{\lambda}/|a|\ensuremath{\gg}1$ and ${R}^{*}/\ensuremath{\lambda}\ensuremath{\gg}1$, we estimate the value of the fourth virial coefficient at unitarity for R∗=0 and we find that it is close to the value obtained in recent experiments. We also derive the equations of state for the pressure, density, and entropy, as well as the spectral function at high temperatures.
We study the virial expansion for three-dimensional Bose and Fermi gases at finite temperature using an approximation that only considers two-body processes and is valid for high temperatures and low densities. The first virial coefficients are computed and the second is exact. The results are obtained for the full range of values of the scattering length, and the unitary limit is recovered as a particular case. A weak coupling expansion is performed and the free case is also obtained as a proper limit. The influence of an anisotropic harmonic trap is considered using the local density approximation (LDA), analytical results are obtained, and the special case of the isotropic trap is discussed in detail.
Using a leading-order semiclassical approximation, we calculate the third- and fourth-order virial coefficients of nonrelativistic spin-1/2 fermions in a harmonic trapping potential in arbitrary spatial dimensions, and as functions of temperature, trapping frequency, and coupling strength. Our simple, analytic results for the interaction-induced changes $\mathrm{\ensuremath{\Delta}}{b}_{3}$ and $\mathrm{\ensuremath{\Delta}}{b}_{4}$ agree qualitatively, and in some regimes quantitatively, with previous numerical calculations for the unitary limit of three-dimensional Fermi gases.