pith. sign in

arxiv: 2604.08714 · v1 · submitted 2026-04-09 · ❄️ cond-mat.quant-gas · physics.atom-ph

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

Lauro Tomio , S. Sabari , Arnaldo Gammal , R. K. Kumar This is my paper

Pith reviewed 2026-05-10 16:42 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas physics.atom-ph
keywords binary Bose-Einstein condensateimmiscible-miscible quenchingvorticessound wavesKolmogorov scalingkinetic energy spectratwo-dimensional superfluid
0
0 comments X

The pith

Sudden reduction in interspecies repulsion drives vortex production and dominant sound waves in two-dimensional binary Bose-Einstein condensates.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines how a sudden drop in the interspecies scattering length converts two-dimensional binary rubidium condensates from immiscible to miscible, triggering instabilities. Numerical evolution shows the transition produces large vortices together with sound waves, and sound waves become the leading feature once the system settles. Spectra of the compressible and incompressible kinetic energy follow a Kolmogorov k to the minus five-thirds power law near the onset of the instability yet display a bottleneck before small-scale dissipation. In the final miscible state both vorticity and sound activity remain steady, while the strength of these excitations depends linearly on the miscibility parameter fixed by the starting density configuration.

Core claim

The central claim is that immiscible-to-miscible quenching instabilities in two-dimensional binary condensates are driven by the simultaneous production of large vortices and sound waves (phonons), with sound waves prevailing in the long-term evolution. The compressible and incompressible kinetic-energy spectra exhibit classical Kolmogorov scaling at the onset of the instabilities, followed by a bottleneck effect before the ultraviolet dissipation range. In the asymptotic miscible regime vorticity and sound-wave production remain practically stable, and a linear relation holds between the miscibility parameter and the initial quenching configuration.

What carries the argument

The instantaneous reduction of the interspecies scattering length a12 that triggers the transition, together with the decomposition of kinetic energy into compressible and incompressible components whose spectra are examined versus wave number k.

If this is right

  • Large vortices appear rapidly and coexist with propagating sound waves.
  • Sound-wave excitations dominate the long-time evolution after the quench.
  • Compressible and incompressible kinetic-energy spectra follow k to the minus five-thirds scaling at the onset of instability.
  • A bottleneck develops in the spectra before the ultraviolet dissipation range is reached.
  • Vorticity and sound-wave production stabilize in the final miscible regime and show a linear dependence on the initial configuration.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The bottleneck feature suggests that higher-resolution simulations would be needed to reach the dissipation range.
  • The same quenching protocol could be applied to other multicomponent superfluids to generate controlled quantum turbulence.
  • The linear relation implies that the initial overlap between the two components largely determines the final level of excitations.

Load-bearing premise

The two-dimensional circular-box confinement and instantaneous reduction in a12 produce dynamics representative of real experimental systems without significant three-dimensional effects or numerical discretization artifacts.

What would settle it

A quasi-two-dimensional experiment that measures the time-dependent vortex number and sound-wave amplitude after a rapid lowering of a12 and tests whether sound waves remain dominant at late times.

read the original abstract

Immiscible to miscible quenching transitions (IMQT) in homogeneous Bose-Einstein condensate are investigated, considering rubidium isotopes $^{85}$Rb and $^{87}$Rb confined in a two-dimensional (2D) circular box, under two different initial configurations. These IMQT instabilities, triggered by sudden reductions in the two-body interspecies scattering length $a_{12}$, are explored under two distinct initialconditions, highlighting the critical role of nonlinear dynamics in their evolution. The numerical simulations indicate that the instability dynamics are primarily driven by the production of large vortices and the propagation of sound waves (phonons), with sound wave excitations prevailing in the long-term evolution. The compressible and incompressible parts of the kinetic energy spectra, in terms of the wave number $k$, are confronted with the classical Kolmogorov scaling, $k^{-5/3}$ for turbulence, which is observed in the onset of instabilities. Before reaching the ultraviolet dissipation region at small scales, the IMQT spectra exhibit a bottleneck effect, indicating a clear departure from classical scaling behavior. In the time asymptotic miscible regime, it is observed that the vorticity and sound-wave production remain practically stable. In this regime, for both cases investigated, a linear relation is also recognized between the miscibility parameter and the initial IMQT configuration.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The manuscript investigates immiscible-to-miscible quenching transitions (IMQT) in two-dimensional binary Bose-Einstein condensates of ^{85}Rb and ^{87}Rb confined in a circular box. Through numerical simulations of the Gross-Pitaevskii equations with sudden reductions in the interspecies scattering length a_{12} under two initial configurations, it concludes that the dynamics are primarily driven by the production of large vortices and the propagation of sound waves, with sound excitations dominating the long-term evolution. The compressible and incompressible kinetic energy spectra are compared to the Kolmogorov k^{-5/3} scaling, revealing a bottleneck effect before dissipation, and stability of vorticity and sound production is reported in the asymptotic miscible regime, along with a linear relation between the miscibility parameter and the initial configuration.

