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

Mathematical and Computational Applications

Comparison of Splitting Methods for Deterministic/Stochastic Gross–Pitaevskii Equation

Jürgen Geiser, Amirbahador Nasari

Published Aug 20, 2019• 7 citations
Publisher

Abstract

In this paper, we discuss the different splitting approaches to numerically solve the Gross–Pitaevskii equation (GPE). The models are motivated from spinor Bose–Einstein condensate (BEC). This system is formed of coupled mean-field equations, which are based on coupled Gross–Pitaevskii equations. We consider conservative finite-difference schemes and spectral methods for the spatial discretisation. Furthermore, we apply implicit or explicit time-integrators and combine these schemes with different splitting approaches. The numerical solutions are compared based on the conservation of the L 2 -norm with the analytical solutions. The advantages of the novel splitting methods for large time-domains are based on the asymptotic conservation of the solution of the soliton’s applications. Furthermore, we have the benefit of larger local time-steps and therefore obtain faster numerical schemes.

DOI
10.3390/mca24030076

Related papers

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

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

A time-splitting spectral method for computing dynamics of spinor F=1 Bose–Einstein condensates

Hanquan Wang

We propose a time-splitting spectral method for the generalized Gross–Pitaevskii equations, which describe the dynamics of spinor F=1 Bose–Einstein condensates at a very low temperature. The new numerical method is explicit, unconditionally stable, and of spectral accuracy in space. Moreover, it conserves the position densities at the discretized level. We apply the method for studying both the dynamic generation of vortices and vortex lattice dynamics for the spinor F=1 Bose–Einstein condensates held in an Ioffe-Pritchard magnetic field.

International Journal of Computer Mathematics•Jun 1, 2007•32 citations

Numerical solution to the time-dependent Gross-Pitaevskii equation

Tsogbayar Tsednee, B. Tsednee, Tsookhuu Khinayat

In this work we employ the split-step technique combined with a Legendre pseudospectral representation to solve various time-dependent GrossPitaevskii equations (GPE). Our findings based on the numerical accuracy of this approach applied for one-dimensional (1D) and two-dimensional (2D) problems show that it can provide accurate and stable solutions. Moreover, this approach has been applied to study the dynamics of the Bose-Einstein condensate which is modeled with the GPE. The breathing of condensate with a repulsive and attractive interactions trapped in 1D and 2D harmonic potentials has been simulated as well.

arXiv:2211.03964•Journal of Applied Science and Engineering A•Dec 31, 2022•4 citations

Computation of the dynamics of rotating 2D/3D Gross-Pitaevskii equations based on the HPC pseudo-spectral solver BEC2HPC

Jérémie Gaidamour, Xavier Antoine

The aim of this paper is to present the extension of the HPC pseudo-spectral solver BEC2HPC to compute the dynamics of 2D/3D rotating Bose-Einstein condensates modeled by the Gross-Pitaevskii equations with a rotation term. Numerical examples are provided to show the efficiency of the solver for large-scale simulations. NEW VERSION PROGRAM SUMMARY Program Title: BEC2HPC CPC Library link to program files: https://doi.org/10.17632/mdzpw4dr4t.2 Licensing provisions: GPLv2 Programming language: C++, Python Journal reference of previous version: Comput. Phys. Commun. 265 (2021) 108007, https://doi.org/10.1016/j.cpc.2021.108007 Does the new version supersede the previous version?: Yes Reasons for the new version: The first version of BEC2HPC provided a stationary-state solver, while this new version also enables the computation of the dynamics of rotating Gross–Pitaevskii equations. Summary of revisions: We have implemented ADI time-splitting and relaxation pseudo-spectral schemes in BEC2HPC, following the same parallel design and domain decomposition strategy as for the stationary-state solver. Nature of problem: This extension of the software BEC2HPC computes the dynamics of rotating Bose–Einstein condensates (BEC) modeled by the Gross–Pitaevskii equation (GPE). The available numerical methods are accurate and efficient for BEC with large rotation speeds and high nonlinearities. The parallel implementation allows performing large-scale simulations of 2D or 3D problems on parallel computing platforms. Solution method: Four time discretization methods are implemented for the dynamics of the GPE: a first-, second- and a fourth-order splitting schemes, and a second-order relaxation scheme. The spatial discretization is based on a Fourier pseudo-spectral method. This solver is implemented in distributed memory using MPI and a decomposition of the computational domain. Additional comments including restrictions and unusual features: The parallel algorithms are implemented in C++ and MPI, with a Python interface available for defining the physics of the problem. A pure Python version of sequential algorithms is also provided for rapid prototyping using the same interface. Results can be exported and visualized with external tools (HDF5 and VTK formats for ParaView or PyVista).

Computer Physics Communications•Jun 19, 2026•0 citations

A Fourth-Order Time-Splitting Laguerre--Hermite Pseudospectral Method for Bose--Einstein Condensates

Weizhu Bao, Jie Shen

