Source author record

S. K. Adhikari

S. 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

56works
13topics
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

56 published item(s)

preprint2022arXiv

Deep inelastic collision of two-dimensional anisotropic dipolar condensate solitons

The possibility of generating stable anisotropic solitons in dipolar Bose-Einstein condensates (BECs) leads to a different scenario not possible in a nondipolar BEC with contact interaction. We study the statics and dynamics of anisotropic bright solitons in quasi-two-dimensional BECs consisting of polarized dipolar atoms.We study the collision dynamics of two such solitons at different velocities for different angles between the polarization and collision directions. The collision is found to be quasi elastic at large velocities. At small velocities the collision is inelastic leading to the formation of a coalesced soliton in an excited scissors mode, monopole mode or quadrupole mode. Also, at small velocities, after collision, a large change of direction of motion of the solitons is possible. The investigation is performed through a numerical solution of the underlying mean-field Gross-Pitaevskii equation.

preprint2022arXiv

Effect of quartic-quintic beyond-mean-field interactions on a self-bound dipolar droplet

We study the effect of beyond-mean-field quantum-fluctuation (QF) Lee-Huang-Yang (LHY) and three-body interactions, with quartic and quintic nonlinearities, respectively, on the formation of a stable self-repulsive (positive scattering length $a$) and a self-attractive (negative $a$) self-bound dipolar Bose-Einstein condensate (BEC) droplet in free space under the action of two-body contact and dipolar interactions. Previous studies of dipolar droplets considered either the LHY interaction or the three-body interaction, as either of these interactions could stabilize a dipolar BEC droplet against collapse. We find that the effect of three-body recombination on the formation of a dipolar droplet could be quite large and for a complete description of the problem both the QF LHY and three-body interactions should be considered simultaneously, where appropriate. In the self-repulsive case for small $a$ and in the self-attractive case, no appropriate LHY interaction is known and only three-body interaction should be used, otherwise both beyond-mean-field interactions should be used. We consider a numerical solution of a highly-nonlinear beyond-mean-field model as well as a variational approximation to it in this investigation and present results for size, shape and energy of a dipolar droplet of polarized $^{164}$Dy atoms. The shape is filament-like, along the polarization direction, and could be long, for a large number of atoms $N$, short for small $N$, thin for negative $a$ and small positive $a$, and fat for large positive $a$.

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

Supersolid-like solitons in spin-orbit coupled spin-2 condensate

We study supersolid-like crystalline structures emerging in the stationary states of a quasi-two-dimensional spin-orbit (SO)-coupled spin-2 condensate in the ferromagnetic, cyclic, and antiferro-magnetic phases by solving a mean-field model.Interplay of different strengths of SO coupling and interatomic interactions gives rise to a variety of non-trivial density patterns in the emergent solutions. For small SO-coupling strengths $γ$ ($γ\approx 0.5$), the ground state is an axisymmetric multi-ring soliton for polar, cyclic and weakly-ferromagnetic interactions, whereas for stronger-ferromagnetic interactions a circularly-asymmetric soliton emerges as the ground state.Depending on the values of interaction parameters, with an increase in SO-coupling strength, a stripe phase may also emerge as the ground state for polar and cyclic interactions. For intermediate values of SO-coupling strength ($γ\approx 1$), in addition to these solitons, one could have a quasi-degenerate triangular-lattice soliton in all magnetic phases. On further increasing the SO-coupling strength ($γ\gtrapprox 4$), a square-lattice and a superstripe soliton emerge as quasi-degenerate states. The emergence of all these solitons can be inferred from a study of solutions of the single-particle Hamiltonian.

preprint2022arXiv

Supersolid-like square- and honeycomb-lattice crystallization of droplets in a dipolar condensate

We demonstrate a supersolid-like spatially-periodic square- and honeycomb-lattice crystallization of droplets, in addition to the commonly-studied triangular-lattice crystallization, in a cylindrically-symmetric quasi-two-dimensional trapped dipolar condensate, using a beyond-mean-field model including a quantum-fluctuation Lee-Huang-Yang-type interaction. These three types of crystallization of droplets may appear for the same atomic interactions and the same trap frequencies. The energy $E$ of all three crystallization as a function of number $N$ of atoms satisfy the universal scaling relation $E\sim N^{0.4}$ indicating that all three arrangements of the droplets should be energetically probable processes of phenomenological interest. The state of square-lattice crystallization may have the central site occupied or unoccupied, corresponding to a parity-symmetric or parity-antisymmetric state, respectively. The state of square-lattice crystallization with the occupied central site and the state of triangular-lattice crystallization, for a fixed $N$, constitute two quasi-degenerate ground states while the other states are low-lying excited states. This makes the square-lattice crystallization with the occupied central site an ideal candidate for future experimental observation.

preprint2021arXiv

Multi-ring, stripe, and super-lattice solitons in a spin-orbit coupled spin-1 condensate

We demonstrate exotic stable quasi-two-dimensional solitons in a Rashba or a Dresselhaus spin-orbit (SO) coupled hyperfine spin-1 ($F=1$) trap-less antiferromagnetic Bose-Einstein condensate using the mean-field Gross-Pitaevskii equation. For weak SO coupling, the solitons are of the $(-1,0,+1)$ or $(+1,0,-1)$ type with intrinsic vorticity, for Rashba or Dresselhaus SO coupling, where the numbers in the parentheses denote angular momentum in spin components $F_z= +1,0,-1$, respectively. For intermediate SO coupling, the solitons have multi-ring structure maintaining the above-mentioned vortices in respective components. For larger SO coupling, superlattice solitons with a square-lattice structure in total density are found in addition to stripe solitons with stripe pattern in component densities with no periodic modulation in total density.

preprint2020arXiv

Spontaneous spatial order in two-dimensional ferromagnetic spin-orbit coupled uniform spin-1 condensate solitons

