pith. sign in

arxiv: 1907.01194 · v1 · pith:EUWDDMFVnew · submitted 2019-07-02 · 🧮 math.NA · cs.NA

Ground states and their characterization of spin-F Bose-Einstein condensates

Tonghua Tian , Yongyong Cai , Xinming Wu , Zaiwen Wen This is my paper

Pith reviewed 2026-05-25 11:14 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Bose-Einstein condensatesground statesmanifold optimizationFourier pseudospectralspin-F systemsnumerical methodsNewton method
0
0 comments X

The pith

An optimization method on a constrained manifold computes ground states of spin-F Bose-Einstein condensates for any integer spin.

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

The paper formulates ground-state computation for spin-F Bose-Einstein condensates as minimization of an energy functional subject to two quadratic constraints. It discretizes the functional and constraints with Fourier pseudospectral schemes, then treats the discrete problem as optimization on the manifold defined by those constraints. Three retractions are constructed to return points to the manifold after each step, which permits an adaptive regularized Newton method accelerated by cascadic multigrid. This combination yields the first practical algorithm that works for arbitrary integer spin, including the spin-3 case, and is tested on one-, two-, and three-dimensional systems with varying interactions and lattice potentials.

Core claim

The discretized energy minimization problem for spin-F BECs is solved by viewing it as an optimization on the manifold of the two quadratic constraints and designing three retractions that enable Newton and multigrid iterations to converge to the ground states, providing the first applicable method for any integer spin.

What carries the argument

Three types of retractions that map points back onto the manifold defined by the two quadratic constraints after each optimization step.

If this is right

  • Ground states become computable for spin-3 BECs in one, two, and three dimensions under diverse interaction strengths and optical lattices.
  • The same discretization and retraction framework applies without change to spin-1 and spin-2 cases.
  • Cascadic multigrid acceleration reduces iteration counts for large-scale three-dimensional problems.
  • Numerical output directly exhibits physical phenomena such as phase separation or magnetization patterns.

Where Pith is reading between the lines

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

  • The retraction construction could be reused for other energy functionals that share the same two quadratic constraints.
  • Extending the method from stationary states to imaginary-time propagation would allow study of spinor dynamics.
  • Systematic comparison of the three retraction types on the same problem could identify which one is most robust for high-spin cases.

Load-bearing premise

The three designed retractions correctly return points to the manifold after each step and the resulting sequence converges to the physically relevant ground state rather than a local minimum or saddle.

What would settle it

A spin-3 computation in which the final state violates one of the two quadratic constraints or yields an energy higher than an independently verified lower bound.

read the original abstract

The computation of the ground states of spin-$F$ Bose-Einstein condensates (BECs) can be formulated as an energy minimization problem with two quadratic constraints. We discretize the energy functional and constraints using the Fourier pseudospectral schemes and view the discretized problem as an optimization problem on manifold. Three different types of retractions to the manifold are designed. They enable us to apply various optimization methods on manifold to solve the problem. Specifically, an adaptive regularized Newton method is used together with a cascadic multigrid technique to accelerate the convergence. According to our limited knowledege, our method is the first applicable algorithm for BECs with an arbitrary integer spin, including the complicated spin-3 BECs. Extensive numerical results on ground states of spin-1, spin-2 and spin-3 BECs with diverse interaction and optical lattice potential in one/two/three dimensions are reported to show the efficiency of our method and to demonstrate some interesting physical phenomena.

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 paper discretizes the energy minimization problem for spin-F Bose-Einstein condensates (with two quadratic constraints) via Fourier pseudospectral schemes, reformulates the discrete problem as manifold-constrained optimization, designs three retractions onto the manifold, and solves it with an adaptive regularized Newton method accelerated by cascadic multigrid. Numerical experiments are presented for spin-1/2/3 cases in 1D/2D/3D with various interactions and potentials; the central claim is that the approach is the first applicable algorithm for arbitrary integer F, including the complicated spin-3 case.

Significance. If the retractions are shown to enforce both constraints exactly and the iterates reliably reach global minimizers, the method would supply a practical, extensible computational framework for ground-state computation in higher-spin BECs where analytic characterizations are unavailable, complementing existing techniques limited to lower F.

