Virial expansion for a strongly correlated Fermi gas

Observation of the Fermionic Joule-Thomson Effect

We report the first observation of the quantum Joule-Thomson (JT) effect in ideal and unitary Fermi gases. We study the temperature dynamics of these systems while they undergo an energy-per-particle conserving rarefaction. For scale-invariant systems, whose equations of state satisfy the relation U∝PV, this rarefaction conserves the specific enthalpy, which makes it thermodynamically equivalent to a JT throttling process. We observe JT heating in an ideal Fermi gas, a direct consequence of Pauli blocking. In a unitary Fermi gas, we observe that the JT heating is marginal in the temperature range 0.2≲T/T_{F}≲0.8 as the repulsive quantum-statistical effect is lessened by the attractive interparticle interactions.

Three-dimensional non-approximate Coulomb interaction between two trapped quantum particles

The two trapped quantum particles interacting problem is generalized to three dimensions, and the exact Coulomb potential is used. The system is solved by expanding the wavefunction in terms of the isotropic harmonic oscillator eigenfunctions and Hydrogen atom eigenfunctions independently, showing that each one results in a prime approximation for different domains of the normalized coupling constant of the relative interactions, suggesting that the combination of the basis is enough to build a well-suited base for the non-approximate problem. The results are compared to previous works that use a model of approximate finite-rage soft-core interaction model of the problem to give insights into the many-body states of strongly correlated ultracold bosons and fermions. We conclude that the proposed three-dimensional approach facilitates the distinction between bosons and fermions while the solutions given by the expansions better define the behavior of the particles for repulsive potentials. In addition, we discuss the substantial differences between our work and the previous approximate model.

A quantum engine in the BEC-BCS crossover

Heat engines convert thermal energy into mechanical work both in the classical and quantum regimes1. However, quantum theory offers genuine non-classical forms of energy, different from heat, which so far have not been exploited in cyclic engines. Here we experimentally realize a quantum many-body engine fuelled by the energy difference between fermionic and bosonic ensembles of ultracold particles that follows from the Pauli exclusion principle2. We employ a harmonically trapped superfluid gas of 6Li atoms close to a magnetic Feshbach resonance3 that allows us to effectively change the quantum statistics from Bose-Einstein to Fermi-Dirac, by tuning the gas between a Bose-Einstein condensate of bosonic molecules and a unitary Fermi gas (and back) through a magnetic field4-10. The quantum nature of such a Pauli engine is revealed by contrasting it with an engine in the classical thermal regime and with a purely interaction-driven device. We obtain a work output of several 106 vibrational quanta per cycle with an efficiency of up to 25%. Our findings establish quantum statistics as a useful thermodynamic resource for work production.

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.

Tunneling Transport of Unitary Fermions across the Superfluid Transition

We investigate the transport of a Fermi gas with unitarity-limited interactions across the superfluid phase transition, probing its response to a direct current (dc) drive through a tunnel junction. As the superfluid critical temperature is crossed from below, we observe the evolution from a highly nonlinear to an Ohmic conduction characteristic, associated with the critical breakdown of the Josephson dc current induced by pair condensate depletion. Moreover, we reveal a large and dominant anomalous contribution to resistive currents, which reaches its maximum at the lowest attained temperature, fostered by the tunnel coupling between the condensate and phononic Bogoliubov-Anderson excitations. Increasing the temperature, while the zeroing of supercurrents marks the transition to the normal phase, the conductance drops considerably but remains much larger than that of a normal, uncorrelated Fermi gas tunneling through the same junction. We attribute such enhanced transport to incoherent tunneling of sound modes, which remain weakly damped in the collisional hydrodynamic fluid of unpaired fermions at unitarity.

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.

Few-Body Perspective of a Quantum Anomaly in Two-Dimensional Fermi Gases

A quantum anomaly manifests itself in the deviation of the breathing mode frequency from the scale invariant value of 2ω in two-dimensional harmonically trapped Fermi gases, where ω is the trapping frequency. Its recent experimental observation with cold atoms reveals an unexpected role played by the effective range of interactions, which requires a quantitative theoretical understanding. Here we provide accurate, benchmark results on a quantum anomaly from a few-body perspective. We consider the breathing mode of a few trapped interacting fermions in two dimensions up to six particles and present the mode frequency as a function of scattering length for a wide range of effective range. We show that the maximum quantum anomaly gradually reduces as the effective range increases while the maximum position shifts towards the weak-coupling limit. We extrapolate our few-body results to the many-body limit and find a good agreement with the experimental measurements. Our results may also be directly applicable to a few-fermion system prepared in microtraps and optical tweezers.

Resummation of Diagrammatic Series with Zero Convergence Radius for Strongly Correlated Fermions

We demonstrate that a summing up series of Feynman diagrams can yield unbiased accurate results for strongly correlated fermions even when the convergence radius vanishes. We consider the unitary Fermi gas, a model of nonrelativistic fermions in three-dimensional continuous space. Diagrams are built from partially dressed or fully dressed propagators of single particles and pairs. The series is resummed by a conformal-Borel transformation that incorporates the large-order behavior and the analytic structure in the Borel plane, which are found by the instanton approach. We report highly accurate numerical results for the equation of state in the normal unpolarized regime, and reconcile experimental data with the theoretically conjectured fourth virial coefficient.

Review of pseudogaps in strongly interacting Fermi gases

A central challenge in modern condensed matter physics is developing the tools for understanding nontrivial yet unordered states of matter. One important idea to emerge in this context is that of a 'pseudogap': the fact that under appropriate circumstances the normal state displays a suppression of the single particle spectral density near the Fermi level, reminiscent of the gaps seen in ordered states of matter. While these concepts arose in a solid state context, they are now being explored in cold gases. This article reviews the current experimental and theoretical understanding of the normal state of strongly interacting Fermi gases, with particular focus on the phenomonology which is traditionally associated with the pseudogap.

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.

