Numerical simulation of the Gross-Pitaevskii equation via vortex tracking

Abstract

This paper deals with the numerical simulation of the Gross-Pitaevskii (GP) equation, for which a well-known feature is the appearance of quantized vortices with core size of the order of a small parameter $\epsilon$. Without a magnetic field and with suitable initial conditions, these vortices interact, in the singular limit $\epsilon \to 0$, through an explicit Hamiltonian dynamics. Using this analytical framework, we develop and analyze a numerical strategy based on the reduced-order Hamiltonian system to efficiently simulate the infinite-dimensional GP equation for small, but finite, $\epsilon$. This method allows us to avoid numerical stability issues in solving the GP equation, where small values of $\epsilon$ typically require very fine meshes and time steps. We also provide a mathematical justification of our method in terms of rigorous error estimates of the error in the supercurrent, together with numerical illustrations.

References in corpus (3)

• Emergence of order from turbulence in an isolated planar superfluid
• Dynamics of nearly parallel vortex filaments for the Gross-Pitaevskii equation
• A numerical study of vortex nucleation in 2D rotating Bose-Einstein condensates

Cited by in corpus (1)

• Simulating Quantum Turbulence with Matrix Product States