Significance. Should the numerical findings prove robust upon detailed validation, the work would offer valuable insights into the mechanisms of quenching-induced instabilities and the emergence of turbulent-like scaling in quantum fluids. The identification of vortex and phonon roles, along with the bottleneck phenomenon, could guide theoretical and experimental studies of non-equilibrium dynamics in multicomponent BECs. The linear relation in the miscible phase adds a quantitative aspect that might be testable.

major comments (3)
  1. [Numerical simulations and methods] The central claims on vortex production, phonon propagation, and long-time sound dominance rest on 2D GPE time evolution, but the manuscript provides no details on spatial discretization (grid size or method), time-stepping scheme, or convergence tests with respect to resolution and domain size. This is load-bearing because discretization artifacts in a circular box can contaminate the incompressible/compressible decomposition and the reported k^{-5/3} onset.
  2. [Initial conditions and quench protocol] The instantaneous reduction in a_{12} is taken as the trigger for IMQT instabilities, yet no analysis is given of how a finite-duration experimental ramp would modify the initial instability spectrum or the subsequent Kolmogorov scaling and bottleneck. This assumption directly affects the attribution of driving mechanisms.
  3. [Model and confinement] The study is performed entirely in 2D with a circular box, without justification for neglecting out-of-plane modes or comparison to 3D effects (e.g., vortex tilting or Kelvin waves) that could shift the relative weights of vortex versus sound contributions at long times. This undermines the claim that the observed dynamics are representative.
minor comments (2)
  1. [Abstract] The abstract refers to 'two different initial configurations' without specifying them, which obscures the generality of the reported linear relation between the miscibility parameter and configuration.
  2. [Results and discussion] Explicit definitions or a table for the pre- and post-quench values of a_{12}, the miscibility parameter, and the box radius would improve clarity when discussing the linear relation in the asymptotic regime.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough review and valuable suggestions. We address each of the major comments below and indicate the revisions we will make to the manuscript.

read point-by-point responses

  1. Referee: [Numerical simulations and methods] The central claims on vortex production, phonon propagation, and long-time sound dominance rest on 2D GPE time evolution, but the manuscript provides no details on spatial discretization (grid size or method), time-stepping scheme, or convergence tests with respect to resolution and domain size. This is load-bearing because discretization artifacts in a circular box can contaminate the incompressible/compressible decomposition and the reported k^{-5/3} onset.

    Authors: We agree that the numerical methods were not described in sufficient detail. In the revised manuscript, we will add a dedicated subsection on the numerical implementation. This will specify the spatial discretization (Fourier pseudospectral method on a grid appropriate for the circular domain), the time-stepping scheme (split-step Fourier method), and the convergence tests with respect to grid resolution and domain size that confirm the robustness of the vortex production, phonon propagation, and kinetic energy spectra. revision: yes

  2. Referee: [Initial conditions and quench protocol] The instantaneous reduction in a_{12} is taken as the trigger for IMQT instabilities, yet no analysis is given of how a finite-duration experimental ramp would modify the initial instability spectrum or the subsequent Kolmogorov scaling and bottleneck. This assumption directly affects the attribution of driving mechanisms.

    Authors: We acknowledge that the instantaneous quench is an idealization. In the revised manuscript, we will include a discussion of the quench protocol, noting that the sudden reduction approximates fast experimental ramps and that the instability spectrum at onset is primarily determined by the change in a_{12}. We will argue that the Kolmogorov scaling and bottleneck are expected to persist for sufficiently rapid ramps, while slower ramps may alter the initial growth rates. A detailed comparison with finite ramps is beyond the current scope but will be mentioned as an avenue for future research. revision: partial

  3. Referee: [Model and confinement] The study is performed entirely in 2D with a circular box, without justification for neglecting out-of-plane modes or comparison to 3D effects (e.g., vortex tilting or Kelvin waves) that could shift the relative weights of vortex versus sound contributions at long times. This undermines the claim that the observed dynamics are representative.

    Authors: Our investigation is restricted to two dimensions to isolate the in-plane vortex and sound dynamics in a controlled manner. We will revise the manuscript to provide justification for the 2D model, explaining that it corresponds to quasi-2D BECs with strong axial confinement where out-of-plane modes are energetically suppressed. We will also note the limitations regarding 3D effects and how they might influence long-time behavior, while emphasizing that the 2D results offer clear insights into the mechanisms that can guide 3D studies. revision: yes

