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:cond-mat/0205172 · Journal of Physics B Atomic Molecular and Optical Physics

Bose-Einstein condensation dynamics from the numerical solution of the Gross-Pitaevskii equation

Sadhan K. Adhikari, Paulsamy Muruganandam

Published Jun 10, 2002• 137 citations
arXivHTMLPDF

Abstract

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.

DOI
10.1088/0953-4075/35/12/317

Related papers

Bose-Einstein Condensation Dynamics from the Numerical Solution of the Gross-Pitaevskii Equation by the Pseudospectral Method

Paulsamy Muruganandam, Sadhan K. Adhikari

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates of weakly interacting alkali atoms described by a nonlinear Gross-Pitaevskii (GP) equation. We suggest a pseudospectral method involving Laguerre polynomials to solve the time dependent GP for a spherically symmetric trap potential. The radial wavefunction and energy values have been calculated for different nonlinearities. Further, we study the effect of suddenly changing the interatomic scattering length or harmonic oscillator trap potential in the condensate. We also investigate the frequency of oscillation due to the variation in the strength of nonlinearity.

arXiv:cond-mat/0310403•Oct 17, 2003•0 citations

Numerical study of the coupled time-dependent Gross-Pitaevskii equation: Application to Bose-Einstein condensation

Sadhan K. Adhikari

We present a numerical study of the coupled time-dependent Gross-Pitaevskii equation, which describes the Bose-Einstein condensate of several types of trapped bosons at ultralow temperature with both attractive and repulsive interatomic interactions. The same approach is used to study both stationary and time-evolution problems. We consider up to four types of atoms in the study of stationary problems. We consider the time-evolution problems where the frequencies of the traps or the atomic scattering lengths are suddenly changed in a stable preformed condensate. We also study the effect of periodically varying these frequencies or scattering lengths on a preformed condensate. These changes introduce oscillations in the condensate, which are studied in detail. Good convergence is obtained in all cases studied.

arXiv:cond-mat/0103565•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Apr 20, 2001•26 citations

Numerical study of the spherically symmetric Gross-Pitaevskii equation in two space dimensions

Sadhan K. Adhikari

We present a numerical study of the time-dependent and time-independent Gross-Pitaevskii (GP) equation in two space dimensions, which describes the Bose-Einstein condensate of trapped bosons at ultralow temperature with both attractive and repulsive interatomic interactions. Both time-dependent and time-independent GP equations are used to study the stationary problems. In addition the time-dependent approach is used to study some evolution problems of the condensate. Specifically, we study the evolution problem where the trap energy is suddenly changed in a stable preformed condensate. In this case the system oscillates with increasing amplitude and does not remain limited between two stable configurations. Good convergence is obtained in all cases studied.

arXiv:cond-mat/0005450•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Aug 1, 2000•79 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

Stationary and Dynamical Solutions of the Gross-Pitaevskii Equation for a Bose-Einstein Condensate in a PT symmetric Double Well

Holger Cartarius, Dennis Dast, Daniel Haag +3 more

We investigate the Gross-Pitaevskii equation for a Bose-Einstein condensate in a PT symmetric double-well potential by means of the time-dependent variational principle and numerically exact solutions. A one-dimensional and a fully three-dimensional setup are used. Stationary states are determined and the propagation of wave function is investigated using the time-dependent Gross-Pitaevskii equation. Due to the nonlinearity of the Gross-Pitaevskii equation the potential dependson the wave function and its solutions decide whether or not the Hamiltonian itself is PT symmetric. Stationary solutions with real energy eigenvalues fulfilling exact PT symmetry are found as well as PT broken eigenstates with complex energies. The latter describe decaying or growing probability amplitudes and are not true stationary solutions of the time-dependent Gross-Pitaevskii equation. However, they still provide qualitative information about the time evolution of the wave functions.

Acta Polytechnica•Jan 3, 2013•7 citations

Bose–Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

Paulsamy Muruganandam, Sadhan K Adhikari

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross–Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose–Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x , y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge–Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank–Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

arXiv:cond-mat/0210177•Journal of Physics B Atomic Molecular and Optical Physics•Jun 5, 2003•123 citations
We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.
We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates of weakly interacting alkali atoms described by a nonlinear Gross-Pitaevskii (GP) equation. We suggest a pseudospectral method involving Laguerre polynomials to solve the time dependent GP for a spherically symmetric trap potential. The radial wavefunction and energy values have been calculated for different nonlinearities. Further, we study the effect of suddenly changing the interatomic scattering length or harmonic oscillator trap potential in the condensate. We also investigate the frequency of oscillation due to the variation in the strength of nonlinearity.
We present a numerical study of the coupled time-dependent Gross-Pitaevskii equation, which describes the Bose-Einstein condensate of several types of trapped bosons at ultralow temperature with both attractive and repulsive interatomic interactions. The same approach is used to study both stationary and time-evolution problems. We consider up to four types of atoms in the study of stationary problems. We consider the time-evolution problems where the frequencies of the traps or the atomic scattering lengths are suddenly changed in a stable preformed condensate. We also study the effect of periodically varying these frequencies or scattering lengths on a preformed condensate. These changes introduce oscillations in the condensate, which are studied in detail. Good convergence is obtained in all cases studied.
We present a numerical study of the time-dependent and time-independent Gross-Pitaevskii (GP) equation in two space dimensions, which describes the Bose-Einstein condensate of trapped bosons at ultralow temperature with both attractive and repulsive interatomic interactions. Both time-dependent and time-independent GP equations are used to study the stationary problems. In addition the time-dependent approach is used to study some evolution problems of the condensate. Specifically, we study the evolution problem where the trap energy is suddenly changed in a stable preformed condensate. In this case the system oscillates with increasing amplitude and does not remain limited between two stable configurations. Good convergence is obtained in all cases studied.
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.
We investigate the Gross-Pitaevskii equation for a Bose-Einstein condensate in a PT symmetric double-well potential by means of the time-dependent variational principle and numerically exact solutions. A one-dimensional and a fully three-dimensional setup are used. Stationary states are determined and the propagation of wave function is investigated using the time-dependent Gross-Pitaevskii equation. Due to the nonlinearity of the Gross-Pitaevskii equation the potential dependson the wave function and its solutions decide whether or not the Hamiltonian itself is PT symmetric. Stationary solutions with real energy eigenvalues fulfilling exact PT symmetry are found as well as PT broken eigenstates with complex energies. The latter describe decaying or growing probability amplitudes and are not true stationary solutions of the time-dependent Gross-Pitaevskii equation. However, they still provide qualitative information about the time evolution of the wave functions.
We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross–Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose–Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x , y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge–Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank–Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.