Structure and symmetry of the Gross-Pitaevskii ground-state manifold

Patrick Henning; Qinglin Tang; Zixu Feng

arxiv: 2603.28174 · v2 · submitted 2026-03-30 · 🧮 math.NA · cs.NA· math.OC

Structure and symmetry of the Gross-Pitaevskii ground-state manifold

Zixu Feng , Patrick Henning , Qinglin Tang This is my paper

Pith reviewed 2026-05-14 02:21 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.OC

keywords Gross-Pitaevskii equationground-state manifoldMorse-Bott conditionsymmetry orbitsRiemannian gradient methodslocal convergencenonlinear Schrödinger equation

0 comments

The pith

The Morse-Bott condition partitions the Gross-Pitaevskii ground-state manifold into finitely many symmetry orbits and marks the exact threshold for local linear convergence of preconditioned Riemannian gradient methods.

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

This paper studies the geometric structure of ground states for the Gross-Pitaevskii energy functional when multiple local minimizers exist due to symmetries. It shows that the Morse-Bott non-degeneracy condition at these minimizers ensures the ground-state set breaks into finitely many embedded submanifolds, each an orbit under phase shifts and spatial rotations. The same condition is necessary and sufficient for local Q-linear convergence of preconditioned Riemannian gradient methods, while its absence produces only sublinear rates. A sympathetic reader cares because the result directly ties the symmetry properties of the physical model to the practical performance of numerical solvers used to compute Bose-Einstein condensates.

Core claim

When local minimizers of the Gross-Pitaevskii energy are non-unique, the Morse-Bott condition ensures that the ground-state set partitions into finitely many embedded submanifolds, each coinciding with an orbit generated by the intrinsic symmetries of phase shifts and spatial rotations. This condition holds if and only if the minimizers decompose into finitely many symmetry orbits and the preconditioned Riemannian gradient method exhibits local Q-linear convergence nearby; otherwise the method converges only sublinearly. These results establish the Morse-Bott condition as the precise separator between linear and sublinear convergence while giving a symmetry-based structural description of 2D

What carries the argument

The Morse-Bott non-degeneracy condition, which requires that the Hessian of the energy functional at each local minimizer has kernel exactly spanned by the infinitesimal generators of the symmetry group actions.

If this is right

The ground-state set decomposes into finitely many embedded submanifolds, each generated by phase and rotation symmetries.
Preconditioned Riemannian gradient methods achieve local Q-linear convergence exactly when the Morse-Bott condition holds at the ground states.
When the condition fails, the same methods converge only sublinearly near the ground-state set.
Minimizers on the ground-state set form finitely many symmetry orbits if and only if the optimization method converges linearly nearby.

Where Pith is reading between the lines

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

Convergence rate monitoring could serve as a practical diagnostic to verify whether a computed ground state satisfies the Morse-Bott condition.
The orbit decomposition may suggest symmetry-reduced formulations that lower the effective dimension of the optimization problem.
Similar Morse-Bott analysis could apply to other nonlinear Schrödinger models whose energy functionals share continuous symmetry groups.

Load-bearing premise

The energy functional is twice continuously differentiable and its Hessian at local minimizers has kernel dimension exactly matching the dimension of the symmetry group.

What would settle it

A direct spectral computation of the Hessian operator at an explicit ground-state minimizer that reveals kernel dimension strictly larger than the number of independent symmetry generators, or a numerical run of the preconditioned Riemannian gradient method near that minimizer that exhibits slower-than-linear convergence.

read the original abstract

The structure and degeneracy of ground states of the Gross-Pitaevskii energy functional play a central role in both analysis and computation, yet a precise characterization of the ground-state manifold in the presence of symmetries remains a fundamental challenge. In this paper, we establish sharp theoretical results describing the geometric structure of local minimizers and its implications for optimization algorithms. We show that when local minimizers are non-unique, the Morse-Bott condition provides a natural and sufficient criterion under which the ground-state set partitions into finitely many embedded submanifolds, each coinciding with an orbit generated by intrinsic symmetries of the energy functional, namely phase shifts and spatial rotations. This yields a structural characterization of the ground-state manifold in terms of these symmetries. Building on this insight, we characterize the local convergence behavior of general preconditioned Riemannian gradient methods (P-RG). Under the Morse-Bott condition, we derive sharp local $Q$-linear convergence estimates and prove that the condition holds if and only if the energy sequence generated by P-RG converges locally $Q$-linearly. In particular, on the ground-state set, the Morse-Bott condition is satisfied if and only if the minimizers decompose into finitely many symmetry orbits and P-RG exhibits local linear convergence nearby. When the condition fails, we establish a local sublinear convergence rate. Taken together, these results show that the Morse-Bott condition is the exact threshold separating linear from sublinear convergence, while determining the symmetry-induced structure of the ground-state manifold, connecting geometry, symmetry, and convergence behavior in a unified framework.

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.

