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

Virial expansion of attractively interacting Fermi gases in one, two, and three dimensions, up to fifth order

Y. Hou, Joaquín E. Drut

Published Sep 11, 2020• 12 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.102.033319

Related papers

Fourth- and Fifth-Order Virial Coefficients from Weak Coupling to Unitarity

Y. Hou, Joaquín E. Drut

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.

arXiv:2004.08685•Physical Review Letters•Jul 28, 2020•18 citations

Stochastic and Semi-Classical Approaches to the Quantum Virial Expansion

C. R. Shill

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.

•Jan 9, 2020•0 citations

Virial coefficients of trapped and untrapped three-component fermions with three-body forces in arbitrary spatial dimensions

A. J. Czejdo, Joaquín E. Drut, Y. Hou +2 more

Using a coarse temporal lattice approximation, we calculate the first few terms of the virial expansion of a three-species fermion system with a three-body contact interaction in $d$ spatial dimensions, both in homogeneous space and in a harmonic trapping potential of frequency $\ensuremath{\omega}$. Using the three-body problem to renormalize, we report analytic results for the change in the fourth- and fifth-order virial coefficients $\mathrm{\ensuremath{\Delta}}{b}_{4}$ and $\mathrm{\ensuremath{\Delta}}{b}_{5}$ as functions of $\mathrm{\ensuremath{\Delta}}{b}_{3}$. Additionally, we argue that in the $\ensuremath{\omega}\ensuremath{\rightarrow}0$ limit the relationship ${b}_{n}^{\text{T}}={n}^{\ensuremath{-}d/2}{b}_{n}$ holds between the trapped (T) and the homogeneous coefficients for arbitrary temperature and coupling strength (not merely in scale-invariant regimes). Finally, we point out an exact, universal (coupling- and frequency-independent) relationship between $\mathrm{\ensuremath{\Delta}}{b}_{3}^{\text{T}}$ in one dimension with three-body forces and $\mathrm{\ensuremath{\Delta}}{b}_{2}^{\text{T}}$ in two dimensions with two-body forces.

arXiv:2003.04958•Physical Review A•Jun 29, 2020•13 citations

Leading- and next-to-leading-order semiclassical approximation to the first seven virial coefficients of spin-1/2 fermions across spatial dimensions

Y. Hou, A. J. Czejdo, Jared De Chant +2 more

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.

arXiv:1907.10120•Physical Review A•Dec 12, 2019•9 citations

High-temperature thermodynamics of a strongly correlated Fermi gas

Xia Liu

We briefly review current research progress in the field of strongly correlated atomic Fermi gases. A deep understanding of strongly correlated quantum gases is plagued by the lack of either exact solutions or a controllable small interaction parameter. In this work, we present a systematic virial expansion method to study the high-temperature thermodynamics of a strongly correlated Fermi gas,where the controllable parameter is given by the fugacity,exp (μkBT),with μ being the chemical potential. A practical way is proposed to obtain the expansion coefficients for both homogeneous and harmonically trapped Fermi gases. We calculate, for the first time accurately, the third-order virial coefficient. The resulting equation of states at high temperatures is compared with a recent thermodynamic measurement and with Monte Carlo simulations.

•Jan 1, 2010•0 citations

Few-Body Correlations in Many-Body Physics

Marcus Barth

In this thesis, cold gases as well as the electron gas are analyzed with two expansion techniques that have, in essence, the same aim: To describe parameter regimes where the few-body aspect dominates over the collective many-body behavior.The Operator Product Expansion (OPE) of quantum field theory is used to extract high-frequency properties of the correlation function that determines the radio-frequency (RF) response of a two-dimensional Fermi gas with zero-range interactions.Exact sum rules that relate the zeroth and first moment of the radio-frequency spectrum to the particle number and the contact of the gas, are derived from these properties.A quantum cluster (virial-) expansion is employed to study the properties of the spectral function of the two-dimensional Fermi gas at small to moderate degeneracy.From the spectral function, observables like the density of states and the momentum distribution are calculated and are found to be in excellent agreement with Tan's universal relations.In the second part of the thesis, the cluster expansion is applied to both bosonic and fermionic three-dimensional quantum gases.Virial coefficients, contacts and momentum distribution are calculated up to third order in the fugacity.This includes three-particle correlations exactly.The results are are in excellent agreement with known universal relations.In the third part of this thesis, an OPE is applied to the problem of the interacting electron gas in both two and three dimensions.A power-law decay of the momentum distribution and the static structure factor at high momenta is found, compatible with earlier findings by J. Kimball.

