Abstract:This paper studies the numerical approximation of the ground state of the Gross-Pitaevskii (GP) eigenvalue problem with a fully discretized Sobolev gradient flow induced by the $H^1$ norm. For the spatial discretization, we consider the finite element method with quadrature using $P^k$ basis on a simplicial mesh and $Q^k$ basis on a rectangular mesh. We prove the global convergence to a critical point of the discrete GP energy, and establish a local exponential convergence to the ground state under the assumption that the linearized discrete Schrödinger operator has a positive spectral gap. We also show that for the $P^1$ finite element discretization with quadrature on an unstructured shape regular simplicial mesh, the eigengap satisfies a mesh-independent lower bound, which implies a mesh-independent local convergence rate for the proposed discrete gradient flow. Numerical experiments with discretization by high order $Q^k$ spectral element methods in two and three dimensions are provided to validate the efficiency of the proposed method.

What problem does this paper attempt to address?

The core problem that this paper attempts to solve is: **How to effectively solve the ground - state of the Gross - Pitaevskii (GP) eigenvalue problem by means of the fully - discretized Sobolev gradient flow method, and prove the global convergence and local exponential convergence of its numerical approximation.** Specifically, the paper focuses on solving the minimization problem of the GP energy functional given the external potential field $ V(x) $ and the scattering length $ a $. The GP energy functional is defined as: \[ E_{GP}(\phi)=\int_{\mathbb{R}^3}\left(|\nabla\phi(x)|^{2}+V(x)|\phi(x)|^{2}+4\pi a|\phi(x)|^{4}\right)dx \] In order to find the ground - state, it is necessary to minimize the above - mentioned energy functional under the normalization condition. For large - scale computations, calculations are usually carried out on a finite domain, and it is assumed that the wave function has a compact support, so it can be approximated as: \[ E_{L}(\phi)=\int_{\Omega}\left(\frac{1}{2}|\nabla\phi(x)|^{2}+\frac{1}{2}V(x)|\phi(x)|^{2}+\frac{\beta}{4}|\phi(x)|^{4}\right)dx \] where $ \Omega =[-L, L]^{d} $ and satisfies the normalization condition $ \int_{\Omega}|\phi(x)|^{2}dx = 1 $. ### Main contributions 1. **Proposed a fully - discretized Sobolev gradient descent method**: - This method performs Riemannian gradient descent on the sphere based on the Riemannian metric induced by the modified $ H^{1} $-norm. 2. **Proved global convergence and local exponential convergence**: - It is proved that the fully - discretized Sobolev gradient descent method globally converges to a critical point of the discrete GP energy under the modified $ H^{1} $-metric. - Under the assumption that the discrete linearized Schrödinger operator has a positive spectral gap, the result of local exponential convergence to the ground - state is established. 3. **Numerical experiment verification**: - Numerical experiments in two - dimensional and three - dimensional spaces using the high - order $ Q_{k} $ spectral element method are provided to verify the effectiveness and efficiency of the proposed method. ### Related work - Self - consistent field iteration (SCF) is a commonly used iterative technique for solving nonlinear eigenvalue problems, but it may diverge unless a good initial guess is provided. - Projected Sobolev gradient flow is a method based on an optimization perspective and has been applied in some works, but its convergence analysis is not yet complete. ### Main theoretical results of the paper - Proved the global convergence and local exponential convergence rate. - For the $ P_{1} $ finite - element method, the positive discrete eigenvalue gap of the discrete eigenvalue problem is proved, which guarantees a mesh - independent local convergence rate. Through these studies, the paper provides an effective and theoretically reliable numerical method for solving the GP eigenvalue problem.

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

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

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

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

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

An extended Fourier pseudospectral method for the Gross-Pitaevskii equation with low regularity potential

Finite volumes for the Gross-Pitaevskii equation

The dependency of spectral gaps on the convergence of the inverse iteration for a nonlinear eigenvector problem

A Multigrid Discretization of Discontinuous Galerkin Method for the Stokes Eigenvalue Problem

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

Numerical simulation of the Gross-Pitaevskii equation via vortex tracking

Localized Orthogonal Decomposition Methods vs. Classical FEM for the Gross-Pitaevskii Equation

A BGK-Based Discontinuous Galerkin Method for the Navier-Stokes Equations on Arbitrary Grids

Properties of ground states of attractive Gross–Pitaevskii equations with multi-well potentials

A Discontinuous Galerkin Method Based On A Bgk Scheme For The Navier-Stokes Equations On Arbitrary Grids

A high-order local discontinuous Galerkin method for the $p$-Laplace equation

Continuity of a spatial gradient of a weak solution to a very singular parabolic equation involving the one-Laplacian

Generalized weak Galerkin methods for Stokes equations

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