Dynamics of one-dimensional Bose-Josephson Junction in a Box Trap: From Coherent Oscillations to Many-Body Dephasing and Dynamical Freezing

Using the multiconfigurational time-dependent Hartree method for bosons, this study reveals how a one-dimensional Bose-Josephson junction in a box trap transitions from coherent oscillations to many-body dephasing and dynamical freezing as interaction strength and initial population imbalance vary, establishing a unified framework for understanding the interplay between coherence and correlation-induced fragmentation.

The Gist

Imagine you have a long, narrow hallway (a "box trap") with a low wall right in the middle. You fill this hallway with 10 tiny, invisible marbles that are actually ultracold atoms. These atoms are special: they are so cold and so alike that they act like a single, giant wave rather than individual particles. This is called a Bose-Einstein Condensate.

The paper you shared is a study of what happens when we start with more marbles on the left side of the wall than the right, and then suddenly knock the wall down (or lower it) to let them mix. The researchers wanted to see how these atoms behave as they try to balance themselves out.

Here is the story of their journey, broken down into three main "moods" or regimes, depending on how much the atoms "hate" or "like" each other (their interaction strength).

The Setup: The "Popcorn" Experiment

Think of the atoms as a bag of popcorn.

• The Box: A long, flat box with a divider in the middle.
• The Imbalance: We put 7 popcorn kernels on the left and 3 on the right.
• The Quench: At time zero, we remove the divider. The popcorn wants to spread out evenly.

The researchers used a super-powerful computer simulation (called MCTDHB) to watch this happen in slow motion, looking at how the "popcorn" moves, spreads, and interacts.

---

Regime 1: The "Polite Dancers" (Weak Interactions)

The Vibe: The atoms barely notice each other. They are like polite dancers who don't bump into one another.

• What happens: If you start with a small imbalance (6 on left, 4 on right), the atoms flow back and forth across the middle like a pendulum. They swing left, then right, then left again. This is called Josephson Oscillation. It's a perfect, rhythmic dance.
• The Twist: If you start with a huge imbalance (9 on left, 1 on right), the dance gets messy. The atoms rush to the right, but because there are so many of them, they get in each other's way. The rhythm slows down, the swings get smaller, and eventually, they just settle into a calm, balanced state (5 on each side).
• The Analogy: Imagine a crowd of people walking through a doorway. If the crowd is small, they flow smoothly back and forth. If the crowd is huge, they bump into each other, slow down, and eventually just stop moving once the room is full.

Regime 2: The "Confused Choir" (Intermediate Interactions)

The Vibe: The atoms are starting to feel each other's presence. They are like a choir where everyone is trying to sing their own song.

• What happens: Here, the atoms interact enough to cause Dephasing.
  • Small Imbalance: They still dance, but the rhythm is a bit wobbly.
  • Medium Imbalance: This is the most interesting part. The atoms start to "collapse and revive." Imagine the choir singing a chord. Suddenly, everyone gets out of sync (the sound collapses into silence). Then, by pure luck, they accidentally get back in sync for a split second (the sound revives). Then they get out of sync again.
  • Large Imbalance: The chaos wins. The atoms rush to the other side, mix up completely, and settle down. They lose their individual rhythm and become a "soup" of particles. They have reached Equilibration.
• The Analogy: Think of a group of runners starting a race at different times. At first, they are spread out. Then, they bunch up, then spread out again, then bunch up. Eventually, they just run in a chaotic pack, and you can't tell who is who anymore. The "order" of the race is lost to "chaos."

Regime 3: The "Frozen Statue Garden" (Strong Interactions)

The Vibe: The atoms absolutely hate being close to each other. They are like magnets with the same pole facing each other—they push away violently.

• What happens: This is the "Dynamical Freezing."
  • Small Imbalance: Even though the wall is gone, the atoms refuse to move. They are so busy pushing away from their neighbors that they get stuck in place. They arrange themselves in a perfect line, like soldiers standing at attention, with equal spacing. They look like a row of individual particles (almost like fermions) rather than a wave. They are frozen.
  • Large Imbalance: They move a tiny bit at first, but then they get stuck again. The system tries to balance, but the "pushing" is so strong that they can't flow freely.
