Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Mathematical Models and Methods in Applied Sciences

High-order compact splitting spectral methods for the rotating spin-1 Bose–Einstein condensates in a magnetic field

Xin Liu, Xiangyu Meng, Qinglin Tang, Yong Zhang

Published Apr 11, 2025• 4 citations
Publisher

Abstract

We propose high-order compact splitting Fourier spectral methods to compute the dynamics of rotating spin-1 Bose–Einstein condensates (BECs) under an Ioffe-Pritchard magnetic field. We split the Hamiltonian into a linear part, which consists of the Laplace, rotation and Zeeman energy terms, and a nonlinear part, which includes all the remaining terms. This “compact” splitting, involving only two operators, significantly simplifies the construction of high-order schemes. For the linear subproblem, we factorize the linear operator as a product of simple sub-operators using the exact classical-quantum correspondence and semigroups decomposition, and such sub-operators are well approximated by Fourier spectral method and integrated exactly in phase space as usual. Importantly and surprisingly, the splitting coefficients, originally determined by a nonlinear equation, can actually be computed exactly by solving a linear system, ending up with explicit formulas. For the nonlinear subproblem, we derive an exact formula of the spin vector, which is not conserved in presence of the magnetic field, and transform it into a linear problem. Then we design explicit high-order schemes using Magnus integrators. Our scheme achieves spectral accuracy in space and high-order precision in time with near-optimal efficiency. It is explicit, easy to implement and unconditionally stable. In addition, we derive some properties of our numerical scheme and conduct a comprehensive investigation, including accuracy confirmation, efficiency test, property verification, interaction of quantized vortices and dynamics under honeycomb potential.

DOI
10.1142/s0218202525500253

Related papers

An Efficient High-Order Compact Splitting Spectral Method for Dipolar Bose-Einstein Condensates with Arbitrary-Angle Rotation

Xin Liu, Yong Zhang

Abstract. We propose high-order compact splitting Fourier spectral methods to compute the dynamics of dipolar Bose–Einstein condensates (BECs) with arbitrary-angle rotation. The Hamiltonian is split into a linear part that consists of the Laplace and rotation and a nonlinear part that includes all the remaining terms. For the linear subproblem, we factorize the linear operator into a product of simple suboperators using the exact classical-quantum correspondence and semigroups decomposition, and these suboperators are integrated exactly in either physical or phase space. Importantly, we first derive explicit formulas for the splitting coefficients originally obtained by numerically solving a nonlinear system. The nonlinear part is solved exactly in physical space, and the convolution-type nonlocal potential with a singular kernel is computed via the Far-field Smooth Approximation (FSA) method with spectral accuracy. This “compact” splitting, involving only two operators, significantly simplifies the design of high-order schemes. Our method achieves high-order temporal precision and spectral spatial accuracy with near-optimal efficiency, and it is explicit, easy to implement and unconditionally stable. In addition, we derive some dynamical properties and conduct a comprehensive investigation, including accuracy confirmation, property verification, interaction of quantized vortices, and dynamics under honeycomb potential.

SIAM Journal on Scientific Computing•Sep 17, 2025•1 citation

An efficient compact splitting Fourier spectral methods for computing the dynamics of rotating spin-orbit coupled spin-2 Bose-Einstein condenstates

Xin Liu, Ziqing Xie, Yongjun Yuan +2 more

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.

arXiv:2603.25282•Mar 26, 2026•0 citations

An Efficient Time-Splitting Method to Simulate the Dynamics of Spin-Orbit Coupled Spin-1 Bose-Einstein Condensates

Xiaolan Chen, Yongjun Yuan, Yong Zhang

