arXiv:2601.08990

Nonlinear Inverse Iterations for Spin-Orbit Coupled Quantum Gases

Patrick Henning, Laura Huynh

Published Jan 13, 2026• 0 citations

Abstract

This work concerns the computation of ground states of two-component spin-orbit coupled Bose-Einstein condensates (SO-coupled BECs), modelled by a coupled nonlinear eigenvalue problem of Gross-Pitaevskii type. In regimes with strongly modulated densities, including supersolid-like phases, conventional numerical approaches such as generalized inverse iterations or gradient descent may converge slowly. To address this issue, we apply the J-method [E. Jarlebring, S. Kvaal, W. Michiels. SIAM J. Sci. Comput. 36-4, 2014], which allows nonlinear inverse iterations to be accelerated through spectral shifting. A central analytical challenge arises from global phase invariance, which causes the target eigenvalue of the J-operator to be multiple and prevents the direct application of existing convergence results based on simple eigenvalues. Under a quasi-isolation assumption, we establish that this eigenvalue has algebraic and geometric multiplicity exactly two. By introducing an auxiliary phase-fixed iteration, we prove local linear convergence rates for the shifted J-method, determined by spectral gaps, in a quotient metric modulo global phase rotations. The particle densities converge at the same rate. Numerical experiments confirm the predicted linear convergence for fixed shifts and indicate superlinear convergence for adaptively chosen shifts.

DOI: 10.48550/arxiv.2601.08990

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 MathematicsOct 30, 20260 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.17950Feb 20, 20260 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.5815SIAM Journal on Applied MathematicsJan 1, 201545 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, 20140 citations

Qualitative and Quantitative Analysis of Riemannian Optimization Methods for Ground States of Rotating Multicomponent Bose-Einstein Condensates

Hermann, Martin, Tatjana Stykel, Mahima Yadav

We develop and analyze Riemannian optimization methods for computing ground states of rotating multicomponent Bose-Einstein condensates, defined as minimizers of the Gross-Pitaevskii energy functional. To resolve the non-uniqueness of ground states induced by phase invariance, we work on a quotient manifold endowed with a general Riemannian metric. By introducing an auxiliary phase-aligned iteration and employing fixed-point convergence theory, we establish a unified local convergence framework for Riemannian gradient descent methods and derive explicit convergence rates. Specializing this framework to two metrics tailored to the energy landscape, we study the energy-adaptive and Lagrangian-based Riemannian gradient descent methods. While monotone energy decay and global convergence are established only for the former, a quantified local convergence analysis is provided for both methods. Numerical experiments confirm the theoretical results and demonstrate that the Lagrangian-based method, which incorporates second-order information on the energy functional and mass constraints, achieves faster local convergence than the energy-adaptive scheme.

arXiv:2512.05939Dec 5, 20250 citations

Numerical Simulations on Ground States for Rotating Two-Component Bose-Einstein Condensates

荣华刘

静态相耦合的Gross-Pitaevskii方程组的解描述了双原子的玻色爱因斯坦凝聚体在极低温度下的基态现象。我们提出一种十分有效的数值方法——梯度法来求解此基态解。我们提出的梯度法在数值上既保持总模量守恒又能使总能量递减；我们严格地证明我们提出的梯度法是一种获得能量函数在给定限制性条件下的最小值(也即基态解)的十分有效的方法。我们通过大量的例子来显示该方法的优点并且应用该方法去研究处于旋转状态下的双原子–玻色爱因斯坦凝聚体在极低温度下所呈现的复杂涡旋现象。 The ground states of rotating two-component Bose-Einstein condensates (BEC) at extremely low temperature are solutions of time-independent coupled Gross-Pitaevskii equations. To compute the ground states, we propose an efficient numerical method—gradient flows with discrete nor-malization. In linear cases and under properly chosen initial data, we can prove that the gradient flows converge into the ground state which has the lowest energy. We show that the method is quite efficient and apply the method to study complicated vortex structure in the ground state solutions of rotating two-component BEC at extremely low temperature.

Advances in Applied MathematicsJan 1, 20170 citations