Efficient continuation methods for spin-1 Bose-Einstein condensates in a magnetic field

We study linear stability analysis for spin-1 Bose-Einstein condensates (BEC). We show that all bounded solutions of this physical system are neutrally stable. In particular, all steady-state solutions of the physical system, and the associated discrete steady-state solutions are neutrally stable. Next, we consider the physical system without the affect of magnetic field. By exploiting the physical properties of both ferromagnetic and antiferromagnetic cases, we develop efficient multi-level pseudo-arclength continuation algorithms combined with a spectral collocation method for these two cases, respectively. When the magnetic field is imposed on the physical system, an additional multi-level continuation algorithm is described for the ferromagnetic case. Extensive numerical results for spin-1 BEC in a magnetic field, and in optical lattices are reported.