Ultracold Fermi Gases Research Articles

Second-order Sobolev gradient flows for computing ground state of ultracold Fermi gases

Based on density functional theory, the ground state for ultracold Fermi gases with dipole–dipole interaction is a functional minimization problem with orthonormality constraints. Many methods such as direct minimization, self-consistent iteration method and first-order gradient flow had been proposed to solve such problem. In this paper, we introduce the second-order Sobolev gradient flows, solve them to the steady state and get the ground state of ultracold Fermi gases numerically. Two contributions have been made: firstly we extend the usual first-order gradient flows to the second-order gradient flows; secondly, we extend the usual L2 gradient to the Sobolev gradient and get the Sobolev gradient flows. We prove that the second-order Sobolev gradient flows have the properties of orthonormality preserving and ‘modified’ energy diminishing, which are desirable for the computation of the ground state solution. Several numerical techniques have been used to discretize the time-dependent second-order Sobolev gradient flows. These include explicit leap-frog method and stabilized semi-implicit method for temporal discretization and Fourier pseudo-spectral method for spatial variables. Extensive numerical examples for the ground state are reported to demonstrate the power of the numerical methods.

Spectral properties and observables in ultracold Fermi gases

Spectral properties and observables in ultracold Fermi gases

Probing molecular spectral functions and unconventional pairing using Raman spectroscopy

An impurity interacting with an ultracold Fermi gas can form either a polaron state or a dressed molecular state, the molaron, in which the impurity forms a bound state with one gas particle. This molaron state features rich physics, including a negative effective mass around unitarity and a first-order transition to the polaron state. However, these features have remained so far experimentally inaccessible. In this work we show theoretically how the molaron state can be directly prepared experimentally even in its excited states using Raman spectroscopy techniques. Initializing the system in the ultrastrong coupling limit, where the binding energy of the molaron is much larger than the Fermi energy, our protocol maps out the momentum-dependent spectral function of the molecule. Using a diagrammatic approach we furthermore show that the molecular spectral function serves as a direct precursor of the elusive Fulde-Ferell-Larkin-Ovchinnikov phase, which is realized for a finite density of fermionic impurity particles. Our results pave the way to a systematic understanding of how composite particles form in quantum many-body environments and provide a basis to develop new schemes for the observation of exotic phases of quantum many-body systems. Published by the American Physical Society 2024

Neural-network quantum states for ultra-cold Fermi gases

Ultra-cold Fermi gases exhibit a rich array of quantum mechanical properties, including the transition from a fermionic superfluid Bardeen-Cooper-Schrieffer (BCS) state to a bosonic superfluid Bose-Einstein condensate (BEC). While these properties can be precisely probed experimentally, accurately describing them poses significant theoretical challenges due to strong pairing correlations and the non-perturbative nature of particle interactions. In this work, we introduce a Pfaffian-Jastrow neural-network quantum state featuring a message-passing architecture to efficiently capture pairing and backflow correlations. We benchmark our approach on existing Slater-Jastrow frameworks and state-of-the-art diffusion Monte Carlo methods, demonstrating a performance advantage and the scalability of our scheme. We show that transfer learning stabilizes the training process in the presence of strong, short-ranged interactions, and allows for an effective exploration of the BCS-BEC crossover region. Our findings highlight the potential of neural-network quantum states as a promising strategy for investigating ultra-cold Fermi gases.

Nonlinear localization of ultracold atomic Fermi gas in moiré optical lattices

Nonlinear localization of ultracold atomic Fermi gas in moiré optical lattices

Magnetic impurity states of a Fermi system in a local magnetic field

