OpenMP Fortran programs for solving the time-dependent dipolar Gross-Pitaevskii equation

• DOI: 10.1016/j.cpc.2023.108669
• Publication Date: 2023-01-25T08:42:37Z
• Authors: Luis E. Young-S.; Paulsamy Murugana...; Antun Balaž; Sadhan K. Adhikari
• Keywords: Fortran; Gross–Pitaevskii equation; Code (set theory)

In this paper we present Open Multi-Processing (OpenMP) Fortran 90/95 versions of previously published numerical programs for solving the dipolar Gross-Pitaevskii (GP) equation including the contact interaction in one, two and three spatial dimensions. The atoms are considered to be polarized along the z axis and we consider different cases, e.g., stationary and non-stationary solutions of the GP equation for a dipolar Bose-Einstein condensate (BEC) in one dimension (along x and z axes), two dimensions (in x-y and x-z planes), and three dimensions. The algorithm used is the split-step semi-implicit Crank-Nicolson scheme for imaginary- and real-time propagation to obtain stationary states and BEC dynamics, respectively, as in the previous version [R. Kishor Kumar et al., Comput. Phys. Commun. 195, 117 (2015)]. These OpenMP versions have significantly reduced execution time in multicore processors.