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.
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 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.
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.
[EN] In this work, we have studied a dilute and ultracold bosonic gas of weakly interacting atoms, \ni.e. a Bose-Einstein Condensate (BEC) by means of the Modi ed Gross-Pitaevskii equation \n(MGPE) taking into account beyond-mean- eld corrections due to quantum \nuctuations. We \nhave considered the case where the cloud of atoms is strongly con ned transversally (within a \ncigar-shaped trap) and therefore, the three-dimensional (3D) MGPE is reduced to an e ective \none-dimensional (1D) form by averaging over the transverse coordinates. Based on this approach, \nwe have performed numerical simulations for analyzing the behavior of the BEC under the \nin \nuence of quantum \nuctuations, which give rise to an additional nonlinear quartic term. After \nreviewing the mathematical theory of Bose-Einstein condensation based on the GPE and the \nnumerical methods used to integrate it we discuss the computation of the ground state and the \nsimulation of its dynamic evolution in a waveguide (free expansion).
The dynamics and interaction of quantized vortices in Bose-Einstein condensates (BECs) are investigated by using the two-dimensional Gross-Pitaevskii equation (GPE) with/without an angular momentum rotation term. If all vortices have the same winding number, they would rotate around the trap center but never collide. In contrast, if the winding numbers are different, their interaction highly depends on the initial distance between vortex centers. The analytical results are presented to describe the dynamics of the vortex centers when β=0. While if β≠0, there is no analytical result but some conclusive numerical findings are provided for the further understanding of vortex interaction in BECs. Finally, the dynamic laws describing the relation of vortex interaction in nonrotating and rotating BECs are presented.
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 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.