In this paper, an efficient time-splitting Fourier spectral method is proposed to simulate the dynamics of spin-orbit coupled spin-1 Bose-Einstein condensates (SOC spin-1 BECs). We split the Hamiltonian into a linear part, which consists of the Laplace and SOC terms, and a nonlinear part that includes all the remaining terms. The linear subproblem is integrated analytically in phase space by solving an ordinary differential system of constant coefficient matrix. While, for the nonlinear subproblem, it is proved the coefficient matrix is actually time-independent in physical space, therefore, the nonlinear subproblem can be integrated exactly. Based on such two-step splitting, we construct high-order schemes to simulate the dynamics. 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 for SOC spin-1 BECs. Extensive numerical results are presented to confirm the accuracy and efficiency, illustrate the dynamical properties at discrete level, and show interesting physics of SOC spin-1 BECs, such as the SOC effects and different wave patterns.

Communications in Computational Physics•Nov 7, 2025•1 citation

A Generalized-Laguerre–Fourier–Hermite Pseudospectral Method for Computing the Dynamics of Rotating Bose–Einstein Condensates

Weizhu Bao, Hai-Liang Li, Jie Shen

A time-splitting generalized-Laguerre–Fourier–Hermite pseudospectral method is proposed for computing the dynamics of rotating Bose–Einstein condensates (BECs) in two and three dimensions. The new numerical method is based on the following: (i) the use of a time-splitting technique for decoupling the nonlinearity; (ii) the adoption of polar coordinates in two dimensions and cylindrical coordinates in three dimensions such that the angular rotation term becomes constant coefficient; and (iii) the construction of eigenfunctions for the linear operator by properly scaling the generalized-Laguerre, Fourier, and Hermite functions. The new method enjoys the following properties: (i) it is explicit, time reversible, and time transverse invariant; (ii) it conserves the position density and is spectrally accurate in space and second-order or fourth-order accurate in time; and (iii) it solves the problem in the original whole space instead of in a truncated artificial computational domain. The method is also extended to solve the coupled Gross–Pitaevskii equations for the dynamics of rotating two-component and spin-1 BECs. Extensive numerical results for the dynamics of BECs are reported to demonstrate the accuracy and efficiency of the new method for rotating BECs.

SIAM Journal on Scientific Computing•Jan 1, 2009•81 citations

A time-splitting spectral method for computing dynamics of spinor F=1 Bose–Einstein condensates

Hanquan Wang

We propose a time-splitting spectral method for the generalized Gross–Pitaevskii equations, which describe the dynamics of spinor F=1 Bose–Einstein condensates at a very low temperature. The new numerical method is explicit, unconditionally stable, and of spectral accuracy in space. Moreover, it conserves the position densities at the discretized level. We apply the method for studying both the dynamic generation of vortices and vortex lattice dynamics for the spinor F=1 Bose–Einstein condensates held in an Ioffe-Pritchard magnetic field.

International Journal of Computer Mathematics•Jun 1, 2007•32 citations

A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose--Einstein Condensates via Rotating Lagrangian Coordinates

Weizhu Bao, Daniel Marahrens, Qinglin Tang +1 more

We propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose--Einstein condensates (BECs) in a rotational frame with or without long-range dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.

