Newton-based alternating methods for the ground state of a class of multi-component Bose-Einstein condensates

Abstract

The computation of the ground states of special multi-component Bose-Einstein condensates (BECs) can be formulated as an energy functional minimization problem with spherical constraints. It leads to a nonconvex quartic-quadratic optimization problem after suitable discretizations. First, we generalize the Newton-based methods for single-component BECs to the alternating minimization scheme for multi-component BECs. Second, the global convergent alternating Newton-Noda iteration (ANNI) is proposed. In particular, we prove the positivity preserving property of ANNI under mild conditions. Finally, our analysis is applied to a class of more general "multi-block" optimization problems with spherical constraints. Numerical experiments are performed to evaluate the performance of proposed methods for different multi-component BECs, including pseudo spin-1/2, anti-ferromagnetic spin-1 and spin-2 BECs. These results support our theory and demonstrate the efficiency of our algorithms.