A fourth-order time-splitting Laguerre--Hermite pseudospectral method is introduced for Bose--Einstein condensates (BECs) in three dimensions with cylindrical symmetry. The method is explicit, time reversible, and time transverse invariant. It conserves the position density and is spectral accurate in space and fourth-order accurate in time. Moreover, the new method has two other important advantages: (i) it reduces a three-dimensional (3-D) problem with cylindrical symmetry to an effective two-dimensional (2-D) problem; (ii) it solves the problem in the whole space instead of in a truncated artificial computational domain. The method is applied to vector Gross--Pitaevskii equations (VGPEs) for multicomponent BECs. Extensive numerical tests are presented for the one-dimensional (1-D) GPE, the 2-D GPE with radial symmetry, the 3-D GPE with cylindrical symmetry, as well as 3-D VGPEs for two-component BECs, to show the efficiency and accuracy of the new numerical method.

arXiv:cond-mat/0310378•SIAM Journal on Scientific Computing•Jan 1, 2005•169 citations
In this paper, we discuss the different splitting approaches to numerically solve the Gross–Pitaevskii equation (GPE). The models are motivated from spinor Bose–Einstein condensate (BEC). This system is formed of coupled mean-field equations, which are based on coupled Gross–Pitaevskii equations. We consider conservative finite-difference schemes and spectral methods for the spatial discretisation. Furthermore, we apply implicit or explicit time-integrators and combine these schemes with different splitting approaches. The numerical solutions are compared based on the conservation of the L 2 -norm with the analytical solutions. The advantages of the novel splitting methods for large time-domains are based on the asymptotic conservation of the solution of the soliton’s applications. Furthermore, we have the benefit of larger local time-steps and therefore obtain faster numerical schemes.
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 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 propose a time-splitting spectral method for the generalized Gross–Pitaevskii equations, which describe the dynamics of spinor F=1 Bose–Einstein condensates at a very low temperature. The new numerical method is explicit, unconditionally stable, and of spectral accuracy in space. Moreover, it conserves the position densities at the discretized level. We apply the method for studying both the dynamic generation of vortices and vortex lattice dynamics for the spinor F=1 Bose–Einstein condensates held in an Ioffe-Pritchard magnetic field.
In this work we employ the split-step technique combined with a Legendre pseudospectral representation to solve various time-dependent GrossPitaevskii equations (GPE). Our findings based on the numerical accuracy of this approach applied for one-dimensional (1D) and two-dimensional (2D) problems show that it can provide accurate and stable solutions. Moreover, this approach has been applied to study the dynamics of the Bose-Einstein condensate which is modeled with the GPE. The breathing of condensate with a repulsive and attractive interactions trapped in 1D and 2D harmonic potentials has been simulated as well.
The aim of this paper is to present the extension of the HPC pseudo-spectral solver BEC2HPC to compute the dynamics of 2D/3D rotating Bose-Einstein condensates modeled by the Gross-Pitaevskii equations with a rotation term. Numerical examples are provided to show the efficiency of the solver for large-scale simulations. NEW VERSION PROGRAM SUMMARY Program Title: BEC2HPC CPC Library link to program files: https://doi.org/10.17632/mdzpw4dr4t.2 Licensing provisions: GPLv2 Programming language: C++, Python Journal reference of previous version: Comput. Phys. Commun. 265 (2021) 108007, https://doi.org/10.1016/j.cpc.2021.108007 Does the new version supersede the previous version?: Yes Reasons for the new version: The first version of BEC2HPC provided a stationary-state solver, while this new version also enables the computation of the dynamics of rotating Gross–Pitaevskii equations. Summary of revisions: We have implemented ADI time-splitting and relaxation pseudo-spectral schemes in BEC2HPC, following the same parallel design and domain decomposition strategy as for the stationary-state solver. Nature of problem: This extension of the software BEC2HPC computes the dynamics of rotating Bose–Einstein condensates (BEC) modeled by the Gross–Pitaevskii equation (GPE). The available numerical methods are accurate and efficient for BEC with large rotation speeds and high nonlinearities. The parallel implementation allows performing large-scale simulations of 2D or 3D problems on parallel computing platforms. Solution method: Four time discretization methods are implemented for the dynamics of the GPE: a first-, second- and a fourth-order splitting schemes, and a second-order relaxation scheme. The spatial discretization is based on a Fourier pseudo-spectral method. This solver is implemented in distributed memory using MPI and a decomposition of the computational domain. Additional comments including restrictions and unusual features: The parallel algorithms are implemented in C++ and MPI, with a Python interface available for defining the physics of the problem. A pure Python version of sequential algorithms is also provided for rapid prototyping using the same interface. Results can be exported and visualized with external tools (HDF5 and VTK formats for ParaView or PyVista).
A fourth-order time-splitting Laguerre--Hermite pseudospectral method is introduced for Bose--Einstein condensates (BECs) in three dimensions with cylindrical symmetry. The method is explicit, time reversible, and time transverse invariant. It conserves the position density and is spectral accurate in space and fourth-order accurate in time. Moreover, the new method has two other important advantages: (i) it reduces a three-dimensional (3-D) problem with cylindrical symmetry to an effective two-dimensional (2-D) problem; (ii) it solves the problem in the whole space instead of in a truncated artificial computational domain. The method is applied to vector Gross--Pitaevskii equations (VGPEs) for multicomponent BECs. Extensive numerical tests are presented for the one-dimensional (1-D) GPE, the 2-D GPE with radial symmetry, the 3-D GPE with cylindrical symmetry, as well as 3-D VGPEs for two-component BECs, to show the efficiency and accuracy of the new numerical method.