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:1212.5341 · Kinetic and Related Models

Mathematical theory and numerical methods for Bose-Einstein condensation

Weizhu Bao, Yongyong Cai

Published Dec 11, 2012• 591 citations
arXivHTMLPDF

Abstract

In this paper, we mainly review recent results on mathematical theory andnumerical methods for Bose-Einstein condensation (BEC),based on the Gross-Pitaevskii equation (GPE). Starting from the simplest case with one-component BEC of the weakly interacting bosons, we study the reduction of GPE to lower dimensions, the ground states of BEC including the existence and uniqueness as well as nonexistence results, and the dynamics of GPE including dynamical laws, well-posedness of the Cauchy problem as well as the finite time blow-up. To compute the ground state, the gradient flow with discrete normalization (or imaginary time) method is reviewed and various full discretization methods are presented and compared. To simulate the dynamics, both finite difference methods and time splitting spectral methods are reviewed, and their error estimates are briefly outlined. When the GPE has symmetric properties, we show how to simplify the numerical methods. Then we compare two widely used scalings, i.e. physical scaling (commonly used) and semiclassical scaling, for BEC in strong repulsive interaction regime (Thomas-Fermi regime), and discuss semiclassical limits of the GPE. Extensions of these results for one-component BEC are thencarried out for rotating BEC by GPE with an angular momentum rotation, dipolar BEC by GPE with long range dipole-dipole interaction, and two-component BEC by coupled GPEs. Finally, as a perspective, we show briefly the mathematical models forspin-1 BEC, Bogoliubov excitation and BEC at finite temperature.

DOI
10.3934/krm.2013.6.1

Related papers

Mathematical Models and Numerical Methods for Bose-Einstein Condensation

Weizhu Bao

The achievement of Bose-Einstein condensation (BEC) in ultracold vapors of alkali atoms has given enormous impulse to the theoretical and experimental study of dilute atomic gases in condensed quantum states inside magnetic traps and optical lattices. This article offers a short survey on mathematical models and theories as well as numerical methods for BEC based on the mean field theory. We start with the Gross-Pitaevskii equation (GPE) in three dimensions (3D) for modeling one-component BEC of the weakly interacting bosons, scale it to obtain a three-parameter model and show how to reduce it to two dimensions (2D) and one dimension (1D) GPEs in certain limiting regimes. Mathematical theories and numerical methods for ground states and dynamics of BEC are provided. Extensions to GPE with an angular momentum rotation term for a rotating BEC, to GPE with long-range anisotropic dipole-dipole interaction for a dipolar BEC and to coupled GPEs for spin-orbit coupled BECs are discussed. Finally, some conclusions are drawn and future research perspectives are discussed.

arXiv:1403.3884•Mar 16, 2014•10 citations

Mathematical Models and Numerical Methods for Spinor Bose-Einstein Condensates.

Weizhu Bao null, Yongyong Cai

In this paper, we systematically review mathematical models, theories and numerical methods for ground states and dynamics of spinor Bose-Einstein condensates (BECs) based on the coupled Gross-Pitaevskii equations (GPEs). We start with a pseudo spin-1/2 BEC system with/without an internal atomic Josephson junction and spin-orbit coupling including (i) existence and uniqueness as well as non-existence of ground states under different parameter regimes, (ii) ground state structures under different limiting parameter regimes, (iii) dynamical properties, and (iv) efficient and accurate numerical methods for computing ground states and dynamics. Then we extend these results to spin-1 BEC and spin-2 BEC. Finally, extensions to dipolar spinor systems and/or general spin-F (F>=3) BEC are discussed.

arXiv:1709.03840•Communications in Computational Physics•Jun 1, 2018•57 citations

Gross–Pitaevskii description of superfluid dynamics at finite temperature: A short review of recent results

Marc Brächet

The Gross–Pitaevskii equation (GPE) describes the dynamics of superflows and Bose–Einstein Condensates (BEC) at very low temperature. When a truncation of Fourier modes is performed, the resulting truncated GPE (TGPE) can also describe the correct thermal behavior of a Bose gas, as long as all relevant modes are highly occupied [M.J. Davis, S.A. Morgan, K. Burnett, Simulations of Bose fields at finite temperature, Phys. Rev. Lett. 87 (16) (2001) 160402]. We review some of our groupʼs recent GPE- and TGPE-based numerical studies of superfluid dynamics and BEC stability. The relations with experiments are discussed.