We demonstrate spontaneous spatial order in stripe and super-lattice solitons in a Rashba spin-orbit (SO) coupled spin-1 uniform quasi-two-dimensional (quasi-2D) ferromagnetic Bose-Einstein condensate. For weak SO coupling, the solitons are of the $(- 1,0,+ 1)$ or $(0,+ 1,+ 2)$ type with intrinsic vorticity, where the numbers in the parentheses denote angular momentum in spin components $F_z=+1,0,-1$, respectively. For intermediate SO coupling, three types of solitons are found: (a) circularly-asymmetric solitons, (b) circularly-symmetric $(- 1,0,+ 1)$- and (c) $(0,+ 1,+ 2)$-type multi-ring solitons maintaining the above-mentioned vortices in respective components. For large SO coupling, quasi-degenerate stripe and super-lattice solitons are found in addition to the circularly-asymmetric solitons. A super-lattice soliton forms a 2D square lattice structure in the total density as in a super-solid; in component densities it may have either (i) a 1D stripe pattern or (ii) a 2D square lattice structure.

preprint2019arXiv

Phase-separated symmetry-breaking vortex-lattice in a binary Bose-Einstein condensate

We study spontaneous-symmetry-breaking circularly-asymmetric phase separation of vortex lattices in a rapidly rotating harmonically-trapped quasi-two-dimensional (quasi-2D) binary Bose-Einstein condensate (BEC) with repulsive inter- and intra-species interactions. The phase separated vortex lattices of the components appear in different regions of space with no overlap between the vortices of the two components, which will permit an efficient experimental observation of such vortices and accurate study of the effect of atomic interaction on such vortex lattice. Such phase separation takes place when the intra-species interaction energies of the two components are equal or nearly equal with relatively strong inter-species repulsion. When the intra-species energies are equal, the two phase-separated vortex lattices have identical semicircular shapes with one being the parity conjugate of the other. When the intra-species energies are nearly equal, the phase separation is also complete but the vortex lattices have different shapes. We demonstrate our claim with a numerical solution of the mean-field Gross-Pitaevskii equation for a rapidly rotating quasi-2D binary BEC.

preprint2016arXiv

Dimensional reduction and localization of a Bose-Einstein condensate in a quasi-1D bichromatic optical lattice

We analyze the localization of a Bose-Einstein condensate (BEC) in a one-dimensional bichromatic quasi-periodic optical-lattice potential by numerically solving the 1D Gross-Pitaevskii equation (1D GPE). We first derive the 1D GPE from the dimensional reduction of the 3D quantum field theory of interacting bosons obtaining two coupled differential equations (for axial wavefuction and space-time dependent transverse width) which reduce to the 1D GPE under strict conditions. Then, by using the 1D GPE we report the suppression of localization in the interacting BEC when the repulsive scattering length between bosonic atoms is sufficiently large.

preprint2016arXiv

Stable and mobile two-dimensional dipolar ring-dark-in-bright Bose-Einstein condensate soliton

We demonstrate robust, stable, mobile two-dimensional (2D) dipolar ring-dark-in-bright (RDB) Bose-Einstein condensate (BEC) solitons for repulsive contact interaction, subject to a harmonic trap along the $y$ direction perpendicular to the polarization direction $z$. Such a RDB soliton has a ring-shaped notch (zero in density) imprinted on a 2D bright soliton free to move in the $x-z$ plane. At medium velocity the head-on collision of two such solitons is found to be quasi elastic with practically no deformation. The possibility of creating the RDB soliton by phase imprinting is demonstrated. The findings are illustrated using numerical simulation employing realistic interaction parameters in a dipolar $^{164}$Dy BEC.

preprint2015arXiv

Analytic models for density of a ground-state spinor condensate

We demonstrate that the ground state of a trapped spin-1 and spin-2 spinor ferromagnetic Bose-Einstein condensate (BEC) can be well approximated by a single decoupled Gross-Pitaevskii (GP) equation. Useful analytic models for the ground-state densities of ferromagnetic BECs are obtained from the Thomas-Fermi approximation (TFA) to this decoupled equation. Similarly, for the ground states of spin-1 anti-ferromagnetic and spin-2 anti-ferromagnetic and cyclic BECs, some of the spin component densities are zero which reduces the coupled GP equation to a simple reduced form. Analytic models for ground state densities are also obtained for anti-ferromagnetic and cyclic BECs from the TFA to the respective reduced GP equations. The analytic densities are illustrated and compared with the full numerical solution of the GP equation with realistic experimental parameters.

preprint2015arXiv

Low temperature HD + $ortho$-/$para$-H$_2$ inelastic scattering of astrophysical interest

State-selected total cross sections and thermal rate coefficients are computed for the HD + $ortho$-/$para$-H$_2$ rotational energy transfer collision at low temperatures: 2 K $\lesssim$ T $\lesssim$ 300 K. A modified H$_2$-H$_2$ potential energy surface (PES) devised by Hinde is used for this pure quantum-mechanical dynamical computation. A comparison of the new results for the HD + $ortho$-/$para$-H$_2$ scattering problem and previous calculations computed with the use of other older PESs is presented and discussed.

preprint2015arXiv

Mobile vector soliton in a spin-orbit coupled spin-$1$ condensate

We study the formation of bound states and three-component bright vector solitons in a quasi-one-dimensional spin-orbit-coupled hyperfine spin $f=1$ Bose-Einstein condensate using numerical solution and variational approximation of a mean-field model. In the antiferromagnetic domain, the solutions are time-reversal symmetric, and the component densities have multi-peak structure. In the ferromagnetic domain, the solutions violate time-reversal symmetry, and the component densities have single-peak structure. The dynamics of the system is not Galelian invariant. From an analysis of Galelian invariance, we establish that the single-peak ferromagnetic vector solitons are true solitons and can move maintaining constant component densities, whereas the antiferromagnetic solitons cannot move with constant component densities.

preprint2015arXiv

Stable matter-wave soliton in the vortex core of a uniform condensate