The ultracold atomic Fermi gas provides a flexible and controllable experimental platform for exploring magnetic impurity states. In this study, we theoretically examine the magnetic impurity states in a Fermi system with spin-1/2 and spin-3/2 under a local magnetic field. Firstly, the spin-1/2 impurity system contains three types of magnetic solutions. One type magnetic solution with higher energy can be regarded as the result of induced magnetization of a trivial solution, while the other two types magnetic solutions with lower energy can be understood as the results of Zeeman splitting of nontrivial solutions with double degeneracy. Notably, the magnetic solution corresponding to the spin state parallel to the Zeeman magnetic field is always the most stable in terms of the energy. Secondly, we extend the idea of studying spin-1/2 impurity to spin-3/2 impurity, where the half-filled state of impurity atoms is calculated. We found that four types of magnetic solutions at zero magnetic field were split into four, six, four, and twelve new magnetic solutions in the presence of a Zeeman magnetic field, respectively. We obtained that the magnetic solution corresponding to the spin state parallel to the magnetic field is the most stable, and its energy is lower than that of the corresponding magnetic solution at zero magnetic field. Therefore, the introduction of a local magnetic field is beneficial for the formation of magnetic impurity states. Finally, the calculated data illustrate that the Anderson impurity model we adopt exhibits invariance under the combined transformations of spin-flip and particle–hole.

Fermionic quantum turbulence: Pushing the limits of high-performance computing.

Ultracold atoms provide a platform for analog quantum computer capable of simulating the quantum turbulence that underlies puzzling phenomena like pulsar glitches in rapidly spinning neutron stars. Unlike other platforms like liquid helium, ultracold atoms have a viable theoretical framework for dynamics, but simulations push the edge of current classical computers. We present the largest simulations of fermionic quantum turbulence to date and explain the computing technology needed, especially improvements in the Eigenvalue soLvers for Petaflop Applications library that enable us to diagonalize matrices of record size (millions by millions). We quantify how dissipation and thermalization proceed in fermionic quantum turbulence by using the internal structure of vortices as a new probe of the local effective temperature. All simulation data and source codes are made available to facilitate rapid scientific progress in the field of ultracold Fermi gases.

A lattice pairing-field approach to ultracold Fermi gases

We develop a pairing-field formalism for ab initio studies of non-relativistic two-component fermions on a (d+1)-dimensional spacetime lattice. More specifically, we focus on theories where the interaction between the two components can be described by the exchange of a corresponding pairing field. The introduction of a pairing field may indeed be convenient for studies of, e.g., the finite-temperature phase structure and critical behavior of, e.g., ultracold atomic Fermi gases. Moreover, such a formalism allows to directly compute the momentum and frequency dependence of the pair propagator, from which the pair-correlation function can be extracted. For a first illustration of the application of our formalism, we compute the density equation of state and the superfluid order parameter for a gas of unpolarized fermions in (0+1) dimensions by employing the complex Langevin approach to surmount the sign problem.

Flat band effects on the ground-state BCS-BEC crossover in atomic Fermi gases in a quasi-two-dimensional Lieb lattice

The ground-state superfluid behavior of ultracold atomic Fermi gases with a short-range attractive interaction in a quasi-two-dimensional Lieb lattice is studied using BCS mean-field theory, within the context of BCS-BEC crossover. We find that the flat band leads to nontrivial exotic effects. As the Fermi level enters the flat band, both the pairing gap and the in-plane superfluid density exhibit an unusual power law as a function of interaction, with strongly enhanced quantum geometric effects, in addition to a dramatic increase of compressibility as the interaction approaches the BCS limit. As the Fermi level crosses the van Hove singularities, the character of pairing changes from particle-like to hole-like or vice versa. We present the computed phase diagram, in which a pair density wave state emerges at high densities with relatively strong interaction strength.

Electromagnetically Induced Transparency Spectra of 6Li Rydberg Atoms

Rydberg atoms possess highly excited valence electrons that are far away from atomic cations. Compared with ground states, Rydberg states are excited states with a high principal quantum number n that exhibit large electric dipole moments and have a variety of applications in quantum information processing. In this communication, we report the measurement of the 6Li Rydberg excitation spectrum by ladder-type electromagnetically induced transparency (EIT) in a vapor cell. The 2p→ns/nd EIT spectra were recorded by sweeping the frequency of an ultraviolet Rydberg pumping laser while keeping the probing laser resonant to the 2s→2p transition. All lasers were locked on an ultrastable optical Fabry-Pérot cavity and measured by an optical frequency comb. Our results provide valuable information to precisely determine quantum defects and enable novel experiments with Rydberg-dressed ultracold Fermi gases.

Recent advances in the theory of the BCS–BEC crossover for fermionic superfluidity

