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

TL;DR

This paper provides OpenMP Fortran 90/95 implementations of numerical programs for solving the dipolar Gross-Pitaevskii equation, enabling faster computations for Bose-Einstein condensates in various dimensions.

Contribution

It introduces parallelized OpenMP versions of existing codes for the dipolar GP equation, significantly reducing execution time on multicore processors.

Findings

Reduced execution time in multicore environments

Applicable to 1D, 2D, and 3D BEC simulations

Supports stationary and non-stationary solutions

Abstract

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.