Universal behavior of pair correlations in a strongly interacting Fermi gas

Precise Photoexcitation Measurement of Tan's Contact in the Entire BCS-BEC Crossover

We study two-body correlations in a spin-balanced ultracold harmonically trapped Fermi gas of ^{6}Li atoms in the crossover from the Bardeen-Cooper-Schrieffer (BCS) to the Bose-Einstein-Condensate (BEC) regime. For this, we precisely measure Tan's contact using a novel method based on photoexcitation of atomic pairs, which was recently proposed by Wang et al. [Photoexcitation measurement of Tan's contact for a strongly interacting Fermi gas, Phys. Rev. A 104, 063309 (2021).PLRAAN2469-992610.1103/PhysRevA.104.063309]. We map out the contact in the entire phase diagram of the BCS-BEC crossover for various temperatures and interaction strengths, probing regions in phase space that have not been investigated yet. Our measurements reach an uncertainty of ≈2% at the lowest temperatures and thus represent a precise quantitative benchmark. By comparison to our data, we localize the regions in phase space where theoretical predictions and interpolations give valid results. In regions where the contact is already well known we find excellent agreement with our measurements. Thus, our results demonstrate that photoinduced loss is a precise probe to measure quantum correlations in a strongly interacting Fermi gas.

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. The ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

Stable recursive auxiliary field quantum Monte Carlo algorithm in the canonical ensemble: Applications to thermometry and the Hubbard model

Many experimentally accessible, finite-sized interacting quantum systems are most appropriately described by the canonical ensemble of statistical mechanics. Conventional numerical simulation methods either approximate them as being coupled to a particle bath or use projective algorithms which may suffer from nonoptimal scaling with system size or large algorithmic prefactors. In this paper, we introduce a highly stable, recursive auxiliary field quantum Monte Carlo approach that can directly simulate systems in the canonical ensemble. We apply the method to the fermion Hubbard model in one and two spatial dimensions in a regime known to exhibit a significant "sign" problem and find improved performance over existing approaches including rapid convergence to ground-state expectation values. The effects of excitations above the ground state are quantified using an estimator-agnostic approach including studying the temperature dependence of the purity and overlap fidelity of the canonical and grand canonical density matrices. As an important application, we show that thermometry approaches often exploited in ultracold atoms that employ an analysis of the velocity distribution in the grand canonical ensemble may be subject to errors leading to an underestimation of extracted temperatures with respect to the Fermi temperature.

Universal pair polaritons in a strongly interacting Fermi gas

Cavity quantum electrodynamics (QED) manipulates the coupling of light with matter, and allows several emitters to couple coherently with one light mode¹. However, even in a many-body system, the light-matter coupling mechanism has so far been restricted to one-body processes. Leveraging cavity QED for the quantum simulation of complex, many-body systems has thus far relied on multi-photon processes, scaling down the light-matter interaction to the low energy and slow time scales of the many-body problem^2-5. Here we report cavity QED experiments using molecular transitions in a strongly interacting Fermi gas, directly coupling cavity photons to pairs of atoms. The interplay of strong light-matter and strong interparticle interactions leads to well-resolved pair polaritons-hybrid excitations coherently mixing photons, atom pairs and molecules. The dependence of the pair-polariton spectrum on interatomic interactions is universal, independent of the transition used, demonstrating a direct mapping between pair correlations in the ground state and the optical spectrum. This represents a magnification of many-body effects by two orders of magnitude in energy. In the dispersive regime, it enables fast, minimally destructive measurements of pair correlations, and opens the way to their measurement at the quantum limit and their coherent manipulation using dynamical, quantized optical fields.

Breakdown of Tan's Relation in Lossy One-Dimensional Bose Gases

In quantum gases with contact repulsion, the distribution of momenta of the atoms typically decays as ∼1/|p|^{4} at large momentum p. Tan's relation connects the amplitude of that 1/|p|^{4} tail to the adiabatic derivative of the energy with respect to the coupling constant or scattering length of the gas. Here it is shown that the relation breaks down in the one-dimensional Bose gas with contact repulsion, for a peculiar class of stationary states. These states exist thanks to the infinite number of conserved quantities in the system, and they are characterized by a rapidity distribution that itself decreases as 1/|p|^{4}. In the momentum distribution, that rapidity tail adds to the usual Tan contact term. Remarkably, atom losses, which are ubiquitous in experiments, do produce such peculiar states. The development of the tail of the rapidity distribution originates from the ghost singularity of the wave function immediately after each loss event. This phenomenon is discussed for arbitrary interaction strengths, and it is supported by exact calculations in the two asymptotic regimes of infinite and weak repulsion.

