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

Abstract 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. Program summary Program title: BEC-GP-SPINOR, consisting of: BEC-GP-SPINOR-OMP package, containing programs spin-SO-imre1d-omp.f90 and spin-SO-imre2d-omp.f90, with util.f90. CPC Library link to program files: https://doi.org/10.17632/j3wr4wn946.1 Licensing provisions: Apache License 2.0 Programming language: OpenMP FORTRAN. The FORTRAN programs are tested with the GNU, Intel, PGI, and Oracle compiler. Nature of problem: The present Open Multi-Processing (OpenMP) FORTRAN programs solve the time-dependent nonlinear partial differential Gross–Pitaevskii (GP) equation for a trapped spinor Bose–Einstein condensate, with or without spin–orbit coupling, in one and two spatial dimensions. Solution method: We employ the split-step Crank–Nicolson rule to discretize the time-dependent GP equation in space and time. The discretized equation is then solved by imaginary- or real-time propagation, employing adequately small space and time steps, to yield the solution of stationary and non-stationary problems, respectively.

[1]  Vladimir Loncar,et al.  OpenMP, OpenMP/MPI, and CUDA/MPI C programs for solving the time-dependent dipolar Gross-Pitaevskii equation , 2016, Comput. Phys. Commun..

[2]  S. K. Adhikari,et al.  Analytic models for the density of a ground-state spinor condensate , 2015, 1507.06883.

[3]  Luis E. Young-S.,et al.  Fortran and C programs for the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap , 2015, Comput. Phys. Commun..

[4]  L. Salasnich,et al.  Effective wave equations for the dynamics of cigar-shaped and disk-shaped Bose condensates , 2002 .

[5]  Victor Galitski,et al.  Spin–orbit coupling in quantum gases , 2013, Nature.

[6]  W. Ketterle,et al.  OPTICAL CONFINEMENT OF A BOSE-EINSTEIN CONDENSATE , 1998 .

[7]  Paulsamy Muruganandam,et al.  C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap , 2012, Comput. Phys. Commun..

[8]  Vladimir Loncar,et al.  C and Fortran OpenMP programs for rotating Bose-Einstein condensates , 2019, Comput. Phys. Commun..

[9]  Weizhu Bao,et al.  Computing Ground States of Spin-1 Bose-Einstein Condensates by the Normalized Gradient Flow , 2007, SIAM J. Sci. Comput..

[10]  W. Ketterle,et al.  Spin domains in ground-state Bose–Einstein condensates , 1998, Nature.

[11]  I. B. Spielman,et al.  Spin–orbit-coupled Bose–Einstein condensates , 2011, Nature.

[12]  Vladimir Loncar,et al.  CUDA programs for solving the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap , 2016, Comput. Phys. Commun..

[13]  Weizhu Bao,et al.  Numerical methods for computing the ground state of spin-1 Bose-Einstein condensates in a uniform magnetic field. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  P. Schuck,et al.  Thomas-Fermi approximation for Bose-Einstein condensates in traps , 2000 .

[15]  S. K. Adhikari,et al.  Vortex-bright solitons in a spin-orbit-coupled spin-1 condensate , 2017 .

[16]  Vladimir Loncar,et al.  OpenMP GNU and Intel Fortran programs for solving the time-dependent Gross-Pitaevskii equation , 2017, Comput. Phys. Commun..

[18]  Paulsamy Muruganandam,et al.  Hybrid OpenMP/MPI programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap , 2016, Comput. Phys. Commun..

[19]  G. Dresselhaus.  Spin-Orbit Coupling Effects in Zinc Blende Structures , 1955 .

[20]  J. Dalibard,et al.  Colloquium: Artificial gauge potentials for neutral atoms , 2010, 1008.5378.

[21]  V. Yukalov.  Dipolar and spinor bosonic systems , 2018, 1804.04873.

[22]  Yun Li,et al.  Bose-Einstein condensation with spin-orbit coupling , 2014, 1410.5526.

[23]  S. Kokkelmans,et al.  Interisotope determination of ultracold rubidium interactions from three high-precision experiments. , 2001, Physical review letters.

[24]  S. K. Adhikari,et al.  Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap , 2009, Comput. Phys. Commun..

[25]  Sadhan K. Adhikari,et al.  OpenMP Fortran and C programs for solving the time-dependent Gross-Pitaevskii equation in an anisotropic trap , 2016, Comput. Phys. Commun..

[26]  T. Kita,et al.  Axisymmetric versus nonaxisymmetric vortices in spinor Bose-Einstein condensates , 2002 .

[27]  Masahito Ueda,et al.  Spinor Bose-Einstein condensates , 2010, Quantum Atom Optics.

[28]  G. Fitzgerald,et al.  'I. , 2019, Australian journal of primary health.

[29]  Sandeep Gautam,et al.  Mobile vector soliton in a spin–orbit coupled spin-1 condensate , 2015, 1505.06193.

[30]  D. Campbell,et al.  Magnetic phases of spin-1 spin–orbit-coupled Bose gases , 2016, Nature Communications.