Desk Editor's Note private letter to a colleague

The paper gives a clean if-and-only-if link between the Morse-Bott condition on the GP energy, the decomposition of ground states into symmetry orbits, and Q-linear convergence of P-RG methods.

read the letter

The main point is that the Morse-Bott non-degeneracy condition on the Gross-Pitaevskii energy is necessary and sufficient for the set of local minimizers to split into finitely many embedded orbits under the natural U(1) phase and SO(d) rotation symmetries, and that the same condition is exactly the threshold for local Q-linear convergence of preconditioned Riemannian gradient methods. When it fails, the methods only achieve sublinear rates. This equivalence is the new structural result. The paper applies standard Morse-Bott theory in the infinite-dimensional setting and shows how the Hessian kernel matches the Lie algebra of the symmetry action, which organizes the manifold and controls the convergence. That connection is useful for anyone running these optimization schemes on quantum-fluid models. The derivations follow from the usual variational and Riemannian assumptions once the functional is C^2 and the group action is smooth, with no visible circularity or post-hoc restrictions in the abstract. The stress-test note found no internal gaps in the kernel identification or rate derivation. A minor soft spot is that the infinite-dimensional embedding and compactness details need to be checked in the full proofs to confirm the submanifolds are properly embedded, but nothing in the outline suggests a load-bearing flaw. This is for readers working on numerical methods for nonlinear Schrödinger equations or geometric optimization on manifolds. Anyone who cares about when symmetry-induced degeneracy produces linear versus sublinear convergence will get concrete value. It deserves a serious referee because the geometric criterion is sharp enough to matter for both theory and implementation.

Referee Report

2 major / 2 minor

Summary. The paper claims that for the Gross-Pitaevskii energy functional, the Morse-Bott non-degeneracy condition at local minimizers is necessary and sufficient for the ground-state set to consist of finitely many embedded submanifolds, each an orbit under the intrinsic symmetries (U(1) phase shifts and SO(d) rotations). It further asserts that this same condition is equivalent to local Q-linear convergence of preconditioned Riemannian gradient (P-RG) methods, with sublinear convergence when the condition fails, thereby providing a sharp geometric threshold linking symmetry-induced degeneracy to algorithmic behavior.

Significance. If the derivations hold, the results supply a clean, if-and-only-if geometric characterization of degenerate ground states that directly governs convergence rates of Riemannian optimization schemes. This unifies symmetry analysis with numerical analysis in a way that could guide both theoretical studies of nonlinear Schrödinger equations and the design of structure-preserving algorithms, especially when ground states are non-unique.

major comments (2)

[Abstract and §4 (Morse-Bott analysis)] The abstract asserts that the Morse-Bott condition holds if and only if minimizers decompose into finitely many symmetry orbits; the central proof must explicitly verify that the Hessian kernel coincides exactly with the Lie algebra of the symmetry action (phase and rotations) and contains no additional directions, without post-hoc restrictions on the potential or domain. This equivalence is load-bearing for both the structural claim and the convergence threshold.
[§5 (convergence analysis)] The claimed sharp Q-linear convergence rate for P-RG under the Morse-Bott condition requires an explicit constant (depending on the smallest positive eigenvalue of the reduced Hessian) and a precise statement of the neighborhood in which the rate holds; the current abstract statement leaves open whether the rate is uniform across all orbits or orbit-dependent.

minor comments (2)

[§3] Define the precise form of the preconditioner in the P-RG method at the first appearance (likely §3) rather than deferring it; this affects readability of the convergence statements.
[Abstract] Clarify whether the finite number of orbits is uniform or depends on the specific symmetry group dimension; add a remark on how the result specializes to d=1 or d=2.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading, positive evaluation, and constructive suggestions. The comments help strengthen the clarity of the Morse-Bott equivalence and the convergence statements. We address each major comment below and have incorporated revisions to make the required verifications and constants explicit.

read point-by-point responses

Referee: [Abstract and §4 (Morse-Bott analysis)] The abstract asserts that the Morse-Bott condition holds if and only if minimizers decompose into finitely many symmetry orbits; the central proof must explicitly verify that the Hessian kernel coincides exactly with the Lie algebra of the symmetry action (phase and rotations) and contains no additional directions, without post-hoc restrictions on the potential or domain. This equivalence is load-bearing for both the structural claim and the convergence threshold.