• The Analogy: Imagine a crowded subway car where everyone is wearing a giant, stiff suit of armor. If you try to move, you bump into the person next to you, and they push back. You end up stuck in your spot, unable to flow. The crowd is "frozen" in place, not because they are cold, but because they are too crowded and aggressive to move.

---

Why Does This Matter?

This paper is like a map of how quantum systems behave when they are out of balance.

1. Coherence: In the beginning, everything is orderly and rhythmic (like the polite dancers).
2. Dephasing: As things get more crowded, the rhythm breaks, and things get chaotic (the confused choir).
3. Freezing: If the pressure gets too high, everything stops moving entirely (the frozen statues).

The researchers found that by simply changing how many atoms are on one side versus the other, and how much they push each other, you can switch between these three totally different worlds.

The Big Takeaway:
Nature is full of surprises. Sometimes, when you push a quantum system hard enough, it doesn't just get chaotic; it gets so chaotic that it freezes solid. This helps scientists understand how quantum computers might work (or fail) and how energy moves in the smallest materials in the universe. It shows that "order" and "chaos" are not just opposites; they are two sides of the same coin, depending on how you squeeze the system.

Technical Summary

1. Problem Statement

The paper addresses a central challenge in quantum many-body physics: understanding how coherent quantum dynamics transition into correlation-dominated behavior in low-dimensional systems. Specifically, the authors investigate a one-dimensional Bose–Josephson Junction (BJJ) confined within a box-shaped potential (hard-wall boundaries with a tunable central barrier).

While mean-field theories (like the Gross-Pitaevskii equation) successfully describe weakly interacting regimes, they fail to capture essential phenomena in strongly interacting or highly imbalanced systems, such as:

• Many-body dephasing: The loss of phase coherence due to interactions.
• Fragmentation: The distribution of the condensate over multiple single-particle orbitals.
• Dynamical freezing: The suppression of transport due to strong correlations.
• Ergodicity vs. Non-ergodicity: Whether the system explores the full Hilbert space or remains trapped in a subspace.

The study aims to map the dynamical phase diagram of a 1D BJJ by varying interaction strength ($\Lambda$) and initial population imbalance ($z_0$), moving beyond simple two-mode models to a full many-body treatment.

2. Methodology

The authors employ a rigorous numerical approach to solve the time-dependent many-body Schrödinger equation:

• System Setup: A system of $N=10$ bosons in a 1D box with hard-wall boundaries. The potential consists of a flat-bottom box with a central Gaussian barrier ($V_1$) creating two weakly coupled reservoirs.
• Quench Protocol: The system is initialized in the ground state of a potential with an offset ($V_{imb}$) creating a population imbalance. At $t=0$, the offset is suddenly removed (quantum quench), and the system evolves under the symmetric box potential.
• Numerical Method: The Multiconfigurational Time-Dependent Hartree method for bosons (MCTDHB) is used, implemented via the MCTDH-X software.
  • Unlike mean-field approaches that assume a single orbital, MCTDHB expands the many-body wavefunction in a basis of $M$ time-dependent single-particle orbitals (natural orbitals).
  • This allows for the dynamic tracking of correlations and fragmentation. The number of orbitals ($M$) is increased until convergence is reached (e.g., $M=4$ for weak interactions, up to $M=10$ for strong interactions).
• Observables: The dynamics are characterized by:
  • Population Imbalance ($z(t)$): To track Josephson oscillations and damping.
  • Fourier Spectrum: To identify oscillation frequencies and dephasing.
  • Natural Orbital Occupations ($n_i$): To quantify fragmentation.
  • Entropies: Coefficient entropy ($S_C$) and Orbital entropy ($S_n$) to measure wavefunction spreading and correlation.
  • Participation Ratio (PR): To measure the effective number of many-body configurations involved.

3. Key Contributions

The paper establishes a unified picture of 1D Josephson dynamics, identifying three distinct regimes governed by the interplay of interaction strength and initial imbalance:

1. Weak Interaction Regime: Transition from coherent oscillations to barrier-controlled damping.
2. Intermediate Interaction Regime: A crossover from oscillations to many-body dephasing (collapse-and-revival) and finally to ergodic equilibration, driven by the initial imbalance.
3. Strong Interaction Regime: The emergence of dynamical freezing, where strong repulsion suppresses tunneling, leading to a near-fermionized state.

4. Detailed Results

A. Weak-Interaction Regime ($\Lambda = 0.01$)

• Small Imbalance: The system exhibits long-lived, coherent Josephson oscillations with negligible damping. The condensate remains non-fragmented (single orbital occupation $n_1 \approx 1$).
• Large Imbalance: Oscillations are damped, and the system equilibrates to $z(t) \approx 0$.
  • Mechanism: This equilibration is driven by the barrier height ($V_1$) rather than strong correlations.
  • Fragmentation: A "two-fold fragmentation" occurs where the population splits between two orbitals ($n_1 \approx n_2 \approx 0.5$). However, unlike the strong interaction case, the many-body wavefunction does not fully relax; entropy measures fluctuate, indicating the system does not reach true ergodicity despite the damping of $z(t)$.

B. Intermediate Interaction Regime ($\Lambda = 0.5$)

Here, interactions play a dominant role, and the dynamics are primarily controlled by the initial population imbalance:

• Weak Imbalance: Coherent oscillations persist with minimal fragmentation (mean-field-like behavior).
• Intermediate Imbalance: The system enters a many-body dephasing regime.
  • Collapse-and-Revival: The population imbalance exhibits repeated cycles of amplitude collapse and revival.
  • Mechanism: This is not irreversible decay but a continuous dephasing and rephasing of multiple many-body modes. The system remains far from ergodic, with coherence redistributed among modes.
• Strong Imbalance: The system undergoes ergodic relaxation.
  • Fast Equilibration: $z(t)$ decays rapidly to zero with no revival.
  • Signatures: The Fourier spectrum becomes broad and featureless. Entropies ($S_C, S_n$) and Participation Ratio saturate to values close to the Gaussian Orthogonal Ensemble (GOE) predictions, indicating the system has explored a large portion of the Hilbert space and reached a chaotic, thermalized state.

C. Strong Interaction Regime ($\Lambda = 5.0$)

In this limit, the system enters a dynamical freezing regime, closely related to the Tonks-Girardeau (fermionization) limit:

• Weak Imbalance: The system is effectively frozen. Tunneling is strongly suppressed, and $z(t)$ remains nearly constant.
  • Structure: The one-body density matrix shows ten well-separated lobes (one per particle), indicating strong spatial localization.
  • Dynamics: All observables (entropy, PR, occupations) remain static. The system is non-ergodic and retains its initial configuration.
• Strong Imbalance:
  • Transient Dynamics: Initial tunneling occurs, followed by rapid equilibration of population.
  • Residual Coherence: Unlike the intermediate regime, the system does not fully thermalize. While $z(t)$ equilibrates, the natural orbital occupations and entropies exhibit persistent, weak oscillations (beating) due to coherent superpositions of discrete eigenstates.
  • Fragmentation: The system remains highly fragmented, but the dynamics are compressed in time, merging dephasing and rephasing processes.

5. Significance

• Unified Framework: The work provides a comprehensive map of how coherence, dephasing, equilibration, and freezing compete in 1D bosonic systems, bridging the gap between mean-field physics and strongly correlated many-body regimes.
• Beyond Mean-Field: It demonstrates that initial conditions (imbalance) can drive a system from coherent oscillations to ergodic chaos or to dynamical freezing, depending on interaction strength.
• Experimental Relevance: The predicted regimes are directly accessible in current ultracold atom experiments using optical box traps and Feshbach resonances to tune interactions. The ability to observe "dynamical freezing" and "many-body dephasing" offers new avenues for studying quantum thermalization and non-equilibrium phases of matter.
• Methodological Validation: The study validates MCTDHB as a powerful tool for capturing complex phenomena like fragmentation and dynamical freezing that are invisible to standard two-mode or mean-field models.

In conclusion, the paper reveals that in low-dimensional systems, the path to equilibrium is not unique; it can be a coherent oscillation, a chaotic thermalization, or a complete dynamical arrest, dictated by the delicate balance between interaction strength and initial state preparation.