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

SIAM Journal on Numerical Analysis

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

Weizhu Bao, Hanquan Wang

Published Jan 1, 2007• 48 citations
Publisher

Abstract

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.

DOI
10.1137/070681624

Related papers

A Normalized Gradient Flow Method for Computing Ground States of Spin-2 Bose–Einstein Condensates

Weizhu Bao, Qinglin Tang, Yongjun Yuan

Abstract. We propose and analyze an efficient and accurate numerical method for computing ground states of spin-2 Bose–Einstein condensates (BECs) by using the normalized gradient flow (NGF). In order to successfully extend the NGF to spin-2 BECs which have five components in the vector wave function but with only two physical constraints on total mass conservation and magnetization conservation, two important techniques are introduced for designing the proposed numerical method. The first one is to systematically investigate the ground state structure and property of spin-2 BECs within a spatially uniform system, which can be used for how to properly choose initial data in the NGF for computing ground states of spin-2 BECs. The second one is to introduce three additional projection conditions based on the relations between the chemical potentials, together with the two existing physical constraints, such that the five projection parameters used in the projection step of the NGF can be uniquely determined. Then a backward-forward Euler finite difference method is adapted to discretize the NGF. We prove rigorously that there exists a unique solution of the nonlinear system for determining the five projection parameters in the full discretization of the NGF under a mild condition on the time step size. Extensive numerical results on ground states of spin-2 BECs with different types of phases and under different potentials are reported to show the efficiency and accuracy of the proposed numerical method and to demonstrate several interesting physical phenomena on ground states of spin-2 BECs.

arXiv:2407.14441•SIAM Journal on Scientific Computing•Jun 6, 2025•2 citations

Computing the Ground State Solution of Bose--Einstein Condensates by a Normalized Gradient Flow

Weizhu Bao, Qiang Du

In this paper, we present a continuous normalized gradient flow (CNGF) and prove its energy diminishing property, which provides a mathematical justification of the imaginary time method used in the physics literature to compute the ground state solution of Bose--Einstein condensates (BEC). We also investigate the energy diminishing property for the discretization of the CNGF. Two numerical methods are proposed for such discretizations: one is the backward Euler centered finite difference (BEFD) method, the other is an explicit time-splitting sine-spectral (TSSP) method. Energy diminishing for BEFD and TSSP for the linear case and monotonicity for BEFD for both linear and nonlinear cases are proven. Comparison between the two methods and existing methods, e.g., Crank--Nicolson finite difference (CNFD) or forward Euler finite difference (FEFD), shows that BEFD and TSSP are much better in terms of preserving the energy diminishing property of the CNGF. Numerical results in one, two, and three dimensions with magnetic trap confinement potential, as well as a potential of a stirrer corresponding to a far-blue detuned Gaussian laser beam, are reported to demonstrate the effectiveness of BEFD and TSSP methods. Furthermore we observe that the CNGF and its BEFD discretization can also be applied directly to compute the first excited state solution in BEC when the initial data is chosen as an odd function.

arXiv:cond-mat/0303241•SIAM Journal on Scientific Computing•Jan 1, 2004•492 citations

Computing Ground States of Spin-1 Bose–Einstein Condensates by the Normalized Gradient Flow

Weizhu Bao, Fong Yin Lim

In this paper, we propose an efficient and accurate numerical method for computing the ground state of spin-1 Bose–Einstein condensates (BECs) by using the normalized gradient flow or imaginary time method. The key idea is to find a third projection or normalization condition based on the relation between the chemical potentials so that the three projection parameters used in the projection step of the normalized gradient flow are uniquely determined by this condition as well as the other two physical conditions given by the conservation of total mass and total magnetization. This allows us to successfully extend the most popular and powerful normalized gradient flow or imaginary time method for computing the ground state of a single-component BEC to compute the ground state of spin-1 BECs. An efficient and accurate discretization scheme, the backward-forward Euler sine-pseudospectral method, is proposed to discretize the normalized gradient flow. Extensive numerical results on ground states of spin-1 BECs with ferromagnetic/antiferromagnetic interaction and harmonic/optical lattice potential in one/three dimensions are reported to demonstrate the efficiency of our new numerical method.

