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

International Journal of Computer Mathematics

High-order numerical methods with mass and energy conservation for spin–orbit-coupled Bose–Einstein condensates

Xiang‐Gui Li, Shu‐Cun Li

Published Jun 15, 2020• 2 citations
Publisher

Abstract

In this work, a high-order numerical method with mass and energy conservation obtained by the operator-compensation (OC) technique is proposed for the coupled Gross–Pitaevskii equations (CGPEs) modelling the spin–orbit-coupled Bose–Einstein condensates. The OC technique is utilized to seek high-order approximations to the first- and the second-order spatial derivatives in the CGPEs, after which a time-dependent semi-discrete system is obtained. Then we obtain a fully discrete scheme by applying the Crank–Nicolson method to discretize this system in time. We theoretically show that the semi-discrete system and the fully discrete scheme conserve both mass and energy on the discrete level. In addition, the stability analysis is discussed for the proposed method. Numerical experiments are given to test accuracy and verify conservative properties. The dynamics of the two-dimensional CGPEs and the effects of the parameters in the CGPEs are simulated.

DOI
10.1080/00207160.2020.1777407

Related papers

Ground States and Dynamics of Spin-Orbit-Coupled Bose--Einstein Condensates

Weizhu Bao, Yongyong Cai

We study analytically and asymptotically, as well as numerically, ground states and dynamics of two-component spin-orbit-coupled Bose--Einstein condensates (BECs) modeled by the coupled Gross--Pitaevskii equations (CGPEs). In fact, due to the appearance of the spin-orbit (SO) coupling in the two-component BEC with a Raman coupling, the ground state structures and dynamical properties become very rich and complicated. For the ground states, we establish existence and nonexistence results under different parameter regimes and obtain their limiting behaviors and/or structures with different combinations of the SO and Raman coupling strengths. For the dynamics, we show that the motion of the center-of-mass is either nonperiodic or with different frequency than the trapping frequency when the external trapping potential is taken as harmonic and/or the initial data is chosen as a stationary state (e.g., ground state) with a shift, which is completely different from the case of a two-component BEC without the SO coupling, and obtain the semiclassical limit of the CGPEs in the linear case via the Wigner transform method. Efficient and accurate numerical methods are proposed for computing the ground states and dynamics, especially for the case of box potentials. Numerical results are reported to demonstrate the efficiency and accuracy of the numerical methods and show the rich phenomenon in the SO-coupled BECs.

arXiv:1407.5815•SIAM Journal on Applied Mathematics•Jan 1, 2015•45 citations

Mathematical and numerical study on the ground states of rotating spin-orbit coupled spin-1 Bose-Einstein condensates

Jing Wang, Wei Yang, Yongjun Yuan +1 more

In this article, we study mathematically and numerically the ground states of three-component rotating spin-orbit coupled (SOC) spin-1 Bose-Einstein condensates modeled by the coupled Gross-Pitaevskii equations (CGPEs). Firstly, we rigorously prove existence result of the ground state and derive some analytical properties, including the virial identity and negativity of SOC energy. Secondly, we propose an efficient and accurate preconditioned nonlinear conjugate gradient (PCG) algorithm to compute the ground states. We truncate the whole space into a bounded rectangular domain and readily apply the Fourier spectral method to approximate the wave function. The PCG method is successfully adapted with appropriate modifications to the adaptive step size control strategy for the one-parameter energy minimization problem and to the choice of preconditioners, achieving great performance in terms of accuracy and efficiency. Lastly, we carry out extensive numerical experiments to verify the existence and property results of the ground states, confirm the spatial spectral accuracy by traversing the most commonly-used initial guesses for each component thanks to its great efficiency, which is also attributed to a utilization of cascadic multigrid and discrete Fast Fourier Transform (FFT). Moreover, we investigate the effects of local interaction, rotation and spin-orbit coupling and external trapping potential on the ground state, and unveil some interesting physical phenomena, such as giant vortex and U-shape vortex line.