Structure Factors of Neutron Matter at Finite Temperature

We compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature.

Heuristic machinery for thermodynamic studies of SU(N) fermions with neural networks

The power of machine learning (ML) provides the possibility of analyzing experimental measurements with a high sensitivity. However, it still remains challenging to probe the subtle effects directly related to physical observables and to understand physics behind from ordinary experimental data using ML. Here, we introduce a heuristic machinery by using machine learning analysis. We use our machinery to guide the thermodynamic studies in the density profile of ultracold fermions interacting within SU(N) spin symmetry prepared in a quantum simulator. Although such spin symmetry should manifest itself in a many-body wavefunction, it is elusive how the momentum distribution of fermions, the most ordinary measurement, reveals the effect of spin symmetry. Using a fully trained convolutional neural network (NN) with a remarkably high accuracy of ~94% for detection of the spin multiplicity, we investigate how the accuracy depends on various less-pronounced effects with filtered experimental images. Guided by our machinery, we directly measure a thermodynamic compressibility from density fluctuations within the single image. Our machine learning framework shows a potential to validate theoretical descriptions of SU(N) Fermi liquids, and to identify less-pronounced effects even for highly complex quantum matter with minimal prior understanding.

Tan's two-body contact across the superfluid transition of a planar Bose gas

Tan’s contact is a quantity that unifies many different properties of a low-temperature gas with short-range interactions, from its momentum distribution to its spatial two-body correlation function. Here, we use a Ramsey interferometric method to realize experimentally the thermodynamic definition of the two-body contact, i.e., the change of the internal energy in a small modification of the scattering length. Our measurements are performed on a uniform two-dimensional Bose gas of ⁸⁷Rb atoms across the Berezinskii–Kosterlitz–Thouless superfluid transition. They connect well to the theoretical predictions in the limiting cases of a strongly degenerate fluid and of a normal gas. They also provide the variation of this key quantity in the critical region, where further theoretical efforts are needed to account for our findings.

Here the authors use Ramsey interferometry to study Tan’s contact in uniform two-dimensional Bose gas of 87Rb atoms across the Berezinskii–Kosterlitz–Thouless superfluid transition. They find that the two-body contact is continuous across the critical point.

Contact in the Unitary Fermi Gas across the Superfluid Phase Transition

A quantity known as the contact is a fundamental thermodynamic property of quantum many-body systems with short-range interactions. Determination of the temperature dependence of the contact for the unitary Fermi gas of infinite scattering length has been a major challenge, with different calculations yielding qualitatively different results. Here we use finite-temperature auxiliary-field quantum Monte Carlo (AFMC) methods on the lattice within the canonical ensemble to calculate the temperature dependence of the contact for the homogeneous spin-balanced unitary Fermi gas. We extrapolate to the continuum limit for 40, 66, and 114 particles, eliminating systematic errors due to finite-range effects. We observe a dramatic decrease in the contact as the superfluid critical temperature is approached from below, followed by a gradual weak decrease as the temperature increases in the normal phase. Our theoretical results are in excellent agreement with the most recent precision ultracold atomic gas experiments. We also present results for the energy as a function of temperature in the continuum limit.

High-Frequency Sound in a Unitary Fermi Gas

We present an experimental and theoretical study of the phonon mode in a unitary Fermi gas. Using two-photon Bragg spectroscopy, we measure excitation spectra at a momentum of approximately half the Fermi momentum, both above and below the superfluid critical temperature T_{c}. Below T_{c}, the dominant excitation is the Bogoliubov-Anderson (BA) phonon mode, driven by gradients in the phase of the superfluid order parameter. The temperature dependence of the BA phonon is consistent with a theoretical model based on the quasiparticle random phase approximation in which the dominant damping mechanism is via collisions with thermally excited quasiparticles. As the temperature is increased above T_{c}, the phonon evolves into a strongly damped collisional mode, accompanied by an abrupt increase in spectral width. Our study reveals strong similarities between sound propagation in the unitary Fermi gas and bosonic liquid helium.

Spectral Response and Contact of the Unitary Fermi Gas

We measure radio frequency (rf) spectra of the homogeneous unitary Fermi gas at temperatures ranging from the Boltzmann regime through quantum degeneracy and across the superfluid transition. For all temperatures, a single spectral peak is observed. Its position smoothly evolves from the bare atomic resonance in the Boltzmann regime to a frequency corresponding to nearly one Fermi energy at the lowest temperatures. At high temperatures, the peak width reflects the scattering rate of the atoms, while at low temperatures, the width is set by the size of fermion pairs. Above the superfluid transition, and approaching the quantum critical regime, the width increases linearly with temperature, indicating non-Fermi-liquid behavior. From the wings of the rf spectra, we obtain the contact, quantifying the strength of short-range pair correlations. We find that the contact rapidly increases as the gas is cooled below the superfluid transition.