Authors: We agree that the kernel identification is central. In the proof of Theorem 4.2 we compute the second variation of the Gross-Pitaevskii energy directly from the Euler-Lagrange equation satisfied by any critical point. The resulting quadratic form on the tangent space is shown to vanish if and only if the variation lies in the span of the infinitesimal generators of the U(1) phase action and the SO(d) rotation action; any orthogonal component produces a strictly positive contribution controlled by the Morse-Bott gap. The argument uses only the standard regularity and growth assumptions on the potential stated in Section 2 and makes no further restrictions on its form or on the domain. We have added an explicit remark after the proof that isolates this kernel equality and confirms its validity for general potentials. revision: yes
Referee: [§5 (convergence analysis)] The claimed sharp Q-linear convergence rate for P-RG under the Morse-Bott condition requires an explicit constant (depending on the smallest positive eigenvalue of the reduced Hessian) and a precise statement of the neighborhood in which the rate holds; the current abstract statement leaves open whether the rate is uniform across all orbits or orbit-dependent.

Authors: We accept the request for sharper quantitative statements. In the revised Theorem 5.3 we now state that the contraction factor is at most 1 - c λ_min, where λ_min > 0 is the smallest positive eigenvalue of the reduced Hessian on the orthogonal complement to the symmetry tangent space and c > 0 depends only on the preconditioner norm. The neighborhood is a tubular neighborhood of radius r around each orbit, with r chosen small enough that the local chart and the Morse-Bott gap dominate the higher-order terms; both λ_min and r are orbit-dependent in general. The abstract has been updated to indicate that the linear rate holds locally near each orbit and is therefore not necessarily uniform across the ground-state set. revision: yes

Circularity Check

0 steps flagged

No significant circularity in Morse-Bott derivation

full rationale

The paper applies the standard Morse-Bott non-degeneracy condition (Hessian kernel exactly equal to the tangent space of the symmetry orbit) to the Gross-Pitaevskii energy on the manifold with U(1) phase and SO(d) rotational symmetries. The claimed equivalences (ground-state set as finite union of embedded orbits, and iff local Q-linear convergence of P-RG) are derived from the definition of the condition plus standard results in infinite-dimensional Morse theory and Riemannian gradient convergence analysis. No self-definitional steps, fitted parameters renamed as predictions, or load-bearing self-citations appear; the central claims follow directly from variational calculus without reducing to the inputs by construction. The derivation is self-contained against external benchmarks in differential geometry.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Analysis rests on standard smoothness and variational assumptions for the Gross-Pitaevskii functional together with the definition of the Morse-Bott condition; no free parameters, new entities, or ad-hoc axioms are introduced in the abstract.

axioms (1)

domain assumption The Gross-Pitaevskii energy functional is sufficiently smooth and the Hessian at critical points has kernel exactly spanned by the infinitesimal generators of the symmetry group.
Required for the Morse-Bott condition to be well-defined and for the orbit decomposition to hold.

pith-pipeline@v0.9.0 · 5589 in / 1339 out tokens · 58789 ms · 2026-05-14T02:21:14.126676+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

48 extracted references · 48 canonical work pages

[1]

Y. Ai, P. Henning, M. Yadav, and S. Yuan, Riemannian conju gate Sobolev gradients and their application to compute ground states of BECs, J. Co mput. Appl. Math., 473 (2026), article 116866

work page 2026
[2]

Altmann, P

R. Altmann, P. Henning, and D. Peterseim, The J-method for the Gross-Pitaevskii eigenvalue problem, Numer. Math., 148 (2021), pp. 575–610

work page 2021
[3]

Altmann, M

R. Altmann, M. Hermann, D. Peterseim, and T. Stykel, Riem annian optimisation meth- ods for ground states of multicomponent Bose-Einstein cond ensates, arXiv:2411.09617

work page arXiv
[4]

M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Observation of Bose-Einstein condensation in a dilute atom ic vapor, Sci., 269 (1995), pp. 198–201

work page 1995
[5]

Antoine and R

X. Antoine and R. Duboscq, Robust and eﬃcient preconditi oned Krylov spectral solvers for computing the ground states of fast rotating and strongl y interacting Bose-Einstein condensates, J. Comput. Phys., 258 (2014), pp. 509–523

work page 2014
[6]

Antoine, A

