Institute of Physics Belgrade

Permanent URI for this community

https://pub.ipb.ac.rs/handle/123456789/1

Browse

Browsing Institute of Physics Belgrade by Author "Adhikari, Sadhan K."

Now showing 1 - 4 of 4
  • No Thumbnail Available
    Item
    C and Fortran OpenMP programs for rotating Bose–Einstein condensates
    Kishor Kumar, Ramavarmaraja; Lončar, Vladimir; Muruganandam, Paulsamy; Adhikari, Sadhan K.; Balaž, Antun
    We present OpenMP versions of C and Fortran programs for solving the Gross–Pitaevskii equation for a rotating trapped Bose–Einstein condensate (BEC) in two (2D) and three (3D) spatial dimensions. The programs can be used to generate vortex lattices and study dynamics of rotating BECs. We use the split-step Crank–Nicolson algorithm for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively. The simulation input parameters for the C programs are provided via input files, while for the Fortran programs they are given at the beginning of each program and therefore their change requires recompilation of the corresponding program. The programs propagate the condensate wave function and calculate several relevant physical quantities, such as the energy, the chemical potential, and the root-mean-square sizes. The imaginary-time propagation starts with an analytic wave function with one vortex at the trap center, modulated by a random phase at different space points. Nevertheless, the converged wave function for a rapidly rotating BEC with a large number of vortices is most efficiently calculated using the pre-calculated converged wave function of a rotating BEC containing a smaller number of vortices as the initial state rather than using an analytic wave function with one vortex as the initial state. These pre-calculated initial states exhibit rapid convergence for fast-rotating condensates to states containing multiple vortices with an appropriate phase structure. This is illustrated here by calculating vortex lattices with up to 61 vortices in 2D and 3D. Outputs of the programs include calculated physical quantities, as well as the wave function and different density profiles (full density, integrated densities in lower dimensions, and density cross-sections). The provided real-time propagation programs can be used to study the dynamics of a rotating BEC using the imaginary-time stationary wave function as the initial state. We also study the efficiency of parallelization of the present OpenMP C and Fortran programs with different compilers. Program summary Program title: BEC-GP-ROT-OMP, consisting of: (1) BEC-GP-ROT-OMP-C package, containing programs (i) bec-gp-rot-2d-th and (ii) bec-gp-rot-3d-th; (2) BEC-GP-ROT-OMP-F package, containing programs (i) bec-gp-rot-2d-th and (ii) bec-gp-rot-3d-th. Program files doi: http://dx.doi.org/10.17632/cw7tkn22v2.2 Licensing provisions: Apache License 2.0 Programming language: OpenMP C; OpenMP Fortran. The C programs are tested with the GNU, Intel, PGI, Oracle, and Clang compiler, and the Fortran programs are tested with the GNU, Intel, PGI, and Oracle compiler. Nature of problem: The present Open Multi-Processing (OpenMP) C and Fortran programs solve the time-dependent nonlinear partial differential Gross–Pitaevskii (GP) equation for a trapped rotating Bose–Einstein condensate in two (2D) and three (3D) spatial dimensions in a fully anisotropic traps. Solution method: We employ the split-step Crank–Nicolson algorithm 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.
  • No Thumbnail Available
    Item
    OpenMP Fortran programs for solving the time-dependent dipolar Gross-Pitaevskii equation
    Young-S., Luis E.; Muruganandam, Paulsamy; Balaž, Antun; Adhikari, Sadhan K.
    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 (Kishor Kumar et al., 2015 [3]). These OpenMP versions have significantly reduced execution time in multicore processors.
  • No Thumbnail Available
    Item
    OpenMP solver for rotating spin-1 spin–orbit- and Rabi-coupled Bose–Einstein condensates
    Muruganandam, Paulsamy; Balaž, Antun; Adhikari, Sadhan K.
    We present OpenMP version of a Fortran program for solving the Gross–Pitaevskii equation for a harmonically trapped three-component rotating spin-1 spinor Bose–Einstein condensate (BEC) in two spatial dimensions with or without spin–orbit (SO) and Rabi couplings. The program uses either Rashba or Dresselhaus SO coupling. We use the split-step Crank–Nicolson discretization scheme for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively.
  • No Thumbnail Available
    Item
    Spin-1 spin–orbit- and Rabi-coupled Bose–Einstein condensate solver
    Ravisankar, Rajamanickam; Vudragović, Dušan; Muruganandam, Paulsamy; Balaž, Antun; Adhikari, Sadhan K.
    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.