Contact and Sum Rules in a Near-Uniform Fermi Gas at Unitarity

We present an experimental study of the high-energy excitation spectra of unitary Fermi gases. Using focused beam Bragg spectroscopy, we locally probe atoms in the central region of a harmonically trapped cloud where the density is nearly uniform, enabling measurements of the dynamic structure factor for a range of temperatures both below and above the superfluid transition. Applying sum rules to the measured Bragg spectra, we resolve the characteristic behavior of the universal contact parameter, C, across the superfluid transition. We also employ a recent theoretical result for the kinetic (second-moment) sum rule to obtain the internal energy of gases at unitarity.

Contact and Momentum Distribution of the Unitary Fermi Gas

A key quantity in strongly interacting resonant Fermi gases is the contact C, which characterizes numerous properties such as the momentum distribution at large momenta or the pair correlation function at short distances. The temperature dependence of C was measured at unitarity, where existing theoretical predictions differ substantially even at the qualitative level. We report accurate data for the contact and the momentum distribution of the unitary gas in the normal phase, obtained by bold diagrammatic Monte Carlo and Borel resummation. Our results agree with experimental data within error bars and provide crucial benchmarks for the development of advanced theoretical treatments and precision measurements.

Dynamics of Three-Body Correlations in Quenched Unitary Bose Gases

We investigate dynamical three-body correlations in the Bose gas during the earliest stages of evolution after a quench to the unitary regime. The development of few-body correlations is theoretically observed by determining the two- and three-body contacts. We find that the growth of three-body correlations is gradual compared to two-body correlations. The three-body contact oscillates coherently, and we identify this as a signature of Efimov trimers. We show that the growth of three-body correlations depends nontrivially on parameters derived from both the density and Efimov physics. These results demonstrate the violation of scaling invariance of unitary bosonic systems via the appearance of log-periodic modulation of three-body correlations.

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.

Quasiparticle Energy in a Strongly Interacting Homogeneous Bose-Einstein Condensate

Using two-photon Bragg spectroscopy, we study the energy of particlelike excitations in a strongly interacting homogeneous Bose-Einstein condensate, and observe dramatic deviations from Bogoliubov theory. In particular, at large scattering length a the shift of the excitation resonance from the free-particle energy changes sign from positive to negative. For an excitation with wave number q, this sign change occurs at a≈4/(πq), in agreement with the Feynman energy relation and the static structure factor expressed in terms of the two-body contact. For a≳3/q we also see a breakdown of this theory, and better agreement with calculations based on the Wilson operator product expansion. Neither theory explains our observations across all interaction regimes, inviting further theoretical efforts.

Connecting Few-Body Inelastic Decay to Quantum Correlations in a Many-Body System: A Weakly Coupled Impurity in a Resonant Fermi Gas

We study three-body recombination in an ultracold Bose-Fermi mixture. We first show theoretically that, for weak interspecies coupling, the loss rate is proportional to Tan's contact. Second, using a ^{7}Li/^{6}Li mixture we probe the recombination rate in both the thermal and dual superfluid regimes. We find excellent agreement with our model in the BEC-BCS crossover. At unitarity where the fermion-fermion scattering length diverges, we show that the loss rate is proportional to n_{f}^{4/3}, where n_{f} is the fermionic density. This unusual exponent signals nontrivial two-body correlations in the system. Our results demonstrate that few-body losses can be used as a quantitative probe of quantum correlations in many-body ensembles.

Universal High-Momentum Asymptote and Thermodynamic Relations in a Spinless Fermi Gas with a Resonant p-Wave Interaction

We investigate universal relations in a spinless Fermi gas near a p-wave Feshbach resonance, and show that the momentum distribution n_{k} has an asymptote proportional to k^{-2} with the proportionality constant-the p-wave contact-scaling with the number of closed-channel molecules. We prove the adiabatic sweep theorem for a p-wave resonance which reveals the thermodynamic implication of the p-wave contact. In contrast to the unitary Fermi gas in which Tan's contact is universal, the p-wave contact depends on the short-range details of the interaction.

Zero-temperature equation of state of mass-imbalanced resonant Fermi gases

We calculate the zero-temperature equation of state of mass-imbalanced resonant Fermi gases in an ab initio fashion, by implementing the recent proposal of imaginary-valued mass difference to bypass the sign problem in lattice Monte Carlo calculations. The fully nonperturbative results thus obtained are analytically continued to real mass-imbalance to yield the physical equation of state, providing predictions for upcoming experiments with mass-imbalanced atomic Fermi gases. In addition, we present an exact relation for the rate of change of the equation of state at small mass imbalances, showing that it is fully determined by the energy of the mass-balanced system.

