pith. sign in

arxiv: 2512.16756 · v3 · submitted 2025-12-18 · ❄️ cond-mat.quant-gas · physics.atom-ph

Structure of the mean-field yrast spectrum of a two-component Bose gas in a ring: role of interaction asymmetry

Pith reviewed 2026-05-16 21:13 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas physics.atom-ph
keywords two-component Bose gasyrast spectruminteraction asymmetryGross-Pitaevskii equationsring geometryplane-wave statessoliton statespersistent currents
0
0 comments X

The pith

Interaction asymmetry in a two-component Bose gas on a ring changes how plane-wave states become the yrast states.

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

The paper examines how breaking the SU(2) symmetry through unequal intra- and inter-component interactions modifies the mean-field yrast spectrum of a two-component Bose gas in a ring. It finds that the critical curves marking the appearance of plane-wave yrast states and the mechanisms producing them depend on whether intra-component repulsion is weaker or stronger than inter-component repulsion. When intra-component interactions are weaker, plane-wave states replace soliton states by continuous evolution, though under tighter conditions than in the symmetric case. When intra-component interactions are stronger, plane-wave states instead overtake soliton states through branch crossings and gain markedly higher stability. These changes directly affect the range of angular momenta supporting stable persistent currents.

Core claim

Numerical solutions of the coupled Gross-Pitaevskii equations show that plane-wave yrast states emerge either by continuous replacement of soliton branches when the intra-component interaction strength is smaller than the inter-component strength, or by discrete branch crossings when the intra-component strength is larger; in the latter regime the plane-wave states acquire significantly enhanced stability.

What carries the argument

Coupled Gross-Pitaevskii equations solved numerically for propagating soliton states, from which the energy versus angular-momentum curves and the critical interaction-asymmetry values separating different yrast regimes are extracted.

If this is right

  • Plane-wave yrast states exist only for more restricted ranges of angular momentum when intra-component interactions are weaker.
  • Branch crossings allow plane-wave states to appear at fractional angular momenta even when soliton states would otherwise remain lower in energy.
  • Persistent-current stability is markedly higher when intra-component interactions exceed inter-component interactions.
  • The non-analytic features of the yrast spectrum survive asymmetry but change their locations and character.

Where Pith is reading between the lines

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

  • The enhanced stability in the stronger-intra-component regime may allow persistent currents to survive at higher temperatures or larger particle numbers than in the symmetric case.
  • The same asymmetry dependence should appear in the spectrum of collective modes around the yrast states.
  • Pushing the asymmetry to extreme values recovers the single-component limit, providing a continuous connection between the two known cases.

Load-bearing premise

The mean-field Gross-Pitaevskii equations remain accurate across the full range of interaction asymmetries examined.

What would settle it

A measurement of the yrast energy spectrum in a ring-trapped two-component Bose gas as the ratio of intra- to inter-component scattering lengths is varied, checking for the predicted shift from continuous replacement to branch-crossing transitions at the calculated critical asymmetry values.

Figures

Figures reproduced from arXiv: 2512.16756 by Eugene Zaremba, Guan-Hua Huang, Hui Tang, Shizhong Zhang, Zhigang Wu.