Circularity Check

0 steps flagged

No circularity: results follow from direct numerical evolution of the governing equations

full rationale

The paper performs direct time-dependent simulations of the 2D Gross-Pitaevskii equations for a binary BEC under instantaneous interspecies scattering-length quench in a circular box. All reported features—vortex production, sound-wave (phonon) dominance at late times, onset of k^{-5/3} scaling in the compressible/incompressible kinetic-energy spectra, bottleneck effect, and the observed linear relation between miscibility parameter and initial configuration—are extracted by post-processing the simulated density and velocity fields. No parameters are fitted to data and then re-predicted, no self-citations supply load-bearing uniqueness theorems or ansatzes, and no quantity is defined in terms of itself. The derivation chain is therefore self-contained against the initial-value problem and the numerical outputs.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard mean-field modeling of binary BECs with chosen quench parameters and box geometry; no new entities are postulated.

free parameters (2)
  • initial interspecies scattering length a12 and its quench depth
    Specific values chosen to trigger the immiscible-to-miscible transition; control parameters set by hand for the two configurations.
  • circular box radius and trap parameters
    Dimensions of the 2D confinement chosen to define the homogeneous system.
axioms (1)
  • domain assumption The binary condensate is accurately described by the two-dimensional coupled Gross-Pitaevskii equations in the mean-field regime.
    Invoked implicitly as the basis for all numerical evolution; standard for dilute ultracold gases at the temperatures considered.