The BCS–BEC crossover realized experimentally with ultra-cold Fermi gases may be considered as one of the important scientific achievements occurred during the last several years. The flexibility for operating on these systems on the experimental side and the full control of the relevant system degrees of freedom on the theoretical side make quite stringent at a fundamental level the comparison between the experimental data and the corresponding theoretical calculations. Here, we briefly survey recent theoretical advances resting on a diagrammatic approach at equilibrium that improves in a systematic way on the widely used t-matrix approach, yielding a quite good comparison between theory and experiments for several physical quantities of interest. It is proposed that the physical phenomena underlying this theoretical approach may also be relevant to the superconducting phase of condensed-matter materials which cannot be described by the standard BCS theory.

Scattering mechanism of work done on an ultracold Fermi gas by a time-dependent interaction

We discuss the work done on repulsive ultracold Fermi gas via collisions tuned by a time-dependent magnetic ramp. If the interaction is strengthened during binary collisions, the kinetic energy of atoms increases with each such inelastic event. The net energy transfer via collisions is compared with the transfer via the mean field. It is found that for a weak interaction, with a scattering length shorter than the inverse mean momentum of atoms, the mean field dominates, while for a strong interaction, the majority of energy enters the gas via collisions.

Quantum Monte Carlo study of the role of p -wave interactions in ultracold repulsive Fermi gases

Single-component ultracold atomic Fermi gases are usually described using noninteracting many-fermion models. However, recent experiments reached a regime where $p$-wave interactions among identical fermionic atoms are important. In this paper, we employ variational and fixed-node diffusion Monte Carlo simulations to investigate the ground-state properties of single-component Fermi gases with short-range repulsive interactions. We determine the zero-temperature equation of state, and elucidate the roles played by the $p$-wave scattering volume and the $p$-wave effective range. A comparison against recently derived second-order perturbative results shows good agreement in a broad range of interaction strength. We also compute the quasiparticle effective mass, and we confirm the perturbative prediction of a linear contribution in the $p$-wave scattering volume, while we find significant deviations from the beyond-mean-field perturbative result, already for moderate interaction strengths. Finally, we determine ground-state energies for two-component unpolarized Fermi gases with both interspecies and intraspecies hard-sphere interactions, finding remarkable agreement with a recently derived fourth-order expansion that includes $p$-wave contributions.

Toward precision Fermi-liquid theory in two dimensions

The ultra-cold and weakly-coupled Fermi gas in two spatial dimensions is studied in an effective field theory framework. It has long been observed that universal corrections to the energy density to two orders in the interaction strength do not agree with Monte Carlo simulations in the weak-coupling regime. Here, universal corrections to three orders in the interaction strength are obtained for the first time, and are shown to provide agreement between theory and simulation. Special consideration is given to the scale ambiguity associated with the non-trivial renormalization of the singular contact interactions. The isotropic superfluid gap is obtained to next-to-leading order, and nonuniversal contributions to the energy density due to effective range effects, p-wave interactions and three-body forces are computed. Results are compared with precise Monte Carlo simulations of the energy density and the contact in the weakly-coupled attractive and repulsive Fermi liquid regimes. In addition, the known all-orders sum of ladder and ring diagrams is compared with Monte Carlo simulations at weak coupling and beyond.

Itinerant ferromagnetism in dilute SU(N) Fermi gases

We present exact analytic results for the energy of a SU(N) repulsive Fermi gas as a function of the spin-channel occupation at second order in the gas parameter. This is an extension of previous results that now incorporates the degree of polarization of the system. Therefore, the magnetic properties of the gas can be obtained, free from numerical uncertainties. For spin 1/2 we find that second-order corrections change the itinerant ferromagnetic transition from continuous to first-order. Instead, for spin larger than 1/2 the phase transition is always of first-order type. The transition critical density reduces when the spin increases, making the phase transition more accessible to experiments with ultracold dilute Fermi gases. Estimations for Fermi gases of Yb and Sr with spin 5/2 and 9/2, respectively, are reported.

Disordered structures in ultracold spin-imbalanced Fermi gas

