New field-theoretic method for the virial expansion

Toward an Automated-Algebra Framework for High Orders in the Virial Expansion of Quantum Matter

The virial expansion provides a non-perturbative view into the thermodynamics of quantum many-body systems in dilute regimes. While powerful, the expansion is challenging as calculating its coefficients at each order n requires analyzing (if not solving) the quantum n-body problem. In this work, we present a comprehensive review of automated algebra methods, which we developed to calculate high-order virial coefficients. The methods are computational but non-stochastic, thus avoiding statistical effects; they are also for the most part analytic, not numerical, and amenable to massively parallel computer architectures. We show formalism and results for coefficients characterizing the thermodynamics (pressure, density, energy, static susceptibilities) of homogeneous and harmonically trapped systems and explain how to generalize them to other observables such as the momentum distribution, Tan contact, and the structure factor.

High Temperature Virial Expansion to Universal Quench Dynamics

High temperature virial expansion is a powerful tool in equilibrium statistical mechanics. In this Letter we generalize the high temperature virial expansion approach to treat far-from-equilibrium quench dynamics. As an application of our framework, we study the dynamics of a Bose gas quenched from noninteracting to unitarity, and we compare our theoretical results with unexplained experimental results by the Cambridge group [Eigen et al., Nature 563, 221 (2018)]. We show that, during the quench dynamics, the momentum distribution decreases for the low-momentum part with k<k^{*}, and increases for high-momentum part with k>k^{*}, where k^{*} is a characteristic momentum scale separating the low- and the high-momentum regimes. We determine the universal value of k^{*}λ that agrees perfectly with the experiment, with λ being the thermal de Broglie wavelength. We also find a jump of the halfway relaxation time across k^{*}λ and the nonmonotonic behavior of energy distribution, both of which agree with the experiment. Finally, we address the issue whether the longtime steady state thermalizes or not, and we find that this state reaches a partial thermalization, namely, it thermalizes for the low-energy part with kλ≲1 but does not thermalize for the very high momentum tail with kλ≫1. Our framework can also be applied to quench dynamics in other systems.

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

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.

Visualizing the Efimov Correlation in Bose Polarons

The Bose polaron is a quasiparticle of an impurity dressed by surrounding bosons. In few-body physics, it is known that two identical bosons and a third distinguishable particle can form a sequence of Efimov bound states in the vicinity of interspecies scattering resonance. On the other hand, in the Bose polaron system with an impurity atom embedded in many bosons, no signature of Efimov physics has been reported in the existing spectroscopy measurements to date. In this Letter, we propose that a large mass imbalance between a light impurity and heavy bosons can help produce visible signatures of Efimov physics in such a spectroscopy measurement. Using the diagrammatic approach in the virial expansion to include three-body effects from pair-wise interactions, we determine the impurity self-energy and its spectral function. Taking the ^{6}Li-^{133}Cs system as a concrete example, we find two visible Efimov branches in the polaron spectrum, as well as their hybridizations with the attractive polaron branch. We also discuss the general scenarios for observing the signature of Efimov physics in polaron systems. This work paves the way for experimentally exploring intriguing few-body correlations in a many-body system in the near future.

Path-Integral Monte Carlo Determination of the Fourth-Order Virial Coefficient for a Unitary Two-Component Fermi Gas with Zero-Range Interactions

The unitary equal-mass Fermi gas with zero-range interactions constitutes a paradigmatic model system that is relevant to atomic, condensed matter, nuclear, particle, and astrophysics. This work determines the fourth-order virial coefficient b_{4} of such a strongly interacting Fermi gas using a customized ab initio path-integral Monte Carlo (PIMC) algorithm. In contrast to earlier theoretical results, which disagreed on the sign and magnitude of b_{4}, our b_{4} agrees within error bars with the experimentally determined value, thereby resolving an ongoing literature debate. Utilizing a trap regulator, our PIMC approach determines the fourth-order virial coefficient by directly sampling the partition function. An on-the-fly antisymmetrization avoids the Thomas collapse and, combined with the use of the exact two-body zero-range propagator, establishes an efficient general means to treat small Fermi systems with zero-range interactions.

Virial expansion coefficients in the harmonic approximation

The virial expansion method is applied within a harmonic approximation to an interacting N-body system of identical fermions. We compute the canonical partition functions for two and three particles to get the two lowest orders in the expansion. The energy spectrum is carefully interpolated to reproduce ground-state properties at low temperature and the noninteracting high-temperature limit of constant virial coefficients. This resembles the smearing of shell effects in finite systems with increasing temperature. Numerical results are discussed for the second and third virial coefficients as functions of dimension, temperature, interaction, and transition temperature between low- and high-energy limits.