pith-pipeline@v0.9.0 · 5545 in / 1371 out tokens · 50707 ms · 2026-05-10T16:42:57.157798+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [1]

    R. K. Kumar, S. Sabari, A. Gammal, and L. Tomio, Rayleigh-Taylor, Kelvin- Helmholtz, and immiscible-to-miscible quenching instabilities in binary Bose- Einstein condensates, Phys. Rev. A112, 033312 (2025)

    work page 2025

  2. [2]

    A. N. Kolmogorov, The local structure of turbulence in incompressible viscous fluid for very large Reynolds’ numbers, Dokl. Akad. Nauk SSSR30, 301 (1941)

    work page 1941

  3. [3]

    Frisch, Turbulence: The Legacy of A

    U. Frisch, Turbulence: The Legacy of A. N. Kolmogorov (Cambridge University Press, Cambridge, 1995)

    work page 1995

  4. [4]

    R. J. Donnelly and C. E. Swanson, Quantum turbulence, J. Fluid Mech.173, 387 (1986)

    work page 1986

  5. [5]

    F. A. Bayocboc Jr., J. Dziarmaga, and W. H. Zurek, Biased dynamics of the miscible-immiscible quantum phase transition in a binary Bose-Einstein condensate, Phys. Rev. B109, 064501 (2024)

    work page 2024

  6. [6]

    Mukherjee, S

    K. Mukherjee, S. I. Mistakidis, P. G. Kevrekidis, and P. Schmelcher, Quench induced vortex-bright-soliton formation in binary Bose-Einstein condensates, J. Phys. B: At. Mol. Opt. Phys.53, 055302 (2020)

    work page 2020

  7. [7]

    Y. Eto, M. Takahashi, M. Kunimi, H. Saito, and T. Hirano, Nonequilibrium dynamics induced by miscible–immiscible transition in binary Bose–Einstein condensates, New J. Phys.18, 073029 (2016)

    work page 2016

  8. [8]

    Kasamatsu and M

    K. Kasamatsu and M. Tsubota, Multiple Domain Formation Induced by Modula- tion Instability in Two-Component Bose-Einstein Condensates, Phys. Rev. Lett. 93, 100402 (2004)

    work page 2004

  9. [9]

    Thiruvalluvar, E

    R.T. Thiruvalluvar, E. Wamba, S. Sabari, K. Porsezian, Impact of higher- order nonlinearity on modulational instability in two-component Bose-Einstein condensates, Phys. Rev. E99, 032202 (2019)

    work page 2019

  10. [10]

    A. S. Bradley, R. K. Kumar, S. Pal, and X. Yu, Spectral analysis for compressible quantum fluids, Phys. Rev. A106, 043322 (2022)

    work page 2022

  11. [11]

    M. T. Reeves, T. P. Billam, B. P. Anderson, and A. S. Bradley, Inverse energy cascade in forced two-dimensional quantum turbulence, Phys. Rev. Lett.110, 13 104501 (2013)

    work page 2013

  12. [12]

    E. A. L. Henn, J. A. Seman, G. Roati, K. M. F. Magalh˜ aes, and V. S. Bagnato, Emergence of turbulence in an oscillating Bose-Einstein condensate, Phys. Rev. Lett.103, 045301 (2009)

    work page 2009

  13. [13]

    M. C. Tsatsos, P. E. S. Tavares, A. Cidrim, A. R. Fritsch, M. A. Caracanhas, F. E. A. dos Santos, C. F. Barenghi, and V. S. Bagnato, Quantum turbulence in trapped atomic Bose–Einstein condensates, Phys. Rep.622, 1 (2016)

    work page 2016

  14. [14]

    M. A. Moreno-Armijos, A. R. Fritsch, A. D. Garc´ ıa-Orozco, S. Sab, G. Telles, Y. Zhu, L. Madeira, S. Nazarenko, V. I. Yukalov, and V. S. Bagnato, Observation of Relaxation Stages in a Nonequilibrium Closed Quantum System: Decaying Turbulence in a Trapped Superfluid, Phys. Rev. Lett.134, 023401 (2025)

    work page 2025

  15. [15]

    Madeira, A

    L. Madeira, A. D. Garc´ ıa-Orozco, M. A. Moreno-Armijos, A. R. Fritsch, and V. S. Bagnato, Universal scaling in far-from-equilibrium quantum systems: An equivalent differential approach, Proc. Natl. Acad. Sci. U.S.A. 121, e2404828121 (2024)

    work page 2024

  16. [16]

    A. N. da Silva, R. K. Kumar, A. S. Bradley, and L. Tomio, Vortex generation in stirred binary Bose-Einstein condensates, Phys. Rev. A107, 033314 (2023)

    work page 2023

  17. [17]

    Tomio, A

    L. Tomio, A. N. da Silva, S. Sabari, R. K. Kumar, Dynamical vortex produc- tion and quantum turbulence in perturbed Bose-Einstein condensates, Few-Body Systems65, 13 (2024)

    work page 2024

  18. [18]

    Sabari, R

    S. Sabari, R. K. Kumar, L. Tomio, Vortex dynamics and turbulence in dipolar Bose-Einstein condensates, Phys. Rev. A109, 023313 (2024)

    work page 2024

  19. [19]

    R. K. Kumar, L. Tomio, and A. Gammal, Vortex patterns in rotating dipolar Bose–Einstein condensate mixtures with squared optical lattices, J. Phys. B: At. Mol. Opt. Phys.52, 025302 (2019)

    work page 2019

  20. [20]

    Timmermans, P

    E. Timmermans, P. Tommasini, M. Hussein, and A. Kerman, Feshbach resonances in atomic Bose-Einstein condensates, Phys. Rep.315, 199 (1999)

    work page 1999

  21. [21]

    R. K. Kumar, P. Muruganandam, L. Tomio, and A. Gammal, Miscibility in cou- pled dipolar and non-dipolar Bose-Einstein condensates, J. Phys. Commun.1, 035012 (2017)

    work page 2017

  22. [22]

    V. S. L’vov, S. V. Nazarenko, and O. Rudenko, Bottleneck crossover between classical and quantum superfluid turbulence, Phys. Rev. B76, 024520 (2007). 14

    work page 2007