We demonstrate a stable, mobile, dipolar or nondipolar three-dimensional matter-wave soliton in the vortex core of a uniform nondipolar condensate. All intra- and inter-species contact interactions can be repulsive for a strongly dipolar soliton. For a weakly dipolar or nondipolar soliton, the intra-species contact interaction in the soliton should be attractive for the formation of a compact soliton. The soliton can propagate with a constant velocity along the vortex core without any deformation. Two such solitons undergo a quasi-elastic collision at medium velocities. We illustrate the findings using realistic interactions in a mean-field model of binary $^{87}$Rb-$^{85}$Rb and $^{87}$Rb-$^{164}$Dy systems.

preprint2015arXiv

Stable spatial and spatiotemporal optical soliton in the core of an optical vortex

We demonstrate a robust, stable, mobile, two-dimensional (2D) spatial and three-dimensional (3D) spatiotemporal optical soliton in the core of an optical vortex, while all nonlinearities are of the cubic (Kerr) type. The 3D soliton can propagate with a constant velocity along the vortex core without any deformation. Stability of the soliton under a small perturbation is established numerically. Two such solitons moving along the vortex core can undergo a quasi-elastic collision at medium velocities. Possibilities of forming such a 2D spatial soliton in the core of a vortical beam are discussed.

preprint2015arXiv

Vector solitons in a spin-orbit coupled spin-$2$ Bose-Einstein condensate

Five-component minimum-energy bound states and mobile vector solitons of a spin-orbit-coupled quasi-one-dimensional hyperfine-spin-2 Bose-Einstein condensate are studied using the numerical solution and variational approximation of a mean-field model. Two distinct types of solutions with single-peak and multi-peak density distribution of the components are identified in different domains of interaction parameters. From an analysis of Galilean invariance and time-reversal symmetry of the Hamiltonian, we establish that vector solitons with multi-peak density distribution preserve time-reversal symmetry, but cannot propagate maintaining the shape of individual components. However, those with single-peak density distribution violate time-reversal symmetry of the Hamiltonian, but can propagate with a constant velocity maintaining the shape of individual components.

preprint2014arXiv

Bright dipolar Bose-Einstein condensate soliton mobile in a direction perpendicular to polarization

We demonstrate stable, robust, bright dipolar Bose-Einstein condensate (BEC) solitons, moving in a direction perpendicular to the polarization direction, formed due to dipolar interaction for repulsive contact interaction. At medium velocity the head on collision of two such solitons is found to be quasi elastic with practically no deformation. Upon small perturbation the solitons are found to exhibit sustained breathing oscillation. The findings are illustrated by numerical simulation using the 3D mean-field Gross-Pitaevskii equation and a reduced 2D model in three and two spatial dimensions employing realistic interaction parameters for a dipolar $^{164}$Dy BEC.

preprint2014arXiv

Demixing and symmetry breaking in binary dipolar Bose-Einstein condensate solitons

We demonstrate fully demixed (separated) robust and stable bright binary dipolar Bose-Einstein condensate soliton in a quasi-one-dimensional setting formed due to dipolar interactions for repulsive contact interactions. For large repulsive interspecies contact interaction the first species may spatially separate from the second species thus forming a demixed configuration, which can be spatially-symmetric or symmetry-broken. In the spatially-symmetric case, one of the the species occupies the central region, whereas the other species separates into two equal parts and stay predominantly out of this central region. In the symmetry-broken case, the two species stay side by side. Stability phase diagrams for the binary solitons are obtained. The results are illustrated with realistic values of parameters in the binary 164Dy-168Er and 164Dy-162Dy mixtures. The demixed solitons are really soliton molecules formed of two types of atoms. A proposal for creating dipolar solitons in experiments is also presented.

preprint2014arXiv

Localization of a spin-orbit coupled Bose-Einstein condensate in a bichromatic optical lattice

We study the localization of a noninteracting and weakly interacting Bose-Einstein condensate with spin-orbit coupling loaded in a quasiperiodic bichromatic optical lattice potential using the numerical solution and variational approximation of a binary mean-field Gross-Pitaevskii equation with two pseudo-spin components. We confirm the existence of the stationary localized states in the presence of the spin-orbit and Rabi couplings for an equal distribution of atoms in the two components. We find that the interaction between the spin-orbit and Rabi couplings favors the localization or delocalization of the BEC depending on the the phase difference between the components. We also studied the oscillation dynamics of the localized states for an initial population imbalance between the two components.

preprint2014arXiv

Phase separation in a spin-orbit coupled Bose-Einstein condensate

We study a spin-orbit (SO) coupled hyperfine spin-1 Bose-Einstein condensate (BEC) in a quasi-one-dimensional trap. For a SO-coupled BEC in a one-dimensional box, we show that in the absence of the Rabi term, any non-zero value of SO coupling will result in a phase separation among the components for a ferromagnetic BEC, like $^{87}$Rb. On the other hand, SO coupling favors miscibility in a polar BEC, like $^{23}$Na. In the presence of a harmonic trap, which favors miscibility, a ferromagnetic BEC phase separates, provided the SO-coupling strength and number of atoms are greater than some critical value. The Rabi term favors miscibility irrespective of the nature of the spin interaction: ferromagnetic or polar.

preprint2014arXiv

Self trapping of a dipolar Bose-Einstein condensate in a double well

We study the Josephson oscillation and self trapping dynamics of a cigar-shaped dipolar Bose-Einstein condensate of $^{52}$Cr atoms polarized along the symmetry axis of an axially-symmetric double-well potential using the numerical solution of a mean-field model, for dominating repulsive contact interaction (large positive scattering length a) over an anisotropic dipolar interaction. Josephson-type oscillation emerges for small and very large number of atoms, whereas self trapping is noted for an intermediate number of atoms. The dipolar interaction pushes the system away from self trapping towards Josephson oscillation. We consider a simple two-mode description for a qualitative understanding of the dynamics.

preprint2014arXiv