arXiv:2602.17950•Feb 20, 2026•0 citations

GROUND STATES AND DYNAMICS OF SPIN-ORBIT-COUPLED

Bose-Einstein Condensates, Yongyong Cai

We study analytically and asymptotically as well as numerically ground states and dynamics of two-component spin-orbit-coupled Bose-Einstein condensates (BECs) modeled by the coupled Gross-Pitaevskii equations (CGPEs). In fact, due to the appearance of the spin-orbit (SO) coupling in the two-component BEC with a Raman coupling, the ground state structures and dynam- ical properties become very rich and complicated. For the ground states, we establish the existence and non-existence results under different parameter regimes, and obtain their limiting behaviors and/or structures with different combinations of the SO and Raman coupling strengths. For the dynamics, we show that the motion of the center-of-mass is either non-periodic or with different frequency to the trapping frequency when the external trapping potential is taken as harmonic and the initial data is chosen as a stationary state (e.g. ground state) with a shift, which is completely different from the case of a two-component BEC without the SO coupling, and obtain the semiclas- sical limit of the CGPEs in the linear case via the Wigner transform method. Efficient and accurate numerical methods are proposed for computing the ground states and dynamics, especially for the case of box potentials. Numerical results are reported to demonstrate the efficiency and accuracy of the numerical methods and show the rich phenomenon in the SO-coupled BECs.

•Jan 1, 2014•0 citations

A Mass and Magnetization Conservative and Energy-Diminishing Numerical Method for Computing Ground State of Spin-1 Bose–Einstein Condensates

Weizhu Bao, Hanquan Wang

In this paper, a mass (or normalization) and magnetization conservative and energy-diminishing numerical method is presented for computing the ground state of spin-1 (or $F=1$ spinor) Bose–Einstein condensates (BECs). We begin with the coupled Gross–Pitaevskii equations, and the ground state is defined as the minimizer of the energy functional under two constraints on the mass and magnetization. By constructing a continuous normalized gradient flow (CNGF) which is mass and magnetization conservative and energy-diminishing, the ground state can be computed as the steady state solution of the CNGF. The CNGF is then discretized by the Crank–Nicolson finite difference method with a proper way to deal with the nonlinear terms, and we prove that the discretization is mass and magnetization conservative and energy-diminishing in the discretized level. Numerical results of the ground state and their energy of spin-1 BECs are reported to demonstrate the efficiency of the numerical method.

SIAM Journal on Numerical Analysis•Jan 1, 2007•48 citations

Semi-implicit finite-difference methods to study the spin-orbit- and coherently-coupled spinor Bose–Einstein condensates

Paramjeet Banger, Pardeep Kaur, Sandeep Gautam

We develop time-splitting finite-difference methods, using an implicit Backward–Euler and a semi-implicit Crank–Nicolson discretization schemes, to study the spin-orbit-coupled (SO-coupled) spinor Bose–Einstein condensates with coherent coupling in quasi-one- and quasi-two-dimensional traps. The split equations involving kinetic energy and spin-orbit-coupling operators are solved using either a time implicit Backward–Euler or a semi-implicit Crank–Nicolson method. We explicitly develop the methods for pseudospin-1/2, spin-1 and spin-2 condensates. The results for ground states obtained with time-splitting Backward–Euler and Crank–Nicolson methods are in excellent agreement with time-splitting Fourier spectral method, which is one of the popular methods to solve the mean-field models for SO-coupled spinor condensates. We confirm the emergence of different phases in SO-coupled pseudospin-1/2, spin-1 and spin-2 condensates with coherent coupling.

arXiv:2106.07762•International Journal of Modern Physics C•Sep 22, 2021•0 citations

Semi-implicit finite-difference methods to study the spin-orbit- and coherently-coupled spinor Bose–Einstein condensates

Paramjeet Banger, Pardeep Kaur, Sandeep Gautam

