Source author record

Sadhan K. Adhikari

Sadhan K. Adhikari appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
16topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

CUDA programs for solving the time-dependent dipolar Gross-Pitaevskii equation in an anisotropic trap

In this paper we present new versions of previously published numerical programs for solving the dipolar Gross-Pitaevskii (GP) equation including the contact interaction in two and three spatial dimensions in imaginary and in real time, yielding both stationary and non-stationary solutions. New versions of programs were developed using CUDA toolkit and can make use of Nvidia GPU devices. The algorithm used is the same split-step semi-implicit Crank-Nicolson method as in the previous version (R. Kishor Kumar et al., Comput. Phys. Commun. 195, 117 (2015)), which is here implemented as a series of CUDA kernels that compute the solution on the GPU. In addition, the Fast Fourier Transform (FFT) library used in the previous version is replaced by cuFFT library, which works on CUDA-enabled GPUs. We present speedup test results obtained using new versions of programs and demonstrate an average speedup of 12 to 25, depending on the program and input size.

preprint2022arXiv

OpenMP, OpenMP/MPI, and CUDA/MPI C programs for solving the time-dependent dipolar Gross-Pitaevskii equation

We present new versions of the previously published C and CUDA programs for solving the dipolar Gross-Pitaevskii equation in one, two, and three spatial dimensions, which calculate stationary and non-stationary solutions by propagation in imaginary or real time. Presented programs are improved and parallelized versions of previous programs, divided into three packages according to the type of parallelization. First package contains improved and threaded version of sequential C programs using OpenMP. Second package additionally parallelizes three-dimensional variants of the OpenMP programs using MPI, allowing them to be run on distributed-memory systems. Finally, previous three-dimensional CUDA-parallelized programs are further parallelized using MPI, similarly as the OpenMP programs. We also present speedup test results obtained using new versions of programs in comparison with the previous sequential C and parallel CUDA programs. The improvements to the sequential version yield a speedup of 1.1 to 1.9, depending on the program. OpenMP parallelization yields further speedup of 2 to 12 on a 16-core workstation, while OpenMP/MPI version demonstrates a speedup of 11.5 to 16.5 on a computer cluster with 32 nodes used. CUDA/MPI version shows a speedup of 9 to 10 on a computer cluster with 32 nodes.

preprint2016arXiv

Hybrid OpenMP/MPI programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

We present hybrid OpenMP/MPI (Open Multi-Processing/Message Passing Interface) parallelized versions of earlier published C programs (D. Vudragovic et al., Comput. Phys. Commun. 183, 2021 (2012)) for calculating both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation in three spatial dimensions. The GP equation describes the properties of dilute Bose-Einstein condensates at ultra-cold temperatures. Hybrid versions of programs use the same algorithms as the C ones, involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method, but consider only a fully-anisotropic three-dimensional GP equation, where algorithmic complexity for large grid sizes necessitates parallelization in order to reduce execution time and/or memory requirements per node. Since distributed memory approach is required to address the latter, we combine MPI programing paradigm with existing OpenMP codes, thus creating fully flexible parallelism within a combined distributed/shared memory model, suitable for different modern computer architectures. The two presented C/OpenMP/MPI programs for real- and imaginary-time propagation are optimized and accompanied by a customizable makefile. We present typical scalability results for the provided OpenMP/MPI codes and demonstrate almost linear speedup until inter-process communication time starts to dominate over calculation time per iteration. Such a scalability study is necessary for large grid sizes in order to determine optimal number of MPI nodes and OpenMP threads per node.

preprint2016arXiv

OpenMP Fortran and C programs for solving the time-dependent Gross-Pitaevskii equation in an anisotropic trap

We present new version of previously published Fortran and C programs for solving the Gross-Pitaevskii equation for a Bose-Einstein condensate with contact interaction in one, two and three spatial dimensions in imaginary and real time, yielding both stationary and non-stationary solutions. To reduce the execution time on multicore processors, new versions of parallelized programs are developed using Open Multi-Processing (OpenMP) interface. The input in the previous versions of programs was the mathematical quantity nonlinearity for dimensionless form of Gross-Pitaevskii equation, whereas in the present programs the inputs are quantities of experimental interest, such as, number of atoms, scattering length, oscillator length for the trap, etc. New output files for some integrated one- and two-dimensional densities of experimental interest are given. We also present speedup test results for the new programs.

preprint2012arXiv

C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