Comptes Rendus Physique•Nov 1, 2012•12 citations

A splitting compact finite difference method for computing the dynamics of dipolar Bose–Einstein condensate

Hanquan Wang

We numerically study the nonlocal Gross–Pitaevskii equation (NGPE) which describes the dynamics of Bose–Einstein condensates (BEC) with dipole–dipole interaction at extremely low temperature. In preparation for the numerics, first we reformulate the dimensionless NGPE into a Schrödinger–Poisson system. Then, we discretize the three-dimensional Schrödinger–Poisson system in space by a sixth-order compact finite difference method and in time by a splitting technique. By means of three-dimensional discrete fast Sine transform, we develop a fast solver for the resulting discretized system. Finally, we present numerical examples in three dimensions to demonstrate the power of the numerical methods and to discuss some physics of dipolar BEC. The merits of the proposed method for the NGPE are that it is fast and unconditionally stable. Moreover, the method is of spectral-like accuracy in space, and conserves the particle number and the energy of the system in the discretized level.

International Journal of Computer Mathematics•Dec 28, 2016•2 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

Gross–Pitaevskii–Poisson Equations for Dipolar Bose–Einstein Condensate with Anisotropic Confinement

Weizhu Bao, Naoufel Ben Abdallah, Yongyong Cai

Ground states and dynamical properties of a dipolar Bose–Einstein condensate are analyzed based on the Gross–Pitaevskii–Poisson system (GPPS) and its dimension reduction models under an anisotropic confining potential. We begin with the three-dimensional (3D) GPPS and review its quasi-two-dimensional (2D) approximate equations when the trap is strongly confined in the z-direction and quasi-one-dimensional (1D) approximate equations when the trap is strongly confined in the x- and y-directions. In fact, in the quasi-2D equations, a fractional Poisson equation with the operator $(-\Delta)^{1/2}$ is involved which brings significant difficulties into the analysis. Existence and uniqueness as well as nonexistence of the ground state under different parameter regimes are established for the quasi-2D and quasi-1D equations. Well-posedness of the Cauchy problem for both types of equations and finite time blow-up in two dimensions are analyzed. Finally, we rigorously prove the convergence with linear convergence rate for the solutions of the 3D GPPS and its quasi-2D and quasi-1D approximate equations in the weak interaction regime.