major comments (2)
  1. [Retraction design] The retraction construction (described after the manifold formulation): the manuscript must supply an algebraic or numerical verification that each of the three retractions maps exactly onto the manifold defined by both quadratic constraints simultaneously, to machine precision, for the F=3 case; without this, the higher-dimensional manifold admits the possibility of drift that would invalidate the optimization steps.
  2. [Adaptive regularized Newton and cascadic multigrid] The optimization and multigrid sections: no convergence theorem, a priori error bound, or comparison against known analytic ground states (e.g., for spin-1) is provided to guarantee that the computed stationary points are global minimizers rather than local minima or saddles; this verification is load-bearing for the claim of applicability to arbitrary F.
minor comments (2)
  1. [Abstract] Abstract: 'knowledege' is a typographical error.
  2. [Discretization] The manuscript should include a brief statement of how the two quadratic constraints are discretized and preserved at the discrete level before the retraction step is introduced.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and will revise the paper accordingly.

read point-by-point responses

  1. Referee: [Retraction design] The retraction construction (described after the manifold formulation): the manuscript must supply an algebraic or numerical verification that each of the three retractions maps exactly onto the manifold defined by both quadratic constraints simultaneously, to machine precision, for the F=3 case; without this, the higher-dimensional manifold admits the possibility of drift that would invalidate the optimization steps.

    Authors: We agree that explicit verification strengthens the presentation. The retractions were constructed to satisfy both quadratic constraints exactly (one via explicit rescaling that solves the two-norm equations simultaneously, the others via Householder-style reflections and a projected Newton step that preserve the manifold by design). In the revision we will add the algebraic verification for general F together with numerical checks confirming satisfaction to machine precision on F=3 examples. revision: yes

  2. Referee: [Adaptive regularized Newton and cascadic multigrid] The optimization and multigrid sections: no convergence theorem, a priori error bound, or comparison against known analytic ground states (e.g., for spin-1) is provided to guarantee that the computed stationary points are global minimizers rather than local minima or saddles; this verification is load-bearing for the claim of applicability to arbitrary F.

    Authors: A rigorous global-convergence theorem is not provided because the energy is non-convex on a non-Euclidean manifold; such analysis lies outside the scope of this computational-methods paper. Applicability to arbitrary F is instead demonstrated through extensive numerical experiments, including recovery of known spin-1 ground states. In revision we will add direct comparisons to analytic spin-1 solutions and additional multi-start tests for F=2,3 to further support robustness against local minima. revision: partial

Circularity Check

0 steps flagged

No significant circularity: algorithmic construction is self-contained

full rationale

The paper presents a discretization via Fourier pseudospectral schemes followed by manifold optimization using three explicitly designed retractions to enforce the two quadratic constraints. These steps constitute an original algorithmic procedure built from standard numerical tools; no parameter is fitted to data and then relabeled as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz or renaming reduces the central claim to its own inputs. The claim of being the first applicable method for arbitrary integer spin is a statement of novelty rather than a derived result that loops back on itself. The derivation chain therefore remains independent of the outputs it produces.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The approach rests on standard assumptions of spectral methods and manifold optimization rather than new physical postulates or fitted constants.

axioms (2)
  • domain assumption Fourier pseudospectral schemes produce a consistent discretization of the energy functional and the two quadratic constraints for the spinor wave function.
    Invoked when the continuous problem is replaced by its discrete counterpart.
  • domain assumption The discretized problem lies on a manifold that admits the three designed retractions as valid maps preserving the constraints.
    Required to apply manifold optimization methods.

pith-pipeline@v0.9.0 · 5708 in / 1587 out tokens · 33173 ms · 2026-05-25T11:14:16.130006+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

