Universal Threshold Laws for Dipolar Evaporation of Molecular Fermi Gases

Technical Summary: Universal Threshold Laws for Dipolar Evaporation of Molecular Fermi Gases

Problem Statement
Evaporative cooling is the primary method for achieving bulk quantum degeneracy in ultracold gases. However, for single-component Fermi gases, a fundamental bottleneck arises at deep degeneracy due to Pauli blocking. As the Fermi sea fills, the available phase space for elastic collisions—required to rethermalize the gas during the removal of energetic molecules—is suppressed. This suppression limits the efficiency of evaporative cooling, preventing the system from reaching lower temperatures. While collisional shielding techniques have successfully suppressed short-range chemistry and sticky collisions in dipolar molecular gases (e.g., KRb), they do not inherently solve the problem of Fermi-suppressed rethermalization. The central question addressed is how the competition between rethermalization and particle loss evolves toward deep degeneracy, and whether the effective dimensionality of the gas (3D vs. 2D) fundamentally alters this limit.

Methodology
The authors employ a combination of analytical variational methods and ab initio numerical calculations to investigate the low-temperature behavior of shielded single-component dipolar Fermi gases.

1. Analytical Framework:

   • The system is modeled as $N$ fermionic dipoles in a harmonic trap, considering both 3D and quasi-2D regimes (where confinement along the $z$-axis is strong enough to occupy the motional ground state but not strictly 2D).
   • The rethermalization rate ($\gamma$) is calculated using the Boltzmann equation and a variational approach. A Fermi suppression factor ($\eta_F$) is defined to quantify the reduction in thermalization relative to a classical gas.
   • The analysis focuses on the Wigner threshold regime, where collision energies are sufficiently low that scattering properties are governed by universal threshold laws.
   • The thermalization-to-loss ratio ($\xi = \gamma/\gamma_L$) is derived by combining the rethermalization rate with the loss rate ($\gamma_L$), which is determined by the loss cross-section ($\sigma_L$).

2. Numerical Verification:

   • To validate the universal scaling laws, the authors perform ab initio close-coupling calculations for specific molecular species: KRb and NaK.
   • These calculations cover different shielding techniques: static electric field shielding (for both KRb and NaK) and microwave shielding (for NaK).
   • The numerical methods utilize a universal loss short-range boundary condition and adaptive step-size propagation, including up to $\ell=9$ collisional partial waves to ensure convergence.
   • Theoretical predictions are benchmarked against experimental data from the JILA KRb experiment.

Key Contributions and Results

• Universal Scaling of Fermi Suppression:
  The paper establishes that the Fermi suppression of elastic rethermalization ($\eta_F$) follows distinct power laws dependent on dimensionality in the deep degeneracy limit ($T/T_F \ll 1$):

  • In 3D, $\eta_F \propto (T/T_F)^2$.
  • In 2D, $\eta_F \propto T/T_F$.
    This indicates that rethermalization is less suppressed in 2D than in 3D at very low temperatures. The difference arises from the density of states (quadratic in 3D, linear in 2D) rather than microscopic scattering details.

• Universal Scaling of the Thermalization-to-Loss Ratio:
  The efficiency of evaporative cooling is governed by the ratio $\xi = \gamma/\gamma_L$. The authors find that while $\xi$ vanishes as $T \to 0$ in both dimensions (rendering evaporation eventually inefficient), the rate of degradation differs:

  • In 3D, $\xi_{3D} \propto T^{5/2}$. This results from the quadratic suppression of rethermalization combined with a loss rate that saturates (due to finite Fermi momentum) and a classical collision rate scaling as $\sqrt{T}$.
  • In 2D, $\xi_{2D} \propto T^2$. The linear suppression of rethermalization and linear scaling of the 2D collision rate lead to a slower degradation of cooling efficiency.

• Validation of Threshold Laws:

  • Ab initio calculations for shielded KRb and NaK confirm that both elastic and loss cross-sections adhere to the predicted threshold laws ($\sigma_{el} \sim E_c^0$ in 3D, $\sigma_{el} \sim E_c^{1/2}$ in 2D; $\sigma_L \sim E_c^{1/2}$ in both).
  • The calculations demonstrate that while the bare two-body loss cross-section follows threshold laws, the effective in-medium loss cross-section saturates at low temperatures due to the finite Fermi momentum, contrasting with the vanishing effective elastic cross-section.
  • The theoretical curves for effective cross-sections and the thermalization-to-loss ratio show strong agreement with experimental measurements from the JILA KRb experiment.

Significance and Claims
The paper claims to identify universal threshold laws that govern the low-temperature limit of evaporative cooling for single-component dipolar Fermi gases. The primary significance lies in demonstrating that dimensionality sets the fundamental scaling exponents for cooling efficiency, independent of the specific shielding technique or molecular species used.

The authors conclude that while evaporative cooling inevitably becomes inefficient as $T \to 0$ in any geometry due to the interplay of Fermi suppression and finite loss rates, 2D geometries are better suited for reaching deeper Fermi degeneracy. The slower degradation of the thermalization-to-loss ratio in 2D ($\propto T^2$) compared to 3D ($\propto T^{5/2}$) suggests that 2D systems can sustain efficient evaporation to lower temperatures. For example, the authors estimate that the onset of inefficient cooling occurs at approximately 3.7 nK in 3D but extends to 440 pK in 2D for KRb.

The work does not propose new experimental protocols but provides a theoretical framework explaining why 2D geometries have been successful in recent experiments and predicts the ultimate limits of cooling for future deep-degeneracy experiments. The authors note that their treatment assumes a dilute gas and neglects Hartree-Fock mean-field effects, which they expect to have only a minor impact on the universal scaling results.