pith. sign in

arxiv: 2506.03486 · v2 · submitted 2025-06-04 · ❄️ cond-mat.quant-gas

Explicit Symplectic Integrators for Massive Point Vortex Dynamics in Binary Mixture of Bose--Einstein Condensates

Tomoki Ohsawa This is my paper

Pith reviewed 2026-05-19 11:58 UTC · model grok-4.3

classification ❄️ cond-mat.quant-gas
keywords symplectic integratorspoint vortex dynamicsBose-Einstein condensatesbinary mixturesHamiltonian systemsnumerical methodsangular momentum conservation
0
0 comments X

The pith

Explicit even-order symplectic integrators simulate massive point vortex dynamics in binary Bose-Einstein condensate mixtures while exactly preserving angular momentum.

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

The paper develops explicit symplectic integrators that reach arbitrary even orders of accuracy for the massive point vortex model of a binary Bose-Einstein condensate mixture. These methods are designed for the small-mass regime, where one component carries only a tiny fraction of the total mass and solutions become highly oscillatory unless initial momenta meet specific conditions. Standard Runge-Kutta schemes suffer large Hamiltonian drift over long times in this regime, while the new integrators keep the Hamiltonian nearly constant without drift and conserve angular momentum exactly. An asymptotic expansion of the modified Hamiltonian supplies a concrete error estimate for the second-order case.

Core claim

We construct explicit integrators of arbitrary even orders of accuracy for massive point vortex dynamics in binary mixture of Bose-Einstein condensates proposed by Richaud et al. The integrators are symplectic and preserve the angular momentum of the system exactly. In the small-mass regime the solution behaviors change significantly depending on the initial momenta and become highly oscillatory unless certain conditions are met. The standard Runge-Kutta method performs very poorly in preserving the Hamiltonian with significant drift, especially for highly oscillatory solutions, whereas our integrators nearly preserve the Hamiltonian without drifts. We also give an estimate of the error in t

What carries the argument

Explicit even-order symplectic integrators constructed directly from the Hamiltonian structure of the massive point vortex model, which enforce exact angular-momentum conservation and produce bounded Hamiltonian error.

If this is right

  • Long-time simulations of vortex motion in the binary mixture become feasible without artificial energy drift.
  • Higher even-order versions improve accuracy while retaining exact angular-momentum conservation and near-conservation of the Hamiltonian.
  • The error estimate derived from the modified Hamiltonian allows quantitative prediction of accumulated deviation over many oscillation periods.
  • The construction supplies a template for building structure-preserving integrators for other Hamiltonian vortex systems with similar small-mass features.

Where Pith is reading between the lines

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

  • The same integrator construction may extend naturally to related multi-component superfluid models that admit a Hamiltonian point-vortex description.
  • Numerical tests on concrete initial conditions with and without the special momentum conditions could directly confirm the predicted change in oscillation behavior.
  • Because the methods are explicit, they remain computationally competitive with standard schemes while adding geometric fidelity.

Load-bearing premise

The massive point vortex model proposed by Richaud et al. accurately captures the dynamics of the binary mixture in the small-mass regime.

What would settle it

Run a long-time integration of a highly oscillatory initial condition with both the second-order symplectic integrator and a high-accuracy reference solver, then verify whether angular momentum remains conserved to machine precision and the Hamiltonian exhibits no secular drift only for the symplectic scheme.

Figures

Figures reproduced from arXiv: 2506.03486 by Tomoki Ohsawa.