X. Antoine, A. Levitt, and Q. Tang, Eﬃcient spectral comp utation of the stationary states of rotating Bose-Einstein condensates by the precon ditioned nonlinear conjugate gradient method, J. Comput. Phys., 343 (2017), pp. 92–109. 24

work page 2017
[7]

Bao and Q

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

work page 2004
[8]

W. Bao, I. Chern, and F. Lim, Eﬃcient and spectrally accur ate numerical methods for computing ground and ﬁrst excited states in Bose-Einste in condensates, J. Comput. Phys., 219 (2006), pp. 836–854

work page 2006
[9]

Bao and Y

W. Bao and Y. Cai, Mathematical theory and numerical meth ods for Bose-Einstein condensation, Kinet. Relat. Models, 6 (2013), pp. 1–135

work page 2013
[10]

C. F. Barenghi, L. Skrbek, and K. R. Sreenivasan, Introd uction to quantum turbulence, PNAS, 111 (2014), pp. 4647–4652

work page 2014
[11]

Bott, Nondegenerate critical manifolds, Ann

R. Bott, Nondegenerate critical manifolds, Ann. of Mat h., 60 (1954), pp. 248–261

work page 1954
[12]

Carusotto and C

I. Carusotto and C. Ciuti, Quantum ﬂuids of light, Rev. M od. Phys., 85 (2013), pp. 299– 366

work page 2013
[13]

Cazenave, Semilinear Schr¨ odinger Equations, Cour ant Lect

T. Cazenave, Semilinear Schr¨ odinger Equations, Cour ant Lect. Notes Math., 10, Amer. Math. Soc., Providence, R.I., 2003

work page 2003
[14]

H. Chen, G. Dong, W. Liu, and Z. Xie, Second-order ﬂows fo r computing the ground states of rotating Bose-Einstein condensates, J. Comput. P hys., 475 (2023), article 111872

work page 2023
[15]

Z. Chen, J. Lu, Y. Lu, and X. Zhang, On the convergence of S obolev gradient ﬂow for the Gross-Pitaevskii eigenvalue problem, SIAM J. Numer. An al., 62 (2024), pp. 667–691

work page 2024
[16]

Chiofalo, S

M. Chiofalo, S. Succi, and M. Tosi, Ground state of trapp ed interacting Bose-Einstein condensates by an explicit imaginary-time algorithm, Phys . Rev. E, 62 (2000), pp. 7438– 7444

work page 2000
[17]

Danaila and P

I. Danaila and P. Kazemi, A new Sobolev gradient method f or direct minimization of the Gross-Pitaevskii energy with rotation, SIAM J. Sci. Comput ., 32 (2010), pp. 2447–2467

work page 2010
[18]

Danaila and B

I. Danaila and B. Protas, Computation of ground states o f the Gross-Pitaevskii func- tional via Riemannian optimization, SIAM J. Sci. Comput., 3 9 (2017), pp. B1102–B1129

work page 2017
[19]

K. B. Davis, M. Mewes, and M. R. Andrews, Bose-Einstein c ondensation in a gas of sodium atoms, Phys. Rev. Lett., 75 (1995), pp. 3969–3973

work page 1995
[20]

C. M. Dion and E. Canc´ es, Ground state of the time-indep endent Gross-Pitaevskii equation, Comput. Phys. Commun., 177 (2007), pp. 787–798

work page 2007
[21]

Dong and Y

L. Dong and Y. V. Kartashov, Rotating multidimensional quantum droplets, Phys. Rev. Lett., 126 (2021), article 244101

work page 2021
[22]

P. M. Feehan and M. Maridakis, /suppress Lojasiewicz-Simon gradient inequalities for analytic and Morse-Bott functions on Banach spaces, J. Reine Angew. Math ., 765 (2020), pp. 35–67

work page 2020
[23]

Feng and Q

Z. Feng and Q. Tang, On preconditioned Riemannian gradi ent methods for minimizing the Gross-Pitaevskii energy functional: algorithms, glob al convergence and optimal local convergence rate, arXiv:2510.13516. 25

work page arXiv
[24]

J. J. Garc ´ ıa. Ripoll and V. M. P´ erez-Garc ´ ıa, Optimizing Schr¨ odinger functionals using Sobolev gradients: Applications to quantum mechanics and n onlinear optics, SIAM J. Sci. Comput., 23 (2001), pp. 1316–1334

work page 2001
[25]

Henning and D

P. Henning and D. Peterseim, Sobolev gradient ﬂow for th e Gross-Pitaevskii eigenvalue problem: global convergence and computational eﬃciency, S IAM J. Numer. Anal., 58 (2020), pp. 1744–1772