Stable and mobile excited two-dimensional dipolar Bose-Einstein condensate solitons

We demonstrate robust, stable, mobile excited states of quasi-two-dimensional (quasi-2D) dipolar Bose-Einstein condensate (BEC) solitons for repulsive contact interaction with a harmonic trap along the $x$ direction perpendicular to the polarization direction $z$. Such a soliton can freely move in the y-z plane. A rich variety of such excitations is considered: one quanta of excitation for movement along (i) y axis or (ii) z axis or (ii) both. A proposal for creating these excited solitonic states in a laboratory by phase imprinting is also discussed. We also consider excited states of quasi-2D dipolar BEC soliton where the sign of the dipolar interaction is reversed by a rotating orienting field. In this sign-changed case the soliton moves freely in the x-y plane under the action of a harmonic trap in the $z$ direction. At medium velocity the head-on collision of two such solitons is found to be quasi elastic with practically no deformation. The findings are illustrated using numerical simulation in three and two spatial dimensions employing realistic interaction parameters for a dipolar $^{164}$Dy BEC.

preprint2014arXiv

Statics and dynamics of a binary dipolar Bose-Einstein condensate soliton

We study the statics and dynamics of a binary dipolar Bose-Einstein condensate soliton for repulsive inter- and intraspecies contact interactions with the two components subject to different spatial symmetries $-$ distinct quasi-one-dimensional and quasi-two-dimensional shapes - using numerical solution and variational approximation of a three-dimensional mean-field model. The results are illustrated with realistic values of parameters in the binary 164Dy-168Er mixture. The possibility of forming robust dipolar solitons of very large number of atoms make them of great experimental interest. The existence of the solitons is illustrated in terms of stability phase diagrams. Exotic shapes of these solitons are illustrated in isodensity plots. The variational results for statics (size and chemical potential) and dynamics (small oscillation) of the binary soliton compare well with the numerical results. A way of preparing and studying these solitons in laboratory is suggested.

preprint2013arXiv

Stability and collapse of fermions in a binary dipolar boson-fermion 164Dy-161Dy mixture

We suggest a time-dependent mean-field hydrodynamic model for a binary dipolar boson-fermion mixture to study the stability and collapse of fermions in the $^{164}$Dy-$^{161}$Dy mixture. The condition of stability of the dipolar mixture is illustrated in terms of phase diagrams. A collapse is induced in a disk-shaped stable binary mixture by jumping the interspecies contact interaction from repulsive to attractive by the Feshbach resonance technique. The subsequent dynamics is studied by solving the time-dependent mean-field model including three-body loss due to molecule formation in boson-fermion and boson-boson channels. Collapse and fragmentation in the fermions after subsequent explosions are illustrated. The anisotropic dipolar interaction leads to anisotropic fermionic density distribution during collapse. The present study is carried out in three-dimensional space using realistic values of dipolar and contact interactions.

preprint2013arXiv

Stability of trapped degenerate dipolar Bose and Fermi gases

Trapped degenerate dipolar Bose and Fermi gases of cylindrical symmetry with the polarization vector along the symmetry axis are only stable for the strength of dipolar interaction below a critical value. In the case of bosons, the stability of such a dipolar Bose-Einstein condensate (BEC) is investigated for different strengths of contact and dipolar interactions using variational approximation and numerical solution of a mean-field model. In the disk shape, with the polarization vector perpendicular to the plane of the disk, the atoms experience an overall dipolar repulsion and this fact should contribute to the stability. However, a complete numerical solution of the dynamics leads to the collapse of a strongly disk-shaped dipolar BEC due to the long-range anisotropic dipolar interaction. In the case of fermions, the stability of a trapped single-component degenerate dipolar Fermi gas is studied including the Hartree-Fock exchange and Brueckner-Goldstone correlation energies in the local density approximation valid for a large number of atoms. Estimates for the maximum allowed number of polar Bose and Fermi molecules in BEC and degenerate Fermi gas are given.

preprint2012arXiv

A comparative study of the low energy HD+o-/p-H2 rotational excitation/de-excitation collisions and elastic scattering

The Diep and Johnson (DJ) H$_2$-H$_2$ potential energy surface (PES) obtained from the first principles [P. Diep, K. Johnson, J. Chem. Phys. 113, 3480 (2000); 114, 222 (2000)], has been adjusted through appropriate rotation of the three-dimensional coordinate system and applied to low-temperature ($T<300$ K) HD+$o$-/$p$-H$_2$ collisions of astrophysical interest. A non-reactive quantum mechanical close-coupling method is used to carry out the computation for the total rotational state-to-state cross sections $σ_{j_1j_2\rightarrow j'_1j'_2}(ε)$ and corresponding thermal rate coefficients $k_{j_1j_2\rightarrow j'_1j'_2}(T)$. A rather satisfactory agreement has been obtained between our results computed with the modified DJ PES and with the newer H$_4$ PES [A.I. Boothroyd, P.G. Martin, W.J. Keogh, M.J. Peterson, J. Chem. Phys. 116, 666 (2002)], which is also applied in this work. A comparative study with previous results is presented and discussed. Significant differences have been obtained for few specific rotational transitions in the H$_2$/HD molecules between our results and previous calculations. The low temperature data for $k_{j_1j_2\rightarrow j'_1j'_2}(T)$ calculated in this work can be used in a future application such as a new computation of the HD cooling function of primordial gas, which is important in the astrophysics of the early Universe.

preprint2012arXiv

A dipolar droplet bound in a trapped Bose-Einstein condensate

We study the statics and dynamics of a dipolar Bose-Einstein condensate (BEC) droplet bound by inter-species contact interaction in a trapped non-dipolar BEC. Our findings are demonstrated in terms of stability plots of a dipolar 164Dy droplet bound in a trapped non-dipolar 87Rb BEC with a variable number of 164Dy atoms and the inter-species scattering length. A trapped non-dipolar BEC of a fixed number of atoms can only bind a dipolar droplet containing atoms less than a critical number for the inter-species scattering length between two critical values. The shape and size (statics) as well as the small breathing oscillation (dynamics) of the dipolar BEC droplet are studied using the numerical and variational solutions of a mean-field model. We also suggest an experimental procedure for achieving such a 164Dy droplet by relaxing the trap on the 164Dy BEC in a trapped binary 87Rb-164Dy mixture.

