Interference of a chain of Bose condensates in the Pitaevskii-Gross approximation

This paper utilizes the Pitaevskii-Gross equation to theoretically model the interference and spatial density spectrum of a freely expanding chain of Bose condensates, successfully reproducing the experimental peak positions and their dependence on interparticle interactions while noting discrepancies in peak heights.

The Gist

The Big Picture: A Dance of Quantum Waves

Imagine you have a long line of dancers (these are Bose-Einstein Condensates, or BECs). In a special laboratory, these dancers are trapped in a grid of invisible "light cages" (an optical lattice). Each dancer is standing perfectly still in their own cage, but they are all holding hands with the same invisible rhythm.

The scientists in this paper wanted to see what happens when they suddenly turn off the lights and let the dancers run free into an open field.

The Experiment: The "Talbot" Magic Trick

When the lights go out, the dancers start to run outward. Because they are quantum particles, they don't just act like little balls; they act like ripples in a pond. As these ripples spread out, they crash into each other and create a pattern of high and low points, known as interference fringes.

There is a magical moment called the Talbot time.

• The Analogy: Imagine you have a row of flashlights shining through a picket fence. If you stand at a specific distance away, the shadow of the fence looks exactly like the fence itself, even though the light has traveled a long way. This is called "self-imaging."
• In the Lab: When the dancers run for exactly this "Talbot time," they naturally re-form their original line-up. It's like a magic trick where the scattered dancers suddenly snap back into their original formation.

The Twist: What if the Dancers are Out of Sync?

In a perfect world, all dancers would start with the exact same rhythm (phase). But in the real world, things get messy. The scientists introduced randomness. Some dancers started slightly ahead, some slightly behind.

• The Result: If you look at the pattern of ripples in a single experiment, it looks chaotic. The lines are wobbly and in different places every time you run the experiment.
• The Surprise: Even though the picture looks different every time, the math behind it is surprisingly stable. When the scientists analyzed the "spectrum" (which is like looking at the frequency of the waves rather than the picture itself), they found a hidden order.

The Two Types of Peaks

The scientists discovered that the "spectrum" of these ripples has two distinct types of peaks (high points on a graph), which act like two different voices in a choir:

1. The "Order" Voice (Coherence Peaks): These are sharp, narrow spikes. They appear because the dancers are still somewhat connected. Even if they are a little out of sync, they remember they were once a line. This is the "Talbot effect" voice.
2. The "Chaos" Voice (Fluctuation Peaks): These are wider, fuzzier bumps. They appear because of the random jitters and phase differences. This voice tells the scientists how "messy" the starting conditions were.

By listening to the balance between these two voices, the scientists can measure exactly how much the dancers were out of sync.

The Math: The "Gross-Pitaevskii" Recipe

To predict this behavior, the authors used a complex recipe called the Gross-Pitaevskii equation.

• The Analogy: Think of this equation as a super-accurate weather forecast for quantum particles. It doesn't just say "it will rain"; it calculates exactly how the wind (interactions between particles) pushes the clouds (the wave functions) around.
• The Interaction: The paper notes that the particles bump into each other (interact). This is like if the dancers were wearing bulky coats; when they run, they push against each other, changing the shape of the ripples. The scientists' math accounted for this "bumping," and it successfully predicted where the peaks in the spectrum would shift.

The Conclusion: A Good Match, With a Few Glitches

The scientists compared their computer simulations (the math recipe) with the actual photos taken in the lab.

• The Win: The math perfectly predicted where the peaks would appear and how much they would shift due to the particles bumping into each other.
• The Glitch: The math sometimes got the height of the peaks wrong. It's like predicting exactly where a wave will crash on the beach, but guessing the wave is 10% taller or shorter than it actually is. The authors suspect this is because their model is a bit too simple and doesn't account for every tiny detail of how the particles lose energy.

Why Does This Matter?

This research is a powerful tool for thermometry (measuring temperature).

• The Analogy: Imagine trying to measure the temperature of a pot of soup without a thermometer, just by watching how the steam rises.
• The Application: By looking at how "messy" the interference pattern is (the ratio of the two types of peaks), scientists can calculate the temperature of the quantum gas with incredible precision, even at temperatures so cold that normal thermometers don't work.

In short: The paper shows that even when a quantum system gets messy and chaotic, there is still a hidden mathematical order. By understanding the "noise," we can learn exactly how the system is behaving.

Technical Summary

1. Problem Statement

The paper addresses the dynamics of a long chain of Bose-Einstein condensates (BECs) released from an optical lattice. Specifically, it investigates the interference patterns formed during free expansion and the resulting spatial density spectrum under two distinct conditions:

1. Coherent Phase: All condensates share identical phases.
2. Fluctuating Phase: Adjacent condensates possess random phase differences (partial or complete incoherence).

