Interaction Quench Dynamics and Stability of Quantum Vortices in Rotating Bose-Einstein Condensates

Using an exact quantum many-body approach, this study reveals that interaction quenches in rotating two-dimensional Bose-Einstein condensates induce distinct dynamical regimes ranging from complete vortex revival to chaotic fragmentation, depending on the initial vortex configuration and the interplay between interaction strength and angular velocity.

The Gist

Imagine a giant, invisible dance floor made of super-cold atoms. This is a Bose-Einstein Condensate (BEC), a state of matter where thousands of atoms stop acting like individual particles and start moving as a single, synchronized team.

Now, imagine you spin this dance floor. Because of the laws of quantum physics, the atoms can't just spin smoothly like a solid disk. Instead, they have to form tiny, swirling whirlpools called quantum vortices. Think of these like little tornadoes or whirlpools in a bathtub, but made of pure energy and matter.

This paper is about what happens when you suddenly change the "stickiness" (interaction strength) of these atoms while they are spinning. The scientists used a super-powerful computer simulation to watch this happen, because doing it in a real lab is incredibly difficult.

Here is the breakdown of their discovery, using simple analogies:

1. The Old Way vs. The New Way

For a long time, scientists used a "Mean-Field" theory to predict how these atoms behave.

• The Old Way (Mean-Field): Imagine predicting the behavior of a crowd by looking at the average person. You assume everyone is identical and follows the same simple rules. It's like saying, "If the crowd spins, they will all just form a perfect circle."
• The New Way (Many-Body): The authors used a method called MCTDHB. This is like looking at every single person in the crowd individually and seeing how they talk to, push, and pull on their neighbors. It captures the messy, complex reality where atoms have their own "personalities" and relationships.

The Result: The old way often got it wrong. It predicted that the atoms would stay in a perfect, calm state. The new way showed that the atoms were actually chaotic, messy, and full of surprises.

2. The "Quench" (The Sudden Change)

The researchers set up the atoms in a spinning state with a specific number of vortices (whirlpools). Then, they performed a "Quench."

• The Analogy: Imagine you are driving a car on a highway with the cruise control set to a high speed. Suddenly, you slam the brakes and the engine cuts out. The car doesn't just stop instantly; it skids, swerves, and the passengers lurch forward.
• In the Lab: They suddenly made the atoms much less "sticky" (reduced the interaction strength). This caused the delicate balance of the spinning vortices to break, sending the system into a chaotic, non-equilibrium state.

3. What Happened Next? (The Four Scenarios)

The scientists watched what happened to the vortices after the "brakes" were slammed. The outcome depended entirely on how many vortices they started with:

• Scenario A: The Single Vortex (The Breathing Giant)

  • What happened: One central vortex started to expand and contract rhythmically, like a giant lung breathing in and out.
  • The Analogy: It's like a single balloon inflating and deflating perfectly in time with a heartbeat. The whole cloud of atoms pulsed in sync.

• Scenario B: Two Vortices (The Dance of Two)

  • What happened: The two vortices got distorted and the surrounding cloud of atoms split into two pieces. These pieces spun in the opposite direction of the original spin!
  • The Analogy: Imagine two ice skaters holding hands and spinning. Suddenly, they let go. Instead of flying apart, they start spinning in the opposite direction, and the crowd around them splits into two groups chasing them. They eventually come back together, but not quite perfectly.

• Scenario C: Three Vortices (The Triangle Breakup)

  • What happened: The three vortices formed a triangle. When the change happened, the cloud split into three pieces. They spun wildly, merged, and then the vortices "revived" (reappeared) in a new shape.
  • The Analogy: Think of a three-legged stool. If you pull the legs apart, the seat breaks into three pieces. But then, magically, the pieces slide back together to reform the stool, though it's slightly twisted.

• Scenario D: Many Vortices (The Chaos)

  • What happened: When they had a lot of vortices (like eight), the system went completely crazy. The vortices distorted, the cloud split into eight pieces, and they never settled back down.
  • The Analogy: This is like a bowl of spaghetti that you stir too fast. The noodles tangle, break, and fly everywhere. There is no rhythm, no pattern, just pure chaos. The "revival" never happens; the system stays messy forever.

4. The Big Takeaway

The most important lesson from this paper is that simplifying the world doesn't always work.

When you have a few atoms, they behave somewhat predictably. But when you have many atoms interacting strongly, the "average" view (Mean-Field) fails completely. The atoms create complex patterns, split apart, and behave in ways that only a detailed, "many-body" view can predict.

Why does this matter?
Understanding how these quantum whirlpools behave helps us build better quantum computers and simulate complex materials. It teaches us that in the quantum world, the whole is often much more chaotic and interesting than the sum of its parts.

In short: The scientists spun a cloud of atoms, suddenly changed the rules, and watched the vortices either breathe, dance, or spiral into chaos, proving that the quantum world is far more complex than our simple models suggest.

Technical Summary

1. Problem Statement

The paper addresses the non-equilibrium dynamics of quantum vortices in two-dimensional (2D) rotating Bose-Einstein Condensates (BECs) following a sudden change (quench) in interaction strength. While the ground-state properties of vortices are well-understood within the Mean-Field Gross-Pitaevskii (GP) framework for weakly interacting systems, the GP equation fails to capture essential physics in regimes of strong interactions and ultra-fast rotation. Specifically, mean-field theory neglects quantum correlations and the phenomenon of fragmentation (where the many-body wavefunction spreads across multiple single-particle orbitals rather than occupying a single condensate state).

The authors aim to investigate:

• How vortex structures form and stabilize under the interplay of rotation and strong interactions.
• The stability and dynamical evolution of these vortices when the interaction strength is suddenly reduced (an interaction quench).
• The limitations of mean-field theory compared to exact many-body solutions in describing these out-of-equilibrium processes.

