Yongjun Yuan

Xin Liu, Ziqing Xie, Yongjun Yuan et al.

This paper investigates the dynamics of spin-2 Bose-Einstein condensates (BECs) with rotation and spin-orbit coupling (SOC). In order to better simulate the dynamics, we present an efficient high-order compact splitting Fourier spectral method. This method splits the Hamiltonian into a linear part, which consists of the Laplace, rotation and SOC terms, and a nonlinear part that includes all the remaining terms. The wave function is well approximated by the Fourier spectral method and is numerically accessed with discrete Fast Fourier transform (FFT). For linear subproblem, the handling of rotation term and SOC term poses a major challenge. Using a function mapping based on rotation, we can integrate the linear subproblem exactly and explicitly. This mapping we propose not only helps eliminate the rotation term, but also prevents the SOC term from evolving into a time-dependent form. The nonlinear subproblem is integrated analytically in physical space. Such "compact" splitting involves only two operators and facilitates the design of high-order splitting schemes. Our method is spectrally accurate in space and high order in time. It is efficient, explicit, unconditionally stable and simple to implement. In addition, we derive some dynamical properties and carry out a systematic study, including accuracy and efficiency tests, dynamical property verification, the SOC effects and dynamics of vortex lattice.

Wei Liu, Tingfeng Wang, Yongjun Yuan et al.

This work focuses on the numerical computation of defocusing action ground states for rotating nonlinear Schrödinger equations (RNLS) using a direct gradient flow (DGF) method. We address theoretical gaps in the existing literature concerning the stability and convergence of this DGF scheme. Firstly, we prove the unconditional stability of the DGF scheme, demonstrating that the action functional is monotonically non-increasing along the discrete flow for arbitrary time step sizes. Secondly, we establish a rigorous convergence analysis, proving global convergence under minor assumptions and local exponential convergence to the action ground state under a reasonable non-degeneracy condition. The analysis relies on the uniform boundedness of sublevel sets of the action functional and introduces a tailored $H^1$-distance between phase-shift equivalence classes to handle complex-valued ground states with quantized vortices. A novel analytical framework is also developed to establish the exponential convergence rate. Numerical experiments are presented to validate the theoretical findings, demonstrating both the global migration towards a neighborhood of the ground state and subsequent exponential convergence.