Source author record

Paulsamy Muruganandam

Paulsamy Muruganandam 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

15works
14topics
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

15 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

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

Many of the static and dynamic properties of an atomic Bose-Einstein condensate (BEC) are usually studied by solving the mean-field Gross-Pitaevskii (GP) equation, which is a nonlinear partial differential equation for short-range atomic interaction. More recently, BEC of atoms with long-range dipolar atomic interaction are used in theoretical and experimental studies. For dipolar atomic interaction, the GP equation is a partial integro-differential equation, requiring complex algorithm for its numerical solution. Here we present numerical algorithms for both stationary and non-stationary solutions of the full three-dimensional (3D) GP equation for a dipolar BEC, including the contact interaction. We also consider the simplified one- (1D) and two-dimensional (2D) GP equations satisfied by cigar- and disk-shaped dipolar BECs. We employ the split-step Crank-Nicolson method with real- and imaginary-time propagations, respectively, for the numerical solution of the GP equation for dynamic and static properties of a dipolar BEC. The atoms are considered to be polarized along the z axis and we consider ten different cases, e.g., stationary and non-stationary solutions of the GP equation for a dipolar BEC in 1D (along x and z axes), 2D (in x-y and x-z planes), and 3D, and we provide working codes in Fortran 90/95 and C for these ten cases (twenty programs in all). We present numerical results for energy, chemical potential, root-mean-square sizes and density of the dipolar BECs and, where available, compare them with results of other authors and of variational and Thomas-Fermi approximations.

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.

preprint2021arXiv

Influence of Rashba spin-orbit and Rabi couplings on the miscibility and ground state phases of binary Bose-Einstein condensates

We study the miscibility properties and ground state phases of two-component spin-orbit (SO) coupled Bose-Einstein condensates (BECs) in a harmonic trap with strong axial confinement. By numerically solving the coupled Gross-Pitaevskii equations in the two-dimensional setting, we analyze the SO-coupled BECs for two possible permutations of the intra- and interspecies interactions, namely (i) weak intra- and weak interspecies interactions (W-W) and (ii) weak intra- and strong interspecies interactions (W-S). Considering the density overlap integral as a miscibility order parameter, we investigate the miscible-immiscible transition by varying the coupling parameters. We obtain various ground state phases, including plane wave, half quantum vortex, elongated plane wave, and different stripe wave patterns for W-W interactions. For finite Rabi coupling, an increase in SO coupling strength leads to the transition from the fully miscible to the partially miscible state. We also characterize different ground states in the coupling parameter space using the root mean square sizes of the condensate. The spin density vector for the ground state phases exhibits density, quadrupole and dipole like spin polarizations. For the W-S interaction, in addition to that observed in the W-W case, we witness semi vortex, mixed mode, and shell-like immiscible phases. We notice a wide variety of spin polarizations, such as density, dipole, quadrupole, symbiotic, necklace, and stripe-like patterns for the W-S case. A detailed investigation in the coupling parameter space indicates immiscible to miscible state phase transition upon varying the Rabi coupling for a fixed Rashba SO coupling. The critical Rabi coupling for the immiscible-miscible phase transition decreases upon increasing the SO coupling strength.

preprint2021arXiv

Spiral wave chimera-like transient dynamics in three-dimensional grid of diffusive ecological systems

In the present article, we demonstrate the emergence and existence of the spiral wave chimera-like transient pattern in coupled ecological systems, composed of prey-predator patches, where the patches are connected in a three-dimensional medium through local diffusion. We explore the transition scenarios among the several collective dynamical behaviors together with transient spiral wave chimera-like states and investigate the long time behavior of these states. The transition from the transient spiral chimera-like pattern to the long time synchronized or desynchronized pattern appears through the deformation of the incoherent region of the spiral core. We discuss the transient dynamics under the influence of the species diffusion at different time instants. By calculating the instantaneous strength of incoherence of the populations, we estimate the duration of the transient dynamics characterized by the persistence of the chimera-like spatial coexistence of coherent and incoherent patterns over the spatial domain. We generalize our observations on the transient dynamics in three-dimensional grid of diffusive ecological systems by considering two different prey-predator systems.

preprint2020arXiv

Quenching dynamics of the bright solitons and other localized states in spin-orbit coupled Bose-Einstein condensates