2. Methodology

The study employs a first-principles (ab initio) quantum many-body approach using the Multiconfigurational Time-Dependent Hartree method for Bosons (MCTDHB), implemented via the MCTDH-X software package.

• System Model: A 2D rotating BEC of $N=100$ atoms confined in a hard-wall disk potential. The atoms interact via short-range contact potentials modeled as narrow Gaussians.
• Hamiltonian: The system is described by the time-dependent Schrödinger equation in a rotating frame, including kinetic energy, the Coriolis term (rotation), the trapping potential, and pairwise interactions.
• Wavefunction Expansion: The many-body wavefunction is expanded as an adaptive superposition of time-dependent permanents constructed from $M$ self-consistent orbitals.
  • For $M=1$, the method reduces to the standard Mean-Field (GP) approximation.
  • For $M=4$ (used in this study), the method captures beyond-mean-field effects and quantum correlations.
  • The number of orbitals ($M=4$) was chosen to ensure convergence of relevant observables while managing computational cost.
• Protocol:
  1. Preparation: The system is initialized in the ground state for specific combinations of interaction strength ($g$) and rotation frequency ($\Omega$) to create distinct vortex configurations (1, 2, 3, or multiple vortices).
  2. Quench: The interaction strength is abruptly reduced to $g_f = g_i / 100$.
  3. Dynamics: Real-time propagation is performed to observe the post-quench evolution of density, vortex geometry, and fragmentation.
• Observables: The study analyzes the one-body density, the degree of fragmentation ($F = 1 - n_1$, where $n_1$ is the largest natural orbital occupation), and the time evolution of natural orbitals.

3. Key Contributions

• Beyond-Mean-Field Validation: The paper demonstrates that mean-field theory fails to predict the correct vortex geometry and formation in strongly interacting regimes. For example, in single and double vortex cases, mean-field theory predicts no vortex or incorrect density profiles, whereas the many-body approach correctly identifies the vortex structures due to significant fragmentation.
• Universal Out-of-Equilibrium Response: The authors identify a universal dynamical response to interaction quenches where the density cloud fragments into pieces equal to the number of initial vortices, rotating in the opposite direction to the imprinting rotation.
• Correlation between Fragmentation and Vortex Dynamics: A direct link is established between the oscillatory behavior of dynamical fragmentation and the timescales of vortex revival and distortion.

4. Key Results

A. Pre-Quench State: Vortex Formation and Fragmentation

• Interplay of Parameters: Vortex geometry is highly sensitive to the competition between rotation ($\Omega$) and interaction strength ($g$).
  • Single Vortex: Achieved at moderate $g$ and low $\Omega$.
  • Multiple Vortices: Higher $\Omega$ and $g$ lead to triangular, pentagonal, diamond, or square arrangements.
• Failure of Mean-Field: In the single, double, and triple vortex cases, the mean-field approach ($M=1$) fails to reproduce the vortex structures seen in the many-body simulations ($M=4$). This discrepancy correlates with high fragmentation ($F$) in the many-body state.
• Fragmentation Behavior: Fragmentation is non-monotonic with respect to rotation frequency, showing peaks and minima depending on the specific vortex configuration.

B. Post-Quench Dynamics

The dynamics following the quench reveal two distinct timescales:

1. Vortex Structure Dynamics: Breathing, distortion, merging, and revival.
2. Cloud Density Dynamics: Outward splitting and rotation of the density cloud.

Specific Regimes Observed:

• Single Vortex: Exhibits a breathing mode (quasi-periodic expansion and contraction of the core). The vortex undergoes a complete revival, returning to its initial state. The breathing frequency is synchronized with oscillations in dynamical fragmentation.
• Double and Triple Vortices:
  • The vortex structures distort and eventually merge.
  • The surrounding density cloud splits into 2 or 3 fragments, rotating counter-clockwise (opposite to the system's rotation).
  • Pseudo-Revival: The split clouds interact and merge, leading to a partial or "pseudo" revival of the vortex structure. This occurs on a timescale linked to fragmentation oscillations.
• Multiple Vortices (e.g., 8 vortices):
  • The system exhibits chaotic many-body dynamics.
  • The density cloud splits into 8 fragments.
  • Unlike simpler cases, no revival is observed; the vortex structure remains distorted and chaotic over long times.
  • The dynamical fragmentation becomes aperiodic, reflecting the lack of a well-defined revival period.

C. Comparison with Mean-Field Dynamics

• In the multiple-vortex case, while mean-field theory predicts the correct number of vortices, it fails to capture the instability and rearrangement.
• Mean-field vortices remain stable and retain their octagonal arrangement, whereas the many-body vortices undergo significant deformation, reorientation, and chaotic evolution.

5. Significance and Conclusion

• Beyond Mean-Field Physics: The study highlights that for finite-sized systems with strong interactions, quantum correlations (fragmentation) are not negligible and fundamentally alter vortex stability and dynamics.
• Universal Features: The splitting of the density cloud into $N_v$ fragments (where $N_v$ is the vortex number) rotating in the opposite direction appears to be a universal feature of interaction quenches in rotating BECs.
• Timescale Synchronization: The synchronization between vortex revival/distortion and orbital fragmentation provides a new diagnostic tool for probing many-body correlations in non-equilibrium quantum fluids.
• Future Directions: The authors suggest that future work should explore asymmetric traps (elongated geometries) where fragmentation induces spatial density splitting, and investigate attractive Bose gases and long-range interactions, which remain underexplored in the context of vortex dynamics.

In summary, this work provides a rigorous, numerically exact characterization of how quantum vortices behave when pushed out of equilibrium, demonstrating that many-body effects are crucial for understanding the stability and revival of topological defects in ultracold quantum fluids.