work page 2020
[26]

Henning, The dependency of spectral gaps on the conve rgence of the inverse iteration for a nonlinear eigenvector problem, Math

P. Henning, The dependency of spectral gaps on the conve rgence of the inverse iteration for a nonlinear eigenvector problem, Math. Mod. Meth. Appl. S., 33 (2023), pp. 1517– 1544

work page 2023
[27]

Henning and M

P. Henning and M. Yadav, Convergence of a Riemannian gra dient method for the Gross- Pitaevskii energy functional in a rotating frame, ESAIM Mat h. Model. Numer. Anal., 59 (2025), pp 1145–1175

work page 2025
[28]

Henning and M

P. Henning and M. Yadav, On discrete ground states of rot ating Bose–Einstein conden- sates, Math. Comp., 94 (2025), pp 1–32

work page 2025
[29]

Henning and E

P. Henning and E. Jarlebring, The Gross-Pitaevskii equ ation and eigenvector nonlinear- ities: numerical methods and algorithms, SIAM Rev., 67 (202 5), pp. 256–317

work page
[30]

Hermann, T

M. Hermann, T. Stykel, and M. Yadav, Qualitative and Qua ntitative Analysis of Rieman- nian Optimization Methods for Ground States of Rotating Mul ticomponent Bose-Einstein Condensates, arXiv:2512.05939

work page arXiv
[31]

W. Hu, R. Barkana, and A. Gruzinov, Fuzzy cold dark matte r: the wave properties of ultralight particles, Phys. Rev. Lett., 85 (2000), pp. 1158 –1161

work page 2000
[32]

Jarlebring, S

E. Jarlebring, S. Kvaal, and W. Michiels, An inverse ite ration method for eigenvalue problems with eigenvector nonlinearities, SIAM J. Sci. Com put., 36 (2014), pp. A1978– A2001

work page 2014
[33]

Kazemi and M

P. Kazemi and M. Eckart, Minimizing the Gross-Pitaevsk ii energy functional with the Sobolev gradient-analytical and numerical results, In t. J. Comput. Meth., 7 (2010), pp. 453–475

work page 2010
[34]

Klaers, J

J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Bose-E instein condensation of photons in an optical microcavity, Nat., 468 (2010), pp. 545-548

work page 2010
[35]

E. H. Lieb and R. Seiringer, Derivation of the Gross-Pit aevskii equation for rotating Bose gases, Commun. Math. Phys., 264 (2006), pp. 505–537

work page 2006
[36]

Lang, Real and Functional Analysis, 3rd ed., Graduat e Texts in Mathematics, vol

S. Lang, Real and Functional Analysis, 3rd ed., Graduat e Texts in Mathematics, vol. 142, Springer-Verlag, New York, 1993

work page 1993
[37]

Liu and Y

W. Liu and Y. Cai, Normalized gradient ﬂow with Lagrange multiplier for computing ground states of Bose-Einstein condensates, SIAM J. Sci. Co mput., 43 (2021), pp. B219– B242

work page 2021
[38]

B. T. Polyak, Introduction to Optimization. New York: O ptimization Software, Inc., 1987

work page 1987
[39]

Rebjock and N

Q. Rebjock and N. Boumal, Fast convergence to non-isola ted minima: four equivalent conditions for C 2 functions, Math. Program., 213 (2025), pp 151–199. 26

work page 2025
[40]

Rupp, On the /suppress Lojasiewicz-Simon gradient inequality on submanifolds, J

F. Rupp, On the /suppress Lojasiewicz-Simon gradient inequality on submanifolds, J. Funct. Anal., 279 (2020), article 108708

work page 2020
[41]

Shamriz, Z

E. Shamriz, Z. Chen, and B. A. Malomed, Suppression of th e quasi-two-dimensional quantum collapse in the attraction ﬁeld by the Lee-Huang-Ya ng eﬀect, Phys. Rev. A., 101 (2020), article 063628

work page 2020
[42]

M. N. Tengstrand, P. St¨ urmer, E. ¨O. Karabulut, and S. M. Reimann, Rotating binary Bose-Einstein condensates and vortex clusters in quantum d roplets, Phys. Rev. Lett., 123 (2019), article 160405

work page 2019
[43]

Y. Wu, C. Liu, and Y. Cai, Normalized ﬂows based on Sobole v gradients for computing ground states of spinor Bose-Einstein condensates, J. Comp ut. Phys., 538 (2025), article 114153

work page 2025
[44]

