Immiscible to miscible quenching instabilities in two-dimensional binary Bose-Einstein condensates

This paper investigates immiscible-to-miscible quenching instabilities in two-dimensional binary Bose-Einstein condensates of rubidium isotopes, revealing that the transition dynamics are driven by vortex production and sound waves, exhibit a bottleneck effect deviating from classical Kolmogorov turbulence scaling, and ultimately stabilize into a linear relationship between the miscibility parameter and the initial configuration.

The Gist

Imagine you have a magical, super-cold bowl of soup made of two different types of "quantum ingredients" (specifically, two kinds of Rubidium atoms). In this super-cold state, called a Bose-Einstein Condensate (BEC), the atoms stop acting like individual particles and start behaving like a single, giant wave.

Usually, if you mix oil and water, they separate because they don't like each other. In this quantum soup, the atoms can also be "immiscible" (they hate each other and separate) or "miscible" (they love each other and mix perfectly).

This paper is about a scientific experiment where the researchers suddenly forced these two ingredients to switch from "hating each other" to "loving each other." They call this a "quenching transition."

Here is the story of what happened, explained simply:

1. The Setup: The Quantum Salad

The researchers put their two types of atoms into a flat, circular bowl (a 2D box).

• Scenario A (The Tennis Ball): They arranged the atoms in three separate zones, like the stripes on a tennis ball.
• Scenario B (The Axial Split): They split the atoms down the middle, like a slice of pizza with two halves.

In both cases, the atoms were initially separated because they were repelling each other.

2. The Trigger: The "Snap"

At the start of the experiment, the researchers suddenly changed the rules. They turned down the "repulsion dial" between the two types of atoms.

• Think of it like suddenly removing the force field that was keeping two magnets apart.
• Because the atoms were now allowed to mix, they rushed toward each other. This sudden rush created chaos, much like dropping a rock into a calm pond, but on a quantum scale.

3. The Chaos: Swirls and Ripples

When the atoms rushed to mix, two main things happened:

• Vortices (The Swirls): Just like water going down a drain, the atoms started spinning into tiny tornadoes. These are called vortices. The researchers found that the stronger the initial "push" to mix, the faster these tornadoes formed.
• Sound Waves (The Ripples): The atoms also created ripples, which are essentially sound waves (phonons) traveling through the soup.

4. The Race: Swirls vs. Ripples

The researchers wanted to know: Which one wins in the long run? The spinning tornadoes or the sound ripples?

• The Result: In the beginning, both happened at once. But as time went on, the sound waves (ripples) took over. The energy of the system became mostly about these ripples, while the spinning tornadoes settled into a stable, quiet background.
• The Analogy: Imagine a party where everyone is dancing wildly (vortices) and shouting (sound waves). After a while, the wild dancing stops, but the music (sound waves) keeps playing, and the room settles into a steady rhythm.

5. The "Traffic Jam" in the Energy

The researchers looked at the "energy spectrum," which is like a map of how much energy is in the different sizes of waves.

• The Classic Rule: In normal fluids (like water in a river), energy follows a famous rule called Kolmogorov scaling. It's like a traffic flow where energy moves smoothly from big waves to small waves.
• The Quantum Twist: In this quantum soup, the energy flow got stuck. Before the energy could break down into tiny, invisible waves, it hit a "bottleneck." It was like a traffic jam where cars (energy) piled up before reaching the exit. This is a unique feature of quantum fluids that doesn't happen in normal water.

6. The Big Discovery: A Simple Math Rule

The most exciting part of the paper is a simple pattern the researchers found.

• They noticed that the bigger the initial "push" to mix (the bigger the change they made to the repulsion dial), the faster the atoms oscillated (shook back and forth) once they settled down.
• It's a straight line: More initial chaos = Faster final shaking.
• They even found a formula that predicts exactly how fast the atoms will shake based on how hard they pushed them to mix in the first place.

Why Does This Matter?

This isn't just about mixing atoms; it's about understanding Quantum Turbulence.

• Scientists want to understand how turbulence works in the universe, from black holes to superfluids.
• By studying how these atoms mix and create "tornadoes" and "ripples," we learn how energy moves in the most extreme, cold environments in the universe.
• The "bottleneck" effect they found suggests that quantum fluids behave differently than the water in your sink, offering new clues about the fundamental laws of physics.

In a nutshell: The researchers pushed two quantum ingredients to mix, watched them create a chaotic storm of tiny tornadoes and sound waves, and discovered that the sound waves eventually win the race. They also found a simple math rule that predicts exactly how wild the party will get based on how hard you start the music.

Technical Summary

1. Problem Statement

The study investigates the dynamical instabilities and turbulence generated in a two-dimensional (2D) binary Bose-Einstein condensate (BEC) mixture of Rubidium isotopes ($^{85}\text{Rb}$ and $^{87}\text{Rb}$) confined within a uniform circular box. Specifically, the research focuses on Immiscible-to-Miscible Quenching Transitions (IMQT).

The core problem is to understand how a sudden reduction in the interspecies scattering length ($a_{12}$) drives the system from a phase-separated (immiscible) state to a mixed (miscible) state. The authors aim to:

• Distinguish dynamics driven by nonlinear interaction variations (quenching) from those driven by external linear forces.
• Analyze the emergence of quantum turbulence, vortex generation, and sound wave (phonon) production.
• Investigate the scaling laws of the kinetic energy spectra (compressible vs. incompressible) and compare them with classical Kolmogorov turbulence scaling ($k^{-5/3}$).
• Determine how different initial spatial configurations and the magnitude of the quench affect the evolution and asymptotic behavior of the system.

2. Methodology

