A numerical scheme for the ground state of rotating spin-1 Bose–Einstein condensates

We study the existence of nontrivial solution branches of three-coupled Gross–Pitaevskii equations (CGPEs), which are used as the mathematical model for rotating spin-1 Bose–Einstein condensates (BEC). The Lyapunov–Schmidt reduction is exploited to test the branching of nontrivial solution curves from the trivial one in some neighborhoods of bifurcation points. A multilevel continuation method is proposed for computing the ground state solution of rotating spin-1 BEC. By properly choosing the constraint conditions associated with the components of the parameter variable, the proposed algorithm can effectively compute the ground states of spin-1 87Rb\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{87}Rb$$\end{document} and 23Na\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{23}Na$$\end{document} under rapid rotation. Extensive numerical results demonstrate the efficiency of the proposed algorithm. In particular, the affect of the magnetization on the CGPEs is investigated.