X. Wu, Z. Wen, and W. Bao, A regularized newton method for computing ground states of Bose-Einstein condensates, J. Sci. Comput., 73 (2017), p p. 303–329

work page 2017
[45]

Zhang and F

T. Zhang and F. Xue, A new preconditioned nonlinear conj ugate gradient method in real arithmetic for computing the ground states of rotation al Bose-Einstein condensate, SIAM J. Sci. Comput., 46 (2024), pp. A1764–A1792

work page 2024
[46]

Zhang, Exponential convergence of Sobolev gradient descent for a class of nonlinear eigenproblems

Z. Zhang, Exponential convergence of Sobolev gradient descent for a class of nonlinear eigenproblems. Commun. Math. Sci., 20 (2022), pp. 377–403

work page 2022
[47]

Zhang, P

C. Zhang, P. Henning, M. Yadav, and W. Chen, Convergence analysis of Sobolev Gradient ﬂows for the rotating Gross-Pitaevskii energy functional, arXiv:2510.15604

work page arXiv
[48]

Zhuang and J

Q. Zhuang and J. Shen, Eﬃcient SA V approach for imaginar y time gradient ﬂows with applications to one- and multi-component Bose-Einstein Co ndensates, J. Comput. Phys., 396 (2019), pp. 72–88. A Proof of Lzφ ∈ H 1 0 (D) In this appendix, for the local minimizer φg, we prove that Lzφg ∈ H 1 0 (D). To this end, we establish the following lemma. Lemma A.1. ...

work page 2019

[1] [1]

Y. Ai, P. Henning, M. Yadav, and S. Yuan, Riemannian conju gate Sobolev gradients and their application to compute ground states of BECs, J. Co mput. Appl. Math., 473 (2026), article 116866

work page 2026

[2] [2]

Altmann, P

R. Altmann, P. Henning, and D. Peterseim, The J-method for the Gross-Pitaevskii eigenvalue problem, Numer. Math., 148 (2021), pp. 575–610

work page 2021

[3] [3]

Altmann, M

R. Altmann, M. Hermann, D. Peterseim, and T. Stykel, Riem annian optimisation meth- ods for ground states of multicomponent Bose-Einstein cond ensates, arXiv:2411.09617

work page arXiv

[4] [4]

M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman, and E. A. Cornell, Observation of Bose-Einstein condensation in a dilute atom ic vapor, Sci., 269 (1995), pp. 198–201

work page 1995

[5] [5]

Antoine and R

X. Antoine and R. Duboscq, Robust and eﬃcient preconditi oned Krylov spectral solvers for computing the ground states of fast rotating and strongl y interacting Bose-Einstein condensates, J. Comput. Phys., 258 (2014), pp. 509–523

work page 2014

[6] [6]

Antoine, A

X. Antoine, A. Levitt, and Q. Tang, Eﬃcient spectral comp utation of the stationary states of rotating Bose-Einstein condensates by the precon ditioned nonlinear conjugate gradient method, J. Comput. Phys., 343 (2017), pp. 92–109. 24

work page 2017

[7] [7]

Bao and Q

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

work page 2004

[8] [8]

W. Bao, I. Chern, and F. Lim, Eﬃcient and spectrally accur ate numerical methods for computing ground and ﬁrst excited states in Bose-Einste in condensates, J. Comput. Phys., 219 (2006), pp. 836–854

work page 2006

[9] [9]

Bao and Y

W. Bao and Y. Cai, Mathematical theory and numerical meth ods for Bose-Einstein condensation, Kinet. Relat. Models, 6 (2013), pp. 1–135

work page 2013

[10] [10]

C. F. Barenghi, L. Skrbek, and K. R. Sreenivasan, Introd uction to quantum turbulence, PNAS, 111 (2014), pp. 4647–4652

work page 2014

[11] [11]

Bott, Nondegenerate critical manifolds, Ann

R. Bott, Nondegenerate critical manifolds, Ann. of Mat h., 60 (1954), pp. 248–261

work page 1954

[12] [12]

Carusotto and C

I. Carusotto and C. Ciuti, Quantum ﬂuids of light, Rev. M od. Phys., 85 (2013), pp. 299– 366

work page 2013

[13] [13]

Cazenave, Semilinear Schr¨ odinger Equations, Cour ant Lect

T. Cazenave, Semilinear Schr¨ odinger Equations, Cour ant Lect. Notes Math., 10, Amer. Math. Soc., Providence, R.I., 2003

work page 2003

[14] [14]

H. Chen, G. Dong, W. Liu, and Z. Xie, Second-order ﬂows fo r computing the ground states of rotating Bose-Einstein condensates, J. Comput. P hys., 475 (2023), article 111872

