Anomalous fluctuations of Bose-Einstein condensates in optical lattices

Through a combined experimental and theoretical study of ultracold $^{87}$Rb atoms in optical lattices, this paper reveals strongly anomalous condensate particle number fluctuations that scale super-linearly with the total atom number, a phenomenon attributed to the interplay of 2D/3D crossover geometry and interactions.

The Gist

The Big Picture: A Dance of Atoms

Imagine you have a massive ballroom filled with thousands of dancers (atoms). Usually, when the music is loud and fast (high temperature), everyone is dancing wildly, bumping into each other, and moving in random directions. This is a normal gas.

But if you turn the music down to a slow, hypnotic beat and cool the room down to near absolute zero, something magical happens. Suddenly, all the dancers stop their chaotic moves and start moving in perfect unison. They become a single, giant "super-dancer." In physics, this is called a Bose-Einstein Condensate (BEC).

This paper is about studying how "jittery" or "wobbly" this super-dancer is when you put them in a very specific, tricky room: an optical lattice.

The Setting: The "Tube City"

Usually, scientists study these super-dancers in a smooth, open ballroom (a continuous system). But in this experiment, the researchers built a city of tubes using lasers.

• The Analogy: Imagine the ballroom floor is covered with a honeycomb grid of invisible, vertical tubes. The atoms are trapped inside these tubes.
• The Twist: The tubes are arranged in a 2D flat sheet (like a pancake), but the atoms can still wiggle a little bit up and down (3D). It's a "crossover" between a flat world and a 3D world.

The Mystery: The "Wobble" (Fluctuations)

In physics, fluctuations are just tiny, random changes in numbers. If you count how many atoms are in the "super-dancer" state versus the "chaotic" state, the number isn't perfectly fixed; it jitters.

For decades, physicists have argued about how much this number should jitter:

1. The "Grand Canyon" Theory: Some theories said the jitter should be huge (as big as the total number of atoms).
2. The "Gaussian" Theory: Others said the jitter should be small and predictable, like rolling dice.

What the researchers found:
They discovered that in their "Tube City," the jitter is wildly abnormal. It's not huge like the Grand Canyon, but it's much bigger and stranger than the simple dice roll.

They found that as they added more atoms (making the dance floor bigger), the "wobble" didn't grow in a normal way. Instead, it grew in a strange, super-linear fashion.

• The Math Metaphor: If you double the number of dancers, the wobble doesn't just double; it explodes much faster. They found a specific "magic number" (an exponent) of about 0.62 to 0.74 that describes this weird growth.

How They Did It: The "Freeze-Frame" Camera

To see this, they used a Matter-Wave Microscope.

• The Problem: Atoms move too fast to count them individually in real-time.
• The Solution: They used a laser trick to "freeze" the atoms in place (like taking a high-speed photo). Then, they used a special lens to magnify the image so they could see the density of atoms in every single tube.
• The Result: They took thousands of these "photos" at different temperatures and counted how many atoms were in the "super-dance" vs. the "chaos" for every single shot.

The Theory: The "Master Plan"

On the computer side, the team built a complex simulation. They didn't just use one rulebook; they combined two different ways of thinking:

1. The Wave Theory: Treating the atoms like ripples in a pond (Bogoliubov theory).
2. The Traffic Light Theory: Treating the atoms like cars stopping and starting at lights (Master Equation).

By mixing these two, they could predict the "wobble" perfectly. Their computer model predicted a magic number of 0.74, and their real-world experiment measured 0.62. These numbers are very close, proving their theory works.

Why Does This Matter?

You might ask, "Who cares about atoms jittering in a tube?"

1. Understanding the Rules of Reality: It helps us understand how the universe behaves when things get very small and very cold. It shows that the shape of the room (the geometry) changes the fundamental rules of how matter behaves.
2. Better Sensors: These "wobbly" atoms are incredibly sensitive. If we can control these fluctuations, we can build super-precise sensors (like atomic clocks or gravity detectors) that are far more accurate than anything we have today.
3. Quantum Computers: Understanding how these atoms interact and fluctuate is a stepping stone toward building quantum computers, which rely on controlling these exact kinds of quantum weirdness.

The Takeaway

The researchers took a bunch of ultra-cold atoms, trapped them in a laser-made honeycomb of tubes, and watched how they danced. They found that the atoms didn't just jitter randomly; they jittered in a strange, powerful, and predictable pattern that depends on the shape of the room they are in.

It's like discovering that in a specific type of room, a crowd of people doesn't just shuffle around randomly, but they all start swaying in a giant, rhythmic wave that gets stronger the more people you add. This discovery helps us understand the deep, hidden rules of the quantum world.

Technical Summary

1. Problem Statement and Motivation

Fluctuations in physical systems are fundamental to understanding phase transitions and the thermodynamic properties of matter. While the fluctuations of the condensate particle number ($N_{BEC}$) in continuous atomic Bose-Einstein condensates (BECs) have been extensively studied, there is a significant gap in knowledge regarding lattice systems.

• Theoretical Context: In continuous systems, the scaling of condensate fluctuations ($\delta N_{BEC}^2$) with total atom number ($N$) depends heavily on dimensionality ($d$) and the dispersion relation ($\epsilon \propto k^\sigma$).
  • For ideal gases, "normal" fluctuations scale as $\delta N^2 \propto N$.
  • "Anomalous" fluctuations scale as $\delta N^2 \propto N^{1+\gamma}$.
  • Previous theoretical predictions for 3D harmonic traps suggested $\gamma \approx 1/3$, while 2D systems suggest $\gamma = 1$.