Figure 1
Figure 1. Figure 1: shows the time evolution of the (x, y)- coordinates of the single massive vortex, computed by the 4th-order symplectic method we shall construct below; see (20) and (21) below with n = 4. The solution with (r(0), p(0)) ∈ K exhibits only small fluctuations that are barely visible on the plot in panel (a), and seems to be dominated by the slow dynamics. On the other hand, the solution with (r(0), p(0)) ∈ K/ … view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Phase portraits on the ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Time evolution of errors in Hamiltonian [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Trajectories of massive vortex dipole; [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Time evolution of relative errors in [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
read the original abstract

We construct explicit integrators of arbitrary even orders of accuracy for massive point vortex dynamics in binary mixture of Bose--Einstein condensates proposed by Richaud et al. The integrators are symplectic and preserve the angular momentum of the system exactly. Our main focus is the small-mass regime in which the minor component of the binary mixture comprises a very small fraction of the total mass. The solution behaviors in this regime change significantly depending on the initial momenta: they are highly oscillatory unless the momenta satisfy certain conditions. The standard Runge--Kutta method performs very poorly in preserving the Hamiltonian showing a significant drift in the long run, especially for highly oscillatory solutions. On the other hand, our integrators nearly preserve the Hamiltonian without drifts. We also give an estimate of the error in the Hamiltonian by finding an asymptotic expansion of the modified Hamiltonian for our second-order integrator.

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 / 2 minor

Summary. The manuscript constructs explicit symplectic integrators of arbitrary even orders for the massive point vortex dynamics in a binary Bose-Einstein condensate mixture, following the Hamiltonian model of Richaud et al. The integrators are symplectic, preserve angular momentum exactly, and are analyzed in the small-mass regime where solutions become highly oscillatory unless initial momenta satisfy special conditions. Standard Runge-Kutta methods exhibit Hamiltonian drift, while the proposed methods nearly preserve the Hamiltonian; an asymptotic expansion of the modified Hamiltonian is derived for the second-order case to estimate the error.

Significance. If the central claims hold, the work supplies a valuable class of structure-preserving integrators for stiff, oscillatory vortex dynamics in multi-component quantum gases. Exact angular-momentum conservation and drift-free Hamiltonian preservation over long times address a practical need in numerical studies of BEC mixtures. The explicit arbitrary-even-order construction and the second-order error estimate are notable strengths that enhance reproducibility and theoretical grounding.

major comments (2)
  1. [§3] §3 (higher-order integrators): The central claim that standard composition methods (Yoshida or triple-jump) applied to the split flows produce explicit maps of arbitrary even order while preserving angular momentum exactly is load-bearing for the main result. In the small-mass regime the inter-component coupling terms couple vortex positions and momenta, and the manuscript must demonstrate explicitly that each subflow remains closed-form integrable; otherwise the “arbitrary even order” and “exact preservation” assertions do not follow.
  2. [§2.2] §2.2 (small-mass regime): The statement that solutions are “highly oscillatory unless the momenta satisfy certain conditions” is used to motivate the integrators, yet the precise conditions on initial momenta that keep the split flows integrable are not stated as an explicit hypothesis. This leaves open whether the claimed explicitness holds for generic initial data in the regime of interest.
minor comments (2)
  1. [Figure 2] Figure 2 caption: the time axis label is missing units; please add “(in units of …)” for clarity.
  2. [Eq. (7)] Notation: the symbol for the inter-component coupling strength is introduced in Eq. (7) but reused without redefinition in the integrator section; a brief reminder would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments. We address the major points below and have revised the manuscript to improve clarity on the issues raised.

read point-by-point responses

  1. Referee: [§3] §3 (higher-order integrators): The central claim that standard composition methods (Yoshida or triple-jump) applied to the split flows produce explicit maps of arbitrary even order while preserving angular momentum exactly is load-bearing for the main result. In the small-mass regime the inter-component coupling terms couple vortex positions and momenta, and the manuscript must demonstrate explicitly that each subflow remains closed-form integrable; otherwise the “arbitrary even order” and “exact preservation” assertions do not follow.

    Authors: We thank the referee for this observation. The Hamiltonian splitting is constructed so that each subflow admits an explicit closed-form solution, including the inter-component coupling terms under the small-mass approximation; this is implicit in the derivation of the second-order integrator and extends directly to the composition methods. Angular-momentum preservation holds for each subflow by rotational invariance and is therefore inherited by the composed map. To strengthen the presentation we will add an explicit verification of closed-form integrability for every subflow in the revised §3. revision: yes

  2. Referee: [§2.2] §2.2 (small-mass regime): The statement that solutions are “highly oscillatory unless the momenta satisfy certain conditions” is used to motivate the integrators, yet the precise conditions on initial momenta that keep the split flows integrable are not stated as an explicit hypothesis. This leaves open whether the claimed explicitness holds for generic initial data in the regime of interest.

    Authors: We agree that an explicit statement of the conditions improves rigor. In the revised manuscript we have added a precise hypothesis in §2.2 that specifies the initial-momentum conditions under which the split flows remain integrable and the solutions avoid the highly oscillatory regime. This makes the domain of applicability of the explicit integrators clear. revision: yes

Circularity Check

0 steps flagged

Derivation self-contained: integrators built from external Hamiltonian model

full rationale

The paper takes the massive point vortex Hamiltonian proposed by Richaud et al. as an external input and applies standard splitting techniques to construct explicit symplectic integrators of even order that preserve angular momentum. No step reduces a claimed prediction or uniqueness result to a fitted parameter, self-citation chain, or definitional tautology within the present work. The central construction relies on the given Hamiltonian structure and known composition methods (Yoshida-type or similar), which are independent of the paper's own outputs. This is the normal case of an applied numerical analysis paper whose claims rest on external model assumptions rather than internal circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work rests on the domain assumption that the binary-mixture dynamics are faithfully represented by the massive point vortex Hamiltonian of Richaud et al.; no free parameters, new entities, or ad-hoc axioms are introduced in the abstract.

axioms (1)
  • domain assumption The binary mixture of Bose-Einstein condensates in the small-mass regime is accurately described by the massive point vortex model proposed by Richaud et al.
    The integrators are built directly on this prior model for the vortex dynamics.

pith-pipeline@v0.9.0 · 5677 in / 1045 out tokens · 41387 ms · 2026-05-19T11:58:58.363665+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

37 extracted references · 37 canonical work pages

  1. [1]

    Richaud, V

    A. Richaud, V. Penna, R. Mayol, and M. Guilleumas, Phys. Rev. A 101, 013630 (2020)

    work page 2020

  2. [2]

    Richaud, V

    A. Richaud, V. Penna, and A. L. Fetter, Phys. Rev. A 103, 023311 (2021)

    work page 2021

  3. [3]

    Feynman, in Progress in Low Temperature Physics , Vol

    R. Feynman, in Progress in Low Temperature Physics , Vol. 1, edited by C. Gorter (Elsevier, 1955) pp. 17–53

    work page 1955

  4. [4]

    Onsager, Il Nuovo Cimento (1943-1954) 6, 279 (1949)

    L. Onsager, Il Nuovo Cimento (1943-1954) 6, 279 (1949)

    work page 1943

  5. [5]

    Kim and A

    J.-k. Kim and A. L. Fetter, Phys. Rev. A 70, 043624 (2004)

    work page 2004

  6. [6]

    V. M. P´ erez-Garc´ ıa, H. Michinel, J. I. Cirac, M. Lewen- stein, and P. Zoller, Phys. Rev. Lett. 77, 5320 (1996)

    work page 1996

  7. [7]

    Sanz-Serna and M

    J. Sanz-Serna and M. Calvo, Numerical Hamiltonian Problems, Dover Books on Mathematics (Dover Publi- cations, 2018)

    work page 2018

  8. [8]

    Leimkuhler and S

    B. Leimkuhler and S. Reich, Simulating Hamiltonian dy- namics, Cambridge Monographs on Applied and Com- putational Mathematics, Vol. 14 (Cambridge University Press, Cambridge, 2004)

    work page 2004

  9. [9]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner, Geometric Numer- ical Integration: Structure-Preserving Algorithms for Or- dinary Differential Equations , 2nd ed. (Springer, Berlin, Heidelberg, 2006)

    work page 2006

  10. [10]

    Ohsawa and A

    T. Ohsawa and A. Richaud, arXiv:2503.19222

    work page arXiv

  11. [11]

    Hairer, S

    E. Hairer, S. Nørsett, and G. Wanner, Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems, Solving Ordinary Differential Equations II: Stiff and Differential-algebraic Problems (Springer, 1993)

    work page 1993

  12. [12]

    J. B. Sturgeon and B. B. Laird, The Journal of Chemi- cal Physics, The Journal of Chemical Physics 112, 3474 (2000)

    work page 2000

  13. [13]

    Blanes, Physical Review E 65, 056703 (2002)

    S. Blanes, Physical Review E 65, 056703 (2002)

    work page 2002

  14. [14]

    Y. K. Wu, E. Forest, and D. S. Robin, Physical Review E 68, 046502 (2003)

    work page 2003

  15. [15]

    R. I. McLachlan and G. R. W. Quispel, BIT Numerical Mathematics 44, 515 (2004)

    work page 2004

  16. [16]

    S. A. Chin, Physical Review E 80, 037701 (2009)

    work page 2009

  17. [17]

    Tao, Journal of Computational Physics 327, 245 (2016)

    M. Tao, Journal of Computational Physics 327, 245 (2016)

    work page 2016

  18. [18]

    Y. Wang, W. Sun, F. Liu, and X. Wu, The Astrophysical Journal 907, 66 (2021)

    work page 2021

  19. [19]

    Y. Wang, W. Sun, F. Liu, and X. Wu, The Astrophysical Journal 909, 22 (2021)

    work page 2021

  20. [20]

    Y. Wang, W. Sun, F. Liu, and X. Wu, The Astrophysical Journal Supplement Series 254, 8 (2021)

    work page 2021

  21. [21]

    X. Wu, Y. Wang, W. Sun, and F. Liu, The Astrophysical Journal 914, 63 (2021)

    work page 2021

  22. [22]

    X. Wu, Y. Wang, W. Sun, F.-Y. Liu, and W.-B. Han, The Astrophysical Journal 940, 166 (2022)

    work page 2022

  23. [23]

    Pihajoki, Celestial Mechanics and Dynamical Astron- omy 121, 211 (2015)

    P. Pihajoki, Celestial Mechanics and Dynamical Astron- omy 121, 211 (2015)

    work page 2015

  24. [24]

    Tao, Physical Review E 94, 043303 (2016)

    M. Tao, Physical Review E 94, 043303 (2016)

    work page 2016

  25. [25]

    Jayawardana and T

    B. Jayawardana and T. Ohsawa, Mathematics of Com- putation 92, 251 (2023)

    work page 2023

  26. [26]

    Ohsawa, SIAM Journal on Numerical Analysis 61, 1293 (2023)

    T. Ohsawa, SIAM Journal on Numerical Analysis 61, 1293 (2023)

    work page 2023

  27. [27]

    Strang, SIAM Journal on Numerical Analysis , SIAM Journal on Numerical Analysis 5, 506 (1968)

    G. Strang, SIAM Journal on Numerical Analysis , SIAM Journal on Numerical Analysis 5, 506 (1968)

    work page 1968

  28. [28]

    Creutz and A

    M. Creutz and A. Gocksch, Phys. Rev. Lett.63, 9 (1989)

    work page 1989

  29. [29]

    Forest, AIP Conference Proceedings, AIP Conference Proceedings 184, 1106 (1989)

    E. Forest, AIP Conference Proceedings, AIP Conference Proceedings 184, 1106 (1989)

    work page 1989

  30. [30]

    Suzuki, Physics Letters A 146, 319 (1990)

    M. Suzuki, Physics Letters A 146, 319 (1990)

    work page 1990

  31. [31]

    Yoshida, Physics Letters A 150, 262 (1990)

    H. Yoshida, Physics Letters A 150, 262 (1990)

    work page 1990

  32. [32]

    Tang, Computers & Mathematics with Applica- tions 27, 31 (1994)

    Y.-F. Tang, Computers & Mathematics with Applica- tions 27, 31 (1994)

    work page 1994

  33. [33]

    Benettin and A

    G. Benettin and A. Giorgilli, Journal of Statistical Physics 74, 1117 (1994)

    work page 1994

  34. [34]

    Reich, SIAM Journal on Numerical Analysis , SIAM Journal on Numerical Analysis 36, 1549 (1999)

    S. Reich, SIAM Journal on Numerical Analysis , SIAM Journal on Numerical Analysis 36, 1549 (1999)

    work page 1999

  35. [35]

    D’Ambroise, W

    J. D’Ambroise, W. Wang, C. Ticknor, R. Carretero- Gonz´ alez, and P. G. Kevrekidis, Phys. Rev. E 111, 034216 (2025)

    work page 2025

  36. [36]

    Bellettini, A

    A. Bellettini, A. Richaud, and V. Penna, The European Physical Journal Plus 138, 676 (2023)

    work page 2023

  37. [37]

    Richaud, P

    A. Richaud, P. Massignan, V. Penna, and A. L. Fetter, Phys. Rev. A 106, 063307 (2022)

    work page 2022