work page 2023

[15] [15]

Z. Chen, J. Lu, Y. Lu, and X. Zhang, On the convergence of S obolev gradient ﬂow for the Gross-Pitaevskii eigenvalue problem, SIAM J. Numer. An al., 62 (2024), pp. 667–691

work page 2024

[16] [16]

Chiofalo, S

M. Chiofalo, S. Succi, and M. Tosi, Ground state of trapp ed interacting Bose-Einstein condensates by an explicit imaginary-time algorithm, Phys . Rev. E, 62 (2000), pp. 7438– 7444

work page 2000

[17] [17]

Danaila and P

I. Danaila and P. Kazemi, A new Sobolev gradient method f or direct minimization of the Gross-Pitaevskii energy with rotation, SIAM J. Sci. Comput ., 32 (2010), pp. 2447–2467

work page 2010

[18] [18]

Danaila and B

I. Danaila and B. Protas, Computation of ground states o f the Gross-Pitaevskii func- tional via Riemannian optimization, SIAM J. Sci. Comput., 3 9 (2017), pp. B1102–B1129

work page 2017

[19] [19]

K. B. Davis, M. Mewes, and M. R. Andrews, Bose-Einstein c ondensation in a gas of sodium atoms, Phys. Rev. Lett., 75 (1995), pp. 3969–3973

work page 1995

[20] [20]

C. M. Dion and E. Canc´ es, Ground state of the time-indep endent Gross-Pitaevskii equation, Comput. Phys. Commun., 177 (2007), pp. 787–798

work page 2007

[21] [21]

Dong and Y

L. Dong and Y. V. Kartashov, Rotating multidimensional quantum droplets, Phys. Rev. Lett., 126 (2021), article 244101

work page 2021

[22] [22]

P. M. Feehan and M. Maridakis, /suppress Lojasiewicz-Simon gradient inequalities for analytic and Morse-Bott functions on Banach spaces, J. Reine Angew. Math ., 765 (2020), pp. 35–67

work page 2020

[23] [23]

Feng and Q

Z. Feng and Q. Tang, On preconditioned Riemannian gradi ent methods for minimizing the Gross-Pitaevskii energy functional: algorithms, glob al convergence and optimal local convergence rate, arXiv:2510.13516. 25

work page arXiv

[24] [24]

J. J. Garc ´ ıa. Ripoll and V. M. P´ erez-Garc ´ ıa, Optimizing Schr¨ odinger functionals using Sobolev gradients: Applications to quantum mechanics and n onlinear optics, SIAM J. Sci. Comput., 23 (2001), pp. 1316–1334

work page 2001

[25] [25]

Henning and D

P. Henning and D. Peterseim, Sobolev gradient ﬂow for th e Gross-Pitaevskii eigenvalue problem: global convergence and computational eﬃciency, S IAM J. Numer. Anal., 58 (2020), pp. 1744–1772

work page 2020

[26] [26]

Henning, The dependency of spectral gaps on the conve rgence of the inverse iteration for a nonlinear eigenvector problem, Math

P. Henning, The dependency of spectral gaps on the conve rgence of the inverse iteration for a nonlinear eigenvector problem, Math. Mod. Meth. Appl. S., 33 (2023), pp. 1517– 1544

work page 2023

[27] [27]

Henning and M

P. Henning and M. Yadav, Convergence of a Riemannian gra dient method for the Gross- Pitaevskii energy functional in a rotating frame, ESAIM Mat h. Model. Numer. Anal., 59 (2025), pp 1145–1175

work page 2025

[28] [28]

Henning and M

P. Henning and M. Yadav, On discrete ground states of rot ating Bose–Einstein conden- sates, Math. Comp., 94 (2025), pp 1–32

work page 2025

[29] [29]

Henning and E

P. Henning and E. Jarlebring, The Gross-Pitaevskii equ ation and eigenvector nonlinear- ities: numerical methods and algorithms, SIAM Rev., 67 (202 5), pp. 256–317

work page

[30] [30]

Hermann, T

M. Hermann, T. Stykel, and M. Yadav, Qualitative and Qua ntitative Analysis of Rieman- nian Optimization Methods for Ground States of Rotating Mul ticomponent Bose-Einstein Condensates, arXiv:2512.05939

work page arXiv

[31] [31]

W. Hu, R. Barkana, and A. Gruzinov, Fuzzy cold dark matte r: the wave properties of ultralight particles, Phys. Rev. Lett., 85 (2000), pp. 1158 –1161

