In this paper, we mainly review recent results on mathematical theory and numerical 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 then carried 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 for spin-1 BEC, Bogoliubov excitation and BEC at finite temperature.Comment: 135 pages and 10 figure

We study the numerical solution of the time-dependent Gross-Pitaevskii equation (GPE) describing a Bose-Einstein condensate (BEC) at zero or very low temperature. In preparation for the numerics we scale the 3d GrossPitaevskii equation and obtain a four-parameter model. Identifying 'extreme parameter regimes', the model is accessible to analytical perturbation theory, which justifies formal procedures well known in the physical literature: reduction to 2d and 1d GPEs, approximation of ground state solutions of the GPE and geometrical optics approximations. Then we use a time-splitting spectral method to discretize the time-dependent GPE. Again, perturbation theory is used to understand the discretization scheme and to choose the spatial/temporal grid in dependence of the perturbation parameter. Extensive numerical examples in 1d, 2d and 3d for weak/strong interactions, defocusing/focusing nonlinearity, and zero/nonzero initial phase data are presented to demonstrate the power of the numerical method and to discuss the physics of Bose-Einstein condensation. *

In this paper, we prove the energy diminishing of a normalized gradient flow which provides a mathematical justification of the imaginary time method used in physical literatures to compute the ground state solution of Bose-Einstein condensates (BEC). We also investigate the energy diminishing property for the discretization of the normalized gradient flow. Two numerical methods are proposed for such discretizations: one is the backward Euler centered finite difference (BEFD), the other one is an explicit time-splitting sine-spectral (TSSP) method. Energy diminishing for BEFD and TSSP for 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 energy diminishing property of the normalized gradient flow. Numerical results in 1d, 2d and 3d 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 normalized gradient flow can also be applied directly to compute the first excited state solution in BEC when the initial data is chosen as an odd function.

No abstract

In this paper, we begin with the nonlinear Schrödinger/Gross-Pitaevskii equation (NLSE/GPE) for modeling Bose-Einstein condensation (BEC) and nonlinear optics as well as other applications, and discuss their dynamical properties ranging from time reversible, time transverse invariant, mass and energy conservation, dispersion relation to soliton solutions. Then, we review and compare different numerical methods for solving the NLSE/GPE including finite difference time domain methods and time-splitting spectral method, and discuss different absorbing boundary conditions. In addition, these numerical methods are extended to the NLSE/GPE with damping terms and/or an angular momentum rotation term as well as coupled NLSEs/GPEs. Finally, applications to simulate a quantized vortex lattice dynamics in a rotating BEC are reported.