Spectral method for the time-dependent Gross-Pitaevskii equation with a harmonic trap

We study the numerical resolution of the time-dependent Gross-Pitaevskii equation, a non-linear Schroedinger equation used to simulate the dynamics of Bose-Einstein condensates. Considering condensates trapped in harmonic potentials, we present an efficient algorithm by making use of a spectral Galerkin method, using a basis set of harmonic oscillator functions, and the Gauss-Hermite quadrature. We apply this algorithm to the simulation of condensate breathing and scissors modes.