preprint2012arXiv

Dipolar Bose-Einstein condensate for large scattering length

A uniform dilute Bose gas of known density has a universal behavior as the atomic scattering length tends to infinity at unitarity while most of its properties are determined by a universal parameter $ξ$ relating the energies of the noninteracting and unitary gases. The usual mean-field equation is not valid in this limit and beyond mean-field corrections become important. We use a dynamical model including such corrections to investigate a trapped disk-shaped dipolar Bose-Einstein condensate (BEC) and a dipolar BEC vortex for large scattering length. We study the sensitivity of our results on the parameter $ξ$ and discuss the possibility of extracting the value of this parameter from experimental observables.

preprint2012arXiv

Dipolar Bose-Einstein condensate in a ring or in a shell

We study properties of a trapped dipolar Bose-Einstein condensate (BEC) in a circular ring or a spherical shell using the mean-field Gross-Pitaevskii equation. In the case of the ring-shaped trap we consider different orientations of the ring with respect to the polarization direction of the dipoles. In the presence of long-range anisotropic dipolar and short-range contact interactions, the anisotropic density distribution of the dipolar BEC in both traps is discussed in detail. The stability condition of the dipolar BEC in both traps is illustrated in phase plot of dipolar and contact interactions. We also study and discuss the properties of a vortex dipolar BEC in these traps.

preprint2012arXiv

Mixing, demixing, and structure formation in a binary dipolar Bose-Einstein condensate

We study static properties of disk-shaped binary dipolar Bose-Einstein condensates of $^{168}$Er-$^{164}$Dy and $^{52}$Cr-$^{164}$Dy mixtures under the action of inter- and intra-species contact and dipolar interactions and demonstrate the effect of dipolar interaction using the mean-field approach. Throughout this study we use realistic values of inter- and intra-species dipolar interactions and the intra-species scattering lengths and consider the inter-species scattering length as a parameter. The stability of the binary mixture is illustrated through phase plots involving number of atoms of the species. The binary system always becomes unstable as the number of atoms increases beyond a certain limit. As the inter-species scattering length increases corresponding to more repulsion, an overlapping mixed state of the two species changes to a separated demixed configuration. During transition from a mixed to a demixed configuration as the inter-species scattering length is increased for parameters just below the stability line, the binary condensate shows special structures in density in the form of red-blood-cell-like biconcave and Saturn-ring-like shapes, which are direct manifestations of dipolar interaction.

preprint2012arXiv

Two-dimensional dipolar Bose-Einstein condensate bright and vortex solitons on one-dimensional optical lattice

By solving the three-dimensional Gross-Pitaevskii equation we generate two-dimensional axially-symmetric and anisotropic dipolar Bose-Einstein condensate bright solitons, for repulsive atomic interaction, stabilized by only a weak one-dimensional optical lattice (OL) aligned along and perpendicular, respectively, to the dipole polarization direction. In the former case vortex solitons can also be created. We show that it is possible to make a stable array of small interacting axially-symmetric dipolar solitons put on alternate OL sites. Further, we demonstrate the elastic nature of the collision of two such solitons.

preprint2012arXiv

Ultracold collisions between two light indistinguishable diatomic molecules: elastic and rotational energy transfer in HD+HD

A close coupling quantum-mechanical calculation is performed for rotational energy transfer in a HD+HD collision at very low energy, down to the ultracold temperatures: $T \sim 10^{-8}$ K. A global six-dimensional H$_2$-H$_2$ potential energy surface is adopted from a previous work [Boothroyd {\it et al.}, J. Chem. Phys., {\bf 116}, 666 (2002).] State-resolved integral cross sections $σ_{ij\rightarrow i'j'}(\varepsilon_{kin})$ of different quantum-mechanical rotational transitions $ij\rightarrow i'j'$ in the HD molecules and corresponding state-resolved thermal rate coefficients $k_{ij\rightarrow i'j'}(T)$ have been computed. Additionally, for comparison, H$_2$+H$_2$ calculations for a few selected rotational transitions have also been performed. The hydrogen and deuterated hydrogen molecules are treated as rigid rotors in this work. A pronounced isotope effect is identified in the cross sections of these collisions at low and ultracold temperatures.

preprint2011arXiv

Anisotropic sound and shock waves in dipolar Bose-Einstein condensate

We study the propagation of anisotropic sound and shock waves in dipolar Bose-Einstein condensate in three dimensions (3D) as well as in quasi-two (2D, disk shape) and quasi-one (1D, cigar shape) dimensions using the mean-field approach. In 3D, the propagation of sound and shock waves are distinct in directions parallel and perpendicular to dipole axis with the appearance of instability above a critical value corresponding to attraction. Similar instability appears in 1D and not in 2D. The numerical anisotropic Mach angle agrees with theoretical prediction. The numerical sound velocity in all cases agrees with that calculated from Bogoliubov theory. A movie of the anisotropic wave propagation in a dipolar condensate is made available as supplementary material.

preprint2011arXiv

Dynamics of gap solitons in a dipolar Bose-Einstein condensate on a three-dimensional optical lattice

We suggest and study the stable disk- and cigar-shaped gap solitons of a dipolar Bose-Einstein condensate of $^{52}$Cr atoms localized in the lowest band gap by three optical-lattice (OL) potentials along orthogonal directions. The one-dimensional version of these solitons of experimental interest confined by an OL along the dipole moment direction and harmonic traps in transverse directions is also considered. Important dynamics of (i) breathing oscillation of a gap soliton upon perturbation and (ii) dragging of a gap soliton by a moving lattice along axial $z$ direction demonstrates the stability of gap solitons. A movie clip of dragging of three-dimensional gap soliton is included.