work page 2000

[32] [32]

Jarlebring, S

E. Jarlebring, S. Kvaal, and W. Michiels, An inverse ite ration method for eigenvalue problems with eigenvector nonlinearities, SIAM J. Sci. Com put., 36 (2014), pp. A1978– A2001

work page 2014

[33] [33]

Kazemi and M

P. Kazemi and M. Eckart, Minimizing the Gross-Pitaevsk ii energy functional with the Sobolev gradient-analytical and numerical results, In t. J. Comput. Meth., 7 (2010), pp. 453–475

work page 2010

[34] [34]

Klaers, J

J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Bose-E instein condensation of photons in an optical microcavity, Nat., 468 (2010), pp. 545-548

work page 2010

[35] [35]

E. H. Lieb and R. Seiringer, Derivation of the Gross-Pit aevskii equation for rotating Bose gases, Commun. Math. Phys., 264 (2006), pp. 505–537

work page 2006

[36] [36]

Lang, Real and Functional Analysis, 3rd ed., Graduat e Texts in Mathematics, vol

S. Lang, Real and Functional Analysis, 3rd ed., Graduat e Texts in Mathematics, vol. 142, Springer-Verlag, New York, 1993

work page 1993

[37] [37]

Liu and Y

W. Liu and Y. Cai, Normalized gradient ﬂow with Lagrange multiplier for computing ground states of Bose-Einstein condensates, SIAM J. Sci. Co mput., 43 (2021), pp. B219– B242

work page 2021

[38] [38]

B. T. Polyak, Introduction to Optimization. New York: O ptimization Software, Inc., 1987

work page 1987

[39] [39]

Rebjock and N

Q. Rebjock and N. Boumal, Fast convergence to non-isola ted minima: four equivalent conditions for C 2 functions, Math. Program., 213 (2025), pp 151–199. 26

work page 2025

[40] [40]

Rupp, On the /suppress Lojasiewicz-Simon gradient inequality on submanifolds, J

F. Rupp, On the /suppress Lojasiewicz-Simon gradient inequality on submanifolds, J. Funct. Anal., 279 (2020), article 108708

work page 2020

[41] [41]

Shamriz, Z

E. Shamriz, Z. Chen, and B. A. Malomed, Suppression of th e quasi-two-dimensional quantum collapse in the attraction ﬁeld by the Lee-Huang-Ya ng eﬀect, Phys. Rev. A., 101 (2020), article 063628

work page 2020

[42] [42]

M. N. Tengstrand, P. St¨ urmer, E. ¨O. Karabulut, and S. M. Reimann, Rotating binary Bose-Einstein condensates and vortex clusters in quantum d roplets, Phys. Rev. Lett., 123 (2019), article 160405

work page 2019

[43] [43]

Y. Wu, C. Liu, and Y. Cai, Normalized ﬂows based on Sobole v gradients for computing ground states of spinor Bose-Einstein condensates, J. Comp ut. Phys., 538 (2025), article 114153

work page 2025

[44] [44]

X. Wu, Z. Wen, and W. Bao, A regularized newton method for computing ground states of Bose-Einstein condensates, J. Sci. Comput., 73 (2017), p p. 303–329

work page 2017

[45] [45]

Zhang and F

T. Zhang and F. Xue, A new preconditioned nonlinear conj ugate gradient method in real arithmetic for computing the ground states of rotation al Bose-Einstein condensate, SIAM J. Sci. Comput., 46 (2024), pp. A1764–A1792

work page 2024

[46] [46]

Zhang, Exponential convergence of Sobolev gradient descent for a class of nonlinear eigenproblems

Z. Zhang, Exponential convergence of Sobolev gradient descent for a class of nonlinear eigenproblems. Commun. Math. Sci., 20 (2022), pp. 377–403

work page 2022

[47] [47]

Zhang, P

C. Zhang, P. Henning, M. Yadav, and W. Chen, Convergence analysis of Sobolev Gradient ﬂows for the rotating Gross-Pitaevskii energy functional, arXiv:2510.15604

work page arXiv

[48] [48]

Zhuang and J

Q. Zhuang and J. Shen, Eﬃcient SA V approach for imaginar y time gradient ﬂows with applications to one- and multi-component Bose-Einstein Co ndensates, J. Comput. Phys., 396 (2019), pp. 72–88. A Proof of Lzφ ∈ H 1 0 (D) In this appendix, for the local minimizer φg, we prove that Lzφg ∈ H 1 0 (D). To this end, we establish the following lemma. Lemma A.1. ...

work page 2019