We investigate properties of spin-imbalanced ultracold Fermi gas in a large range of spin polarizations at low temperatures. We present results of microscopic calculations based on mean-field and density functional theory approaches, with no symmetry constraints. At low polarization values we predict the structure of the system as consisting of several spin-polarized droplets. As the polarization increases, the system self-organizes into a disordered structures similar to liquid crystals, and energetically they can compete with ordered structures such as grid-like domain walls. At higher polarizations the system starts to develop regularities that, in principle, can be called supersolid, where periodic density modulation and pairing correlations coexist. The robustness of the results has been checked with respect to temperature effects, dimensionality, and the presence of a trapping potential. Dynamical stability has also been investigated.

Spectrum, Lifshitz transitions and orbital current in frustrated fermionic ladders with a uniform flux

Ultracold Fermi gases with synthetic gauge field represent an excellent platform to study the combined effect of lattice frustration and an effective magnetic flux close to one flux quantum per particle. The minimal theoretical model to accomplish this task is a system of spinless noninteracting fermions on a triangular two-chain flux ladder. In this paper we consider this model close to half-filling, with interchain hopping amplitudes alternating along the zigzag bonds of the ladder and being small. In such setting qualitative changes in the ground state properties of the model, most notably the flux-induced topological Lifshitz transitions are shifted towards the values of the flux per a diatomic plaquette ($f$) close to the flux quantum ($f\sim 1/2$). Geometrical frustration of the lattice breaks the $k \to \pi - k$ particle-hole spectral symmetry present in non-frustrated ladders, and leads to splitting of the degeneracies at metal-insulator transitions. A remarkable feature of a translationally invariant triangular ladder is the appearance of isolated Dirac node at the Brillouin zone boundary, rendering the ground state at $f=1/2$ semi-metallic. We calculate the flux dependence of the orbital current and identify Lifshitz critical points by the singularities in this dependence at a constant chemical potential $\mu$ and constant particle density $\rho$. The absence of the particle-hole symmetry in the triangular ladder leads to the asymmetry between particle ($\rho>1$) and hole ($\rho<1$) doping and to qualitatively different results for the phase diagram, Lifshitz points and the flux dependence of the orbital current in the settings $\rho={\rm const}$ and $\mu={\rm const}$.

Rydberg electromagnetically induced transparency in 40K ultracold Fermi gases

Rydberg electromagnetically induced transparency in 40K ultracold Fermi gases

Quantized Nonlinear Transport with Ultracold Atoms

In this letter, we propose how to measure the quantized nonlinear transport using two-dimensional ultracold atomic Fermi gases in a harmonic trap. This scheme requires successively applying two optical pulses in the left and lower half-planes and then measuring the number of extra atoms in the first quadrant. In ideal situations, this nonlinear density response to two successive pulses is quantized, and the quantization value probes the Euler characteristic of the local Fermi sea at the trap center. We investigate the practical effects in experiments, including finite pulse duration, finite edge width of pulses, and finite temperature, which can lead to deviation from quantization. We propose a method to reduce the deviation by averaging measurements performed at the first and third quadrants, inspired by symmetry considerations. With this method, the quantized nonlinear response can be observed reasonably well with experimental conditions readily achieved with ultracold atoms.

Feasibility of a Fulde-Ferrell-Larkin-Ovchinnikov superfluid Fermi atomic gas

We theoretically explore a promising route to achieve the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state in a spin-imbalanced ultracold Fermi gas. In the current stage of cold atom physics, search for this exotic Fermi superfluid is facing two serious difficulties: One is the desperate destruction of the FFLO long-range order by FFLO pairing fluctuations, which precludes entering the phase through a second-order transition, even in three dimension. The other is the fierce competition with the phase separation into the BCS (Bardeen-Cooper-Schrieffer) state and the spin-polarized normal state. Including strong FFLO pairing fluctuations within the framework of the strong-coupling theory developed by Nozi\`eres and Schmitt-Rink, we show that the anisotropy of Fermi surface introduced by an optical lattice makes the FFLO state stable against the paring fluctuations. This stabilized FFLO state is also found to be able to overcome the competition with the phase separation under a certain condition. Since the realization of unconventional Fermi superfluids is one of the most exciting challenges in cold atom physics, our results would contribute to the further development of this field.