2019/03/18 by Ramavarmaraja Kishor Kumar, R. Kishor Kumar, Vladimir Loncar +5
Physics and Astronomy · #Bose–Einstein condensate #Classical mechanics #Cold Atom Physics and Bose-Einstein Condensates #Fortran #Function (biology) #Gross–Pitaevskii equation #Mechanics #Physics #Physics of Superconductivity and Magnetism #Quantum mechanics #Statistical physics #Strong Light-Matter Interactions #Vortex #Wave function #cond-mat.quant-gas #nlin.PS #physics.comp-ph #quant-ph
paper · pdf · doi:10.1016/j.cpc.2019.03.004
published as Comput. Phys. Commun. 240 (2019) 74 · 12 pages, 7 figures; to download the programs, click 'Other formats' and download the source
openalex publication_date 2019/03/18 · arxiv created 2019/06/14 · arxiv updated 2019/06/17 · openalex created_date 2025/10/10 · openalex updated_date 2026/08/06
We present OpenMP versions of C and Fortran programs for solving the Gross-Pitaevskii equation for a rotating trapped Bose-Einstein condensate (BEC) in two (2D) and three (3D) spatial dimensions. The programs can be used to generate vortex lattices and study dynamics of rotating BECs. We use the split-step Crank-Nicolson algorithm for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively. The programs propagate the condensate wave function and calculate several relevant physical quantities, such as the energy, the chemical potential, and the root-mean-square sizes. The imaginary-time propagation starts with an analytic wave function with one vortex at the trap center, modulated by a random phase at different space points. Nevertheless, the converged wave function for a rapidly rotating BEC with a large number of vortices is most efficiently calculated using the pre-calculated converged wave function of a rotating BEC containing a smaller number of vortices as the initial state rather than using an analytic wave function with one vortex as the initial state. These pre-calculated initial states exhibit rapid convergence for fast-rotating condensates to states containing multiple vortices with an appropriate phase structure. This is illustrated here by calculating vortex lattices with up to 61 vortices in 2D and 3D. Outputs of the programs include calculated physical quantities, as well as the wave function and different density profiles (full density, integrated densities in lower dimensions, and density cross-sections). The provided real-time propagation programs can be used to study the dynamics of a rotating BEC using the imaginary-time stationary wave function as the initial state. We also study the efficiency of parallelization of the present OpenMP C and Fortran programs with different compilers.