Computing ground states of Bose-Einstein condensates with higher order interaction via a regularized density function formulation

DOI10.1137/19M1240393MaRDI QIDQ5204005zbMATH OpenOpenAlexWikidataFDO

Authors Weizhu Bao, Xinran Ruan

Publication date 9 December 2019

Published in SIAM Journal on Scientific Computing (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1908.09096

zbMATH Keywords

accelerated projected gradient methodBose-Einstein condensatedensity function formulationground statehigher-order interaction

Mathematics Subject Classification ID

Smoothness and regularity of solutions to PDEs (35B65) NLS equations (nonlinear Schrödinger equations) (35Q55) Numerical optimization and variational techniques (65K10) Finite difference methods for boundary value problems involving PDEs (65N06) Stability and convergence of numerical methods for boundary value problems involving PDEs (65N12) Numerical methods for eigenvalue problems for boundary value problems involving PDEs (65N25) Computational density functional analysis in statistical mechanics (82M36) Convex programming (90C25) Nonlinear programming (90C30)

Abstract: We propose and analyze a new numerical method for computing the ground state of the modified Gross-Pitaevskii equation for modeling the Bose-Einstein condensate with a higher order interaction by adapting the density function formulation and the accelerated projected gradient method. By reformulating the energy functional $E(phi)$ with $phi$, the wave function, in terms of the density $ho=|phi{|}^2$, the original non-convex minimization problem for defining the ground state is then reformulated to a convex minimization problem. In order to overcome the semi-smoothness of the function $sqrtho$ in the kinetic energy part, a regularization is introduced with a small parameter $0<varepsilonll1$. Convergence of the regularization is established when $varepsilono0$. The regularized convex optimization problem is discretized by the second order finite difference method. The convergence rates in terms of the density and energy of the discretization are established. The accelerated projected gradient method is adapted for solving the discretized optimization problem. Numerical results are reported to demonstrate the efficiency and accuracy of the proposed numerical method. Our results show that the proposed method is much more efficient than the existing methods in the literature, especially in the strong interaction regime.

Recommendations

• A regularized Newton method for computing ground states of Bose-Einstein condensates
• Ground-state solution of Bose--Einstein condensate by directly minimizing the energy functional
• A hybrid discontinuous Galerkin method for computing the ground state solution of Bose-Einstein condensates
• A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction
• Normalized gradient flow with Lagrange multiplier for computing ground states of Bose-Einstein condensates

Cites work

• A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems
• A handbook of -convergence
• A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction
• A quasi-local Gross-Pitaevskii equation for attractive Bose-Einstein condensates
• A regularized Newton method for computing ground states of Bose-Einstein condensates
• An algorithm for a singly constrained class of quadratic programs subject upper and lower bounds
• An O(n) algorithm for quadratic knapsack problems
• Computation of ground states of the Gross-Pitaevskii functional via Riemannian optimization
• Computing the Ground State Solution of Bose--Einstein Condensates by a Normalized Gradient Flow
• Efficient and spectrally accurate numerical methods for computing ground and first excited states in Bose-Einstein condensates
• Efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates
• Efficient spectral computation of the stationary states of rotating Bose-Einstein condensates by preconditioned nonlinear conjugate gradient methods
• Gauss-Seidel-type methods for energy states of a multi-component Bose-Einstein condensate
• Ground States and Dynamics of Multicomponent Bose--Einstein Condensates
• Ground states of Bose-Einstein condensates with higher order interaction
• Ground-state solution of Bose--Einstein condensate by directly minimizing the energy functional
• scientific article; zbMATH DE number 3850830 (Why is no real title available?)
• scientific article; zbMATH DE number 1151368 (Why is no real title available?)
• scientific article; zbMATH DE number 6803232 (Why is no real title available?)
• scientific article; zbMATH DE number 2214428 (Why is no real title available?)
• Mathematical theory and numerical methods for Bose-Einstein condensation
• Robust and efficient preconditioned Krylov spectral solvers for computing the ground states of fast rotating and strongly interacting Bose-Einstein condensates
• Templates for convex cone problems with applications to sparse signal recovery
• The ground state energy and density of interacting bosons in a trap
• The mathematics of the Bose gas and its condensation.
• The Thomas-Fermi theory of atoms, molecules and solids
• The Thomas-Fermi-von Weizsäcker theory of atoms and molecules

Cited in

(12)

• A regularized Newton method for computing ground states of Bose-Einstein condensates
• Ground-state solution of Bose--Einstein condensate by directly minimizing the energy functional
• Ground states of Bose-Einstein condensates with higher order interaction
• Computing the least action ground state of the nonlinear Schrödinger equation by a normalized gradient flow
• Computable error estimates for ground state solution of Bose-Einstein condensates
• Error estimates of local energy regularization for the logarithmic Schrödinger equation
• Newton--Noda Iteration for Computing the Ground States of Nonlinear Schrödinger Equations
• Semi-discretization and full-discretization with improved accuracy for charged-particle dynamics in a strong nonuniform magnetic field
• A new preconditioned nonlinear conjugate gradient method in real arithmetic for computing the ground states of rotational Bose-Einstein condensate
• Newton-based alternating methods for the ground state of a class of multicomponent Bose-Einstein condensates
• An infeasible projection-type algorithm for finding nonnegative ground state solutions of nonlinear Schrödinger equations
• A positive-conservative Fourier optimization method for ground states of Bose--Einstein condensates with higher-order interactions

Describes a project that uses

Uses Software

• TFOCS

This page was built for publication: Computing ground states of Bose-Einstein condensates with higher order interaction via a regularized density function formulation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5204005)