An efficient numerical method for computing dynamics of spin F = 2 Bose–Einstein condensates

In this paper, we extend the efficient time-splitting Fourier pseudospectral method to solve the generalized Gross–Pitaevskii (GP) equations, which model the dynamics of spin F = 2 Bose–Einstein condensates at extremely low temperature. Using the time-splitting technique, we split the generalized GP equations into one linear part and two nonlinear parts: the linear part is solved with the Fourier pseudospectral method; one of nonlinear parts is solved analytically while the other one is reformulated into a matrix formulation and solved by diagonalization. We show that the method keeps well the conservation laws related to generalized GP equations in 1D and 2D. We also show that the method is of second-order in time and spectrally accurate in space through a one-dimensional numerical test. We apply the method to investigate the dynamics of spin F = 2 Bose–Einstein condensates confined in a uniform/nonuniform magnetic field.