Nuclear neutron-proton contact and the photoabsorption cross section

The nuclear neutron-proton contact is introduced, generalizing Tan's work, and evaluated from medium energy nuclear photodisintegration experiments. To this end we reformulate the quasideuteron model of nuclear photodisintegration and establish the bridge between the Levinger constant and the contact. Using experimental evaluations of Levinger's constant, we extract the value of the neutron-proton contact in finite nuclei and in symmetric nuclear matter. Assuming isospin symmetry we propose to evaluate the neutron-neutron contact through the measurement of photonuclear spin correlated neutron-proton pairs.

Local observation of pair condensation in a Fermi gas at unitarity

We present measurements of the local (homogeneous) density-density response function of a Fermi gas at unitarity using spatially resolved Bragg spectroscopy. By analyzing the Bragg response across one axis of the cloud, we extract the response function for a uniform gas which shows a clear signature of the Bose-Einstein condensation of pairs of fermions when the local temperature drops below the superfluid transition temperature. The method we use for local measurement generalizes a scheme for obtaining the local pressure in a harmonically trapped cloud from the line density and can be adapted to provide any homogeneous parameter satisfying the local density approximation.

Energy and contact of the one-dimensional Fermi polaron at zero and finite temperature

We use the T-matrix approach for studying highly polarized homogeneous Fermi gases in one dimension with repulsive or attractive contact interactions. Using this approach, we compute ground state energies and values for the contact parameter that show excellent agreement with exact and other numerical methods at zero temperature, even in the strongly interacting regime. Furthermore, we derive an exact expression for the value of the contact parameter in one dimension at zero temperature. The model is then extended and used for studying the temperature dependence of ground state energies and the contact parameter.

Precise determination of the structure factor and contact in a unitary Fermi gas

We present a high-precision determination of the universal contact parameter in a strongly interacting Fermi gas. In a trapped gas at unitarity, we find the contact to be 3.06±0.08 at a temperature of 0.08 of the Fermi temperature in a harmonic trap. The contact governs the high-momentum (short-range) properties of these systems, and this low-temperature measurement provides a new benchmark for the zero-temperature homogeneous contact. The experimental measurement utilizes Bragg spectroscopy to obtain the dynamic and static structure factors of ultracold Fermi gases at high momentum in the unitarity and molecular Bose-Einstein condensate regimes. We have also performed quantum Monte Carlo calculations of the static properties, extending from the weakly coupled BCS regime to the strongly coupled Bose-Einstein condensate case, that show agreement with experiment at the level of a few percent.

Measurement of the homogeneous contact of a unitary Fermi gas

By selectively probing the center of a trapped gas, we measure the local, or homogeneous, contact of a unitary Fermi gas as a function of temperature. Tan's contact, C, is proportional to the derivative of the energy with respect to the interaction strength and is thus an essential thermodynamic quantity for a gas with short-range correlations. Theoretical predictions for the temperature dependence of C differ substantially, especially near the superfluid transition, T(c), where C is predicted to either sharply decrease, sharply increase, or change very little. For T/T(F)>0.4, our measurements of the homogeneous gas contact show a gradual decrease of C with increasing temperature, as predicted by theory. We observe a sharp decrease in C at T/T(F)=0.16, which may be due to the superfluid phase transition. While a sharp decrease in C below T(c) is predicted by some many-body theories, we find that none of the predictions fully account for the data.

Dynamic spin response of a strongly interacting Fermi gas

We present an experimental investigation of the dynamic spin response of a strongly interacting Fermi gas using Bragg spectroscopy. By varying the detuning of the Bragg lasers, we show that it is possible to measure the response in the spin and density channels separately. At low Bragg energies, the spin response is suppressed due to pairing, whereas the density response is enhanced. These experiments yield the first independent measurements of the spin-parallel and spin-antiparallel dynamic and static structure factors, which provide insight into the different features of density and spin response functions. At high momentum the spin-antiparallel dynamic structure factor displays a universal high frequency tail, proportional to ω(-5/2), where ℏω is the probe energy.

Measurements of Tan's contact in an atomic Bose-Einstein condensate