We study the dynamics of binary Bose-Einstein condensates made of ultracold and dilute alkali-metal atoms in a quasi-one-dimensional setting. Numerically solving the two coupled Gross-Pitaevskii equations which accurately describe the system dynamics, we demonstrate that the spin transport can be controlled by suitably quenching spin-orbit (SO) and Rabi coupling strengths. Moreover, we predict a variety of dynamical features induced by quenching: broken oscillations, breathers-like oscillating patterns, spin-mixing-demixing, miscible-immiscible transition, emerging dark-bright states, dark solitons, and spin-trapping dynamics. We also outline the experimental relevance of the present study in manipulating the spin states in $^{39}$K condensates.

preprint2016arXiv

Dynamics of trapped interacting vortices in Bose-Einstein condensates: Role of breathing degree of freedom

With use of a variational principle, we investigate a role of breathing width degree of freedom in the effective theory of interacting vortices in a trapped single-component Bose-Einstein condensates in 2 dimensions under the strong repulsive cubic nonlinearity. As for the trial function, we choose a product of two vortex functions assuming a pair interaction and employ the amplitude form of each vortex function in the Padé approximation which accommodates a hallmark of the vortex core. We have obtained Lagrange equation for the interacting vortex-core coordinates coupled with the time-derivative of width and also its Hamilton formalism by having recourse to a non-standard Poisson bracket. By solving the Hamilton equation, we find rapid radial breathing oscillations superposed on the slower rotational motion of vortex cores, consistent with numerical solutions of Gross-Pitaevskii equation. In higher-energy states of 2 vortex systems, the breathing width degree of freedom plays role of a kicking in the kicked rotator and generates chaos with a structure of sea-urchin needles. Byproduct of the present variational approach includes: (1) the charge-dependent logarithmic inter-vortex interaction multiplied with a pre-factor which depends on the scalar product of a pair of core-position vectors; (2) the charge-independent short-range repulsive inter-vortex interaction and spring force.

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.

preprint2016arXiv

Three-dimensional vortex structures in a rotating dipolar Bose-Einstein condensate

We study the three-dimensional ground state vortex lattice structures of purely dipolar Bose-Einstein condensate (BEC). By using the mean-field model we obtain a stability diagram for the vortex states in purely dipolar BECs as a function of harmonic trap aspect ratio ($λ$) and dipole-dipole interaction strength ($D$) under rotation. Rotating the condensate within the unstable region leads to collapse, while in the stable region furnishes stable vortex lattices of dipolar BECs. We analyse stable vortex lattice structures by solving three-dimensional time dependent Gross-Pitaevskii equation in imaginary time. Further, the stability of vortex states is examined by evolution in real-time. We also investigate distribution of vortices in fully anisotropic trap by increasing eccentricity of the external trapping potential. Breaking up of the condensate in two parts has been observed with equally shared vortices on each when the trap is sufficiently weak and the rotation frequency is high.

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

Identifying financial crises in real time

Following the thermodynamic formulation of multifractal measure that was shown to be capable of detecting large fluctuations at an early stage, here we propose a new index which permits us to distinguish events like financial crisis in real time . We calculate the partition function from where we obtain thermodynamic quantities analogous to free energy and specific heat. The index is defined as the normalized energy variation and it can be used to study the behavior of stochastic time series, such as financial market daily data. Famous financial market crashes - Black Thursday (1929), Black Monday (1987) and Subprime crisis (2008) - are identified with clear and robust results. The method is also applied to the market fluctuations of 2011. From these results it appears as if the apparent crisis of 2011 is of a different nature from the other three. We also show that the analysis has forecasting capabilities.

preprint2007arXiv

Analytical calculation of the transition to complete phase synchronization in coupled oscillators

Here we present a system of coupled phase oscillators with nearest neighbors coupling, which we study for different boundary conditions. We concentrate at the transition to total synchronization. We are able to develop exact solutions for the value of the coupling parameter when the system becomes completely synchronized, for the case of periodic boundary conditions as well as for an open chain with fixed ends. We compare the results with those calculated numerically.

preprint2007arXiv

Low dimensional behavior in three-dimensional coupled map lattices

The analysis of one-, two-, and three-dimensional coupled map lattices is here developed under a statistical and dynamical perspective. We show that the three-dimensional CML exhibits low dimensional behavior with long range correlation and the power spectrum follows $1/f$ noise. This approach leads to an integrated understanding of the most important properties of these universal models of spatiotemporal chaos. We perform a complete time series analysis of the model and investigate the dependence of the signal properties by change of dimension.

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.