We develop time-splitting finite-difference methods, using an implicit Backward–Euler and a semi-implicit Crank–Nicolson discretization schemes, to study the spin-orbit-coupled (SO-coupled) spinor Bose–Einstein condensates with coherent coupling in quasi-one- and quasi-two-dimensional traps. The split equations involving kinetic energy and spin-orbit-coupling operators are solved using either a time implicit Backward–Euler or a semi-implicit Crank–Nicolson method. We explicitly develop the methods for pseudospin-1/2, spin-1 and spin-2 condensates. The results for ground states obtained with time-splitting Backward–Euler and Crank–Nicolson methods are in excellent agreement with time-splitting Fourier spectral method, which is one of the popular methods to solve the mean-field models for SO-coupled spinor condensates. We confirm the emergence of different phases in SO-coupled pseudospin-1/2, spin-1 and spin-2 condensates with coherent coupling.

•Jan 1, 2022•10 citations
In this work, a high-order numerical method with mass and energy conservation obtained by the operator-compensation (OC) technique is proposed for the coupled Gross–Pitaevskii equations (CGPEs) modelling the spin–orbit-coupled Bose–Einstein condensates. The OC technique is utilized to seek high-order approximations to the first- and the second-order spatial derivatives in the CGPEs, after which a time-dependent semi-discrete system is obtained. Then we obtain a fully discrete scheme by applying the Crank–Nicolson method to discretize this system in time. We theoretically show that the semi-discrete system and the fully discrete scheme conserve both mass and energy on the discrete level. In addition, the stability analysis is discussed for the proposed method. Numerical experiments are given to test accuracy and verify conservative properties. The dynamics of the two-dimensional CGPEs and the effects of the parameters in the CGPEs are simulated.
We study analytically and asymptotically, as well as numerically, ground states and dynamics of two-component spin-orbit-coupled Bose--Einstein condensates (BECs) modeled by the coupled Gross--Pitaevskii equations (CGPEs). In fact, due to the appearance of the spin-orbit (SO) coupling in the two-component BEC with a Raman coupling, the ground state structures and dynamical properties become very rich and complicated. For the ground states, we establish existence and nonexistence results under different parameter regimes and obtain their limiting behaviors and/or structures with different combinations of the SO and Raman coupling strengths. For the dynamics, we show that the motion of the center-of-mass is either nonperiodic or with different frequency than the trapping frequency when the external trapping potential is taken as harmonic and/or the initial data is chosen as a stationary state (e.g., ground state) with a shift, which is completely different from the case of a two-component BEC without the SO coupling, and obtain the semiclassical limit of the CGPEs in the linear case via the Wigner transform method. Efficient and accurate numerical methods are proposed for computing the ground states and dynamics, especially for the case of box potentials. Numerical results are reported to demonstrate the efficiency and accuracy of the numerical methods and show the rich phenomenon in the SO-coupled BECs.
In this article, we study mathematically and numerically the ground states of three-component rotating spin-orbit coupled (SOC) spin-1 Bose-Einstein condensates modeled by the coupled Gross-Pitaevskii equations (CGPEs). Firstly, we rigorously prove existence result of the ground state and derive some analytical properties, including the virial identity and negativity of SOC energy. Secondly, we propose an efficient and accurate preconditioned nonlinear conjugate gradient (PCG) algorithm to compute the ground states. We truncate the whole space into a bounded rectangular domain and readily apply the Fourier spectral method to approximate the wave function. The PCG method is successfully adapted with appropriate modifications to the adaptive step size control strategy for the one-parameter energy minimization problem and to the choice of preconditioners, achieving great performance in terms of accuracy and efficiency. Lastly, we carry out extensive numerical experiments to verify the existence and property results of the ground states, confirm the spatial spectral accuracy by traversing the most commonly-used initial guesses for each component thanks to its great efficiency, which is also attributed to a utilization of cascadic multigrid and discrete Fast Fourier Transform (FFT). Moreover, we investigate the effects of local interaction, rotation and spin-orbit coupling and external trapping potential on the ground state, and unveil some interesting physical phenomena, such as giant vortex and U-shape vortex line.
We study analytically and asymptotically as well as numerically ground states and dynamics of two-component spin-orbit-coupled Bose-Einstein condensates (BECs) modeled by the coupled Gross-Pitaevskii equations (CGPEs). In fact, due to the appearance of the spin-orbit (SO) coupling in the two-component BEC with a Raman coupling, the ground state structures and dynam- ical properties become very rich and complicated. For the ground states, we establish the existence and non-existence results under different parameter regimes, and obtain their limiting behaviors and/or structures with different combinations of the SO and Raman coupling strengths. For the dynamics, we show that the motion of the center-of-mass is either non-periodic or with different frequency to the trapping frequency when the external trapping potential is taken as harmonic and the initial data is chosen as a stationary state (e.g. ground state) with a shift, which is completely different from the case of a two-component BEC without the SO coupling, and obtain the semiclas- sical limit of the CGPEs in the linear case via the Wigner transform method. Efficient and accurate numerical methods are proposed for computing the ground states and dynamics, especially for the case of box potentials. Numerical results are reported to demonstrate the efficiency and accuracy of the numerical methods and show the rich phenomenon in the SO-coupled BECs.
In this paper, a mass (or normalization) and magnetization conservative and energy-diminishing numerical method is presented for computing the ground state of spin-1 (or F=1F=1F=1 spinor) Bose–Einstein condensates (BECs). We begin with the coupled Gross–Pitaevskii equations, and the ground state is defined as the minimizer of the energy functional under two constraints on the mass and magnetization. By constructing a continuous normalized gradient flow (CNGF) which is mass and magnetization conservative and energy-diminishing, the ground state can be computed as the steady state solution of the CNGF. The CNGF is then discretized by the Crank–Nicolson finite difference method with a proper way to deal with the nonlinear terms, and we prove that the discretization is mass and magnetization conservative and energy-diminishing in the discretized level. Numerical results of the ground state and their energy of spin-1 BECs are reported to demonstrate the efficiency of the numerical method.
We develop time-splitting finite-difference methods, using an implicit Backward–Euler and a semi-implicit Crank–Nicolson discretization schemes, to study the spin-orbit-coupled (SO-coupled) spinor Bose–Einstein condensates with coherent coupling in quasi-one- and quasi-two-dimensional traps. The split equations involving kinetic energy and spin-orbit-coupling operators are solved using either a time implicit Backward–Euler or a semi-implicit Crank–Nicolson method. We explicitly develop the methods for pseudospin-1/2, spin-1 and spin-2 condensates. The results for ground states obtained with time-splitting Backward–Euler and Crank–Nicolson methods are in excellent agreement with time-splitting Fourier spectral method, which is one of the popular methods to solve the mean-field models for SO-coupled spinor condensates. We confirm the emergence of different phases in SO-coupled pseudospin-1/2, spin-1 and spin-2 condensates with coherent coupling.
We develop time-splitting finite-difference methods, using an implicit Backward–Euler and a semi-implicit Crank–Nicolson discretization schemes, to study the spin-orbit-coupled (SO-coupled) spinor Bose–Einstein condensates with coherent coupling in quasi-one- and quasi-two-dimensional traps. The split equations involving kinetic energy and spin-orbit-coupling operators are solved using either a time implicit Backward–Euler or a semi-implicit Crank–Nicolson method. We explicitly develop the methods for pseudospin-1/2, spin-1 and spin-2 condensates. The results for ground states obtained with time-splitting Backward–Euler and Crank–Nicolson methods are in excellent agreement with time-splitting Fourier spectral method, which is one of the popular methods to solve the mean-field models for SO-coupled spinor condensates. We confirm the emergence of different phases in SO-coupled pseudospin-1/2, spin-1 and spin-2 condensates with coherent coupling.