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

Numerical Methods for Partial Differential Equations

A simple Dufort‐Frankel‐type scheme for the Gross‐Pitaevskii equation of Bose‐Einstein condensates on different geometries

Ming‐Chih Lai, Chung‐Yin Huang, Te‐Sheng Lin

Published Feb 26, 2004• 22 citations
Publisher

Abstract

Abstract We develop a simple Dufort‐Frankel‐type scheme for solving the time‐dependent Gross‐Pitaevskii equation (GPE). The GPE is a nonlinear Schrödinger equation describing the Bose‐Einstein condensation (BEC) at very low temperature. Three different geometries including 1D spherically symmetric, 2D cylindrically symmetric, and 3D anisotropic Cartesian domains are considered. The present finite difference method is explicit, linearly unconditional stable and is able to handle the coordinate singularities in a natural way. Furthermore, the scheme is time reversible and satisfies a discrete analogue of density conservation law. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004

DOI
10.1002/num.20008

Related papers

Analytical and numerical aspects in solving the controlled 3D Gross-Pitaevskii equation

Renato Fedele, Dušan Jovanović, Sergio M. De Nicola +4 more

The results of recently developed investigations, that have been carried out within the framework of the controlling potential method (CPM), are reviewed. This method allows one to decompose a three dimensional (3D) Gross‐Pitaevskii equation (GPE) into the pair of coupled Schrödinger‐type equations. Under suitable mathematical conditions, the solutions of the 3D controlled GPE can be constructed from the solutions of a 2D linear Schrödinger equation (the transverse component of the GPE) coupled with a 1D nonlinear Schrödinger equation (the longitudinal component of the GPE). Such decomposition allows one to cast the solutions in the form of the product of the solutions of the transverse and the longitudinal components of the GPE. The coupling between these two equations is the functional of both the transverse and the longitudinal profiles. It is shown that the CPM can be used to obtain a new class of three‐dimensional solitary waves solutions of the GPE, which governs the dynamics of Bose‐Einstein condensates. By imposing an external controlling potential, the desired time‐dependent shape of the localized BECs is obtained. The stability of the exact solutions was checked with direct simulations of the time ‐dependent, three‐dimensional GPE. Our simulations show that the localized condensates are stable with respect to perturbed initial conditions.

AIP conference proceedings•Jan 1, 2009•0 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

Mathematical theory and numerical methods for Bose-Einstein condensation

Weizhu Bao, Yongyong Cai

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.

arXiv:1212.5341•Kinetic and Related Models•Dec 11, 2012•591 citations

Numerical solution of the Gross-Pitaevskii equation using an explicit finite-difference scheme: An application to trapped Bose-Einstein condensates

Maria Mercede Cerimele, M. L. Chiofalo, Francesca Pistella +2 more

We present the application of a fast, explicit time-marching scheme for the solution of the Gross-Pitaevskii equation in cylindrical geometry. The scheme is validated on simple analytical tests and demonstrated for two situations of physical interest in experiments on the Bose-Einstein condensation (BEC) of trapped alkali-metal vapors. It is tested by reproducing known results on the free expansion of a BEC after removing a cylindrical trap, and it is then used to address the formation of matter-wave pulses that result from gravity-induced transport of a condensate in an optical potential.

Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jul 1, 2000•146 citations

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

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
Abstract We develop a simple Dufort‐Frankel‐type scheme for solving the time‐dependent Gross‐Pitaevskii equation (GPE). The GPE is a nonlinear Schrödinger equation describing the Bose‐Einstein condensation (BEC) at very low temperature. Three different geometries including 1D spherically symmetric, 2D cylindrically symmetric, and 3D anisotropic Cartesian domains are considered. The present finite difference method is explicit, linearly unconditional stable and is able to handle the coordinate singularities in a natural way. Furthermore, the scheme is time reversible and satisfies a discrete analogue of density conservation law. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004
The results of recently developed investigations, that have been carried out within the framework of the controlling potential method (CPM), are reviewed. This method allows one to decompose a three dimensional (3D) Gross‐Pitaevskii equation (GPE) into the pair of coupled Schrödinger‐type equations. Under suitable mathematical conditions, the solutions of the 3D controlled GPE can be constructed from the solutions of a 2D linear Schrödinger equation (the transverse component of the GPE) coupled with a 1D nonlinear Schrödinger equation (the longitudinal component of the GPE). Such decomposition allows one to cast the solutions in the form of the product of the solutions of the transverse and the longitudinal components of the GPE. The coupling between these two equations is the functional of both the transverse and the longitudinal profiles. It is shown that the CPM can be used to obtain a new class of three‐dimensional solitary waves solutions of the GPE, which governs the dynamics of Bose‐Einstein condensates. By imposing an external controlling potential, the desired time‐dependent shape of the localized BECs is obtained. The stability of the exact solutions was checked with direct simulations of the time ‐dependent, three‐dimensional GPE. Our simulations show that the localized condensates are stable with respect to perturbed initial conditions.
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.
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.
We present the application of a fast, explicit time-marching scheme for the solution of the Gross-Pitaevskii equation in cylindrical geometry. The scheme is validated on simple analytical tests and demonstrated for two situations of physical interest in experiments on the Bose-Einstein condensation (BEC) of trapped alkali-metal vapors. It is tested by reproducing known results on the free expansion of a BEC after removing a cylindrical trap, and it is then used to address the formation of matter-wave pulses that result from gravity-induced transport of a condensate in an optical potential.
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.
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.