While the Talbot effect (self-imaging of the wave function) is well understood for coherent sources, the behavior of the interference spectrum when phases fluctuate is less intuitive. In classical optics, fluctuating phases typically destroy spatial order. However, in this quantum system, a reproducible spatial spectrum persists despite fringe position fluctuations. The authors aim to model this phenomenon, specifically analyzing how interparticle interactions (mean-field effects) influence the spectrum and comparing theoretical predictions with recent experimental data [Phys. Rev. Lett. 122, 090403 (2019)].

2. Methodology

The authors developed a theoretical model based on the Pitaevskii-Gross (Gross-Pitaevskii) equation, which describes the dynamics of the condensate wave function $\Psi(\mathbf{r}, t)$ including interparticle interactions.

Key Modeling Steps:

• Initial Conditions: The system is initialized as a chain of $K$ condensates trapped in an optical lattice with period $d$. The initial wave function is a product of a longitudinal part $\psi(z,0)$ and a radial part $\chi(\rho,0)$.
  • Longitudinal: Modeled as a sum of Gaussian wave packets, each with an arbitrary phase $\phi_j$. The phase differences between neighbors are treated as normally distributed random variables, characterized by a coherence factor $\alpha = \langle \cos(\phi_j - \phi_{j+1}) \rangle$.
  • Radial: Modeled using the Thomas-Fermi approximation due to strong anisotropy.
• Self-Consistency: The axial width $\sigma$ of the condensates is determined self-consistently by minimizing the energy functional, accounting for the expansion caused by repulsive interactions. This results in $\sigma$ being larger than the harmonic oscillator length ($1.3\text{--}1.5 \times l_z$).
• Dimensional Reduction: Due to strong anisotropy (radial expansion is much slower than axial), the radial expansion is neglected. The 3D Gross-Pitaevskii equation is averaged over the radial coordinate to derive a 1D effective equation for $\psi(z,t)$.
• Simulation: The 1D equation is solved numerically to evolve the system to the Talbot time $T_d = md^2/(\pi\hbar)$.
• Analysis: The column density $n_2(x,z)$ is calculated, and its spatial spectrum $\tilde{n}_1(k)$ is obtained via Fourier transform. The results are averaged over multiple realizations of random phases to match experimental averaging.

3. Key Contributions

• Identification of Dual Spectral Peaks: The model successfully distinguishes two types of peaks in the spatial density spectrum:
  1. Narrow Peaks: Associated with the Talbot effect (coherence between condensates). Their positions are determined by the lattice period $d$ and are unaffected by mean-field interactions.
  2. Wide Peaks: Arising from phase fluctuations. Their positions are shifted by interparticle interactions (mean-field effects) to lower momenta.
• Quantitative Reproduction of Interaction Shifts: The study demonstrates that the Pitaevskii-Gross approximation accurately reproduces the experimentally observed shift in the wide spectral peaks caused by interatomic interactions. This confirms that the broadening of the condensate (due to interactions) is a critical factor in the quantitative agreement.
• Thermometry via Coherence Factor: The paper establishes a method to determine the coherence factor $\alpha$ (and thus the temperature) by analyzing the relative heights of the narrow (Talbot) and wide (fluctuation) peaks in the spectrum. This allows for thermometry at temperatures significantly below the critical temperature where conventional bimodal fitting fails.

4. Results

• Spectral Reproducibility: Even when the interference fringes fluctuate in position and intensity due to random phases, the spectrum of the spatial density distribution remains reproducible and exhibits equidistant peaks.
• Phase Dependence:
  • For $\alpha = 1$ (perfect coherence), the spectrum shows narrow peaks at $k = 2\pi l/d$.
  • For $\alpha = 0$ (complete disorder), the spectrum shows wider peaks at $k = \pi l T_d / t$.
  • For intermediate $\alpha$, both peak types coexist.
• Interaction Effects:
  • The model correctly predicts that interparticle interactions shift the wide peaks to lower momenta but leave the narrow Talbot peaks stationary.
  • Without accounting for the interaction-induced broadening of the condensate (using the self-consistent $\sigma$), the predicted shift is 2–3 times smaller than observed.
• Comparison with Experiment:
  • The calculated peak positions and widths show excellent quantitative agreement with experimental data [17].
  • Discrepancy: In some cases (specifically at lower lattice depths $s \leq 18.4$), the calculated peak heights exceed experimental values. The authors attribute this to the Pitaevskii-Gross approximation neglecting condensate depletion and the deviation of the initial wave function from a Gaussian profile at lower lattice depths.

5. Significance

This work provides a robust theoretical framework for understanding the interference of matter waves in the presence of both phase disorder and interparticle interactions. It bridges the gap between the idealized Talbot effect and realistic experimental conditions involving thermal fluctuations and interactions.

The ability to extract the coherence factor $\alpha$ from the spectral shape offers a powerful diagnostic tool for quantum gases, enabling precise temperature measurements in regimes where standard methods fail. Furthermore, the successful modeling of interaction-induced spectral shifts validates the use of the Gross-Pitaevskii equation for describing collective interference phenomena in quasi-2D Bose gases, despite the noted limitations regarding peak amplitudes in specific parameter regimes.