A powerful set of universal relations, centered on a quantity called the contact, connects the strength of short-range two-body correlations to the thermodynamics of a many-body system with zero-range interactions. We report on measurements of the contact, using rf spectroscopy, for an (85)Rb atomic Bose-Einstein condensate (BEC). For bosons, the fact that contact spectroscopy can be used to probe the gas on short time scales is useful given the decreasing stability of BECs with increasing interactions. A complication is the added possibility, for bosons, of three-body interactions. In investigating this issue, we have located an Efimov resonance for (85)Rb atoms with loss measurements and thus determined the three-body interaction parameter. In our contact spectroscopy, in a region of observable beyond-mean-field effects, we find no measurable contribution from three-body physics.

Probing anisotropic superfluidity in atomic Fermi gases with Rashba spin-orbit coupling

Motivated by the prospect of realizing a Fermi gas with a synthetic non-Abelian gauge field, we investigate theoretically a strongly interacting Fermi gas in the presence of a Rashba spin-orbit coupling. As the twofold spin degeneracy is lifted by spin-orbit interaction, bound pairs with mixed singlet and triplet components emerge, leading to an anisotropic superfluid. We calculate the relevant physical quantities, such as the momentum distribution, the single-particle spectral function, and the spin structure factor, that characterize the system.

Universal energy functional for trapped Fermi gases with short range interactions

Alhassid conjectured that the total energy of a harmonically trapped two-component Fermi gas with a short range interaction is a linear functional of the occupation probabilities of single-particle energy eigenstates. We confirm his conjecture and derive the functional explicitly. We show that the functional applies to all smooth (namely, differentiable) potentials having a minimum, not just harmonic traps. We also calculate the occupation probabilities of high energy states.

Momentum distribution and contact of the unitary Fermi gas

We calculate the momentum distribution n(k) of the unitary Fermi gas by using quantum Monte Carlo calculations at finite temperature T/ϵ(F) as well as in the ground state. At large momenta k/k(F), we find that n(k) falls off as C/k⁴, in agreement with the Tan relations. From the asymptotics of n(k), we determine the contact C as a function of T/ϵ(F) and present a comparison with theory. At low T/ϵ(F), we find that C increases with temperature, and we tentatively identify a maximum around T/ϵ(F) ≃ 0.4. Our calculations are performed on lattices of spatial extent up to N(x) = 14 with a particle number per unit volume of ≃ 0.03-0.07.

Temperature dependence of the universal contact parameter in a unitary Fermi gas

The contact I, introduced by Tan, has emerged as a key parameter characterizing universal properties of strongly interacting Fermi gases. For ultracold Fermi gases near a Feshbach resonance, the contact depends upon two quantities: the interaction parameter 1/(k(F)a), where k(F) is the Fermi wave vector and a is the s-wave scattering length, and the temperature T/T(F), where T(F) is the Fermi temperature. We present the first measurements of the temperature dependence of the contact in a unitary Fermi gas using Bragg spectroscopy. The contact is seen to follow the predicted decay with temperature and shows how pair-correlations at high momentum persist well above the superfluid transition temperature.

Universal relations for identical bosons from three-body physics

Systems consisting of identical bosons with a large scattering length satisfy universal relations determined by 2-body physics that are similar to those for fermions with two spin states. They require the momentum distribution to have a large-momentum 1/k(4) tail and the radio-frequency transition rate to have a high-frequency 1/ω(3/2) tail, both of which are proportional to the 2-body contact. Identical bosons also satisfy additional universal relations that are determined by 3-body physics and involve the 3-body contact, which measures the probability of 3 particles being very close together. The coefficients of the 3-body contact in the 1/k(5) tail of the momentum distribution and in the 1/ω(2) tail of the radio-frequency transition rate are log-periodic functions of k and ω that depend on the Efimov parameter.

Speckle imaging of spin fluctuations in a strongly interacting Fermi gas

Spin fluctuations and density fluctuations are studied for a two-component gas of strongly interacting fermions along the Bose-Einstein condensate-BCS crossover. This is done by in situ imaging of dispersive speckle patterns. Compressibility and magnetic susceptibility are determined from the measured fluctuations. This new sensitive method easily resolves a tenfold suppression of spin fluctuations below shot noise due to pairing, and can be applied to novel magnetic phases in optical lattices.

Establishing the presence of coherence in atomic fermi superfluids: spin-flip and spin-preserving Bragg scattering at finite temperatures

We show how in ultracold Fermi gases the difference between the finite temperature T structure factors, called S_(ω,q), associated with spin and density, reflects coherent order at all ω, q, k(F)a, and T. This observation can be exploited in two photon Bragg scattering experiments on gases which are subject to variable attractive interactions. Our calculations incorporate spin and particle number conservation laws which lead to compatibility at general T with two f-sum rules. Because of its generality a measurement of S_(ω,q) can be a qualitative, direct, in situ approach for establishing superfluid order.