FORTRESS II: FORTRAN programs for solving coupled Gross-Pitaevskii equations for spin-orbit coupled spin-2 Bose-Einstein condensate

TL;DR

This paper introduces three parallelized FORTRAN programs that efficiently compute ground states and dynamics of spin-2 Bose-Einstein condensates with spin-orbit coupling by solving coupled Gross-Pitaevskii equations using spectral methods.

Contribution

The paper provides a set of OpenMP parallelized FORTRAN programs for solving coupled Gross-Pitaevskii equations in spin-2 BECs with spin-orbit coupling, including validation and performance analysis.

Findings

Excellent agreement between full and simplified models in density profiles, energies, and chemical potentials.

High speedup and efficiency demonstrated for the parallelized codes.

Versatile implementation for 1D, 2D, and 3D simulations.

Abstract

We provide here a set of three OpenMP parallelized FORTRAN 90/95 programs to compute the ground states and the dynamics of trapped spin-2 Bose-Einstein condensates (BECs) with anisotropic spin-orbit (SO) coupling by solving a set of five coupled Gross-Pitaevskii equations using a time-splitting Fourier spectral method. Depending on the nature of the problem, without any loss of generality, we have employed the Cartesian grid spanning either three-, two-, or one-dimensional space for numerical discretization. To illustrate the veracity of the package, wherever feasible, we have compared the numerical ground state solutions of the full mean-field model with those from the simplified scalar models. The two set of results show excellent agreement, in particular, through the equilibrium density profiles, energies and chemical potentials of the ground-states. We have also presented test results for OpenMP performance parameters like speedup and the efficiency of the three codes.