arXiv:0711.0568•SIAM Journal on Scientific Computing•Jan 1, 2008•90 citations

Numerical methods for computing the ground state of spin-1 Bose-Einstein condensates in a uniform magnetic field

Fong Yin Lim, Weizhu Bao

We propose efficient and accurate numerical methods for computing the ground-state solution of spin-1 Bose-Einstein condensates subjected to a uniform magnetic field. The key idea in designing the numerical method is based on the normalized gradient flow with the introduction of a third normalization condition, together with two physical constraints on the conservation of total mass and conservation of total magnetization. Different treatments of the Zeeman energy terms are found to yield different numerical accuracies and stabilities. Numerical comparison between different numerical schemes is made, and the best scheme is identified. The numerical scheme is then applied to compute the condensate ground state in a harmonic plus optical lattice potential, and the effect of the periodic potential, in particular to the relative population of each hyperfine component, is investigated through comparison to the condensate ground state in a pure harmonic trap.

Physical Review E•Dec 30, 2008•50 citations

The Numerical Study of the Ground States of Spin-1 Bose-Einstein Condensates with Spin-Orbit-Coupling

Yongjun Yuan, Zhiguo Xu, Qinglin Tang +1 more

A projection gradient method for computing the ground state of the spin-orbit-coupled spin-1 Bose-Einstein condensate at extremely low temperatures is proposed. The continuous gradient flows are discretised by a second-order finite difference method in space and the Crank-Nicolson method in time. Our discretisation preserves the total mass conservation and the energy diminishing property. Numerical results show the efficiency of the method.

East Asian Journal on Applied Mathematics•Sep 17, 2018•7 citations

Ground States and Dynamics of Multicomponent Bose--Einstein Condensates

Weizhu Bao

We study numerically the time-independent vector Gross--Pitaevskii equations (VGPEs) for ground states and time-dependent VGPEs with (or without) an external driven field for dynamics describing a multicomponent Bose--Einstein condensate (BEC) at zero or a very low temperature. In preparation for the numerics, we scale the three-dimensional (3d) VGPEs, approximately reduce it to lower dimensions, present a continuous normalized gradient flow (CNGF) to compute ground states of multicomponent BEC, prove energy diminishing of the CNGF, which provides a mathematical justification, and discretize it by the backward Euler finite difference (BEFD), which is monotone in linear and nonlinear cases and preserves energy diminishing property in the linear case. Then we use a time-splitting sine-spectral (TSSP) method to discretize the time-dependent VGPEs with an external driven field for computing dynamics of multicomponent BEC. The merits of the TSSP method for VGPEs are that it is explicit, unconditionally stable, time reversible and time transverse invariant if the VGPEs is, has "good" resolution in the semiclassical regime, is of spectral-order accuracy in space and second-order accuracy in time, and conserves the total particle number in the discretized level. Extensive numerical examples in three dimensions for ground states and dynamics of multicomponent BEC are presented to demonstrate the power of the numerical methods and to discuss the physics of multicomponent BEC.

