Effects of intertube dipole-dipole interactions in nearly integrable one-dimensional $^{162}$Dy gases

This study demonstrates that while intertube dipole-dipole interactions in nearly integrable one-dimensional $^{162}$Dy gases slightly alter both equilibrium properties and rapidity measurements, these opposing effects nearly cancel each other out, resulting in measured rapidity distributions that closely match predictions made without considering such interactions.

The Gist

Imagine you have a giant, invisible grid made of invisible strings. On this grid, you place thousands of tiny, super-cold magnets (atoms of Dysprosium). Because they are so cold, they behave like a single, giant wave rather than individual particles.

In a perfect world, if you line these magnets up in long, separate rows (like lanes on a highway), they would only talk to the neighbors in their own lane. Physicists call these "nearly integrable" systems because they are so orderly that they don't easily forget their initial state or "thermalize" (reach a messy, random equilibrium).

However, in the real world, these magnets are also long-range magnets. Even though they are in separate lanes, a magnet in Lane 1 can still feel a tiny tug from a magnet in Lane 2. This is called inter-tube dipole-dipole interaction.

For a long time, scientists studying these systems ignored the "tugs" between lanes because they thought they were too weak to matter. They assumed the lanes were completely independent. But recently, experiments showed that the real-world data didn't quite match the perfect theoretical predictions.

The Big Question:
Could those tiny "tugs" between the lanes be the reason the math didn't match the experiment?

The Story of the Paper:
This paper is a detective story where the authors (Yicheng Zhang and colleagues) decided to put those "tugs" back into the math to see if they could fix the mismatch.

1. The Setup: The "Magic Angle"

Imagine you are trying to arrange these magnetic atoms. If you tilt your magnetic field at a very specific angle (about 55 degrees), the magnets in the same lane stop pushing or pulling each other. It's like turning off the engine in every car in a single lane.

But here's the catch: even if the cars in the same lane are quiet, the cars in neighboring lanes are still whispering to each other. The previous experiments ignored these whispers. The authors of this paper decided to listen to them.

2. The Two-Step Dance

The experiment happens in two main phases, and the "whispers" affect them in opposite ways:

• Phase A: Getting Ready (State Preparation)
  Imagine you are packing these atoms into their lanes. The "whispers" between lanes act like a weak, invisible repulsive force. It's as if the lanes are slightly pushing each other apart.

  • The Effect: This push makes the atoms spread out a little more and actually makes them slightly cooler (more organized). It's like a gentle breeze that clears the air.
  • Result: The "rapidity distribution" (a fancy way of measuring how fast and in what direction the atoms are moving) becomes slightly narrower and more focused.

• Phase B: The Release (Expansion)
  Next, the scientists turn off the walls holding the atoms in their lanes and let them zoom out to be measured.

  • The Effect: Now, that same "whisper" (the repulsive force between lanes) acts like a spring. As the atoms zoom out, the repulsion between lanes gives them a little extra kick, pushing them apart faster.
  • Result: This extra kick makes the "rapidity distribution" slightly wider and more spread out.

3. The Plot Twist: The Perfect Cancellation

Here is the most surprising part of the story.

When the authors did the math including the "whispers," they found that the two effects almost perfectly canceled each other out!

• The "whispers" made the atoms more focused during the setup.
• But then, the "whispers" made them spread out more during the release.

The Final Result: When you look at the final measurement, the "whispers" between lanes barely changed anything at all. The final picture looked almost exactly the same as if the lanes had been completely silent and independent.

The Conclusion: Why the Mismatch?

So, if the "whispers" didn't cause the difference between the theory and the experiment, what did?

The authors conclude that the "whispers" (inter-tube interactions) are not the culprit. The mismatch must be caused by something else entirely. They suspect it's due to the fact that these systems are "nearly integrable."

Think of it like this: In a perfectly ordered system, the atoms remember everything perfectly. In a "nearly" ordered system, they remember almost everything, but with a few tiny glitches. The authors suggest that the real-world experiment is suffering from these "glitches" (non-thermal effects) that their current models aren't sophisticated enough to capture yet.

Summary in a Nutshell

• The Problem: Theory and experiment didn't match. Scientists thought maybe the atoms were talking to neighbors in other lanes.
• The Investigation: They added the "neighbor talk" to their computer models.
• The Discovery: The "neighbor talk" had two opposite effects that canceled each other out like a tug-of-war where both sides pull with equal strength.
• The Verdict: The "neighbor talk" isn't the problem. The problem is likely that our current understanding of how these "nearly perfect" quantum systems behave is still missing a few pieces of the puzzle.

It's a great example of how science works: you test a hypothesis (maybe the neighbors are talking?), do the hard math, and even if the answer is "no," you learn something valuable about what isn't causing the problem, narrowing down the search for the real answer.

Technical Summary

1. Problem Statement

The paper addresses a discrepancy between theoretical models and experimental observations in recent studies of nearly integrable one-dimensional (1D) dipolar Bose gases composed of Dysprosium-162 ($^{162}\text{Dy}$) atoms.