preprint2011arXiv

Dynamics of quasi-one-dimensional bright and vortex solitons of a dipolar Bose-Einstein condensate with repulsive atomic interaction

By numerical and variational analysis of the three-dimensional Gross-Pitaevskii equation we study the formation and dynamics of bright and vortex-bright solitons in a cigar-shaped dipolar Bose-Einstein condensate for large repulsive atomic interactions. Phase diagram showing the region of stability of the solitons is obtained. We also study the dynamics of breathing oscillation of the solitons as well as the collision dynamics of two solitons at large velocities. Two solitons placed side-by-side at rest coalesce to form a stable bound soliton molecule due to dipolar attraction.

preprint2011arXiv

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

We develop simple numerical algorithms for both stationary and non-stationary solutions of the time-dependent Gross-Pitaevskii (GP) equation describing the properties of Bose-Einstein condensates at ultra low temperatures. In particular, we consider algorithms 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 linear, two-dimensional circularly symmetric, and the three-dimensional spherically-symmetric 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. Numerical results for the chemical potential and root-mean-square size of stationary states are reported using imaginary-time propagation programs for all the cases and compared with previously obtained results. Also presented are numerical results of non-stationary oscillation for different trap symmetries using real-time propagation programs. A set of convenient working codes developed in Fortran 77 are also provided for all these cases (twelve programs in all). In the case of two or three space variables, {Fortran 90/95 versions provide some simplification over the Fortran 77 programs}, and these programs are also included (six programs in all).

preprint2011arXiv

Localization of a Bose-Fermi mixture in a bichromatic optical lattice

We study the localization of a cigar-shaped super-fluid Bose-Fermi mixture in a quasi-periodic bichromatic optical lattice (OL) for inter-species attraction and intra-species repulsion. The mixture is described by the Gross-Pitaevskii equation for the bosons, coupled to a hydrodynamic mean-field equation for fermions at unitarity. We confirm the existence of the symbiotic localized states in the Bose-Fermi mixture and Anderson localization of the Bose component in the interacting Bose-Fermi mixture on a bichromatic OL. The phase diagram in boson and fermion numbers showing the regions of the symbiotic and Anderson localization of the Bose component is presented. Finally, the stability of symbiotic and Anderson localized states is established under small perturbations.

preprint2011arXiv

Localization of collisionally inhomogeneous condensates in a bichromatic optical lattice

By direct numerical simulation and variational solution of the Gross-Pitaevskii equation, we studied the stationary and dynamic characteristics of a cigar-shaped, localized, collisionally inhomogeneous Bose-Einstein condensate trapped in a one-dimensional bichromatic quasi-periodic optical-lattice potential, as used in a recent experiment on the localization of a Bose-Einstein condensate [Roati et al., Nature (London) {\bf 453}, 895 (2008)]. The effective potential characterizing the spatially modulated nonlinearity is obtained. It is found that the collisional inhomogeneity has influence not only on the central region but also on the tail of the Bose-Einstein condensate. The influence depends on the sign and value of the spatially modulated nonlinearity coefficient. We also demonstrate the stability of the stationary localized stat$ performing a standard linear stability analysis. Where possible, the numerical results are shown to be in good agreement with the variational results.

preprint2011arXiv

Matter-wave localization in a weakly perturbed optical lattice

By numerical solution and variational approximation of the Gross-Pitaevskii equation, we studied the localization of a noninteracting and weakly-interacting Bose-Einstein condensate in a weakly perturbed optical lattice in one and three dimensions. The perturbation achieved through a weak delocalizing expulsive or a linear potential as well as a weak localizing harmonic potential removes the periodicity of the optical lattice and leads to localization. We also studied some dynamics of the localized state confirming its stability.

preprint2011arXiv

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

We suggest a simple Gaussian Lagrangian variational scheme for the reduced time-dependent quasi-one- and quasi-two-dimensional Gross-Pitaevskii (GP) equations of a dipolar Bose-Einstein condensate (BEC) in cigar and disk configurations, respectively. The variational approximation for stationary states and breathing oscillation dynamics in reduced dimensions agrees well with the numerical solution of the GP equation even for moderately large short-range and dipolar nonlinearities. The Lagrangian variational scheme also provides much physical insight about soliton formation in dipolar BEC.

preprint2011arXiv

Self-trapping of a binary Bose-Einstein condensate induced by interspecies interaction

The problem of self-trapping of a Bose-Einstein condensate (BEC) and a binary BEC in an optical lattice (OL) and double well (DW) is studied using the mean-field Gross-Pitaevskii equation. For both DW and OL, permanent self-trapping occurs in a window of the repulsive nonlinearity $g$ of the GP equation: $g_{c1}<g<g_{c2}$. In case of OL, the critical nonlinearities $g_{c1}$ and $g_{c2}$ correspond to a window of chemical potentials $μ_{c1}<μ<μ_{c2}$ defining the band gap(s) of the periodic OL. The permanent self-trapped BEC in an OL usually represents a breathing oscillation of a stable stationary gap soliton. The permanent self-trapped BEC in a DW, on the other hand, is a dynamically stabilized state without any stationary counterpart. For a binary BEC with intraspecies nonlinearities outside this window of nonlinearity, a permanent self trapping can be induced by tuning the interspecies interaction such that the effective nonlinearities of the components fall in the above window.

preprint2010arXiv

BCS-BEC crossover in a trapped Fermi super-fluid using a density-functional equation

We derive a generalized time-dependent Galilean-invariant density-functional (DF) equation appropriate to study the stationary and non-stationary properties of a trapped Fermi super-fluid in the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensation (BCS) crossover. This equation is equivalent to a quantum hydrodynamical equation for a Fermi super-fluid. The bulk chemical potential of this equation has the proper (model-independent) dependence on the Fermi-Fermi scattering length in the BCS and BEC limits. We apply this DF equation to the study of stationary density profile and size of a cigar-shaped Fermi super-fluid of $^6$Li atoms and the results are in good agreement with the experiment of Bartenstein {\it et al.} in the BCS-BEC crossover. We also apply the DF equation to the study of axial and radial breathing oscillation and our results for these frequencies are in good agreement with experiments in the BCS-BEC crossover.

