A continuation BSOR-Lanczos--Galerkin method for positive bound states of a multi-component Bose-Einstein condensate (Q2568062)

The authors develop a continuation block successive over-relaxation Lanczos-Galerkin method for the computation of positive bound states of a multi-component Bose-Einstein condensate (BEC). The bifurcation diagram of positive eigenvalues and eigenvectors of a nonlinear algebraic eigenvalue problem, obtained as result of a discretization of a coupled Gross-Pitaevskii equation and its associated energy functional, is traced by the proposed continuation method, different from other methods for the computation of the ground states of a multicomponent BEC and applicable for all possible bound states. They prove the occurrence of a phase separation of \(m\) ground/bound states at a finite value of the repulsive scattering length. For a two-component BEC it is shown that two identical ground/bound states bifurcate into different \(\Pi_{\theta}\)-symmetry ones.