Lellouch-Lüscher relation for ultracold few-atom systems under confinement

This paper derives and numerically validates an analog of the Lellouch-Lüscher relation for few-body bosonic systems, establishing a robust framework to link harmonically trapped state properties to scattering loss rates and enabling precise determination of multi-body scattering rates in confined ultracold experiments.

The Gist

Imagine you are trying to understand how three people bump into each other and stick together to form a pair, leaving one person behind. In the vast, empty space of a park (free space), this is hard to study because the people are moving too fast, there are too many of them, and it's chaotic. You can only guess the rules of the game by watching the crowd.

Now, imagine you put just three of these people into a tiny, invisible, bouncy box (an optical trap). You can control their movements perfectly. You can see exactly how long they stay in the box before one of them gets kicked out.

This paper is about building a mathematical bridge between what happens in that tiny, controlled box and what happens in the chaotic, open world.

Here is the breakdown of the paper's big idea, using simple analogies:

1. The Problem: The "Box" vs. The "World"

For decades, scientists have been great at studying atoms in huge clouds (like a foggy room). They measure how fast atoms disappear (loss rates) to figure out how they interact.

• The Challenge: Recently, scientists learned how to trap just 2 or 3 atoms in a tiny "box" (using lasers called tweezers). This is amazing because it's like watching a single billiard ball instead of a whole pool table.
• The Catch: Because the atoms are in a box, they bounce around differently than they would in free space. The "rules" of the box change the outcome. Scientists needed a way to translate the results from the "box experiment" back to the "real world" rules.

2. The Solution: The "Lellouch-Lüscher" Translator

The authors created a new version of a famous formula (originally used in particle physics) called the Lellouch-Lüscher (LL) relation.

Think of this formula as a universal translator or a currency exchange rate:

• Input: You measure how long a specific "trapped state" lasts (its lifetime) and how much energy it has inside the box.
• Output: The formula instantly tells you the exact speed at which those atoms would collide and react if they were in free space.

It's like measuring how fast a car slows down on a specific test track and using a formula to tell you exactly how fast it would crash on a highway.

3. The Key Discovery: The "Width" of the State

In the paper, the scientists focus on "width" ($\Gamma$).

• The Analogy: Imagine a bell ringing. If the bell rings clearly and for a long time, it has a "narrow" sound. If the bell is muffled and stops quickly, it has a "broad" sound.
• In quantum physics, a trapped atom state that disappears quickly (short life) is "broad." One that stays for a long time is "narrow."
• The paper proves that how "broad" the state is in the trap is directly linked to how fast the atoms would collide in free space.

4. Why This is a Big Deal (The "Super-Scanner")

The authors tested this with computer simulations of Rubidium atoms (a common element in these experiments). They found that this "translator" works perfectly, even when the atoms are interacting very strongly.

Why does this matter to the real world?

• Precision: In a big cloud of gas, atoms are moving at different speeds and hitting each other in different ways (like a mosh pit). It's hard to tell exactly which type of collision caused a reaction.
• Clarity: In a trapped few-atom system, you can isolate a single "dance move" (a specific quantum state). The LL relation allows scientists to measure the rules of that one specific move with incredible precision.
• Future Tech: This helps us understand complex quantum chemistry and could lead to better quantum computers. It turns the "messy" gas experiments into "clean" laboratory measurements.

5. The "N-Body" Extension

The paper doesn't just stop at three atoms. They generalized the math to work for N atoms (4, 5, 10, etc.).

• The Analogy: If you know how to translate the rules for a trio, you can now translate the rules for a whole dance troupe. This allows scientists to figure out how groups of atoms interact, which is crucial for understanding new states of matter.

Summary

The paper says: "Stop guessing how atoms behave in the wild by looking at the whole crowd. Instead, trap a few of them in a box, measure how long they stay, and use our new formula to know exactly how they would behave in the real world."

It transforms the study of the very small from a game of "guessing the crowd's behavior" into a precise science of "measuring individual interactions."

Technical Summary

1. Problem Statement

Trapped ultracold atoms in optical lattices and tweezers offer unprecedented control over few-body quantum systems, enabling the study of fundamental chemical reactions and multi-body interactions. However, a critical theoretical challenge remains: relating observables measured in finite-volume confined systems (traps) to free-space scattering parameters.

In bulk ultracold gases, interactions are characterized by scattering loss rates (e.g., the three-body recombination rate $L_3$). In contrast, trapped systems are characterized by discrete energy levels and lifetimes of specific quantum states. While the Lellouch-Lüscher (LL) relation is a well-established tool in high-energy physics (specifically Lattice QCD) to extract scattering observables from finite-volume energy spectra, an analogous rigorous framework for ultracold few-body bosonic systems in harmonic traps was lacking. Without such a relation, extracting precise, partial-wave-resolved scattering rates from trap experiments is difficult, especially when distinguishing between different decay pathways and thermal averaging effects.

2. Methodology

The authors developed a theoretical framework and performed numerical simulations to derive and validate an LL relation for few-body systems.