preprint2010arXiv

Dimensional reduction of a binary Bose-Einstein condensate in mixed dimensions

We present effective reduced equations for the study of a binary Bose-Einstein condensate (BEC), where the confining potentials of the two BEC components have distinct asymmetry so that the components belong to different space dimensions as in a recent experiment [G. Lamporesi et al., Phys. Rev. Lett. 104, 153202 (2010)]. Starting from a binary three-dimensional (3D) Gross-Pitaevskii equation (GPE) and using a Lagrangian variational approach we derive a binary effective nonlinear Schroedinger equation with components in different reduced dimensions, e. g., the first component in one dimension and the second in two dimensions as appropriate to represent a cigar-shaped BEC coupled to a disk-shaped BEC. We demonstrate that the effective reduced binary equation, which depend on the geometry of the system, is quite reliable when compared with the binary 3D GPE and can be efficiently used to perform numerical simulation and analytical calculation for the investigation of static and dynamic properties of a binary BEC in mixed dimensions.

preprint2010arXiv

Localization of a Bose-Einstein condensate vortex in a bichromatic optical lattice

By numerical simulation of the time-dependent Gross-Pitaevskii equation we show that a weakly interacting or noninteracting Bose-Einstein condensate (BEC) vortex can be localized in a three-dimensional bichromatic quasi-periodic optical-lattice (OL) potential generated by the superposition of two standing-wave polarized laser beams with incommensurate wavelengths. This is a generalization of the localization of a BEC in a one-dimensional bichromatic OL as studied in a recent experiment [Roati et al., Nature 453, 895 (2008)]. We demonstrate the stability of the localized state by considering its time evolution in the form of a stable breathing oscillation in a slightly altered potential for a large period of time. {Finally, we consider the localization of a BEC in a random 1D potential in the form of several identical repulsive spikes arbitrarily distributed in space.

preprint2010arXiv

Localization of a dipolar Bose-Einstein condensate in a bichromatic optical lattice

By numerical simulation and variational analysis of the Gross-Pitaevskii equation we study the localization, with an exponential tail, of a dipolar Bose-Einstein condensate (DBEC) of $^{52}$Cr atoms in a three-dimensional bichromatic optical-lattice (OL) generated by two monochromatic OL of incommensurate wavelengths along three orthogonal directions. For a fixed dipole-dipole interaction, a localized state of a small number of atoms ($\sim 1000$) could be obtained when the short-range interaction is not too attractive or not too repulsive. A phase diagram showing the region of stability of a DBEC with short-range interaction and dipole-dipole interaction is given.

preprint2010arXiv

Matter-wave localization in a random potential

By numerical and variational solution of the Gross-Pitaevskii equation, we studied the localization of a noninteracting and weakly-interacting Bose-Einstein condensate (BEC) in a disordered cold atom lattice and a speckle potential. In the case of a single BEC fragment, the variational analysis produced good results. For a weakly disordered potential, the localized BECs are found to have an exponential tail as in weak Anderson localization. We also investigated the expansion of a noninteracting BEC in these potential. We find that the BEC will be locked in an appropriate localized state after an initial expansion and will execute breathing oscillation around a mean shape when a BEC at equilibrium in a harmonic trap is suddenly released into a disorder potential.

preprint2010arXiv

Spatially-antisymmetric localization of matter wave in a bichromatic optical lattice

By direct numerical simulation of the time-dependent Gross-Pitaevskii equation using the split-step Fourier spectral method we study the double-humped localization of a cigar-shaped Bose-Einstein condensate (BEC) in a one-dimensional bichromatic quasi-periodic optical-lattice potential, as used in a recent experiment on the localization of a BEC [Roati et al., Nature 453, 895 (2008)]. Such states are spatially antisymmetric and are excited modes of Anderson localization. Where possible, we have compared the numerical results with a variational analysis. We also demonstrate the stability of the localized double-humped BEC states under small perturbation.

preprint2010arXiv

Spontaneous symmetry breaking of Bose-Fermi mixtures in double-well potentials

We study the spontaneous symmetry breaking (SSB) of a superfluid Bose-Fermi (BF) mixture in a double-well potential (DWP). The mixture is described by the Gross-Pitaevskii equation (GPE) for the bosons, coupled to an equation for the order parameter of the Fermi superfluid, which is derived from the respective density functional in the unitarity limit (a similar model applies to the BCS regime too). Straightforward SSB in the degenerate Fermi gas loaded into a DWP is impossible, as it requires an attractive self-interaction, while the intrinsic nonlinearity in the Fermi gas is repulsive. Nonetheless, we demonstrate that the symmetry breaking is possible in the mixture with attraction between fermions and bosons, like 40K and 87Rb. Numerical results are represented by dependencies of asymmetry parameters for both components on particle numbers of the mixture, N_F and N_B, and by phase diagrams in the (N_F,N_B) plane, which displays regions of symmetric and asymmetric ground states. The dynamical picture of the SSB, induced by a gradual transformation of the single-well potential into the DWP, is reported too. An analytical approximation is proposed for the case when GPE for the boson wave function may be treated by means of the Thomas-Fermi (TF) approximation. Under a special linear relation between N_F and N_B, the TF approximation allows us to reduce the model to a single equation for the fermionic function, which includes competing repulsive and attractive nonlinear terms. The latter one directly displays the mechanism of the generation of the effective attraction in the Fermi superfluid, mediated by the bosonic component of the mixture.

preprint2010arXiv

Symmetry breaking in a localized interacting binary BEC in a bi-chromatic optical lattice

By direct numerical simulation of the time-dependent Gross-Pitaevskii equation using the split-step Fourier spectral method we study different aspects of the localization of a cigar-shaped interacting binary (two-component) Bose-Einstein condensate (BEC) in a one-dimensional bi-chromatic quasi-periodic optical-lattice potential, as used in a recent experiment on the localization of a BEC [Roati et al., Nature 453, 895 (2008)]. We consider two types of localized states: (i) when both localized components have a maximum of density at the origin x=0, and (ii) when the first component has a maximum of density and the second a minimum of density at x=0. In the non-interacting case the density profiles are symmetric around x=0. We numerically study the breakdown of this symmetry due to inter-species and intra-species interaction acting on the two components. Where possible, we have compared the numerical results with a time-dependent variational analysis. We also demonstrate the stability of the localized symmetry-broken BEC states under small perturbation.

preprint2009arXiv

Localization of a Bose-Einstein condensate in a bichromatic optical lattice

By direct numerical simulation of the time-dependent Gross-Pitaevskii equation we study different aspects of the localization of a non-interacting ideal Bose-Einstein condensate (BEC) in a one-dimensional bichromatic quasi-periodic optical-lattice potential. Such a quasi-periodic potential, used in a recent experiment on the localization of a BEC [Roati et al., Nature 453, 895 (2008)], can be formed by the superposition of two standing-wave polarized laser beams with different wavelengths. We investigate the effect of the variation of optical amplitudes and wavelengths on the localization of a non-interacting BEC. We also simulate the non-linear dynamics when a harmonically trapped BEC is suddenly released into a quasi-periodic potential, {as done experimentally in a laser speckle potential [Billy et al., Nature 453, 891 (2008)]$ We finally study the destruction of the localization in an interacting BEC due to the repulsion generated by a positive scattering length between the bosonic atoms.

preprint2009arXiv

Self-trapping of a Fermi super-fluid in a double-well potential in the BEC-unitarity crossover

We derive a generalized Gross-Pitaevskii density-functional equation appropriate to study the Bose-Einstein condensate (BEC) of dimers formed of singlet spin-half Fermi pairs in the BEC-unitarity crossover while the dimer-dimer scattering length $a$ changes from 0 to $\infty$. Using an effective one-dimensional form of this equation, we study the phenomenon of dynamical self-trapping of a cigar-shaped Fermi super-fluid in the entire BEC-unitarity crossover in a double-well potential. A simple two-mode model is constructed to provide analytical insights. We also discuss the consequence of our study on the self-trapping of an atomic BEC in a double-well potential.

preprint2008arXiv

Gap solitons in a model of a superfluid fermion gas in optical lattices

We consider a dynamical model for a Fermi gas in the Bardeen-Cooper-Schrieffer (BCS) superfluid state, trapped in a combination of a 1D or 2D optical lattice (OL) and a tight parabolic potential acting in the transverse direction(s). The model is based on an equation for the order parameter (wave function), which is derived from the energy density for the weakly coupled BCS superfluid. The equation includes a nonlinear self-repulsive term of power 7/3, which accounts for the Fermi pressure. Reducing the equation to the 1D or 2D form, we construct families of stable 1D and 2D gap solitons (GSs) by means of numerical simulations, which are guided by the variational approximation (VA). The GSs are, chiefly, compact objects trapped in a single cell of the OL potential. In the linear limit, the VA predicts almost exact positions of narrow Bloch bands that separate the semi-infinite and first finite gaps, as well as the first and second finite ones. Families of stable even and odd bound states of 1D GSs are constructed too. We also demonstrate that the GS can be dragged without much distortion by an OL moving at a moderate velocity ($\sim $ 1 mm/s, in physical units). The predicted GSs contain $\sim 10^{3}-10^{4}$ and $\sim 10^{3}$ atoms per 1D and 2D settings, respectively.

preprint2000arXiv

Cooper pair dispersion relation for weak to strong coupling

Cooper pairing in two dimensions is analyzed with a set of renormalized equations to determine its binding energy for any fermion number density and all coupling assuming a generic pairwise residual interfermion interaction. \ Also considered are Cooper pairs (CPs) with nonzero center-of-mass momentum (CMM)--usually neglected in BCS theory--and their binding energy is expanded analytically in powers of the CMM up to quadratic terms. A Fermi-sea-dependent {\it linear} term in the CMM dominates the pair excitation energy in weak coupling (also called the BCS regime) while the more familiar quadratic term prevails in strong coupling (the Bose regime). The crossover, though strictly unrelated to BCS theory {\it per se,} is studied numerically as it is expected to play a central role in a model of superconductivity as a Bose-Einstein condensation of CPs where the transition temperature vanishes for all dimensionality $d\leq 2$ for quadratic dispersion, but is {\it nonzero} for all $d\geq 1$ for linear dispersion.

preprint2000arXiv

Superconductivity as a Bose-Einstein condensation?

Bose-Einstein condensation (BEC) in two dimensions (2D) (e.g., to describe the quasi-2D cuprates) is suggested as the possible mechanism widely believed to underlie superconductivity in general. A crucial role is played by nonzero center-of-mass momentum Cooper pairs (CPs) usually neglected in BCS theory. Also vital is the unique {\it linear} dispersion relation appropriate to weakly-coupled "bosonic" CPs moving in the Fermi sea--rather than in vacuum where the dispersion would be quadratic but only for very strong coupling, and for which BEC is known to be impossible in 2D.

preprint1997arXiv

Application of Renormalization to Potential Scattering

A recently proposed renormalization scheme can be used to deal with nonrelativistic potential scattering exhibiting ultraviolet divergence in momentum space. A numerical application of this scheme is made in the case of potential scattering with $r^{-2}$ divergence for small r, common in molecular and nuclear physics, by the use of cut-offs in momentum and configuration spaces. The cut-off is finally removed in terms of a physical observable and model-independent result is obtained at low energies. The expected variation of the off-shell behavior of the t matrix arising from the renormalization scheme is also discussed.