We present C programming language versions of earlier published Fortran programs (Muruganandam and Adhikari, Comput. Phys. Commun. 180 (2009) 1888) for calculating both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation. The GP equation describes the properties of dilute Bose-Einstein condensates at ultra-cold temperatures. C versions of programs use the same algorithms as the Fortran ones, involving real- and imaginary-time propagation based on a split-step Crank-Nicolson method. In a one-space-variable form of the GP equation, we consider the one-dimensional, two-dimensional, circularly-symmetric, and the three-dimensional spherically-symmetric harmonic-oscillator traps. In the two-space-variable form, we consider the GP equation in two-dimensional anisotropic and three-dimensional axially-symmetric traps. The fully-anisotropic three-dimensional GP equation is also considered. In addition to these twelve programs, for six algorithms that involve two and three space variables, we have also developed threaded (OpenMP parallelized) programs, which allow numerical simulations to use all available CPU cores on a computer. All 18 programs are optimized and accompanied by makefiles for several popular C compilers. We present typical results for scalability of threaded codes and demonstrate almost linear speedup obtained with the new programs, allowing a decrease in execution times by an order of magnitude on modern multi-core computers.

preprint2012arXiv

Study of a degenerate dipolar Fermi gas of 161Dy atoms

We study properties of a single-component (spin polarized) degenerate dipolar Fermi gas of 161Dy atoms using a hydrodynamic description. Under axially-symmetric trapping we suggest reduced one- (1D) and two-dimensional (2D) description of the same for cigar and disk shapes, respectively. In addition to a complete numerical solution of the hydrodynamic model we also consider a variational approximation of the same. For a trapped system under appropriate conditions, the variational approximation as well as the reduced 1D and 2D models are found to yield results for shape, size and chemical potential of the system in agreement with the full numerical solution of the three-dimensional (3D) model. For the uniform system we consider anisotropic sound propagation in 3D. An analytical result for anisotropic sound propagation in uniform dipolar degenerate Fermi gas is found to be in agreement with results of numerical simulation in 3D.

preprint2009arXiv

Quenching of $para$-H$_2$ with an ultra-cold anti-hydrogen atom $\bar{H}_{1s}$

