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:0904.3440 · Physical Review E

Numerical method for evolving the dipolar projected Gross-Pitaevskii equation

P. B. Blakie, Christopher Ticknor, Ashton S. Bradley, A. Martin, Matthew J. Davis, Yuki Kawaguchi

Published Jul 15, 2009• 23 citations
arXivHTMLPDF

Abstract

We describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for an interacting Bose gas in a harmonic-oscillator potential, with the inclusion of a long-range dipolar interaction. The central difficulty in solving this equation is the requirement that the field is restricted to a small set of prescribed modes that constitute the low-energy c -field region of the system. We present a scheme, using a Hermite-polynomial-based spectral representation, which precisely implements this mode restriction and allows an efficient and accurate solution of the dipolar PGPE. We introduce a set of auxiliary oscillator states to perform a Fourier transform necessary to evaluate the dipolar interaction in reciprocal space. We extensively characterize the accuracy of our approach and derive Ehrenfest equations for the evolution of the angular momentum.

DOI
10.1103/physreve.80.016703

Related papers

Numerical method for evolving the projected Gross-Pitaevskii equation

P. B. Blakie

In this paper we describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for a Bose gas in a harmonic oscillator potential. The central difficulty in solving this equation is the requirement that the classical field is restricted to a small set of prescribed modes that constitute the low energy classical region of the system. We present a scheme, using a Hermite-polynomial based spectral representation, that precisely implements this mode restriction and allows an efficient and accurate solution of the PGPE. We show equilibrium and nonequilibrium results from the application of the PGPE to an anisotropic trapped three-dimensional Bose gas.

arXiv:0803.3664•Physical Review E•Aug 8, 2008•39 citations

Numerical method for the stochastic projected Gross-Pitaevskii equation

S. J. Rooney, P. B. Blakie, Ashton S. Bradley

We present a method for solving the stochastic projected Gross-Pitaevskii equation (SPGPE) for a three-dimensional weakly interacting Bose gas in a harmonic-oscillator trapping potential. The SPGPE contains the challenge of both accurately evolving all modes in the low-energy classical region of the system, and evaluating terms from the number-conserving scattering reservoir process. We give an accurate and efficient procedure for evaluating the scattering terms using a Hermite-polynomial based spectral-Galerkin representation, which allows us to precisely implement the low-energy mode restriction. Stochastic integration is performed using the weak semi-implicit Euler method. We extensively characterize the accuracy of our method, finding a faster-than-expected rate of stochastic convergence. Physical consistency of the algorithm is demonstrated by considering thermalization of initially random states.

arXiv:1310.0161•Physical Review E•Jan 10, 2014•29 citations

A simple and efficient numerical method for computing the dynamics of rotating Bose-Einstein condensates via a rotating Lagrangian coordinate

Weizhu Bao, Daniel Marahrens, Qinglin Tang +1 more

We propose a simple, efficient and accurate numerical method for simulating the dynamics of rotating Bose-Einstein condensates (BECs) in a rotational frame with/without a long-range dipole-dipole interaction. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range dipole-dipole interaction, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential and review some dynamical laws related to the 2D and 3D GPE. By introducing a rotating Lagrangian coordinate system, the original GPEs are re-formulated to GPEs without the angular momentum rotation which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable and very efficient in computation. It is spectral order accurate in space and second-order accurate in time, and conserves the mass in the discrete level. Extensive numerical results are reported to demonstrate the efficiency and accuracy of the new numerical method. Finally, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation and center-of-mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without/with the long-range dipole-dipole interaction.

arXiv:1305.1378•May 7, 2013•1 citation

A Simple and Efficient Numerical Method for Computing the Dynamics of Rotating Bose--Einstein Condensates via Rotating Lagrangian Coordinates

Weizhu Bao, Daniel Marahrens, Qinglin Tang +1 more

We propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose--Einstein condensates (BECs) in a rotational frame with or without long-range dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.

SIAM Journal on Scientific Computing•Jan 1, 2013•42 citations

Projected Gross–Pitaevskii Equation for Ring-Shaped Bose–Einstein Condensates

O. O. Prikhodko, Y. M. Bidasyuk

We propose an alternative implementation of the projected Gross–Pitaevskki equation adapted for a numerical modeling of the atomic Bose–Einstein condensate trapped in a toroidally shaped potential. We present an accurate efficient scheme to evaluate the required matrix elements and to calculate tthe ime evolution of the matter wave field. We analyze the stability and accuracy of the developed method for equilibrium and nonequilibrium solutions in a ring-shaped trap with an additional barrier potential corresponding to recent experimental realizations.

arXiv:2003.06921•Ukrainian Journal of Physics•Apr 7, 2021•3 citations

Accurate and Efficient Numerical Methods for Computing Ground States and Dynamics of Dipolar Bose-Einstein Condensates via the Nonuniform FFT

Weizhu Bao, Qinglin Tang, Yong Zhang

We propose efficient and accurate numerical methods for computing the ground state and dynamics of the dipolar Bose-Einstein condensates utilising a newly developed dipole-dipole interaction (DDI) solver that is implemented with the non-uniform fast Fourier transform (NUFFT) algorithm. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with a DDI term and present the corresponding two-dimensional (2D) model under a strongly anisotropic confining potential. Different from existing methods, the NUFFT based DDI solver removes the singularity by adopting the spherical/polar coordinates in the Fourier space in 3D/2D, respectively, thus it can achieve spectral accuracy in space and simultaneously maintain high efficiency by making full use of FFT and NUFFT whenever it is necessary and/or needed. Then, we incorporate this solver into existing successful methods for computing the ground state and dynamics of GPE with a DDI for dipolar BEC. Extensive numerical comparisons with existing methods are carried out for computing the DDI, ground states and dynamics of the dipolar BEC. Numerical results show that our new methods outperform existing methods in terms of both accuracy and efficiency.

