Abstract:This paper investigates the numerical approximation of ground states of rotating Bose-Einstein condensates. This problem requires the minimization of the Gross-Pitaevskii energy $E$ on a Hilbert manifold $\mathbb{S}$. To find a corresponding minimizer $u$, we use a generalized Riemannian gradient method that is based on the concept of Sobolev gradients in combination with an adaptively changing metric on the manifold. By a suitable choice of the metric, global energy dissipation for the arising gradient method can be proved. The energy dissipation property in turn implies global convergence to the density $|u|^2$ of a critical point $u$ of $E$ on $\mathbb{S}$. Furthermore, we present a precise characterization of the local convergence rates in a neighborhood of each ground state $u$ and how these rates depend on the first spectral gap of $E^{\prime\prime}(u)$ restricted to the $L^2$-orthogonal complement of $u$. With this we establish the first convergence results for a Riemannian gradient method to minimize the Gross-Pitaevskii energy functional in a rotating frame. At the same, we refine previous results obtained in the case without rotation. The major complication in our new analysis is the missing isolation of minimizers, which are at most unique up to complex phase shifts. For that, we introduce an auxiliary iteration in the tangent space $T_{\mathrm{i} u} \mathbb{S}$ and apply the Ostrowski theorem to characterize the asymptotic convergence rates through a weighted eigenvalue problem. Afterwards, we link the auxiliary iteration to the original Riemannian gradient method and bound the spectrum of the weighted eigenvalue problem to obtain quantitative convergence rates. Our findings are validated in numerical experiments.

Convergence of a Riemannian gradient method for the Gross-Pitaevskii energy functional in a rotating frame

On discrete ground states of rotating Bose-Einstein condensates

On discrete ground states of rotating Bose–Einstein condensates

Riemannian conjugate Sobolev gradients and their application to compute ground states of BECs

On the Convergence of Sobolev Gradient Flow for the Gross–Pitaevskii Eigenvalue Problem

Convergence Analysis of Gradient Algorithms on Riemannian Manifolds Without Curvature Constraints and Application to Riemannian Mass

Riemannian optimisation methods for ground states of multicomponent Bose-Einstein condensates

An efficient and conservative compact finite difference scheme for the coupled Gross-Pitaevskii equations describing spin-1 Bose-Einstein condensate.

Structure-preserving finite element methods for computing dynamics of rotating Bose-Einstein condensate

Fully discretized Sobolev gradient flow for the Gross-Pitaevskii eigenvalue problem

Uniform $L^\infty$-bounds for energy-conserving higher-order time integrators for the Gross-Pitaevskii equation with rotation

Randomized Gradient Descents on Riemannian Manifolds: Almost Sure Convergence to Global Minima in and beyond Quantum Optimization

A conservative splitting high-order finite difference method for coupled Gross–Pitaevskii equations in 2D

Riemannian Newton methods for energy minimization problems of Kohn-Sham type

A Regularized Newton Method for Computing Ground States of Bose–Einstein Condensates

Mixed finite elements for the Gross-Pitaevskii eigenvalue problem: a priori error analysis and guaranteed lower energy bound

Super-localised wave function approximation of Bose-Einstein condensates

Mixed finite elements for the Gross–Pitaevskii eigenvalue problem: a priori error analysis and guaranteed lower energy bound

Projection Gradient Method for Energy Functional Minimization with A Constraint and Its Application to Computing the Ground State of Spin-Orbit-Coupled Bose-Einstein Condensates