arXiv:1212.4945•SIAM Journal on Mathematical Analysis•Jan 1, 2012•34 citations
In this paper, we mainly review recent results on mathematical theory andnumerical methods for Bose-Einstein condensation (BEC),based on the Gross-Pitaevskii equation (GPE). Starting from the simplest case with one-component BEC of the weakly interacting bosons, we study the reduction of GPE to lower dimensions, the ground states of BEC including the existence and uniqueness as well as nonexistence results, and the dynamics of GPE including dynamical laws, well-posedness of the Cauchy problem as well as the finite time blow-up. To compute the ground state, the gradient flow with discrete normalization (or imaginary time) method is reviewed and various full discretization methods are presented and compared. To simulate the dynamics, both finite difference methods and time splitting spectral methods are reviewed, and their error estimates are briefly outlined. When the GPE has symmetric properties, we show how to simplify the numerical methods. Then we compare two widely used scalings, i.e. physical scaling (commonly used) and semiclassical scaling, for BEC in strong repulsive interaction regime (Thomas-Fermi regime), and discuss semiclassical limits of the GPE. Extensions of these results for one-component BEC are thencarried out for rotating BEC by GPE with an angular momentum rotation, dipolar BEC by GPE with long range dipole-dipole interaction, and two-component BEC by coupled GPEs. Finally, as a perspective, we show briefly the mathematical models forspin-1 BEC, Bogoliubov excitation and BEC at finite temperature.
The achievement of Bose-Einstein condensation (BEC) in ultracold vapors of alkali atoms has given enormous impulse to the theoretical and experimental study of dilute atomic gases in condensed quantum states inside magnetic traps and optical lattices. This article offers a short survey on mathematical models and theories as well as numerical methods for BEC based on the mean field theory. We start with the Gross-Pitaevskii equation (GPE) in three dimensions (3D) for modeling one-component BEC of the weakly interacting bosons, scale it to obtain a three-parameter model and show how to reduce it to two dimensions (2D) and one dimension (1D) GPEs in certain limiting regimes. Mathematical theories and numerical methods for ground states and dynamics of BEC are provided. Extensions to GPE with an angular momentum rotation term for a rotating BEC, to GPE with long-range anisotropic dipole-dipole interaction for a dipolar BEC and to coupled GPEs for spin-orbit coupled BECs are discussed. Finally, some conclusions are drawn and future research perspectives are discussed.
In this paper, we systematically review mathematical models, theories and numerical methods for ground states and dynamics of spinor Bose-Einstein condensates (BECs) based on the coupled Gross-Pitaevskii equations (GPEs). We start with a pseudo spin-1/2 BEC system with/without an internal atomic Josephson junction and spin-orbit coupling including (i) existence and uniqueness as well as non-existence of ground states under different parameter regimes, (ii) ground state structures under different limiting parameter regimes, (iii) dynamical properties, and (iv) efficient and accurate numerical methods for computing ground states and dynamics. Then we extend these results to spin-1 BEC and spin-2 BEC. Finally, extensions to dipolar spinor systems and/or general spin-F (F>=3) BEC are discussed.
The Gross–Pitaevskii equation (GPE) describes the dynamics of superflows and Bose–Einstein Condensates (BEC) at very low temperature. When a truncation of Fourier modes is performed, the resulting truncated GPE (TGPE) can also describe the correct thermal behavior of a Bose gas, as long as all relevant modes are highly occupied [M.J. Davis, S.A. Morgan, K. Burnett, Simulations of Bose fields at finite temperature, Phys. Rev. Lett. 87 (16) (2001) 160402]. We review some of our groupʼs recent GPE- and TGPE-based numerical studies of superfluid dynamics and BEC stability. The relations with experiments are discussed.
We numerically study the nonlocal Gross–Pitaevskii equation (NGPE) which describes the dynamics of Bose–Einstein condensates (BEC) with dipole–dipole interaction at extremely low temperature. In preparation for the numerics, first we reformulate the dimensionless NGPE into a Schrödinger–Poisson system. Then, we discretize the three-dimensional Schrödinger–Poisson system in space by a sixth-order compact finite difference method and in time by a splitting technique. By means of three-dimensional discrete fast Sine transform, we develop a fast solver for the resulting discretized system. Finally, we present numerical examples in three dimensions to demonstrate the power of the numerical methods and to discuss some physics of dipolar BEC. The merits of the proposed method for the NGPE are that it is fast and unconditionally stable. Moreover, the method is of spectral-like accuracy in space, and conserves the particle number and the energy of the system in the discretized level.
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.
Ground states and dynamical properties of a dipolar Bose–Einstein condensate are analyzed based on the Gross–Pitaevskii–Poisson system (GPPS) and its dimension reduction models under an anisotropic confining potential. We begin with the three-dimensional (3D) GPPS and review its quasi-two-dimensional (2D) approximate equations when the trap is strongly confined in the z-direction and quasi-one-dimensional (1D) approximate equations when the trap is strongly confined in the x- and y-directions. In fact, in the quasi-2D equations, a fractional Poisson equation with the operator (−Δ)1/2(-\Delta)^{1/2}(−Δ)1/2 is involved which brings significant difficulties into the analysis. Existence and uniqueness as well as nonexistence of the ground state under different parameter regimes are established for the quasi-2D and quasi-1D equations. Well-posedness of the Cauchy problem for both types of equations and finite time blow-up in two dimensions are analyzed. Finally, we rigorously prove the convergence with linear convergence rate for the solutions of the 3D GPPS and its quasi-2D and quasi-1D approximate equations in the weak interaction regime.