•Jan 1, 2016•2 citations
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 \ensuremathΔb3,\ensuremathΔb4\mathrm{\ensuremath{\Delta}}{b}_{3}, \mathrm{\ensuremath{\Delta}}{b}_{4}\ensuremathΔb3​,\ensuremathΔb4​, and \ensuremathΔb5\mathrm{\ensuremath{\Delta}}{b}_{5}\ensuremathΔb5​ as functions of \ensuremathΔb2\mathrm{\ensuremath{\Delta}}{b}_{2}\ensuremathΔb2​, 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 \ensuremathΔb3\mathrm{\ensuremath{\Delta}}{b}_{3}\ensuremathΔb3​ with previous calculations in all dimensions and we formulate predictions for \ensuremathΔb4\mathrm{\ensuremath{\Delta}}{b}_{4}\ensuremathΔb4​ and \ensuremathΔb5\mathrm{\ensuremath{\Delta}}{b}_{5}\ensuremathΔb5​ in one and two dimensions. We also provide, for a range of couplings, the subspace contributions \ensuremathΔb31,\ensuremathΔb22,\ensuremathΔb41\mathrm{\ensuremath{\Delta}}{b}_{31}, \mathrm{\ensuremath{\Delta}}{b}_{22}, \mathrm{\ensuremath{\Delta}}{b}_{41}\ensuremathΔb31​,\ensuremathΔb22​,\ensuremathΔb41​, and \ensuremathΔb32\mathrm{\ensuremath{\Delta}}{b}_{32}\ensuremathΔb32​, 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.
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.
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.
Using a coarse temporal lattice approximation, we calculate the first few terms of the virial expansion of a three-species fermion system with a three-body contact interaction in ddd spatial dimensions, both in homogeneous space and in a harmonic trapping potential of frequency $\ensuremath{\omega}$. Using the three-body problem to renormalize, we report analytic results for the change in the fourth- and fifth-order virial coefficients $\mathrm{\ensuremath{\Delta}}{b}_{4}$ and $\mathrm{\ensuremath{\Delta}}{b}_{5}$ as functions of $\mathrm{\ensuremath{\Delta}}{b}_{3}$. Additionally, we argue that in the $\ensuremath{\omega}\ensuremath{\rightarrow}0$ limit the relationship ${b}_{n}^{\text{T}}={n}^{\ensuremath{-}d/2}{b}_{n}$ holds between the trapped (T) and the homogeneous coefficients for arbitrary temperature and coupling strength (not merely in scale-invariant regimes). Finally, we point out an exact, universal (coupling- and frequency-independent) relationship between $\mathrm{\ensuremath{\Delta}}{b}_{3}^{\text{T}}$ in one dimension with three-body forces and $\mathrm{\ensuremath{\Delta}}{b}_{2}^{\text{T}}$ in two dimensions with two-body forces.
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.
We briefly review current research progress in the field of strongly correlated atomic Fermi gases. A deep understanding of strongly correlated quantum gases is plagued by the lack of either exact solutions or a controllable small interaction parameter. In this work, we present a systematic virial expansion method to study the high-temperature thermodynamics of a strongly correlated Fermi gas,where the controllable parameter is given by the fugacity,exp (μkBT),with μ being the chemical potential. A practical way is proposed to obtain the expansion coefficients for both homogeneous and harmonically trapped Fermi gases. We calculate, for the first time accurately, the third-order virial coefficient. The resulting equation of states at high temperatures is compared with a recent thermodynamic measurement and with Monte Carlo simulations.
In this thesis, cold gases as well as the electron gas are analyzed with two expansion techniques that have, in essence, the same aim: To describe parameter regimes where the few-body aspect dominates over the collective many-body behavior.The Operator Product Expansion (OPE) of quantum field theory is used to extract high-frequency properties of the correlation function that determines the radio-frequency (RF) response of a two-dimensional Fermi gas with zero-range interactions.Exact sum rules that relate the zeroth and first moment of the radio-frequency spectrum to the particle number and the contact of the gas, are derived from these properties.A quantum cluster (virial-) expansion is employed to study the properties of the spectral function of the two-dimensional Fermi gas at small to moderate degeneracy.From the spectral function, observables like the density of states and the momentum distribution are calculated and are found to be in excellent agreement with Tan's universal relations.In the second part of the thesis, the cluster expansion is applied to both bosonic and fermionic three-dimensional quantum gases.Virial coefficients, contacts and momentum distribution are calculated up to third order in the fugacity.This includes three-particle correlations exactly.The results are are in excellent agreement with known universal relations.In the third part of this thesis, an OPE is applied to the problem of the interacting electron gas in both two and three dimensions.A power-law decay of the momentum distribution and the static structure factor at high momenta is found, compatible with earlier findings by J. Kimball.