Phenomenological model of decaying Bose polarons

This paper introduces a phenomenological model for decaying Bose polarons based on a minimal variational wave function and complex interaction strength, which successfully explains the significant line broadening observed in strong-interaction experiments by accounting for the decay of single-boson correlated states into multi-boson states.

The Gist

The Big Picture: The "Party Crasher" in a Crowd

Imagine a massive, perfectly synchronized dance party. Everyone is moving in perfect unison, swaying to the same rhythm. This is a Bose-Einstein Condensate (BEC), a state of matter where atoms act like a single, giant super-atom.

Now, imagine one person (the Impurity) walks into this dance floor. They don't know the dance. As they try to move, they bump into the dancers, causing ripples and changes in the crowd's movement. The crowd starts to swirl around this new person, forming a protective "cloud" of dancers.

In physics, this new combined entity—the clumsy person plus the swirling crowd—is called a Bose Polaron. It acts like a new, heavier particle with its own unique personality.

The Problem: The "Fuzzy" Signal

Scientists love studying these Polaron "party crashers" because they tell us how particles interact. Usually, when you look at a Polaron, you see a sharp, clear signal, like a distinct note on a piano.

However, when the interaction between the intruder and the crowd gets very strong, something weird happens. Instead of a sharp note, the signal becomes a fuzzy, broad blur. It's like the piano key is being held down while someone is shaking the piano.

For a long time, physicists were confused. Why does the Polaron lose its sharp identity so quickly when the interactions get strong? Is the Polaron even real anymore, or is it just a fleeting moment before it falls apart?

The New Idea: The "One-Step" Dance

The authors of this paper (a team from Austria, Denmark, and Austria) decided to solve this mystery with a new, simpler way of thinking.

The Old Way (The Hard Math):
Previously, scientists tried to calculate exactly how the intruder interacts with every single person in the crowd simultaneously. They tried to track the intruder holding hands with 1 dancer, then 2 dancers, then 3, and so on.

• The Problem: This is like trying to track every single conversation in a stadium. It's computationally impossible to do perfectly, and the math gets messy.

The New Way (The "One-Step" Model):
The authors proposed a simpler idea: The most important thing that happens is the intruder interacting with just one dancer at a time.

Think of it like this: When you walk into a crowded room, you might bump into one person first. That interaction is the loudest and most noticeable. Even if that person eventually introduces you to their friends (a group of 3 or 4), your first impression is defined by that single bump.

The authors realized that experimental tools (like radio waves used to measure the Polaron) mostly "see" this single-dancer interaction. They don't see the complex, deep group hugs happening in the background.

The Secret Ingredient: The "Leaky Bucket"

Here is the clever part. The authors knew that this "single-dancer" state is unstable. It's like a balloon that is slowly leaking air. The Polaron forms, but then it quickly decays into a more complex state (involving many dancers) that the experiment can't easily see.

To model this without doing impossible math, they introduced a "Leaky Bucket" concept (technically called a complex interaction strength).

• Real World: You pour water into a bucket (the Polaron forms).
• The Leak: The bucket has a hole (the decay). The water level drops over time.
• The Math: Instead of calculating exactly where every drop of water falls out of the hole, they just added a "leak factor" to their equation. This factor represents the fact that the Polaron is dying or "decohering" (losing its rhythm) as it tries to interact with the crowd.

What Did They Find?

They tested this "Leaky Bucket" model against two real-world experiments (one from Cambridge and one from Aarhus).

1. The Spectrum (The Sound): When they looked at the "fuzzy" signals in the experiments, their simple model predicted the exact width and position of the blur. It matched perfectly.
2. The Dynamics (The Rhythm): They also looked at how the Polaron behaves over time. The experiments showed the Polaron's rhythm dying out much faster than old theories predicted. Their "Leaky Bucket" model captured this rapid decay perfectly.

Why Does This Matter?

This paper is a breakthrough because it offers a practical toolkit for scientists.

• Before: Scientists were stuck trying to solve an unsolvable equation involving millions of variables, or they were confused why their theories didn't match the "fuzzy" experiments.
• Now: They have a simple, intuitive framework. They can say, "Okay, the Polaron is mostly interacting with one particle, but it's leaking energy into the rest of the crowd."

The Takeaway:
Sometimes, to understand a complex system, you don't need to track every single detail. You just need to understand the most likely interaction and acknowledge that things fall apart. By accepting that the Polaron is a "leaky" state that decays quickly, the authors finally explained why the experimental data looks so blurry.

It's like realizing that the reason a song sounds muffled isn't because the singer is bad, but because the microphone is slightly broken. Once you account for the broken microphone (the decay), the song makes perfect sense.

Technical Summary

1. Problem Statement

The Bose polaron describes a mobile impurity immersed in a Bose-Einstein condensate (BEC). While theoretical models successfully describe weak to moderate interactions, recent experiments (specifically with $^{39}$K atoms) reveal significant discrepancies in the strong interaction regime:

• Spectral Broadening: Instead of sharp quasiparticle peaks, experiments observe broad spectral lines for strong interactions, even in homogeneous box potentials where density broadening is minimized.
• Rapid Decoherence: Non-equilibrium dynamics (measured via Ramsey interferometry) show a much faster loss of coherence than predicted by standard theories.
• Theoretical Gap: Unlike Fermi polarons, Bose polarons lack a Pauli exclusion principle, allowing the impurity to correlate with $N \geq 2$ bosons (forming trimers, tetramers, etc.). Standard variational approaches including up to two bosons (or Bogoliubov modes) fail to capture the full spectral width and decay dynamics, suggesting that the experimentally observed "polaron" is not the true ground state but a broad resonance that decays into lower-energy, multi-boson correlated states.

2. Methodology

The authors propose a phenomenological model that bridges the gap between microscopic complexity and experimental observability.

A. Motivation via Variational Analysis

• The authors first analyzed the system using an elaborate variational wavefunction (Eq. 2) that includes correlations with up to two Bogoliubov modes (and closed-channel molecules).
• Key Finding: The spectral function reveals two branches:
  1. A sharp, low-energy ground state with vanishing spectral weight (hard to detect experimentally).
  2. A broad, high-energy branch with large spectral weight, corresponding to an impurity dressed by at most one boson.
• Conclusion: Experimental probes (RF spectroscopy, Ramsey interferometry) couple primarily to the single-boson dressed state. This state is metastable and decays into the lower-energy multi-boson ground state.

B. The Phenomenological Model

To model this decay without computationally intractable many-body calculations, the authors introduce a complex impurity-boson interaction strength ($g_I$).

• Hamiltonian: They use a single-channel point interaction model where $g_I$ is complex:
  $$\frac{1}{g_I} \simeq \frac{1}{4\pi a} - \frac{\Lambda}{2\pi^2} + i\Gamma + \gamma \frac{n_0^{1/3}|a|}{4\pi|a|}$$
  • The real part recovers the standard scattering length physics.
  • The imaginary part ($i\Gamma + \dots$) introduces damping, representing the decay of the single-boson dressed state into the continuum of multi-boson states (and other loss channels like three-body recombination).
• Approximation: The model truncates the variational ansatz to include only the impurity and zero or one Bogoliubov mode ($C=D=E=0$ in Eq. 2), justified by the large overlap of these states with experimental probes.

C. Computational Techniques

• Spectral Function: Calculated using the Krylov subspace method to efficiently access the low-energy spectrum of the variational Hamiltonian.
• Dynamics: Time evolution of the contrast $|A(t)|$ (Ramsey signal) is simulated using the time-dependent coefficients of the variational ansatz, incorporating the complex interaction to induce decay.

3. Key Contributions

1. Identification of the "Visible" State: The paper establishes that the experimentally observed Bose polaron in strong interactions is likely a broad resonance (single-boson dressed state) rather than the true ground state, explaining the lack of sharp peaks.
2. Complex Interaction Formalism: It introduces a minimal phenomenological framework using a complex scattering length to effectively model the decay of the polaron into multi-boson correlations, bypassing the need for explicit calculation of all $N \geq 2$ decay channels.
3. Unified Description: The model successfully describes both spectral properties (peak position and width) and non-equilibrium dynamics (decoherence rates) using a consistent set of parameters.

4. Results

The model was tested against two recent high-precision experiments:

• Spectral Properties (Cambridge Experiment, Ref. [20]):

  • The model reproduces the broad spectral linewidths observed near the unitary limit.
  • It correctly captures the asymmetric behavior of the attractive and repulsive branches.
  • The extracted polaron width scales as $n^{2/3}$ at unitarity and linearly with scattering length in the weak regime, matching experimental trends.

• Non-Equilibrium Dynamics (Aarhus Experiment, Ref. [19]):

  • Standard ladder approximations (real $g_I$) predict oscillations between repulsive and attractive branches that are too pronounced compared to experimental data.
  • The phenomenological model with complex $g_I$ ($\Gamma \approx 0.05, \gamma \approx 0.5$) accurately reproduces the rapid damping of these oscillations and the overall decay of the Ramsey contrast $|A(t)|$.
  • The model captures the experimental data within uncertainties, whereas the standard variational ansatz (without complex interaction) fails to match the decay rate.

5. Significance

• Intuitive Framework: The work provides a practical tool for analyzing experimental data in the strong-coupling regime, where full microscopic theories are computationally prohibitive.
• Decay Mechanism: It highlights that the "missing" physics in previous theories is the decay of the single-boson dressed state into lower-lying multi-boson states. This decay is the primary source of the observed spectral broadening and decoherence.
• Future Directions: The paper suggests that future theoretical efforts must focus on the microscopic description of few-body correlations (trimers, etc.) within a many-body environment. It also opens avenues for studying Bose polarons in different geometries (lattices), at finite temperatures, and with mass-imbalanced systems.

In summary, the paper resolves the discrepancy between theory and experiment for Bose polarons in strong interactions by identifying the experimentally relevant state as a decaying resonance and modeling this decay via a complex interaction strength.