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

arXiv:2602.17950

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

Jing Wang, Wei Yang, Yongjun Yuan, Yong Zhang

Published Feb 20, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.2602.17950

Related papers

On the Ground States of Rotating Spin-Orbit Coupled Two-Component Bose-Einstein Condensates

Jing Wang, Wei Yang, Yong Zhang

In this paper, we study mathematically and numerically the ground states of rotating spin-orbit coupled two-component Bose-Einstein condensates modeled by the coupled Gross-Pitaevskii equations. Firstly, we rigorously prove the existence and nonexistence results for the ground states 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 confirm the spatial spectral accuracy by traversing the most commonly-used initial guesses for each component. Moreover, we investigate the effects of local interaction, angular velocity, spin-orbit coupling and external trapping potential on the ground states, and predict some interesting physical phenomena, such as vortex lattices.

CSIAM Transactions on Applied Mathematics•Oct 30, 2026•0 citations

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

A Preconditioned Riemannian Conjugate Gradient Method to Compute Ground States of Spin-1 Bose–Einstein Condensates

Wen Tao Huang, Zilin Yang, Qinglin Tang

Abstract. In this paper, we propose a preconditioned Riemannian conjugate gradient method for computing ground states of spin-1 Bose-Einstein condensates, which can be reformulated as an optimization problem with the mass and magnetization constraints. The energy functional and constraints can be discretized using Fourier pseudospectral schemes, thereby transforming the problem into an optimization task on a manifold. We derive vector transports by differentiating three existing retractions and propose an initial step size selection strategy based on the second-order approximation of the energy function for the proposed Riemannian optimization algorithm. In addition, a preconditioner is derived and used. To further accelerate the convergence, we combine the proposed algorithm with a multigrid method. Numerical experiments demonstrate the efficiency and accuracy of the proposed method while confirming the usefulness of our preconditioning technique and the effectiveness of the proposed strategy for initial step size selection.

SIAM Journal on Scientific Computing•Apr 1, 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

An efficient and accurate numerical method for computing the ground states of three-dimensional rotating dipolar Bose-Einstein condensates under strongly anisotropic trap

Qinglin Tang, Hanquan Wang, Shaobo Zhang +1 more

In this article, we propose an efficient and spectrally accurate numerical method to compute the ground states of three-dimensional (3D) rotating dipolar Bose-Einstein condensates (BEC) under strongly anisotropic trapping potentials.The kernel singularity, convolution non-locality and density anisotropy together complicate the dipolar potential evaluation. The fast rotation mechanism not only induces a complicated energy landscape with many local minima, but also creates a large number of vortices in the condensates. Such factors collectively make the ground state computation challenging in terms of convergence, accuracy and efficiency, especially for 3D anisotropic systems. Coupled with Fourier spectral discretization, we proposed a preconditioned conjugate gradient method (PCG) by integrating the anisotropic truncated kernel method (ATKM) for the dipolar potential evaluation. An adaptive step size control strategy is designed and ATKM allows for a spectral accuracy without introducing any extra anisotropy-dependent memory requirement or computational time. Our algorithm is spectrally accurate, highly efficient and memory-economic. Extensive numerical results are presented to confirm the accuracy and efficiency, together with applications to study impacts of the model parameters on critical rotational frequency, energies and chemical potential. Furthermore, these simulations reveal additional novel ground state patterns, such as bent vortices.

arXiv:2603.05194•Mar 5, 2026•0 citations

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

Sirilak Sriburadet, Yin-Tzer Shih, B.‐W. Jeng +2 more

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 [Formula: see text] and [Formula: see text] 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.

Scientific Reports•Nov 23, 2021•1 citation
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.
In this paper, we study mathematically and numerically the ground states of rotating spin-orbit coupled two-component Bose-Einstein condensates modeled by the coupled Gross-Pitaevskii equations. Firstly, we rigorously prove the existence and nonexistence results for the ground states 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 confirm the spatial spectral accuracy by traversing the most commonly-used initial guesses for each component. Moreover, we investigate the effects of local interaction, angular velocity, spin-orbit coupling and external trapping potential on the ground states, and predict some interesting physical phenomena, such as vortex lattices.
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.
Abstract. In this paper, we propose a preconditioned Riemannian conjugate gradient method for computing ground states of spin-1 Bose-Einstein condensates, which can be reformulated as an optimization problem with the mass and magnetization constraints. The energy functional and constraints can be discretized using Fourier pseudospectral schemes, thereby transforming the problem into an optimization task on a manifold. We derive vector transports by differentiating three existing retractions and propose an initial step size selection strategy based on the second-order approximation of the energy function for the proposed Riemannian optimization algorithm. In addition, a preconditioner is derived and used. To further accelerate the convergence, we combine the proposed algorithm with a multigrid method. Numerical experiments demonstrate the efficiency and accuracy of the proposed method while confirming the usefulness of our preconditioning technique and the effectiveness of the proposed strategy for initial step size selection.
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 article, we propose an efficient and spectrally accurate numerical method to compute the ground states of three-dimensional (3D) rotating dipolar Bose-Einstein condensates (BEC) under strongly anisotropic trapping potentials.The kernel singularity, convolution non-locality and density anisotropy together complicate the dipolar potential evaluation. The fast rotation mechanism not only induces a complicated energy landscape with many local minima, but also creates a large number of vortices in the condensates. Such factors collectively make the ground state computation challenging in terms of convergence, accuracy and efficiency, especially for 3D anisotropic systems. Coupled with Fourier spectral discretization, we proposed a preconditioned conjugate gradient method (PCG) by integrating the anisotropic truncated kernel method (ATKM) for the dipolar potential evaluation. An adaptive step size control strategy is designed and ATKM allows for a spectral accuracy without introducing any extra anisotropy-dependent memory requirement or computational time. Our algorithm is spectrally accurate, highly efficient and memory-economic. Extensive numerical results are presented to confirm the accuracy and efficiency, together with applications to study impacts of the model parameters on critical rotational frequency, energies and chemical potential. Furthermore, these simulations reveal additional novel ground state patterns, such as bent vortices.
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 [Formula: see text] and [Formula: see text] 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.