arXiv:1504.02897•Communications in Computational Physics•Jan 1, 2016•32 citations
We describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for an interacting Bose gas in a harmonic-oscillator potential, with the inclusion of a long-range dipolar interaction. The central difficulty in solving this equation is the requirement that the field is restricted to a small set of prescribed modes that constitute the low-energy c -field region of the system. We present a scheme, using a Hermite-polynomial-based spectral representation, which precisely implements this mode restriction and allows an efficient and accurate solution of the dipolar PGPE. We introduce a set of auxiliary oscillator states to perform a Fourier transform necessary to evaluate the dipolar interaction in reciprocal space. We extensively characterize the accuracy of our approach and derive Ehrenfest equations for the evolution of the angular momentum.
In this paper we describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for a Bose gas in a harmonic oscillator potential. The central difficulty in solving this equation is the requirement that the classical field is restricted to a small set of prescribed modes that constitute the low energy classical region of the system. We present a scheme, using a Hermite-polynomial based spectral representation, that precisely implements this mode restriction and allows an efficient and accurate solution of the PGPE. We show equilibrium and nonequilibrium results from the application of the PGPE to an anisotropic trapped three-dimensional Bose gas.
We present a method for solving the stochastic projected Gross-Pitaevskii equation (SPGPE) for a three-dimensional weakly interacting Bose gas in a harmonic-oscillator trapping potential. The SPGPE contains the challenge of both accurately evolving all modes in the low-energy classical region of the system, and evaluating terms from the number-conserving scattering reservoir process. We give an accurate and efficient procedure for evaluating the scattering terms using a Hermite-polynomial based spectral-Galerkin representation, which allows us to precisely implement the low-energy mode restriction. Stochastic integration is performed using the weak semi-implicit Euler method. We extensively characterize the accuracy of our method, finding a faster-than-expected rate of stochastic convergence. Physical consistency of the algorithm is demonstrated by considering thermalization of initially random states.
We propose a simple, efficient and accurate numerical method for simulating the dynamics of rotating Bose-Einstein condensates (BECs) in a rotational frame with/without a long-range dipole-dipole interaction. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range dipole-dipole interaction, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential and review some dynamical laws related to the 2D and 3D GPE. By introducing a rotating Lagrangian coordinate system, the original GPEs are re-formulated to GPEs without the angular momentum rotation which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable and very efficient in computation. It is spectral order accurate in space and second-order accurate in time, and conserves the mass in the discrete level. Extensive numerical results are reported to demonstrate the efficiency and accuracy of the new numerical method. Finally, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation and center-of-mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without/with the long-range dipole-dipole interaction.
We propose a simple, efficient, and accurate numerical method for simulating the dynamics of rotating Bose--Einstein condensates (BECs) in a rotational frame with or without long-range dipole-dipole interaction (DDI). We begin with the three-dimensional (3D) Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and/or long-range DDI, state the two-dimensional (2D) GPE obtained from the 3D GPE via dimension reduction under anisotropic external potential, and review some dynamical laws related to the 2D and 3D GPEs. By introducing a rotating Lagrangian coordinate system, the original GPEs are reformulated to GPEs without the angular momentum rotation, which is replaced by a time-dependent potential in the new coordinate system. We then cast the conserved quantities and dynamical laws in the new rotating Lagrangian coordinates. Based on the new formulation of the GPE for rotating BECs in the rotating Lagrangian coordinates, a time-splitting spectral method is presented for computing the dynamics of rotating BECs. The new numerical method is explicit, simple to implement, unconditionally stable, and very efficient in computation. It is spectral-order accurate in space and second-order accurate in time and conserves the mass on the discrete level. We compare our method with some representative methods in the literature to demonstrate its efficiency and accuracy. In addition, the numerical method is applied to test the dynamical laws of rotating BECs such as the dynamics of condensate width, angular momentum expectation, and center of mass, and to investigate numerically the dynamics and interaction of quantized vortex lattices in rotating BECs without or with the long-range DDI.
We propose an alternative implementation of the projected Gross–Pitaevskki equation adapted for a numerical modeling of the atomic Bose–Einstein condensate trapped in a toroidally shaped potential. We present an accurate efficient scheme to evaluate the required matrix elements and to calculate tthe ime evolution of the matter wave field. We analyze the stability and accuracy of the developed method for equilibrium and nonequilibrium solutions in a ring-shaped trap with an additional barrier potential corresponding to recent experimental realizations.
We propose efficient and accurate numerical methods for computing the ground state and dynamics of the dipolar Bose-Einstein condensates utilising a newly developed dipole-dipole interaction (DDI) solver that is implemented with the non-uniform fast Fourier transform (NUFFT) algorithm. We begin with the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with a DDI term and present the corresponding two-dimensional (2D) model under a strongly anisotropic confining potential. Different from existing methods, the NUFFT based DDI solver removes the singularity by adopting the spherical/polar coordinates in the Fourier space in 3D/2D, respectively, thus it can achieve spectral accuracy in space and simultaneously maintain high efficiency by making full use of FFT and NUFFT whenever it is necessary and/or needed. Then, we incorporate this solver into existing successful methods for computing the ground state and dynamics of GPE with a DDI for dipolar BEC. Extensive numerical comparisons with existing methods are carried out for computing the DDI, ground states and dynamics of the dipolar BEC. Numerical results show that our new methods outperform existing methods in terms of both accuracy and efficiency.