• The Gap: Experimental and theoretical studies for BECs in optical lattices were missing. Specifically, the behavior of fluctuations in a 2D/3D crossover geometry (a triangular array of tubes confined by a weak harmonic trap) had not been investigated.
• Goal: To experimentally and theoretically characterize the scaling of condensate number fluctuations across the phase transition in an optical lattice BEC and determine the anomalous scaling exponent $\gamma$.

2. Methodology

The study employs a combined experimental and theoretical approach using ultracold $^{87}\text{Rb}$ atoms.

A. Experimental Setup

• System: $^{87}\text{Rb}$ atoms prepared in a magnetic trap ($F=2, m_F=2$) and loaded into a triangular optical lattice formed by 1064 nm laser beams.
• Geometry: The system consists of a triangular lattice of tubes within a harmonic confinement. The trapping frequency in the lattice plane ($\omega_{xy}$) is much higher than in the perpendicular direction ($\omega_z$), creating a 2D/3D crossover regime.
  • Lattice depth: $1 E_{rec}$ (recoil energy).
  • Tunneling energy: $J = h \times 12$ Hz.
• Imaging: A matter-wave microscope is used for high-resolution, single-tube resolution imaging.
  • Protocol: The lattice is ramped up to freeze density, followed by a matter-wave expansion pulse (44x magnification) and absorption imaging.
• Data Extraction: A semi-ideal bimodal model is fitted to the density profiles to extract the condensate fraction ($N_{BEC}$), thermal fraction, and temperature ($T$) for each experimental shot.
• Control: Temperatures are varied by changing the end-point of the radio-frequency evaporation ramp.

B. Theoretical Framework

The authors use a hybrid theoretical approach to model the system, as standard Bogoliubov theory breaks down near the critical temperature ($T_c$) for fixed particle numbers.

1. Hamiltonian: The system is described by a Bose-Hubbard Hamiltonian including tunneling ($J$), on-site interaction ($U$), and a harmonic trapping potential.
2. Bogoliubov-De Gennes (BdG) Equations: The field operators are decomposed into a classical condensate field ($\Phi$) and quantum fluctuations ($\tilde{\psi}$). The authors solve the coupled Gross-Pitaevskii equation (GPE) and the Bogoliubov excitation equations self-consistently.
3. Master Equation Approach: To handle the fixed particle number constraint and thermalization correctly (avoiding the "grand-canonical catastrophe"), the non-condensed particles are treated as a reservoir. A master equation is derived for the condensate particle number, incorporating thermal corrections based on the statistical properties of Bogoliubov quasiparticles.
4. Simulation: Numerical simulations are performed for various particle numbers ($N$) and interaction strengths ($U/J$) to extract the scaling exponent.

3. Key Contributions

1. First Experimental Observation in Lattices: This work provides the first experimental measurement of condensate number fluctuations in an optical lattice system, filling a critical gap between continuous gas studies and lattice theory.
2. Hybrid Theoretical Method: The authors successfully combine the Bogoliubov quasiparticle framework with a master equation analysis. This allows for accurate modeling of fluctuations across the phase transition ($T \approx T_c$) where standard approximations fail.
3. Quantification of 2D/3D Crossover: The study explicitly characterizes the fluctuation behavior in a geometry that bridges 2D and 3D physics, demonstrating how trap geometry dictates the nature of quantum fluctuations.

4. Results

A. Condensate Fraction

• The experimental condensate fraction vs. reduced temperature ($T/T_c^0$) matches well with numerical simulations.
• Increasing the interaction strength ($U$) shifts the critical temperature to lower values, consistent with theoretical expectations.

B. Fluctuation Scaling (The Core Finding)

The variance of the condensate number scales with the total atom number as:
$$\delta N_{BEC}^2 \propto N^{1+\gamma}$$

• Theoretical Result: Numerical simulations yield an anomalous scaling exponent of $\gamma_{theo} = 0.74$. This value is largely independent of temperature in the absence of interactions and remains stable for interacting systems at sufficiently high temperatures.
• Experimental Result: Analysis of experimental data yields $\gamma_{exp} = 0.62$.
• Agreement: The experimental and theoretical values match within a ~20% margin, confirming the presence of strongly anomalous fluctuations.
• Comparison to Limits:
  • The observed $\gamma$ (0.62–0.74) lies between the expected values for a purely 2D lattice ($\gamma = 1$) and a purely 3D lattice ($\gamma = 1/3$).
  • This confirms that the specific 2D/3D crossover geometry (triangular tubes in a harmonic trap) drives the system into a unique scaling regime.

C. Temperature Dependence

• The fluctuations exhibit a pronounced peak below the critical temperature ($T \approx 0.65 T_c^0$).
• Increasing interaction strength ($U$) causes the fluctuation peak to shift to lower temperatures and slightly increases the scaling exponent (e.g., $\gamma \approx 0.80$ for higher $U$ in simulations), attributed to interaction-induced depletion and finite-size effects.

5. Significance and Impact

• Fundamental Physics: The study resolves the controversy regarding fluctuation scaling in interacting systems by providing clear experimental evidence of anomalous scaling in a controlled lattice environment. It highlights that the trap geometry is a decisive factor in determining the universality class of fluctuations.
• Quantum Metrology: Understanding and controlling condensate fluctuations is crucial for quantum metrology. Large fluctuations can be detrimental, but in specific regimes, they can be harnessed to generate entangled atom pairs for atom interferometers, improving measurement precision beyond the standard quantum limit.
• Future Directions: The paper opens avenues for studying full counting statistics, local fluctuations, and the effects of dipolar interactions or higher lattice bands. It also suggests that similar hybrid methods could be applied to other open quantum systems.

In summary, Jalali-Mola et al. have successfully demonstrated that Bose-Einstein condensates in optical lattices exhibit strongly anomalous number fluctuations, with a scaling exponent intermediate between 2D and 3D limits, driven by the specific geometry of the trap and inter-particle interactions.