35 extracted references · 35 canonical work pages

  1. [1]

    Absil, C

    P.-A. Absil, C. Baker, and K. Gallivan , Trust-region methods on Riemannian manifolds , Found. Comput. Math., 7 (2007), pp. 303–330

    work page 2007

  2. [2]

    Absil, R

    P.-A. Absil, R. Mahony, and R. Sepulchre , Optimization algorithms on matrix manifolds , Princeton University Press, Princeton, NJ, 2008

    work page 2008

  3. [3]

    Absil and J

    P.-A. Absil and J. Malick , Projection-like retractions on matrix manifolds, SIAM J. Optim., 22 (2012), pp. 135– 158

    work page 2012

  4. [4]

    Anderson , Theory of the weakly interacting Bose gas , Rev

    J. Anderson , Theory of the weakly interacting Bose gas , Rev. Mod. Phys., 76 (2004), pp. 599–639

    work page 2004

  5. [5]

    Anderson, J

    M. Anderson, J. Ensher, M. Mattews, C. Wieman, and E. Cornell , Observation of Bose-Einstein condensa- tion in a dilute atomic vapor , Science, 269 (1995), pp. 198–201. GROUND STATES AND THEIR CHARACTERIZATION OF SPIN-F BOSE-EINSTEIN CONDENSATES 25 Fig. 3. Contour plots for the wave functions of the ground state, i.e. φ2(x,y, 0),φ1(x,y, 0),φ0(x,y, 0),φ−1(...

    work page 1995

  6. [6]

    Bao , Ground states and dynamics of multicomponent Bose-Einstein condensates , Multiscale Model

    W. Bao , Ground states and dynamics of multicomponent Bose-Einstein condensates , Multiscale Model. Simul., 2 (2004), pp. 210–236

    work page 2004

  7. [7]

    Bao, I.-L

    W. Bao, I.-L. Chern, and Y. Zhang , Efficient numerical methods for computing ground states of spin-1 Bose- Einstein condensates based on their characterizations , J. Comput. Phys., 253 (2013), pp. 189–208

    work page 2013

  8. [8]

    Bao and Q

    W. Bao and Q. Du , Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow, SIAM J. Sci. Comput., 25 (2004), pp. 1674–1697

    work page 2004

  9. [9]

    Bao and F

    W. Bao and F. Lim , Computing ground states of spin-1 Bose-Einstein condensates by the normalized gradient flow, SIAM J. Sci. Comput., 30 (2008), pp. 1925–1948

    work page 2008

  10. [10]

    W. Bao, Q. Tang, and Y. Yuan , Computing ground states of spin-2 Bose-Einstein condensates by the normalized gradient flow. preprint

    work page

  11. [11]

    Bao and H

    W. Bao and H. W ang , A mass and magnetization conservative and energy-diminishing numerical method for computing ground state of spin-1 Bose-Einstein condensates , SIAM J. Numer. Anal., 45 (2007), pp. 2177– 2200

    work page 2007

  12. [12]

    Barrett, J

    M. Barrett, J. Sauer, and M. Chapman , All-optical formation of an atomic Bose-Einstein condensate , Phys. Rev. Lett., 87 (2001), p. 010404

    work page 2001

  13. [13]

    Bornemann and P

    F. Bornemann and P. Deuflhard , The cascadic multigrid method for elliptic problems, Numer. Math., 75 (1996), pp. 135–152

    work page 1996

  14. [14]

    Bradley, C

    C. Bradley, C. Sackett, J. Tollett, and R. Hulet , Evidence of Bose-Einstein condensation in an atomic gas with attractive interations, Phys. Rev. Lett., 75 (1995), pp. 1687–1690

    work page 1995

  15. [15]

    Dalfovo, S

    F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari , Theory of Bose-Einstein condensation in trapped gases, Rev. Mod. Phys., 71 (1999), pp. 463–512

    work page 1999

  16. [16]

    Davis, M

    K. Davis, M. Mewes, M. Andrews, N. van Druten, D. Durfee, D. Kurn, and W. Ketterle , Bose-Einstein 26 T. TIAN, Y. CAI, X. WU, AND Z. WEN Fig. 4 . Wave functions of the ground state, i.e. φ3(x), φ2(x), φ1(x), φ0(x), φ−1(x), φ−2(x) and φ−3(x) of a spin-3 BEC for Case I in 1D under different magnetizations M = 0, 0.5, 1.5. Fig. 5 . Contour plots for the wave ...

    work page 1995

  17. [17]

    Fetter , Rotating trapped Bose-Einstein condensates, Rev

    A. Fetter , Rotating trapped Bose-Einstein condensates, Rev. Mod. Phys., 81 (2009), pp. 647–691

    work page 2009

  18. [18]

    Gautam and S

    S. Gautam and S. K. Adhikari , Analytic models for the density of a ground-state spinor condensate , Phys. Rev. A, 92 (2015), p. 023616

    work page 2015

  19. [19]

    G ¨orlitz, T

    A. G ¨orlitz, T. L. Gustavson, A. E. Leanhardt, R. L ¨ow, A. P. Chikkatur, S. Gupta, S. Inouye, D. E. Pritchard, and W. Ketterle , Sodium Bose-Einstein condensates in the F = 2 state in a large-volume optical trap, Phys. Rev. Lett., 90 (2003), p. 090401

    work page 2003

  20. [20]

    Ho , Spinor Bose condensates in optical traps , Phys

    T. Ho , Spinor Bose condensates in optical traps , Phys. Rev. Lett., 81 (1998), pp. 742–745

    work page 1998

  21. [21]

    J. Hu, A. Milzarek, Z. Wen, and Y. Yuan , Adaptive quadratically regularized Newton method for Riemannian optimization, SIAM J. Matrix Anal. Appl., 39 (2018), pp. 1181–1207

    work page 2018

  22. [22]

    C. Law, H. Pu, and N. Bigelow , Quantum spins mixing in spinor Bose-Einstein condensates , Phys. Rev. Lett., 81 (1998), pp. 5257–5261

    work page 1998

  23. [23]

    Leggett , Bose-Einstein condensation in the alkali gases: Some fundamental concepts , Rev

    A. Leggett , Bose-Einstein condensation in the alkali gases: Some fundamental concepts , Rev. Mod. Phys., 73 (2001), pp. 307–356

    work page 2001

  24. [24]

    Miesner, D

    H. Miesner, D. Stamper-Kurn, J. Stenger, S. Inouye, A. Chikkatur, and W. Ketterle , Observation of metastable states in spinor Bose-Einstein condensates , Phys. Rev. Lett., 82 (1999), pp. 2228–2231

    work page 1999

  25. [25]

    Morsch and M

    O. Morsch and M. Oberthaler , Dynamics of Bose-Einstein condensates in optical lattices , Rev. Mod. Phys., GROUND STATES AND THEIR CHARACTERIZATION OF SPIN-F BOSE-EINSTEIN CONDENSATES 27 Fig. 6 . Contour plots for the wave functions of the ground state, i.e. φ3(x,y, 0), φ2(x,y, 0), φ1(x,y, 0), φ0(x,y, 0), φ−1(x,y, 0), φ−2(x,y, 0), φ−3(x,y, 0) of a spin-3 ...

    work page 2006

  26. [26]

    Ohmi and K

    T. Ohmi and K. Machida , Bose-Einstein condensation with internal degrees of freedom in alkali atom gases , J. Phys. Soc. Japan, 67 (1998), pp. 1822–1825

    work page 1998

  27. [27]

    Ozeri, N

    R. Ozeri, N. Katz, J. Steinhauer, and N. Davidson , Colloquium: Bulk Bogoliubov excitations in a Bose-Einstein condensate, Rev. Mod. Phys., 77 (2005), pp. 187–205

    work page 2005

  28. [28]

    Posazhennikova , Colloquium: Weakly interacting, dilute Bose gases in 2D , Rev

    A. Posazhennikova , Colloquium: Weakly interacting, dilute Bose gases in 2D , Rev. Mod. Phys., 78 (2006), pp. 1111–1134

    work page 2006

  29. [29]

    Stamper-Kurn, M

    D. Stamper-Kurn, M. Andrews, A. Chikkatur, S. Inouye, H. Miesner, J. Stenger, and W. Ketterle , Optical confinement of a Bose-Einstein condensate , Phys. Rev. Lett., 80 (1998), pp. 2027–2030

    work page 1998

  30. [30]

    Stamper-Kurn and W

    D. Stamper-Kurn and W. Ketterle , Spinor condensates and light scattering from Bose-Einstein condensates , in Coherent atomic matter waves, 2001, pp. 139–217

    work page 2001

  31. [31]

    Stenger, S

    J. Stenger, S. Inouye, D. Stamper-Kurn, H. Miesner, A. Chikkatur, and W. Ketterle , Spin domains in ground state Bose-Einstein condensates , Nature, 396 (1998), pp. 345–348

    work page 1998

  32. [32]

    W ang, A projection gradient method for computing ground state of spin-2 Bose-Einstein condensates , J

    H. W ang, A projection gradient method for computing ground state of spin-2 Bose-Einstein condensates , J. Com- put. Phys., 274 (2014), pp. 473–488

    work page 2014

  33. [33]

    Yang, L.-H

    W. Yang, L.-H. Zhang, and R. Song , Optimality conditions for the nonlinear programming problems on Rie- mannian manifolds, Pac. J. Optim., 10 (2014), pp. 415–434

    work page 2014

  34. [34]

    Zhang, S

    W. Zhang, S. Yi, and L. You , Mean field ground state of a spin-1 condensate in a magnetic field , New J. Phys., 5 (2003), pp. 77–89

    work page 2003

  35. [35]

    Zhang and L

    W. Zhang and L. You , An effective quasi-one-dimensional description of a spin-1 atomic condensate , Phys. Rev. A, 71 (2005), p. 025603

    work page 2005