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.

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

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

Positivity preserving finite element method for the Gross-Pitaevskii ground state: discrete uniqueness and global convergence

Mixed methods and lower eigenvalue bounds

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

The Lower/Upper Bound Property of Approximate Eigenvalues by Nonconforming Finite Element Methods for Elliptic Operators

The a posteriori error estimates of the Ciarlet-Raviart mixed finite element method for the biharmonic eigenvalue problem

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

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

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

Energy-norm error estimates for finite element discretization of parabolic problems

Error analysis of mixed finite element methods for nonlinear parabolic equations

The Lower Bounds for Eigenvalues of Elliptic Operators --By Nonconforming Finite Element Methods

H 1 -conforming virtual element method for the Laplacian eigenvalue problem in mixed form

A new optimal error analysis of a mixed finite element method for advection–diffusion–reaction Brinkman flow

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

Error estimates for a mixed finite element method for the Maxwell's transmission eigenvalue problem

An equilibrated estimator for mixed finite element discretizations of the curl-curl problem

Asymptotically Exact a Posteriori Error Estimates of Eigenvalues by the Crouzeix-Raviart Element and Enriched Crouzeix-Raviart Element

Raviart-Thomas finite elements of Petrov-Galerkin type

Unconditionally optimal error estimates of a linearized weak Galerkin finite element method for semilinear parabolic equations