Spin-1 spin–orbit- and Rabi-coupled Bose–Einstein condensate solver

We present OpenMP versions of FORTRAN programs for solving the Gross–Pitaevskii equation for a harmonically trapped three-component spin-1 spinor Bose–Einstein condensate (BEC) in one (1D) and two (2D) spatial dimensions with or without spin–orbit (SO) and Rabi couplings. Several different forms of SO coupling are included in the programs. We use the split-step Crank–Nicolson discretization for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively. The imaginary-time propagation programs calculate the lowest-energy stationary state. The real-time propagation programs can be used to study the dynamics. The simulation input parameters are provided at the beginning of each program. The programs propagate the condensate wave function and calculate several relevant physical quantities. Outputs of the programs include the wave function, energy, root-mean-square sizes, different density profiles (linear density for the 1D program, linear and surface densities for the 2D program). The imaginary- or real-time propagation can start with an analytic wave function or a pre-calculated numerical wave function. The imaginary-time propagation usually starts with an analytic wave function, while the real-time propagation is often initiated with the previously calculated converged imaginary-time wave function.