A Regularized Newton Method for Computing Ground States of Bose-Einstein condensates

Abstract

In this paper, we propose a regularized Newton method for computing ground states of Bose-Einstein condensates (BECs), which can be formulated as an energy minimization problem with a spherical constraint. The energy functional and constraint are discretized by either the finite difference, or sine or Fourier pseudospectral discretization schemes and thus the original infinite dimensional nonconvex minimization problem is approximated by a finite dimensional constrained nonconvex minimization problem. Then an initial solution is first constructed by using a feasible gradient type method, which is an explicit scheme and maintains the spherical constraint automatically. To accelerate the convergence of the gradient type method, we approximate the energy functional by its second-order Taylor expansion with a regularized term at each Newton iteration and adopt a cascadic multigrid technique for selecting initial data. It leads to a standard trust-region subproblem and we solve it again by the feasible gradient type method. The convergence of the regularized Newton method is established by adjusting the regularization parameter as the standard trust-region strategy. Extensive numerical experiments on challenging examples, including a BEC in three dimensions with an optical lattice potential and rotating BECs in two dimensions with rapid rotation and strongly repulsive interaction, show that our method is efficient, accurate and robust.

arXiv comment25 pages, 6 figures

References in corpus (12)

• Theory of Bose-Einstein condensation in trapped gases
• Vortices in a Bose-Einstein Condensate
• Mathematical theory and numerical methods for Bose-Einstein condensation
• Vortex Nucleation in a Stirred Bose-Einstein Condensate
• Bosons in a Trap: A Rigorous Derivation of the Gross-Pitaevskii Energy Functional
• Derivation of the Gross-Pitaevskii Equation for Rotating Bose Gases
• Ground state solution of Bose-Einstein condensate by directly minimizing the energy functional
• Bose-Einstein condensate: critical velocities and energy diagrams in the Thomas-Fermi regime
• Efficient numerical methods for computing ground states and dynamics of dipolar Bose-Einstein condensates
• Numerical Solution of the Two-Dimensional Gross-Pitaevskii Equation for Trapped Interacting Atoms
• A numerical approach to the ground and excited states of a Bose-Einstein condensed gas confined in a completely anisotropic trap
• Three-dimensional vortex configurations in a rotating Bose Einstein condensate

Cited by in corpus (6)

• Computation of Ground States of the Gross-Pitaevskii Functional via Riemannian Optimization
• Constrained high-index saddle dynamics for the solution landscape with equality constraints
• Second-order flows for computing the ground states of rotating Bose-Einstein condensates
• Computing the least action ground state of the nonlinear Schrödinger equation by a normalized gradient flow
• Computing the action ground state for the rotating nonlinear Schrödinger equation
• A normalized gradient flow method with attractive-repulsive splitting for computing ground states of Bose-Einstein condensates with higher-order interaction