Stable solvers and block elimination for bordered systems
SIAM Journal on Matrix Analysis and Applications
Algorithm 652: HOMPACK: a suite of codes for globally convergent homotopy algorithms
ACM Transactions on Mathematical Software (TOMS)
Large-Scale Continuation and Numerical Bifurcation for Partial Differential Equations
SIAM Journal on Numerical Analysis
Numerical methods for bifurcations of dynamical equilibria
Numerical methods for bifurcations of dynamical equilibria
Ground-state solution of Bose--Einstein condensate by directly minimizing the energy functional
Journal of Computational Physics
Numerical solution of the Gross--Pitaevskii equation for Bose--Einstein condensation
Journal of Computational Physics
SIAM Journal on Numerical Analysis
Computing the Ground State Solution of Bose--Einstein Condensates by a Normalized Gradient Flow
SIAM Journal on Scientific Computing
Gauss-Seidel-type methods for energy states of a multi-component Bose-Einstein condensate
Journal of Computational Physics
Gauss-Seidel-type methods for energy states of a multi-component Bose-Einstein condensate
Journal of Computational Physics
Journal of Computational Physics
Applied Numerical Mathematics
Journal of Computational Physics
Hi-index | 31.46 |
We develop a continuation block successive over-relaxation (BSOR)-Lanczos-Galerkin method for the computation of positive bound states of time-independent, coupled Gross-Pitaevskii equations (CGPEs) which describe a multi-component Bose-Einstein condensate (BEC). A discretization of the CGPEs leads to a nonlinear algebraic eigenvalue problem (NAEP). The solution curve with respect to some parameter of the NAEP is then followed by the proposed method. For a single-component BEC, we prove that there exists a unique global minimizer (the ground state) which is represented by an ordinary differential equation with the initial value. For a multi-component BEC, we prove that m identical ground/bound states will bifurcate into m different ground/bound states at a finite repulsive inter-component scattering length. Numerical results show that various positive bound states of a two/three-component BEC are solved efficiently and reliably by the continuation BSOR-Lanczos-Galerkin method.