• Theoretical Derivation:

  • The authors utilized the hyperspherical adiabatic representation to solve the three-body Schrödinger equation.
  • They modeled the system as three identical bosons in an isotropic harmonic trap ($V_{ho} = \mu \omega_{ho}^2 R^2 / 2$).
  • They identified that the mechanism governing both the free-space recombination rate ($L_3$) and the trapped state width ($\Gamma_q$) is the same: inelastic transitions occurring at short distances ($R \le R_0 \sim r_{vdW}$).
  • By defining a reactive three-body probability density ($D_q$) within the reaction region, they derived a proportionality between the decay rate of a trapped state and the free-space loss coefficient.
  • They generalized this derivation to $N$-body systems.

• Numerical Simulations:

  • The team performed numerical simulations for $^{85}\text{Rb}$ (strongly attractive interactions) and $^{87}\text{Rb}$ (weakly repulsive interactions).
  • They used realistic interatomic potentials derived from Born-Oppenheimer potentials, incorporating full atomic hyperfine spin structures and over 100 atom-dimer channels.
  • They calculated the time delay ($\tau_D$) to extract the energies ($E_q$) and widths ($\Gamma_q$) of the trapped states, identifying resonances as Lorentzian peaks.
  • They compared these numerical widths against predictions made using the derived LL relation and free-space $L_3(E)$ calculations.

3. Key Contributions

A. Derivation of the Three-Body LL Relation

The paper derives a specific LL relation connecting the free-space three-body loss coefficient $L_3(E_q)$ to the width $\Gamma_q$ of the $q$-th trapped state with energy $E_q$:

$$L_3(E_q) = C \frac{\hbar^4}{m^3} \frac{\Gamma_q}{\omega_{ho} [E_q^2 - (\hbar\omega_{ho})^2]}$$

Where:

• $m$ is the atomic mass.
• $\omega_{ho}$ is the trap frequency.
• $C = 72\sqrt{3}\pi^3$ is a universal constant.
• The relation holds for the lowest partial wave ($J=0, \lambda=0$).

B. Generalization to N-Body Systems

The authors extended the relation to $N$ identical bosons, providing a formula to extract the $N$-body loss coefficient $L_N$ from the width and energy of trapped $N$-body states. This allows for the systematic determination of multi-body scattering rates, which is crucial for understanding many-body physics in optical lattices.

C. Handling Interaction Strengths

The derivation accounts for interaction-induced shifts in energy and spatial extent. The authors demonstrate that while the relation is exact in the weakly interacting limit ($|a|/a_{ho} \ll 1$), it remains robust in the non-perturbative regime ($|a|/a_{ho} \sim 1$) if the actual interacting energy $E_q$ (rather than the non-interacting harmonic oscillator energy) is used in the formula.

D. Thermal Equilibrium Extension

The paper provides an extension of the LL relation for experiments where atoms are loaded into traps at finite temperatures ($k_B T \gg \hbar\omega_{ho}$), deriving expressions for the thermal decay rate in both the threshold and unitary limits.

4. Key Results

• Validation of the Relation: Numerical simulations for $^{85}\text{Rb}$ and $^{87}\text{Rb}$ show excellent agreement between the widths $\Gamma_q$ calculated directly from the time-delay method and those predicted by the LL relation using free-space $L_3(E)$ data.
• Isotope-Specific Behavior:
  • $^{87}\text{Rb}$ (Weak Interaction): The loss rate $L_3(E)$ is relatively constant (threshold behavior) but enhanced by a d-wave shape resonance. Consequently, the trapped state width $\Gamma_q$ shows a strong dependence on energy $E_q$.
  • $^{85}\text{Rb}$ (Strong Interaction): The system enters the unitary regime where $L_3(E) \propto 1/E^2$. This leads to a width $\Gamma_q$ that is nearly independent of the state energy $E_q$, consistent with the LL prediction.
• Partial Wave Resolution: The relation allows for the extraction of scattering rates for a single partial wave ($J=0$). This is a significant advantage over bulk gas experiments where thermal averaging obscures individual partial-wave contributions and resonance details.
• Multi-Body Decay Pathways: The generalized $N$-body relation enables the deconvolution of total decay rates into fundamental few-body rates (e.g., separating 3-body, 4-body, and 5-body contributions in a system of $N$ atoms).

5. Significance

• Theoretical Framework: This work establishes a rigorous bridge between finite-volume trapped experiments and free-space scattering theory, filling a gap in few-body physics similar to the impact of the original LL relations in lattice QCD.
• Experimental Utility: It provides a practical tool for experimentalists using optical tweezers and lattices to determine precise scattering rates and resonance parameters (like Efimov parameters) without the complications of thermal averaging found in bulk gases.
• Precision Metrology: By enabling the measurement of individual partial-wave contributions, the LL relation offers a more sensitive probe for quantum interference and resonance phenomena, potentially leading to better control over chemical reactions and quantum state engineering.
• Scalability: The generalization to $N$-body systems opens new avenues for studying complex many-body interactions and the emergence of few-body physics in larger ensembles, which is essential for the development of quantum simulators and quantum information processors.

In summary, the paper successfully adapts the Lellouch-Lüscher formalism to ultracold atomic physics, providing a powerful, validated method to extract fundamental scattering parameters from the discrete spectra of confined few-body systems.