Collective modes in a unitary Fermi gas across the superfluid phase transition

We provide a joint theoretical and experimental investigation of the temperature dependence of the collective oscillations of first sound nature exhibited by a highly elongated harmonically trapped Fermi gas at unitarity, including the region below the critical temperature for superfluidity. Differently from the lowest axial breathing mode, the hydrodynamic frequencies of the higher-nodal excitations show a temperature dependence, which is calculated starting from Landau two-fluid theory and using the available experimental knowledge of the equation of state. The experimental results agree with high accuracy with the predictions of theory and provide the first evidence for the temperature dependence of the collective frequencies near the superfluid phase transition.

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.

Universal equation of state of a unitary fermionic gas

It is suggested that for a Fermi gas at unitarity, the two-body bond plays a special role. We propose an equation of state using an ansatz relating the interaction part of the l-body cluster to its two-body counterpart. This allows a parameter-free comparison with the recently measured equation of state by the ENS group. The agreement between the two over a range of fugacity (z<5 for a homogeneous gas, and z<10 for the trapped gas) leads us to perform the calculations of more sensitive quantities measured recently by the MIT group.

Few-body physics with ultracold atomic and molecular systems in traps

Few-body physics has played a prominent role in atomic, molecular and nuclear physics since the early days of quantum mechanics. It is now possible-thanks to tremendous progress in cooling, trapping and manipulating ultracold samples-to experimentally study few-body phenomena in trapped atomic and molecular systems with unprecedented control. This review summarizes recent studies of few-body phenomena in trapped atomic and molecular gases, with an emphasis on small trapped systems. We start by introducing the free-space scattering properties and then investigate what happens when two particles, bosons or fermions, are placed in an external confinement. Next, various three-body systems are treated analytically in limiting cases. Our current understanding of larger two-component Fermi systems and Bose systems is reviewed, and connections with the corresponding bulk systems are established. Lastly, future prospects and challenges are discussed. Throughout this review, commonalities with other systems such as nuclei or quantum dots are highlighted.

Density and spin response of a strongly interacting Fermi gas in the attractive and quasirepulsive regime

Recent experimental advances in ultracold Fermi gases allow for exploring response functions under different dynamical conditions. In particular, the issue of obtaining a "quasirepulsive" regime starting from a Fermi gas with an attractive interparticle interaction while avoiding the formation of the two-body bound state is currently debated. Here, we provide a calculation of the density and spin response for a wide range of temperature and coupling both in the attractive and quasirepulsive regime, whereby the system is assumed to evolve nonadiabatically toward the "upper branch" of the Fermi gas. A comparison is made with the available experimental data for these two quantities.

Quantum statistics and thermodynamics in the harmonic approximation

We describe a method to compute thermodynamic quantities in the harmonic approximation for identical bosons and fermions in an external confining field. We use the canonical partition function where only energies and their degeneracies enter. The number of states of given energy and symmetry is found by separating the center-of-mass motion, and by counting the remaining states of given symmetry and excitation energy of the relative motion. The oscillator frequencies that enter the harmonic Hamiltonian can be derived from realistic model parameters, and the method corresponds to an effective interaction approach based on harmonic interactions. To demonstrate the method, we apply it to systems in two dimensions. Numerical calculations are compared to a brute force method, which is considerably more computationally intensive.

Few-particle Green's functions for strongly correlated systems on infinite lattices

We show how few-particle Green's functions can be calculated efficiently for models with nearest-neighbor hopping, for infinite lattices in any dimension. As an example, for one-dimensional spinless fermions with both nearest-neighbor and second-nearest-neighbor interactions, we investigate the ground states for up to 5 fermions. This allows us not only to find the stability region of various bound complexes, but also to infer the phase diagram at small but finite concentrations.

New field-theoretic method for the virial expansion

We develop a graphical method for computing the virial expansion coefficients for a nonrelativistic quantum field theory. As an example we compute the third virial coefficient b3 for unitary fermions, a nonperturbative system. By calculating several graphs and performing an extrapolation, we arrive at b3=-0.2930, within 0.7% of a recent computation b3=-0.290 952 95 by Liu, Hu, and Drummond, which involved summing 10,000 energy levels for three unitary fermions in a harmonic trap.

Universal behavior of pair correlations in a strongly interacting Fermi gas

We show that short-range pair correlations in a strongly interacting Fermi gas follow a simple universal law described by Tan's relations. This is achieved through measurements of the static structure factor which displays a universal scaling proportional to the ratio of Tan's contact to the momentum C/q. Bragg spectroscopy of ultracold 6Li atoms from a periodic optical potential is used to measure the structure factor for a wide range of momenta and interaction strengths, providing broad confirmation of this universal law. We calibrate our Bragg spectra using the f-sum rule, which is found to improve the accuracy of the structure factor measurement.

Pseudogap pairing in ultracold Fermi atoms

The Bose-Einstein condensate to Bardeen-Cooper-Schrieffer crossover in ultracold Fermi gases creates an ideal environment to enrich our knowledge of many-body systems. It is relevant to a wide range of fields from condensed matter to astrophysics. The nature of pairing in strongly interacting Fermi gases can be readily studied. This aids our understanding of related problems in high-Tc superconductors, whose mechanism is still under debate due to the large interaction parameter. Here, we calculate the dynamical properties of a normal, trapped strongly correlated Fermi gas, by developing a quantum cluster expansion. Our calculations for the single-particle spectral function agree with recent rf spectroscopy measurements, and clearly demonstrate pseudogap pairing in the strongly interacting regime.