arXiv:cond-mat/0305309•Multiscale Modeling and Simulation•Jan 1, 2004•160 citations
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.
Abstract. We propose and analyze an efficient and accurate numerical method for computing ground states of spin-2 Bose–Einstein condensates (BECs) by using the normalized gradient flow (NGF). In order to successfully extend the NGF to spin-2 BECs which have five components in the vector wave function but with only two physical constraints on total mass conservation and magnetization conservation, two important techniques are introduced for designing the proposed numerical method. The first one is to systematically investigate the ground state structure and property of spin-2 BECs within a spatially uniform system, which can be used for how to properly choose initial data in the NGF for computing ground states of spin-2 BECs. The second one is to introduce three additional projection conditions based on the relations between the chemical potentials, together with the two existing physical constraints, such that the five projection parameters used in the projection step of the NGF can be uniquely determined. Then a backward-forward Euler finite difference method is adapted to discretize the NGF. We prove rigorously that there exists a unique solution of the nonlinear system for determining the five projection parameters in the full discretization of the NGF under a mild condition on the time step size. Extensive numerical results on ground states of spin-2 BECs with different types of phases and under different potentials are reported to show the efficiency and accuracy of the proposed numerical method and to demonstrate several interesting physical phenomena on ground states of spin-2 BECs.
In this paper, we present a continuous normalized gradient flow (CNGF) and prove its energy diminishing property, which provides a mathematical justification of the imaginary time method used in the physics literature to compute the ground state solution of Bose--Einstein condensates (BEC). We also investigate the energy diminishing property for the discretization of the CNGF. Two numerical methods are proposed for such discretizations: one is the backward Euler centered finite difference (BEFD) method, the other is an explicit time-splitting sine-spectral (TSSP) method. Energy diminishing for BEFD and TSSP for the linear case and monotonicity for BEFD for both linear and nonlinear cases are proven. Comparison between the two methods and existing methods, e.g., Crank--Nicolson finite difference (CNFD) or forward Euler finite difference (FEFD), shows that BEFD and TSSP are much better in terms of preserving the energy diminishing property of the CNGF. Numerical results in one, two, and three dimensions with magnetic trap confinement potential, as well as a potential of a stirrer corresponding to a far-blue detuned Gaussian laser beam, are reported to demonstrate the effectiveness of BEFD and TSSP methods. Furthermore we observe that the CNGF and its BEFD discretization can also be applied directly to compute the first excited state solution in BEC when the initial data is chosen as an odd function.
In this paper, we propose an efficient and accurate numerical method for computing the ground state of spin-1 Bose–Einstein condensates (BECs) by using the normalized gradient flow or imaginary time method. The key idea is to find a third projection or normalization condition based on the relation between the chemical potentials so that the three projection parameters used in the projection step of the normalized gradient flow are uniquely determined by this condition as well as the other two physical conditions given by the conservation of total mass and total magnetization. This allows us to successfully extend the most popular and powerful normalized gradient flow or imaginary time method for computing the ground state of a single-component BEC to compute the ground state of spin-1 BECs. An efficient and accurate discretization scheme, the backward-forward Euler sine-pseudospectral method, is proposed to discretize the normalized gradient flow. Extensive numerical results on ground states of spin-1 BECs with ferromagnetic/antiferromagnetic interaction and harmonic/optical lattice potential in one/three dimensions are reported to demonstrate the efficiency of our new numerical method.
We propose efficient and accurate numerical methods for computing the ground-state solution of spin-1 Bose-Einstein condensates subjected to a uniform magnetic field. The key idea in designing the numerical method is based on the normalized gradient flow with the introduction of a third normalization condition, together with two physical constraints on the conservation of total mass and conservation of total magnetization. Different treatments of the Zeeman energy terms are found to yield different numerical accuracies and stabilities. Numerical comparison between different numerical schemes is made, and the best scheme is identified. The numerical scheme is then applied to compute the condensate ground state in a harmonic plus optical lattice potential, and the effect of the periodic potential, in particular to the relative population of each hyperfine component, is investigated through comparison to the condensate ground state in a pure harmonic trap.
We study numerically the time-independent vector Gross--Pitaevskii equations (VGPEs) for ground states and time-dependent VGPEs with (or without) an external driven field for dynamics describing a multicomponent Bose--Einstein condensate (BEC) at zero or a very low temperature. In preparation for the numerics, we scale the three-dimensional (3d) VGPEs, approximately reduce it to lower dimensions, present a continuous normalized gradient flow (CNGF) to compute ground states of multicomponent BEC, prove energy diminishing of the CNGF, which provides a mathematical justification, and discretize it by the backward Euler finite difference (BEFD), which is monotone in linear and nonlinear cases and preserves energy diminishing property in the linear case. Then we use a time-splitting sine-spectral (TSSP) method to discretize the time-dependent VGPEs with an external driven field for computing dynamics of multicomponent BEC. The merits of the TSSP method for VGPEs are that it is explicit, unconditionally stable, time reversible and time transverse invariant if the VGPEs is, has "good" resolution in the semiclassical regime, is of spectral-order accuracy in space and second-order accuracy in time, and conserves the total particle number in the discretized level. Extensive numerical examples in three dimensions for ground states and dynamics of multicomponent BEC are presented to demonstrate the power of the numerical methods and to discuss the physics of multicomponent BEC.