The authors employ a numerical approach based on the Gross-Pitaevskii (GP) mean-field formalism adapted for a 2D binary system.

• Physical Model:
  • System: A binary mixture of $^{85}\text{Rb}$ (species 1) and $^{87}\text{Rb}$ (species 2) with equal atom numbers ($N = 3.5 \times 10^4$).
  • Confinement: A 2D circular box potential with radius $R$ and height $V_0 \gg \mu$ (chemical potential).
  • Equations: Coupled dimensionless 2D GP equations are solved using the split-step Fourier method on an $800 \times 800$ grid.
  • Parameters: Intraspecies scattering lengths are fixed at $a_{11} = a_{22} = 90a_0$. The interspecies length $a_{12}$ is the control parameter.
• Quenching Protocol:
  • The system is initialized in an immiscible ground state where $a_{12} > a_{ii}$.
  • At $t=0$, $a_{12}$ is suddenly reduced (quenched) to a value where $a_{12} < a_{ii}$, triggering the transition to miscibility.
  • Two distinct initial spatial configurations are tested:
    1. Central (Tripartite): Three non-overlapping spatial domains (Species 2 in the center, Species 1 on the sides). Initial $a_{12} = 100a_0 \to 60a_0$ ($\Delta\delta = 4/9$).
    2. Axial (Bipartite): Two axially separated domains. Initial $a_{12} = 110a_0 \to 60a_0$ ($\Delta\delta = 5/9$).
• Analysis Techniques:
  • Energy Decomposition: The kinetic energy is decomposed into compressible (linked to density fluctuations/sound waves) and incompressible (linked to vorticity/rotational motion) components using Helmholtz decomposition.
  • Spectral Analysis: Power spectra of the kinetic energy are calculated in momentum space ($k$-space) to identify scaling laws.
  • Miscibility Parameter ($\Lambda$): Defined via the overlap integral of the two density profiles to track the transition from immiscible ($\Lambda=0$) to miscible ($\Lambda=1$).

3. Key Contributions

• Enhanced IMQT Magnitude: Unlike previous studies (Ref. [1]) that used smaller quench magnitudes, this study employs significantly larger transitions ($\Delta\delta = 4/9$ and $5/9$). This allows for the observation of faster dynamical responses and clearer separation of time scales.
• Linear Correlation of Frequency and Quench: The authors establish a novel linear relationship between the asymptotic oscillation frequency of the miscibility parameter ($\nu_\Lambda$) and the initial quench magnitude ($\Delta\delta$).
  • Formula: $\Delta\delta \approx \alpha \nu_\Lambda + 0.15$.
  • The slope $\alpha$ depends on the initial geometry ($\alpha \approx 0.37$ for 3 domains, $\alpha \approx 1.17$ for 2 domains), while the intercept ($0.15$) represents a threshold IMQT independent of geometry.
• Bottleneck Effect Identification: The study confirms the presence of a "bottleneck effect" in the energy spectra before reaching the ultraviolet dissipation region, indicating a departure from classical scaling behavior in the intermediate scales.
• Dominance of Compressible Modes: In the long-term asymptotic regime, the compressible kinetic energy (sound waves) is found to dominate the incompressible energy (vortices) by a factor of approximately 5, regardless of the initial configuration.

4. Key Results

• Dynamics and Timescales:
  • Larger IMQT magnitudes lead to significantly faster instability onset and evolution. The system reaches the asymptotic miscible regime much earlier (e.g., $t \approx 10$) compared to previous studies with smaller quenches ($t \approx 20$).
  • The system generates numerous vortex dipoles and turbulent flows immediately after the quench.
• Spectral Scaling (Kolmogorov):
  • During the onset of instabilities (short time intervals), the incompressible and compressible energy spectra exhibit the classical Kolmogorov scaling ($k^{-5/3}$), characteristic of turbulence.
  • This scaling is observed in specific time windows (e.g., $0.5 \le t \le 0.8$ for the central case; $0.7 \le t \le 1.2$ for the axial case).
  • As the system evolves, the spectra deviate from this scaling, showing a bottleneck effect before entering the dissipation region.
• Asymptotic Stability:
  • In the long-term limit ($t > 20$), the vorticity and sound-wave production stabilize.
  • The number of vortices ($N_v$) remains constant (non-zero), while the density overlap $\Lambda$ oscillates with a frequency directly proportional to the initial quench strength.
• Configuration Differences:
  • Central (3 domains): The instability onset is slightly faster, and the spectral bottleneck is more pronounced in the outer species ($^{85}\text{Rb}$).
  • Axial (2 domains): The system shows interference fringes localized near the trap boundaries, distinct from the central case.

5. Significance

This work advances the understanding of quantum turbulence and non-equilibrium dynamics in binary quantum fluids.

• Universal Scaling: It reinforces the connection between quantum turbulence in BECs and classical fluid turbulence, validating the applicability of Kolmogorov scaling in the early stages of IMQT.
• Control Mechanisms: The discovery of a linear relationship between the quench magnitude and the asymptotic oscillation frequency provides a potential control mechanism for tuning the dynamical properties of binary condensates.
• Energy Partitioning: The finding that compressible modes (sound) eventually dominate incompressible modes (vortices) in the asymptotic regime offers insight into the energy cascade and dissipation mechanisms in 2D quantum fluids.
• Experimental Relevance: The use of larger quench magnitudes and the identification of specific time scales for turbulence onset provide actionable parameters for experimentalists working with Feshbach resonances to manipulate scattering lengths in Rubidium mixtures.

In summary, the paper demonstrates that the magnitude of the interaction quench is a critical parameter governing the speed of instability development, the spectral characteristics of turbulence, and the final oscillatory state of the miscible condensate.