In this work we report the results concerning calculations for quantum-mechanical rotational transitions in molecular hydrogen, H$_2$, induced by an ultra-cold ground state anti-hydrogen atom $\bar{H}_{1s}$. The calculations are accomplished using a non-reactive close-coupling quantum-mechanical approach. The H$_2$ molecule is treated as a rigid rotor. The total elastic scattering cross section $σ_{el}(ε)$ at energy $ε$, state-resolved rotational transition cross sections $σ_{jj'}(ε)$ between states $j$ and $j'$ and corresponding thermal rate coefficients $k_{jj'}(T)$ are computed in the temperature range 0.004 K $ \lesssim T \lesssim$ 4 K. Satisfactory agreement with other calculations (variational) has been obtained for $σ_{el}(ε)$.

preprint2009arXiv

Universal scaling in a trapped Fermi super-fluid in the BCS-unitarity crossover

Using numerical simulation based on a density-functional equation for a trapped Fermi super-fluid valid along the BCS-unitarity crossover, we establish robust scaling over many orders of magnitude in the observables of the system as a function of fermion number. This scaling allows to predict the static properties of the system, such as energy, chemical potential, etc., for a large number of fermions, over the crossover, from the knowledge of those for a small number ($\sim 4-10$) of fermions.

preprint2006arXiv

Black soliton in a quasi-one-dimensional trapped fermion-fermion mixture

Employing a time-dependent mean-field-hydrodynamic model we study the generation of black solitons in a degenerate fermion-fermion mixture in a cigar-shaped geometry using variational and numerical solutions. The black soliton is found to be the first stationary vibrational excitation of the system and is considered to be a nonlinear continuation of the vibrational excitation of the harmonic oscillator state. We illustrate the stationary nature of the black soliton, by studying different perturbations on it after its formation.

preprint2006arXiv

Dissipation-managed soliton in a quasi-one-dimensional Bose-Einstein condensate

We use the time-dependent mean-field Gross-Pitaevskii equation to study the formation of a dynamically-stabilized dissipation-managed bright soliton in a quasi-one-dimensional Bose-Einstein condensate (BEC). Because of three-body recombination of bosonic atoms to molecules, atoms are lost (dissipated) from a BEC. Such dissipation leads to the decay of a BEC soliton. We demonstrate by a perturbation procedure that an alimentation of atoms from an external source to the BEC may compensate for the dissipation loss and lead to a dynamically-stabilized soliton. The result of the analytical perturbation method is in excellent agreement with mean-field numerics. It seems possible to obtain such a dynamically-stabilized BEC soliton without dissipation in laboratory.

preprint2003arXiv

Matter-wave interference, Josephson oscillation and its disruption in a Bose-Einstein condensate on an optical lattice

Using the axially-symmetric time-dependent mean-field Gross-Pitaevskii equation we study the Josephson oscillation in a repulsive Bose-Einstein condensate trapped by a harmonic plus an one-dimensional optical-lattice potential to describe the experiments by Cataliotti et. al. [Science 293 (2001) 843, New J. Phys. 5 (2003) 71.1]. After a study of the formation of matter-wave interference upon releasing the condensate from the optical trap, we directly investigate the alternating atomic superfluid Josephson current upon displacing the harmonic trap along the optical axis. The Josephson current is found to be disrupted upon displacing the harmonic trap through a distance greater than a critical distance signaling a superfluid to a classical insulator transition in the condensate.

preprint2002arXiv

Effect of an impulsive force on vortices in a rotating Bose-Einstein condensate

The effects of a sudden increase and decrease of the interatomic interaction and harmonic-oscillator trapping potential on vortices in a quasi two-dimensional rotating Bose-Einstein condensate are investigated using the mean-field Gross-Pitaevskii equation. Upon increasing the strength of interaction suddenly the condensate enters a nonstationary oscillating phase which starts to develop more vortices. The opposite happens if the strength is reduced suddenly. Eventually, the number of vortices attains a final value at large times. Similarly, the number of vortices increases (decreases) upon a sudden reduction (augmentation) in the trapping potential. We also study the decay of vortices when the rotation of the condensate is suddenly stopped. Upon a free expansion of a rotating BEC with vortices the radius of the vortex core increases more rapidly than the radius of the condensate. This makes the counting and detection of multiple vortex easier after a free expansion.

preprint2002arXiv

Low-energy muon-transfer reaction from hydrogen isotopes to helium isotopes

Direct muon transfer in low-energy collisions of the muonic hydrogen H$_μ$ and helium (He$^{++}$) is considered in a three-body quantum-mechanical framework of coordinate-space integro-differential Faddeev-Hahn-type equations within two- and six-state close coupling approximations. The final-state Coulomb interaction is treated without any approximation employing appropriate Coulomb waves in the final state. The present results agree reasonably well with previous semiclassical calculations.

preprint2001arXiv

Positronium-positronium interaction: Resonance, scattering length, and Bose-Einstein condensation

The low-energy scattering of ortho positronium (Ps) by ortho Ps has been studied in a full quantum mechanical coupled-channel approach. In the singlet channel (total spin $s_T=0$) we find S- and P-wave resonances at 3.35 eV (width 0.02 eV) and 5.05 eV (width 0.04 eV), respectively and a binding of 0.44 eV of Ps$_2$. The scattering length for $s_T=0$ is 3.95 Åand for $s_T=2$ is 0.83 Å. The small $s_T=2$ scattering length makes the spin-polarized ortho Ps atoms an almost noninteracting ideal gas which may undergo Bose-Einstein condensation.

preprint2000arXiv

Numerical Solution of the Two-Dimensional Gross-Pitaevskii Equation for Trapped Interacting Atoms

We present a numerical scheme for solving the time-independent nonlinear Gross-Pitaevskii equation in two dimensions describing the Bose-Einstein condensate of trapped interacting neutral atoms at zero temperature. The trap potential is taken to be of the harmonic-oscillator type and the interaction both attractive and repulsive. The Gross-Pitaevskii equation is numerically integrated consistent with the correct boundary conditions at the origin and in the asymptotic region. Rapid convergence is obtained in all cases studied. In the attractive case there is a limit to the maximum number of atoms in the condensate.

preprint2000arXiv

Scattering of positronium by H, He, Ne, and Ar

The low-energy scattering of ortho positronium (Ps) by H, He, Ne, and Ar atoms has been investigated in the coupled-channel framework by using a recently proposed time-reversal-symmetric nonlocal electron-exchange model potential with a single parameter $C$. For H and He we use a three-Ps-state coupled-channel model and for Ar and Ne we use a static-exchange model. The sensitivity of the results is studied with respect to the parameter $C$. Present low-energy cross sections for He, Ne and Ar are in good agreement with experiment.

preprint1999arXiv

Electron exchange model potential: Application to positronium-helium scattering

The formulation of a suitable nonlocal model potential for electron exchange is presented, checked with electron-hydrogen and electron-helium scattering, and applied to the study of elastic and inelastic scattering and ionization of ortho positronium (Ps) by helium. The elastic scattering and the $n = 2$ excitations of Ps are investigated using a three-Ps-state close-coupling approximation. The higher ($n\ge 3$) excitations and ionization of Ps atom are treated in the framework of Born approximation with present exchange. Calculations are reported of phase shifts, and elastic, Ps-excitation, and total cross sections. The present target elastic total cross section agrees well with experimental results at thermal to medium energies.

preprint1993arXiv

Relativistic three-particle dynamical equations: II. Application to the trinucleon system

We calculate the contribution of relativistic dynamics on the neutron-deuteron scattering length and triton binding energy employing five sets trinucleon potential models and four types of three-dimensional relativistic three-body equations suggested in the preceding paper. The relativistic correction to binding energy may vary a lot and even change sign depending on the relativistic formulation employed. The deviations of these observables from those obtained in nonrelativistic models follow the general universal trend of deviations introduced by off- and on-shell variations of two- and three-nucleon potentials in a nonrelativistic model calculation. Consequently, it will be difficult to separate unambiguously the effect of off- and on-shell variations of two- and three- nucleon potentials on low-energy three-nucleon observables from the effect of relativistic dynamics.