Figure 1
Figure 1. Figure 1: FIG. 1. Illustration of the yrast spectrum for a single-component [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The internal part of the mean-field yrast spectrum of an [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. The internal part of the mean-field yrast spectrum of an [PITH_FULL_IMAGE:figures/full_fig_p004_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. The internal part of the mean-field yrast spectrum of an [PITH_FULL_IMAGE:figures/full_fig_p004_4.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The critical curves for the plane-wave yrast state [PITH_FULL_IMAGE:figures/full_fig_p005_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. The upper and middle panels, respectively, show the amplitude and phase of the condensate wave function at [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The critical curves for the plane-wave yrast state [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. The critical curves for plane-wave yrast states [PITH_FULL_IMAGE:figures/full_fig_p007_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. The upper and middle panels, respectively, show the amplitude and phase of the condensate wave function at [PITH_FULL_IMAGE:figures/full_fig_p008_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. The two lowest spectral branches in the vicinity of [PITH_FULL_IMAGE:figures/full_fig_p008_12.png] view at source ↗
read the original abstract

The mean-field yrast spectrum of an SU(2)-symmetric two-component Bose gas confined to a ring geometry is known to exhibit an intricate nonanalytic structure that is absent in single-component systems. In particular, due to the interplay between the species concentration and the atomic interactions, a sequence of plane-wave states can emerge as yrast states at fractional values of the angular momentum per particle. This behavior stands in sharp contrast to the single-component case, where plane-wave states occur only at integer angular momenta. In this paper, we investigate how the structure of the yrast spectrum in a two-component Bose gas is modified by interaction asymmetry. By numerically solving the coupled Gross-Pitaevskii equations for propagating soliton states, we compute the mean-field yrast spectrum and, in particular, determine the critical curves associated with the emergence of various plane-wave yrast states. We find that both the behavior of these critical curves and the mechanisms by which plane-wave yrast states arise depend sensitively on the relative strengths of the inter- and intra-component interactions. When the intra-component interaction is weaker, the plane-wave yrast states replace soliton states through a continuous evolution, as in the SU(2)-symmetric case, although the conditions for their existence become more restrictive. In contrast, when the intra-component interaction is stronger, plane-wave yrast states may emerge by overtaking soliton states via branch crossings, and their stability is significantly enhanced. Our results have important implications for the existence and stability of persistent currents in asymmetric, two-component Bose gases.

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

2 major / 1 minor

Summary. The paper examines the mean-field yrast spectrum of a two-component Bose gas in a ring geometry when the intra- and inter-component interaction strengths are asymmetric. By numerically propagating soliton initial conditions in the coupled Gross-Pitaevskii equations, the authors compute the lowest-energy curves E(L) and identify critical curves separating regions where plane-wave states appear as yrast states at fractional angular momenta. They report that the topology of these curves and the mechanism of plane-wave emergence (continuous replacement of solitons versus overtaking via branch crossings) depend sensitively on whether the intra-component interaction is weaker or stronger than the inter-component interaction, with the latter case also yielding enhanced stability.

Significance. If the numerical results hold under convergence checks, the work extends the known SU(2)-symmetric yrast structure to the asymmetric regime and supplies concrete predictions for the stability of persistent currents in two-component ring condensates. The distinction between continuous evolution and branch-crossing mechanisms is a clear, falsifiable outcome of the mean-field model.

major comments (2)
  1. [Numerical methods / Results] The central distinction between continuous replacement and branch-crossing mechanisms (abstract and results section) rests on the precise topology of the numerically obtained E(L) branches. No spatial discretization (grid points per ring), time-step size, or convergence tests with respect to these parameters are reported, so it is not possible to assess whether small numerical errors could convert an avoided crossing into an apparent crossing or vice versa.
  2. [Results on critical curves] The claim of 'significantly enhanced' stability when intra-component interactions are stronger is stated without quantitative support such as the magnitude of the energy gap at the crossing point or a linear-stability analysis of the plane-wave states (e.g., Bogoliubov-de Gennes spectrum).
minor comments (1)
  1. [Abstract] The abstract refers to 'the range of interaction asymmetries studied' but does not quote the specific ratios g11/g12 and g22/g12 that were scanned; adding these values would help readers reproduce the critical curves.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We are grateful to the referee for their thorough review and valuable suggestions, which have helped us improve the clarity and rigor of our manuscript. Below, we provide detailed responses to each major comment and outline the revisions we will make.

read point-by-point responses
  1. Referee: [Numerical methods / Results] The central distinction between continuous replacement and branch-crossing mechanisms (abstract and results section) rests on the precise topology of the numerically obtained E(L) branches. No spatial discretization (grid points per ring), time-step size, or convergence tests with respect to these parameters are reported, so it is not possible to assess whether small numerical errors could convert an avoided crossing into an apparent crossing or vice versa.

    Authors: We fully agree that providing the numerical discretization parameters and convergence tests is crucial to substantiate the reported topologies of the E(L) branches. In the revised version of the manuscript, we will include these details: we use a spatial grid with 1024 points around the ring, a time step of 0.0001 in the imaginary-time propagation, and we have performed convergence checks by doubling the grid size and halving the time step, confirming that the locations of branch crossings and the continuous replacement mechanisms are unchanged within numerical precision. This ensures that the distinction between the mechanisms is not affected by discretization errors. revision: yes

  2. Referee: [Results on critical curves] The claim of 'significantly enhanced' stability when intra-component interactions are stronger is stated without quantitative support such as the magnitude of the energy gap at the crossing point or a linear-stability analysis of the plane-wave states (e.g., Bogoliubov-de Gennes spectrum).

    Authors: We appreciate this comment and recognize that a quantitative assessment would better support our statement on enhanced stability. In the revised manuscript, we will add quantitative data on the energy gaps at the relevant crossing points for different interaction asymmetries. Additionally, we will include results from a Bogoliubov-de Gennes linear stability analysis applied to the plane-wave states, demonstrating that the lowest excitation frequencies are higher (indicating greater stability) in the regime where intra-component interactions exceed inter-component ones. These additions will provide the requested quantitative support. revision: yes

Circularity Check

0 steps flagged

No circularity in numerical yrast spectrum computation

full rationale

The paper derives its central claims about critical curves and plane-wave yrast emergence mechanisms solely from direct numerical propagation of soliton initial conditions in the coupled Gross-Pitaevskii equations on a ring. No load-bearing step reduces to a self-definition, a fitted parameter renamed as a prediction, or a self-citation chain; the SU(2)-symmetric reference supplies background only and is not invoked to force the asymmetric results. The derivation chain is therefore self-contained against external numerical benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the mean-field Gross-Pitaevskii equations for the asymmetric two-component system and on the numerical accuracy of the soliton-propagation method; no new particles or forces are introduced.

free parameters (1)
  • relative intra- versus inter-component interaction strengths
    Varied continuously as control parameters to trace critical curves; values are chosen by the authors rather than derived.
axioms (1)
  • domain assumption Mean-field Gross-Pitaevskii equations accurately describe the two-component Bose gas in a ring
    All spectrum calculations are performed by solving these equations; the assumption is invoked throughout the abstract.

pith-pipeline@v0.9.0 · 5602 in / 1355 out tokens · 29332 ms · 2026-05-16T21:13:47.888400+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    By numerically solving the coupled Gross-Pitaevskii equations for propagating soliton states, we compute the mean-field yrast spectrum and, in particular, determine the critical curves associated with the emergence of various plane-wave yrast states.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages

  1. [1]

    Byers and C

    N. Byers and C. N. Yang, Theoretical considerations concern- ing quantized magnetic flux in superconducting cylinders, Phys. Rev. Lett.7, 46 (1961)

  2. [2]

    Bloch, Superfluidity in a ring, Phys

    F. Bloch, Superfluidity in a ring, Phys. Rev. A7, 2187 (1973)

  3. [3]

    Mottelson, Yrast spectra of weakly interacting bose-einstein condensates, Phys

    B. Mottelson, Yrast spectra of weakly interacting bose-einstein condensates, Phys. Rev. Lett.83, 2695 (1999)

  4. [4]

    B. S. Deaver and W. M. Fairbank, Experimental evidence for quantized flux in superconducting cylinders, Phys. Rev. Lett.7, 43 (1961)

  5. [5]

    Doll and M

    R. Doll and M. N ¨abauer, Experimental proof of magnetic flux quantization in a superconducting ring, Phys. Rev. Lett.7, 51 (1961)

  6. [6]

    C. Ryu, M. F. Andersen, P. Clad ´e, V . Natarajan, K. Helmer- son, and W. D. Phillips, Observation of persistent flow of a bose-einstein condensate in a toroidal trap, Phys. Rev. Lett.99, 260401 (2007)

  7. [7]

    Ramanathan, K

    A. Ramanathan, K. C. Wright, S. R. Muniz, M. Zelan, W. T. Hill, C. J. Lobb, K. Helmerson, W. D. Phillips, and G. K. Camp- bell, Superflow in a toroidal bose-einstein condensate: An atom circuit with a tunable weak link, Phys. Rev. Lett.106, 130401 (2011)

  8. [8]

    Moulder, S

    S. Moulder, S. Beattie, R. P. Smith, N. Tammuz, and Z. Hadz- ibabic, Quantized supercurrent decay in an annular bose- einstein condensate, Phys. Rev. A86, 013629 (2012)

  9. [9]

    Beattie, S

    S. Beattie, S. Moulder, R. J. Fletcher, and Z. Hadzibabic, Per- sistent currents in spinor condensates, Phys. Rev. Lett.110, 025301 (2013)

  10. [10]

    K. C. Wright, R. B. Blakestad, C. J. Lobb, W. D. Phillips, and G. K. Campbell, Driving phase slips in a superfluid atom circuit with a rotating weak link, Phys. Rev. Lett.110, 025302 (2013)

  11. [11]

    Y . Cai, D. G. Allman, P. Sabharwal, and K. C. Wright, Persis- tent currents in rings of ultracold fermionic atoms, Phys. Rev. Lett.128, 150401 (2022)

  12. [12]

    Del Pace, K

    G. Del Pace, K. Xhani, A. Muzi Falconi, M. Fedrizzi, N. Grani, D. Hernandez Rajkov, M. Inguscio, F. Scazza, W. J. Kwon, and G. Roati, Imprinting persistent currents in tunable fermionic rings, Phys. Rev. X12, 041037 (2022)

  13. [13]

    J. Polo, W. Chetcuti, T. Haug, A. Minguzzi, K. Wright, and L. Amico, Persistent currents in ultracold gases, Physics Re- ports1137, 1 (2025), persistent currents in ultracold gases

  14. [14]

    Kawaguchi and M

    Y . Kawaguchi and M. Ueda, Spinor bose–einstein condensates, Physics Reports520, 253 (2012), spinor Bose–Einstein conden- sates

  15. [15]

    Lin and I

    Y .-J. Lin and I. B. Spielman, Synthetic gauge potentials for ul- tracold neutral atoms, Journal of Physics B: Atomic, Molecular and Optical Physics49, 183001 (2016)

  16. [16]

    Smyrnakis, S

    J. Smyrnakis, S. Bargi, G. M. Kavoulakis, M. Magiropou- los, K. K ¨arkk¨ainen, and S. M. Reimann, Mixtures of bose gases confined in a ring potential, Phys. Rev. Lett.103, 100404 (2009)

  17. [17]

    Anoshkin, Z

    K. Anoshkin, Z. Wu, and E. Zaremba, Persistent currents in a bosonic mixture in the ring geometry, Phys. Rev. A88, 013609 10 (2013)

  18. [18]

    Wu and E

    Z. Wu and E. Zaremba, Mean-field yrast spectrum of a two- component bose gas in ring geometry: Persistent currents at higher angular momentum, Phys. Rev. A88, 063640 (2013)

  19. [19]

    M. Abad, A. Sartori, S. Finazzi, and A. Recati, Persistent cur- rents in two-component condensates in a toroidal trap, Phys. Rev. A89, 053602 (2014)

  20. [20]

    Z. Wu, E. Zaremba, J. Smyrnakis, M. Magiropoulos, N. K. Efremidis, and G. M. Kavoulakis, Mean-field yrast spectrum and persistent currents in a two-component bose gas with inter- action asymmetry, Phys. Rev. A92, 033630 (2015)

  21. [21]

    Pecci, G

    G. Pecci, G. Aupetit-Diallo, M. Albert, P. Vignolo, and A. Min- guzzi, Persistent currents in a strongly interacting multicom- ponent Bose gas on a ring, Comptes Rendus. Physique24, 87 (2023)

  22. [22]

    Smyrnakis, M

    J. Smyrnakis, M. Magiropoulos, N. K. Efremidis, and G. M. Kavoulakis, Persistent currents in a two-component bose–einstein condensate confined in a ring potential, Journal of Physics B: Atomic, Molecular and Optical Physics47, 215302 (2014)

  23. [23]

    Sandin, M

    P. Sandin, M. ¨Ogren, and M. Gulliksson, Numerical solution of the stationary multicomponent nonlinear schr ¨odinger equation with a constraint on the angular momentum, Phys. Rev. E93, 033301 (2016)

  24. [24]

    E. T. D. Matsushita and E. J. V . de Passos, Mixture of two ultra cold bosonic atoms confined in a ring: stability and persistent currents, Journal of Physics Communications2, 035023 (2018)

  25. [25]

    Roussou, J

    A. Roussou, J. Smyrnakis, M. Magiropoulos, N. K. Efremidis, G. M. Kavoulakis, P. Sandin, M. ¨Ogren, and M. Gulliksson, Excitation spectrum of a mixture of two bose gases confined in a ring potential with interaction asymmetry, New Journal of Physics20, 045006 (2018)