SIAM Journal on Scientific Computing•Jan 1, 2013•42 citations
We propose high-order compact splitting Fourier spectral methods to compute the dynamics of rotating spin-1 Bose–Einstein condensates (BECs) under an Ioffe-Pritchard magnetic field. We split the Hamiltonian into a linear part, which consists of the Laplace, rotation and Zeeman energy terms, and a nonlinear part, which includes all the remaining terms. This “compact” splitting, involving only two operators, significantly simplifies the construction of high-order schemes. For the linear subproblem, we factorize the linear operator as a product of simple sub-operators using the exact classical-quantum correspondence and semigroups decomposition, and such sub-operators are well approximated by Fourier spectral method and integrated exactly in phase space as usual. Importantly and surprisingly, the splitting coefficients, originally determined by a nonlinear equation, can actually be computed exactly by solving a linear system, ending up with explicit formulas. For the nonlinear subproblem, we derive an exact formula of the spin vector, which is not conserved in presence of the magnetic field, and transform it into a linear problem. Then we design explicit high-order schemes using Magnus integrators. Our scheme achieves spectral accuracy in space and high-order precision in time with near-optimal efficiency. It is explicit, easy to implement and unconditionally stable. In addition, we derive some properties of our numerical scheme and conduct a comprehensive investigation, including accuracy confirmation, efficiency test, property verification, interaction of quantized vortices and dynamics under honeycomb potential.
Abstract. We propose high-order compact splitting Fourier spectral methods to compute the dynamics of dipolar Bose–Einstein condensates (BECs) with arbitrary-angle rotation. The Hamiltonian is split into a linear part that consists of the Laplace and rotation and a nonlinear part that includes all the remaining terms. For the linear subproblem, we factorize the linear operator into a product of simple suboperators using the exact classical-quantum correspondence and semigroups decomposition, and these suboperators are integrated exactly in either physical or phase space. Importantly, we first derive explicit formulas for the splitting coefficients originally obtained by numerically solving a nonlinear system. The nonlinear part is solved exactly in physical space, and the convolution-type nonlocal potential with a singular kernel is computed via the Far-field Smooth Approximation (FSA) method with spectral accuracy. This “compact” splitting, involving only two operators, significantly simplifies the design of high-order schemes. Our method achieves high-order temporal precision and spectral spatial accuracy with near-optimal efficiency, and it is explicit, easy to implement and unconditionally stable. In addition, we derive some dynamical properties and conduct a comprehensive investigation, including accuracy confirmation, property verification, interaction of quantized vortices, and dynamics under honeycomb potential.
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.
In this paper, an efficient time-splitting Fourier spectral method is proposed to simulate the dynamics of spin-orbit coupled spin-1 Bose-Einstein condensates (SOC spin-1 BECs). We split the Hamiltonian into a linear part, which consists of the Laplace and SOC terms, and a nonlinear part that includes all the remaining terms. The linear subproblem is integrated analytically in phase space by solving an ordinary differential system of constant coefficient matrix. While, for the nonlinear subproblem, it is proved the coefficient matrix is actually time-independent in physical space, therefore, the nonlinear subproblem can be integrated exactly. Based on such two-step splitting, we construct high-order schemes to simulate the dynamics. 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 for SOC spin-1 BECs. Extensive numerical results are presented to confirm the accuracy and efficiency, illustrate the dynamical properties at discrete level, and show interesting physics of SOC spin-1 BECs, such as the SOC effects and different wave patterns.
A time-splitting generalized-Laguerre–Fourier–Hermite pseudospectral method is proposed for computing the dynamics of rotating Bose–Einstein condensates (BECs) in two and three dimensions. The new numerical method is based on the following: (i) the use of a time-splitting technique for decoupling the nonlinearity; (ii) the adoption of polar coordinates in two dimensions and cylindrical coordinates in three dimensions such that the angular rotation term becomes constant coefficient; and (iii) the construction of eigenfunctions for the linear operator by properly scaling the generalized-Laguerre, Fourier, and Hermite functions. The new method enjoys the following properties: (i) it is explicit, time reversible, and time transverse invariant; (ii) it conserves the position density and is spectrally accurate in space and second-order or fourth-order accurate in time; and (iii) it solves the problem in the original whole space instead of in a truncated artificial computational domain. The method is also extended to solve the coupled Gross–Pitaevskii equations for the dynamics of rotating two-component and spin-1 BECs. Extensive numerical results for the dynamics of BECs are reported to demonstrate the accuracy and efficiency of the new method for rotating BECs.
We propose a time-splitting spectral method for the generalized Gross–Pitaevskii equations, which describe the dynamics of spinor F=1 Bose–Einstein condensates at a very low temperature. The new numerical method is explicit, unconditionally stable, and of spectral accuracy in space. Moreover, it conserves the position densities at the discretized level. We apply the method for studying both the dynamic generation of vortices and vortex lattice dynamics for the spinor F=1 Bose–Einstein condensates held in an Ioffe-Pritchard magnetic field.
We propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose--Einstein condensates (BECs) in a rotational frame with or without long-range dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.