Jones-Roberts solitary waves and the onset of rotation in a spherical surface condensate
Pith reviewed 2026-05-19 23:58 UTC · model grok-4.3
The pith
In a spherical-shell Bose-Einstein condensate, Jones-Roberts solitary waves mark the onset of rotation, evolving from vortex dipoles to hybrid equatorial modes whose speed limits the entire family via a Landau critical velocity.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The propagation speed of equatorially confined modes plays the role of a Landau critical velocity, thereby setting the upper limiting angular speed of the entire Jones-Roberts family.
Load-bearing premise
The analysis assumes that the nonlinear excitations remain solitary waves that rotate at strictly constant angular speed on an idealized thin spherical shell without significant radial leakage or damping.
Figures
read the original abstract
The nonlinear excitations underlying the onset of rotation in a dilute Bose-Einstein condensate confined to a thin spherical shell are studied. These excitations correspond to solitary waves rotating about the sphere at constant angular speed: at low speeds they appear as dipoles of singly quantized vortices with opposite circulation, while at higher speeds they evolve into vortex-free Jones-Roberts solitons. With further increase of the angular speed, these excitations hybridize with equatorially confined modes whose azimuthal wave number is set by the sphere radius measured in units of the healing length. The propagation speed of these modes is shown to play the role of a Landau critical velocity, thereby setting the upper limiting angular speed of the entire Jones-Roberts family.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies nonlinear excitations in a dilute Bose-Einstein condensate confined to a thin spherical shell. These are solitary waves rotating at constant angular speed: low-speed states appear as oppositely circulating singly-quantized vortex dipoles that evolve into vortex-free Jones-Roberts solitons at higher speeds. With further increase in angular speed the excitations hybridize with equatorially confined modes whose azimuthal wave number is fixed by the sphere radius in healing-length units; the propagation speed of these modes is identified as a Landau critical velocity that caps the upper angular speed of the entire Jones-Roberts family.
Significance. If the central identification holds, the work supplies a concrete mechanism linking the onset of rotation in spherical-shell condensates to the termination of a solitary-wave branch at a Landau velocity. This supplies a falsifiable prediction for the maximum sustainable angular speed and connects Jones-Roberts physics on curved surfaces to established critical-velocity concepts, which is of interest for both theory and future experiments on shell-trapped BECs.
major comments (2)
- [Section on hybridization and critical velocity (near the discussion of Landau velocity)] The claim that equatorial-mode speed sets the strict upper limit for the Jones-Roberts family is load-bearing for the abstract and conclusion. The manuscript must show either (i) an explicit dispersion relation for small-amplitude equatorial waves (with radius in healing lengths fixing the azimuthal number) or (ii) numerical continuation data demonstrating that no higher-speed branches exist beyond this velocity. Without one of these, the termination-by-hybridization argument remains an assumption rather than a demonstrated result.
- [Numerical methods and results section] The numerical or analytic continuation procedure that tracks the solitary-wave family from the vortex-dipole regime through the Jones-Roberts regime and into the hybridization regime is central. Convergence checks, error estimates on the angular speed at which hybridization occurs, and a direct comparison of that speed with the independently computed equatorial-mode speed are required to substantiate the limiting-velocity identification.
minor comments (2)
- [Abstract] The abstract states the qualitative transitions clearly but does not indicate the numerical or analytic method used to obtain the family; a single sentence on the approach would improve readability.
- [Figure captions] Figures displaying the density and phase profiles at successive angular speeds should include the value of the predicted Landau velocity for direct visual comparison.
Simulated Author's Rebuttal
We thank the referee for their thorough review and valuable suggestions. We have carefully considered the major comments and have made revisions to address the concerns regarding the demonstration of the limiting velocity and the numerical procedures. We believe these changes strengthen the manuscript and clarify the key results.
read point-by-point responses
-
Referee: The claim that equatorial-mode speed sets the strict upper limit for the Jones-Roberts family is load-bearing for the abstract and conclusion. The manuscript must show either (i) an explicit dispersion relation for small-amplitude equatorial waves (with radius in healing lengths fixing the azimuthal number) or (ii) numerical continuation data demonstrating that no higher-speed branches exist beyond this velocity. Without one of these, the termination-by-hybridization argument remains an assumption rather than a demonstrated result.
Authors: We agree that an explicit demonstration is necessary to substantiate the claim. In the revised version, we have included a derivation of the dispersion relation for small-amplitude equatorial waves on the sphere. The azimuthal wave number m is determined by the condition that the wavelength fits the circumference, specifically m ≈ R / ξ where R is the sphere radius and ξ the healing length. The phase velocity from this dispersion matches the speed at which the Jones-Roberts branch terminates in our numerical continuations. This is now presented in a new subsection with supporting analytical expressions and numerical verification. revision: yes
-
Referee: The numerical or analytic continuation procedure that tracks the solitary-wave family from the vortex-dipole regime through the Jones-Roberts regime and into the hybridization regime is central. Convergence checks, error estimates on the angular speed at which hybridization occurs, and a direct comparison of that speed with the independently computed equatorial-mode speed are required to substantiate the limiting-velocity identification.
Authors: We have expanded the numerical methods section to detail the continuation procedure, including the use of Newton-Krylov solvers with adaptive mesh refinement. Convergence has been checked by varying the spatial resolution (from 128x128 to 512x512 grid points) and solver tolerances (down to 10^{-10}). The angular speed at hybridization is determined with an estimated error of less than 0.5%, and we now include a direct comparison showing agreement with the equatorial mode speed to within this margin. A new figure compares the two speeds explicitly. revision: yes
Circularity Check
No significant circularity detected
full rationale
The paper's central claim is that hybridization with equatorially confined modes occurs and that their propagation speed functions as a Landau critical velocity capping the Jones-Roberts family. This is presented as a derived result from studying nonlinear excitations on the spherical shell, not as a self-definition or a parameter fitted to the target quantity itself. No equations or self-citations are visible in the provided text that reduce the upper angular-speed limit directly to an input by construction (e.g., no fitted dispersion relation renamed as a prediction). The thin-shell and constant-angular-speed assumptions are modeling choices whose validity can be checked externally against the underlying Gross-Pitaevskii dynamics; they do not create a tautological loop. The derivation therefore remains self-contained and independent of the result it claims to establish.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The condensate is dilute and can be described by the Gross-Pitaevskii equation on a thin spherical shell.
- domain assumption Nonlinear excitations exist that rotate at constant angular speed.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The propagation speed of these modes is shown to play the role of a Landau critical velocity, thereby setting the upper limiting angular speed of the entire Jones-Roberts family.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
ωBog(l) = ± sqrt[ l(l+1)/R² (l(l+1)/R² + 2) ]
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
-
[1]
Initial condition Our numerical analysis begins with the initial condition for a pair of vortices on the surface of the sphere. In the flat plane, the phase profile of a vortex-antivortex pair, situated at posi- tions (x+,y +) and (x−,y −), takes the well-known form S pl(x,y)=atan2(y−y +,x−x +)−atan2(y−y −,x−x −),(C1) where atan2 is the two-argument arcta...
-
[2]
Spectral methods The order parameterψis represented using spherical har- monicsY l,m(θ, ϕ) with complex coefficientscl,m, ψ(ϕ, θ)= lmaxX l=0 lX m=−l cl,mYl,m(ϕ, θ),(C4) whose indices are truncated to 0≤l≤l max and|m| ≤l max for numerical purposes. Two different numerical libraries have been used in the numerical simulations, both leading to the same resul...
-
[3]
δψ ∂ψ∗ sw ∂ϕ +c.c. # dAϕ,θ =−iΩR 2 Z
Newton-Raphson method Our goal is to find a solution of 0=iΩ ∂ψsw ∂ϕ − 1 R2 ∆ϕ,θψsw +|ψ sw|2ψsw −ψ sw .(C8) In the following we drop the index SW for simplicity. We define the complex functionalFfor the Newton-Raphson method using the previous equation, resulting in F(Reψ,Imψ)=iΩ ∂ψ ∂ϕ − 1 R2 ∆ϕ,θψ+|ψ| 2ψ−ψ .(C9) Given the root of the functionalFbeingf ro...
-
[4]
for more details. In the following we chooseθ ξ =1/R, i.e., a core arc-length radius ofRθ ξ =1, in units of the healing length. By taking the limitR≫1, where one expects the Rankine vortex solution to work better, we can show that lim R→∞ Lz =2ρΓR 2 +O (1) (E4) and lim R→∞ Ekin = ρΓ2 1+4 log (2R) 8π +O R−2 .(E5) Hence, we find that the angular momentum sc...
- [5]
- [6]
-
[7]
Y . Colombe, E. Knyazchyan, O. Morizot, B. Mercier, V . Lorent, and H. Perrin, EPL67, 593 (2004)
work page 2004
-
[8]
B. M. Garraway and H. Perrin, J. Phys. B: At. Mol. Opt. Phys. 49, 172001 (2016)
work page 2016
-
[9]
Y . Guo, E. M. Gutierrez, D. Rey, T. Badr, A. Perrin, L. Longchambon, V . S. Bagnato, H. Perrin, and R. Dubessy, New J. Phys.24, 093040 (2022)
work page 2022
-
[10]
D. Rey, S. Thomas, R. Sharma, T. Badr, L. Longchambon, R. Dubessy, and H. Perrin, J. Appl. Phys.132, 214401 (2022)
work page 2022
-
[11]
D. Aveline, J. Williams, E. Elliott, C. Dutenhoffer, J. R. Kel- logg, J. M. Kohel, N. E. Lay, K. Oudrhiri, R. F. Shotwell, N. Yu, and R. J. Thompson, Nature582, 193–197 (2020)
work page 2020
-
[12]
R. A. Carollo, D. C. Aveline, B. Rhyno, S. Vishveshwara, C. Lannert, J. D. Murphree, E. R. Elliott, J. R. Williams, R. J. Thompson, and N. Lundblad, Nature606, 281–286 (2022)
work page 2022
-
[13]
K. Frye, S. Abend, W. Bartosch, A. Bawamia, D. Becker, H. Blume, C. Braxmaier, S.-W. Chiow, M. A. Efremov, W. Ert- mer, P. Fierlinger, T. Franz, N. Gaaloul, J. Grosse, C. Grzeschik, O. Hellmig, V . A. Henderson, W. Herr, U. Israelsson, J. Kohel, M. Krutzik, C. Kürbis, C. Lämmerzahl, M. List, D. Lüdtke, N. Lundblad, J. P. Marburger, M. Meister, M. Mihm, H....
work page 2021
-
[14]
D. Becker, M. D. Lachmann, S. T. Seidel, H. Ahlers, A. N. Dinkelaker, J. Grosse, O. Hellmig, H. Müntinga, V . Schkolnik, T. Wendrich, A. Wenzlawski, B. Weps, R. Corgier, T. Franz, N. Gaaloul, W. Herr, D. Lüdtke, M. Popp, S. Amri, H. Duncker, M. Erbe, A. Kohfeldt, A. Kubelka-Lange, C. Braxmaier, E. Charron, W. Ertmer, M. Krutzik, C. Lämmerzahl, A. Peters, ...
work page 2018
-
[15]
T. van Zoest, N. Gaaloul, Y . Singh, H. Ahlers, W. Herr, S. T. Seidel, W. Ertmer, E. Rasel, M. Eckart, E. Kajari, S. Arnold, G. Nandi, W. P. Schleich, R. Walser, A. V ogel, K. Seng- stock, K. Bongs, W. Lewoczko-Adamczyk, M. Schiemangk, T. Schuldt, A. Peters, T. Könemann, H. Müntinga, C. Läm- merzahl, H. Dittus, T. Steinmetz, T. W. Hänsch, and J. Reichel, ...
work page 2010
-
[16]
F. Jia, Z. Huang, L. Qiu, R. Zhou, Y . Yan, and D. Wang, Phys. 10 Rev. Lett.129, 243402 (2022)
work page 2022
- [17]
- [18]
-
[19]
S. J. Bereta, L. Madeira, V . S. Bagnato, and M. A. Caracanhas, Am. J. Phys.87, 924 (2019)
work page 2019
- [20]
- [21]
- [22]
- [23]
- [24]
-
[25]
N. S. Móller, F. E. A. dos Santos, V . S. Bagnato, and A. Pelster, New J. Phys.22, 063059 (2020)
work page 2020
-
[26]
K. Padavi ´c, K. Sun, C. Lannert, and S. Vishveshwara, Phys. Rev. A102, 043305 (2020)
work page 2020
- [28]
- [29]
- [30]
- [31]
-
[32]
S. J. Bereta, M. A. Caracanhas, and A. L. Fetter, Phys. Rev. A 103, 053306 (2021)
work page 2021
-
[33]
N.-E. Guenther, P. Massignan, and A. L. Fetter, Phys. Rev. A 96, 063608 (2017)
work page 2017
- [34]
-
[35]
M. A. Caracanhas, P. Massignan, and A. L. Fetter, Phys. Rev. A 105, 023307 (2022)
work page 2022
-
[36]
C. A. Jones and P. H. Roberts, J. Phys. A: Math. Gen.15, 2599 (1982)
work page 1982
- [37]
-
[38]
M. Baker-Rasooli, T. Aladjidi, N. A. Krause, A. S. Bradley, and Q. Glorieux, Phys. Rev. Lett.134, 233401 (2025)
work page 2025
-
[39]
See Supplemental Material at [URL to be inserted] for further details
- [40]
-
[41]
J. Boyd,Chebyshev and F ourier Spectral Methods: Second Re- vised Edition, Dover Books on Mathematics (Dover Publica- tions, 2013)
work page 2013
-
[42]
P. K. Newton,The N-V ortex Problem: Analytical Techniques (Springer New York, New York, NY , 2001)
work page 2001
-
[43]
N. G. Berloff, J. Phys. A: Math. Gen.37, 1617 (2004)
work page 2004
-
[44]
See video material at this URL
-
[45]
C. A. Jones, S. J. Putterman, and P. H. Roberts, J. Phys. A: Math. Gen.19, 2991 (1986)
work page 1986
-
[46]
J. D’Ambroise, R. Carretero-González, P. Schmelcher, and P. G. Kevrekidis, Phys. Rev. A105, 063325 (2022)
work page 2022
-
[47]
N. A. Krause and A. S. Bradley, Phys. Rev. A110, 053302 (2024)
work page 2024
-
[48]
Y . K. Chembo and C. R. Menyuk, Phys. Rev. A87, 053852 (2013)
work page 2013
-
[49]
T. Herr, V . Brasch, J. D. Jost, C. Y . Wang, N. M. Kondratiev, M. L. Gorodetsky, and T. J. Kippenberg, Nature Photonics8, 145 (2014)
work page 2014
-
[50]
J. Tempere, I. F. Silvera, and J. T. Devreese, Phys. Rev. B65, 195418 (2002)
work page 2002
-
[51]
J. Tempere, V . N. Gladilin, I. F. Silvera, and J. T. Devreese, Phys. Rev. B72, 094506 (2005)
work page 2005
-
[52]
V . N. Gladilin, J. Tempere, I. F. Silvera, and J. T. Devreese, Phys. Rev. B74, 104512 (2006)
work page 2006
-
[53]
J. Tempere, V . N. Gladilin, I. F. Silvera, and J. T. Devreese, J Low. Temp. Phys.148, 181 (2007)
work page 2007
-
[54]
J. Tempere, V . N. Gladilin, I. F. Silvera, J. T. Devreese, and V . V . Moshchalkov, Phys. Rev. B79, 134516 (2009)
work page 2009
- [55]
-
[56]
C. D. Brown, Y . Wang, M. Namazi, G. I. Harris, M. T. Uysal, and J. G. E. Harris, Phys. Rev. Lett.130, 216001 (2023)
work page 2023
-
[57]
C. Lannert, T.-C. Wei, and S. Vishveshwara, Phys. Rev. A75, 013611 (2007)
work page 2007
-
[58]
K. Sun, K. Padavi ´c, F. Yang, S. Vishveshwara, and C. Lannert, Phys. Rev. A98, 013609 (2018)
work page 2018
-
[59]
K. Padavi ´c, K. Sun, C. Lannert, and S. Vishveshwara, EPL120, 20004 (2017)
work page 2017
-
[60]
V . P. Ruban, JETP Letters116, 329–334 (2022)
work page 2022
-
[61]
L. Cavicchioli, C. Fort, F. Ancilotto, M. Modugno, F. Minardi, and A. Burchianti, Phys. Rev. Lett.134, 093401 (2025)
work page 2025
- [62]
-
[63]
A. C. White, Phys. Rev. A109, 013301 (2024)
work page 2024
-
[64]
M. A. Wieczorek and M. Meschede, Geochemistry, Geo- physics, Geosystems19, 2574 (2018)
work page 2018
-
[65]
Schaeffer, Geochemistry, Geophysics, Geosystems14, 751 (2013)
N. Schaeffer, Geochemistry, Geophysics, Geosystems14, 751 (2013)
work page 2013
- [66]
-
[67]
A. H. Baker, E. R. Jessup, and T. Manteuffel, SIAM J. Matrix Anal. Appl.26, 962 (2005)
work page 2005
-
[68]
See Wolfram Mathematica notebooks at this URL
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.