• Context: In previous work (Ref. [43]), an array of 1D tubes was created using a 2D optical lattice. Theoretical modeling successfully described most experimental features by accounting for intratube dipole-dipole interactions (DDI) via a modification of the contact interaction, while explicitly ignoring intertube DDI (interactions between atoms in different tubes).
• The Issue: Despite this simplification, the theoretical predictions for the rapidity distributions (which characterize the equilibrium states of integrable systems) did not perfectly match experimental data.
• The Question: Could the neglected intertube DDI be the source of this disagreement? Specifically, how do these long-range interactions affect the state preparation (loading into the lattice) and the subsequent expansion dynamics used to measure rapidity?

2. Methodology

The authors developed a self-consistent theoretical framework to incorporate intertube DDI into the modeling of the experiment described in Ref. [43].

• Modeling Approach:
  • Mean-Field Approximation: The leading effect of the intertube DDI is treated as a modification to the 1D confining potentials ($V_{1D}$). This potential is calculated as a sum of mean-field contributions from all other tubes in the 2D array.
  • Iterative Self-Consistency:
    1. State Preparation: Starting from a 3D Bose-Einstein Condensate (BEC), the system is loaded into a 2D optical lattice. The authors calculate the atom number ($N_\ell$) and temperature ($T_\ell$) for each tube $\ell$ at the "decoupling point" (where tubes become independent). They iteratively update the density distributions to account for the intertube DDI potential until convergence.
    2. Adiabatic Evolution: The model assumes adiabatic changes (constant entropy) as experimental parameters (lattice depth, trap frequency, magnetic field angle) are adjusted to reach the final equilibrium state.
    3. Expansion Dynamics: To simulate the measurement, the authors use Generalized Hydrodynamics (GHD) to model the expansion of the 1D gases after the 1D trap is turned off. The intertube DDI acts as a time-dependent, weakening effective potential during this expansion.
• Key Parameters:
  • Magnetic field angle ($\theta_B$): Varied to tune intratube DDI and contact interactions.
  • Contact interaction strength ($g_{1D}^{vdW}$): Tuned via Feshbach resonance.
  • Specific Focus: The study focuses on the case where $\theta_B = 55^\circ$, a "magic angle" where intratube DDI vanishes, isolating the effect of intertube DDI.

3. Key Contributions

• First Systematic Inclusion of Intertube DDI: This work provides the first rigorous theoretical treatment of intertube dipole-dipole interactions in the context of 1D dipolar gas arrays, moving beyond the approximation of independent tubes.
• Decoupling of Effects: The authors successfully separated the effects of intertube DDI into two distinct phases:
  1. Equilibrium State Preparation: How intertube DDI alters the initial density and temperature profiles.
  2. Expansion Dynamics: How intertube DDI alters the evolution of the rapidity distribution during the time-of-flight expansion.
• Refinement of Numerical Methods: The authors removed the simplification of rounding atom numbers to integers (used in Ref. [43]) and excluded tubes with very low occupancy (<2 atoms) to improve the accuracy of the density and rapidity calculations.

4. Key Results

The study reveals a remarkable cancellation effect that explains why previous models (ignoring intertube DDI) were surprisingly close to experimental results.

• Effect on Equilibrium State (State Preparation):

  • The intertube DDI acts as a weak antitrapping potential in the transverse ($y-z$) and longitudinal ($x$) directions.
  • Cooling Effect: This antitrapping potential effectively lowers the temperature of the 1D gases, particularly in tubes with lower atomic filling (those further from the center of the array).
  • Narrowing: The combination of lower temperature and the antitrapping potential results in a narrowing of both the spatial density distribution and the initial rapidity distribution compared to the model without intertube DDI.

• Effect on Expansion Dynamics (Measurement):

  • During the expansion (when the 1D trap is turned off), the intertube DDI energy is converted into kinetic energy.
  • This conversion accelerates the quasiparticles, causing the rapidity distribution to broaden.
  • This broadening effect is the opposite of the narrowing observed during state preparation.

• The Cancellation Phenomenon:

  • Crucially, the narrowing of the rapidity distribution caused by intertube DDI during state preparation is almost perfectly canceled by the broadening caused by the same interaction during the expansion dynamics.
  • Result: The final measured rapidity distribution (after expansion) with intertube DDI included is virtually indistinguishable from the distribution predicted by the model that ignores intertube DDI.

• Comparison with Experiment:

  • Since including intertube DDI does not significantly change the theoretical prediction, the remaining discrepancy between theory and experiment (observed in Ref. [43]) is not caused by intertube DDI.
  • The authors conclude that the discrepancy likely arises from non-thermal effects due to the near-integrability of the 1D gases (e.g., failure of the system to reach thermal equilibrium during the decoupling or adiabatic processes).

5. Significance

• Validation of Previous Models: The work validates the robustness of the "independent tube" approximation used in previous studies, demonstrating that intertube DDI is not a dominant source of error in these specific experimental configurations.
• Identification of New Physics: By ruling out intertube DDI, the paper shifts the focus to more subtle physics. It suggests that the disagreement between theory and experiment stems from the breakdown of thermalization assumptions in nearly integrable systems, specifically regarding how tubes decouple and evolve.
• Methodological Advance: The study provides a blueprint for modeling complex, multi-component dipolar systems where long-range interactions between different subsystems (tubes) must be accounted for self-consistently.
• Future Directions: The authors propose that future models should account for a distribution of decoupling depths (rather than a single depth) and explore non-thermal descriptions of the adiabatic processes following decoupling to resolve the remaining experimental discrepancies.