Pitaevskii Equation Research Articles

The Casimir-like effect in a one-dimensional Bose gas

The electromagnetic Casimir effect manifests as the interaction between uncharged conducting objects that are placed in a vacuum. More generally, the Casimir-like effect denotes an induced interaction between external bodies in a fluctuating medium. We study the Casimir-like interaction between two impurities embedded in a weakly interacting one-dimensional Bose gas. We develop a theory based on the Gross–Pitaevskii equation that accounts for the effect of quantum fluctuations. At small separations, the induced interaction between the impurities decays exponentially with the distance. This is a classical result that can be understood using the mean-field Gross–Pitaevskii equation. We find that at larger distances, the induced interaction crosses over into a power law dependence due to the quantum fluctuations. We obtain an analytic expression for the interaction that interpolates between the two limiting behaviors. The obtained result does not require any regularization.

Magnetic solitons in binary mixtures of Bose–Einstein condensates

Solitons, that is stable localized perturbations of a medium, are the topological excitations of nonlinear systems. They can be stable and live for a long time and may have promising applications for telecommunication. The basic one is the Tsuzuki dark soliton, which can be described by an analytical solution of the Gross–Pitaevskii equation (GPE). Ultracold Bose–Einstein condensed (BEC) gases are an important example for the investigation of solitons which can be created by phase and density imprinting. New possibilities arise in mixtures of different hyperspin states of ultra-cold gases, where the so-called magnetic solitons (MS), that is localized magnetized regions, can exist. We will see that these MS permit an analytical description. New peculiar phenomena can take place in the presence of a coherent Rabi coupling between the spin states, where two different type of solitons exist—so-called $$2\pi $$ and $$0\pi $$ solitons. $$2\pi $$ solitons, unlike the usual Tsuzuki solitons, have at small velocity a positive effective mass and consequently do not undergo the snake instability. Solitary waves can oscillate in BEC gases along elongated traps. The theoretical description of this motion requires the knowledge of the effective soliton mass and the effective number of particles in the soliton. These quantities are calculated.

Long wave limit for Schrödinger maps

We study long wave limits for general Schrödinger map systems into Kähler manifolds with a constraining potential vanishing on a Lagrangian submanifold. We obtain KdV-type systems set on the tangent space of the submanifold. Our general theory is applied to study the long wave limits of the Gross–Pitaevskii equation and of the Landau–Lifshitz systems for ferromagnetic and anti-ferromagnetic chains.

Semiclassical analysis with new Galilean transformations for a Gross–Pitaevskii system with nonzero conditions at infinity

Recently, a rich variety of the micro-phenomena of the superfluid passing an obstacle has been observed in the binary mixture of rotating Bose--Einstein condensates (BECs). Among such phenomena, the interaction of dark--bright solitons is one of the most important issues. In this work we investigate the semi-classical limit for a coupled system of Gross--Pitaevskii (GP) equations with rotating fields and trap potentials in a two-dimensional exterior domain, where the superfluid is non-vanishing at infinity. We establish a new Galilean type transformation and follow the argument of the modulated energy functional (a Lyapunov type functional) in \cite{ll08,lz05} to control the propagation of mass densities and linear momenta of the solution via a compressible Euler equation with Coriolis force in a semi-classical regime. Moreover, the effect of the rotating field on the superfluid in the region far away from the obstacle is precisely described.

Vortex-lattice in a uniform Bose–Einstein condensate in a box trap

We study numerically the vortex-lattice formation in a rapidly rotating uniform quasi-two-dimensional Bose–Einstein condensate (BEC) in a box trap. We consider two types of boxes: square and circle. In a square-shaped 2D box trap, when the number of generated vortices is the square of an integer, the vortices are found to be arranged in a perfect square lattice, although deviations near the center are found when the number of generated vortices is arbitrary. In case of a circular box trap, the generated vortices in the rapidly rotating BEC lie on concentric closed orbits. Near the center, these orbits have the shape of polygons, whereas near the periphery the orbits are circles. The circular box trap is equivalent to the rotating cylindrical bucket used in early experiment(s) with liquid He II. The number of generated vortices in both cases is in qualitative agreement with Feynman’s universal estimate. The numerical simulation for this study is performed by a solution of the underlying mean-field Gross–Pitaevskii (GP) equation in the rotating frame, where the wave function for the generated vortex lattice is a stationary state. Consequently, the imaginary-time propagation method can be used for a solution of the GP equation, known to lead to an accurate numerical solution. We also demonstrated the dynamical stability of the vortex lattices in real-time propagation upon a small change of the angular frequency of rotation, using the converged imaginary-time wave function as the initial state.

Hybrid matter-wave - microwave solitons on the lattice

We introduce a two-component system which models a pseudospinor Bose–Einstein condensate (BEC), with a microwave field coupling its two components. The feedback of BEC onto the field (the local-field effect) is taken into account by dint of the respective Poisson equation, which is solved using the Green’s function. This gives rise to an effective long-range self-trapping interaction, which may act alone, or be combined with the contact cubic nonlinearity. The system is made discrete by loading the BEC into a deep optical-lattice potential. Numerical solutions of the discrete system demonstrate that onsite-centered fundamental solitons are stable in the cases of attractive or zero contact interactions, while offsite-centered solitons are unstable. In the case of the repulsive onsite nonlinearity, offsite solitons are stable, while their onsite-centered counterparts are stable only at sufficiently small values of the norm, where bistability between the off- and onsite-centered mode takes place. The shape of the onsite-centered solitons is very accurately predicted by a variational approximation (which includes essential technical novelties), and their super-exponentially decaying tails are found by means of direct analytical consideration. Spatially-antisymmetric (“twisted”) solitons are stable at small values of the norm, being unstable at larger norms. In the strongly asymmetric version of the two-component system, which includes the Zeeman splitting, the system is reduced to a single discrete Gross–Pitaevskii equation, by eliminating the small higher-energy component.

Vortex states in rotating two-component dipolar Bose–Einstein condensates

We consider the stationary state properties of pseudo-spin-1/2 rotating dipolar Bose–Einstein condensates (BECs) by numerical simulations of the Gross–Pitaevskii equation. Different vortex structures in each component are studied, depending on the competition between the dipole–dipole interactions (DDIs) and rotational. We also investigate the differences of vortex number in the two components, showing that anisotropic nature of DDIs plays a significant role in vortices formation process.

Weakly bound solitons and two-soliton molecules in dipolar Bose–Einstein condensates

Analytic expressions have been derived for the interaction potential between dipolar bright solitons and the binding energy of a two-soliton molecule. The properties of these localized structures are explored with a focus on their behavior in the weakly bound regime, with a small binding energy. Using the variational approach a coupled system of ordinary differential equations for the parameters of a soliton molecule is obtained for the description of their evolution. Predictions of the model are compared with numerical simulations of the governing nonlocal Gross–Pitaevskii equation and good qualitative agreement between them is demonstrated.

Tools for designing atom interferometers in a microgravity environment

We present a variational model suitable for rapid preliminary design of atom interferometers in a microgravity environment. The model approximates the solution of the 3D rotating--frame Gross--Pitaevskii equation (GPE) as the sum of Nc Gaussian clouds. Each Gaussian cloud is assumed to have time--dependent center positions, widths, and linear and quadratic phase parameters. We applied the Lagrangian Variational Method (LVM) with this trial wave function to derive equations of motion for these parameters that can be adapted to any external potential. We also present a 1D version of this variational model. As an example we apply the model to a 1D atom interferometry scheme for measuring Newton's gravitational constant, G, in a microgravity environment. We show how the LVM model can (1) constrain the experimental parameter space size, (2) show how the value of G can be obtained from the experimental conditions and interference pattern characteristics, and (3) show how to improve the sensitivity of the measurement and construct a preliminary error budget.

W-shaped solitons in inhomogeneous cigar-shaped Bose–Einstein condensates with repulsive interatomic interactions

We analyze the formation and dynamics of W-shaped solitons in an inhomogeneous cigar-shaped Bose–Einstein condensate with repulsive interatomic interactions in the presence of spatially and temporally nonuniform nonlinearity and external potential. The condensate is assumed to have a large number of particles; hence it may be approximated by an effective one-dimensional equation, which is derived from the full three-dimensional Gross–Pitaevskii equation. New classes of W-shaped soliton solutions, which may propagate in the condensate system with the same hump value, are found for the first time by means of the similarity transformation technique. The functional form of these localized nonlinear structures is different from the usual ones reported in the setting of optical fibers. The dynamic behaviors of the W-shaped soliton pulses under different nonlinearity and external potential management that are of particular interest for applications in Bose–Einstein condensates are analyzed.

The numerical study for the ground and excited states of fractional Bose–Einstein condensates

In this paper, we study the ground and first excited states of the fractional Bose–Einstein condensates (BEC) which is modeled by fractional Gross–Pitaevskii (GP) equation. We first introduce the normalized gradient flow method and prove its energy diminishing property. Then the weighted shifted Grünwald–Letnikov difference (WSGD) method is used to discretize the Gross–Pitaevskii equation. The corresponding normalization and energy diminishing property for the semi-discrete scheme are proved. For the time discretization, we use the implicit integration factor (IIF) method which decouples the diffusion and nonlinear terms separately. Finally the numerical methods are applied to compute the ground and first excited states of fractional BEC with harmonic oscillator, harmonic-plus-optical lattice and box potential. Our numerical results show that the ground and excited states in fractional GP equation differ from those of the standard (non-fractional) GP equation.

A NEW LINEAR AND CONSERVATIVE FINITE DIFFERENCE SCHEME FOR THE GROSS–PITAEVSKII EQUATION WITH ANGULAR MOMENTUM ROTATION

A new linear and conservative finite difference scheme which preserves discrete mass and energy is developed for the two-dimensional Gross–Pitaevskii equation with angular momentum rotation. In addition to the energy estimate method and mathematical induction, we use the lifting technique as well as some well-known inequalities to establish the optimal $H^{1}$ -error estimate for the proposed scheme with no restrictions on the grid ratio. Unlike the existing numerical solutions which are of second-order accuracy at the most, the convergence rate of the numerical solution is proved to be of order $O(h^{4}+\unicode[STIX]{x1D70F}^{2})$ with time step $\unicode[STIX]{x1D70F}$ and mesh size $h$ . Numerical experiments have been carried out to show the efficiency and accuracy of our new method.

C and Fortran OpenMP programs for rotating Bose–Einstein condensates

We present OpenMP versions of C and Fortran programs for solving the Gross–Pitaevskii equation for a rotating trapped Bose–Einstein condensate (BEC) in two (2D) and three (3D) spatial dimensions. The programs can be used to generate vortex lattices and study dynamics of rotating BECs. We use the split-step Crank–Nicolson algorithm for imaginary- and real-time propagation to calculate stationary states and BEC dynamics, respectively. The simulation input parameters for the C programs are provided via input files, while for the Fortran programs they are given at the beginning of each program and therefore their change requires recompilation of the corresponding program. The programs propagate the condensate wave function and calculate several relevant physical quantities, such as the energy, the chemical potential, and the root-mean-square sizes. The imaginary-time propagation starts with an analytic wave function with one vortex at the trap center, modulated by a random phase at different space points. Nevertheless, the converged wave function for a rapidly rotating BEC with a large number of vortices is most efficiently calculated using the pre-calculated converged wave function of a rotating BEC containing a smaller number of vortices as the initial state rather than using an analytic wave function with one vortex as the initial state. These pre-calculated initial states exhibit rapid convergence for fast-rotating condensates to states containing multiple vortices with an appropriate phase structure. This is illustrated here by calculating vortex lattices with up to 61 vortices in 2D and 3D. Outputs of the programs include calculated physical quantities, as well as the wave function and different density profiles (full density, integrated densities in lower dimensions, and density cross-sections). The provided real-time propagation programs can be used to study the dynamics of a rotating BEC using the imaginary-time stationary wave function as the initial state. We also study the efficiency of parallelization of the present OpenMP C and Fortran programs with different compilers. Program summaryProgram title: BEC-GP-ROT-OMP, consisting of: (1) BEC-GP-ROT-OMP-C package, containing programs (i) bec-gp-rot-2d-th and (ii) bec-gp-rot-3d-th; (2) BEC-GP-ROT-OMP-F package, containing programs (i) bec-gp-rot-2d-th and (ii) bec-gp-rot-3d-th.Program files doi:http://dx.doi.org/10.17632/cw7tkn22v2.2Licensing provisions: Apache License 2.0Programming language: OpenMP C; OpenMP Fortran. The C programs are tested with the GNU, Intel, PGI, Oracle, and Clang compiler, and the Fortran programs are tested with the GNU, Intel, PGI, and Oracle compiler.Nature of problem: The present Open Multi-Processing (OpenMP) C and Fortran programs solve the time-dependent nonlinear partial differential Gross–Pitaevskii (GP) equation for a trapped rotating Bose–Einstein condensate in two (2D) and three (3D) spatial dimensions in a fully anisotropic traps.Solution method: We employ the split-step Crank–Nicolson algorithm to discretize the time-dependent GP equation in space and time. The discretized equation is then solved by imaginary- or real-time propagation, employing adequately small space and time steps, to yield the solution of stationary and non-stationary problems, respectively.

Stability analysis and continuation for the coupled Gross–Pitaevskii equations

We study linear stability analysis for the coupled Gross–Pitaevskii equations (CGPEs), which are used as the governing equations for rotating two-component Bose–Einstein condensates (BEC), and (ultrarapidly) rotating two-component Rydberg-dressed BEC. The CGPEs describe the coexistence of superfluidity and crystallization of BEC simultaneously. We show numerically that the linear stability of discrete steady state solutions is related to the angular velocity. Next, we describe some continuation algorithms for computing the ground state solutions of the CGPEs. Specifically, one of them can efficiently handle the case when long-range intra-component interactions of the two components are different. For ultrarapid rotation our numerical experiments demonstrate a competition between the supersolid crystal structure and the rotating-induced vortex lattice which gives rise to new phases.

Parafermionic behavior of Bose-Einstein condensates

Bright solitons of 7Li atoms in a quasi one-dimensional optical trap, formed in a stable Bose–Einstein condensate in which the interactions have been magnetically tuned from repulsive to attractive, have been seen to exhibit repulsive interactions among them when set in motion by offsetting the optical potential. Solving first the Gross–Pitaevskii equation for the special conditions corresponding to the experiment, we show then that this system can be described in terms of generalized parafermionic oscillators, the order of the parafermions being related to the strength of the interaction among the atoms and being a measure of the bosonic behavior vs. the fermionic behavior of the system.

Derivation of the 1d nonlinear Schrödinger equation from the 3d quantum many-body dynamics of strongly confined bosons

We consider the dynamics of N interacting bosons initially exhibiting Bose–Einstein condensation. Due to an external trapping potential, the bosons are strongly confined in two spatial directions, with the transverse extension of the trap being of order ε. The non-negative interaction potential is scaled such that its scattering length is positive and of order (N/ε2)−1 and the range of the interaction scales as (N/ε2)−β for β ∈ (0, 1). We prove that in the simultaneous limit N → ∞ and ε → 0, the condensation is preserved by the dynamics and the time evolution is asymptotically described by a cubic defocusing nonlinear Schrödinger equation in one dimension, where the strength of the nonlinearity depends on the interaction and on the confining potential. This is the first derivation of a lower-dimensional effective evolution equation for singular potentials scaling with β≥12 and lays the foundations for the derivation of the physically relevant one-dimensional Gross–Pitaevskii equation (β = 1). For our analysis, we adapt an approach by Pickl to the problem with strong confinement.

Collision of parallel vortex dipoles in a Bose–Einstein condensate

Understanding how two vortex dipoles interact with each other microscopically is helpful for understanding the dynamics of quantum turbulence in a two-dimensional Bose–Einstein condensate (BEC). We investigate the collision of two vortex dipoles propagating in opposite directions along parallel lines by numerically solving the two-dimensional Gross–Pitaevskii equation. By varying the impact parameter, we identify three general dynamical modes of the vortex dipoles in a homogeneous BEC, which are the vortex recombination mode, the encircling mode and the flyby mode. The existence of these three modes has been corroborated by numerical computations of the ordinary differential equations describing the vortex interactions. Furthermore, we show the corresponding parameter ranges for these modes.

Self-trapped quantum balls in binary Bose–Einstein condensates

We study the formation of a stable self-trapped spherical quantum ball in a binary Bose–Einstein condensate (BEC) with two-body inter-species attraction and intra-species repulsion employing the beyond-mean-field Lee–Huang–Yang and the three-body interactions. We find that either of these interactions or a combination of them can stabilize the binary BEC quantum ball with very similar stationary results, and for a complete description of the problem both the terms should be considered. These interactions lead to higher-order nonlinearities, e.g. quartic and quintic, respectively, in a nonlinear dynamical equation compared to the cubic nonlinearity of the two-body contact interaction in the mean-field Gross–Pitaevskii equation. The higher-order nonlinearity makes the energy infinitely large at the center of the binary ball and thus avoids its collapse. In addition to the formation of stationary binary balls, we also study a collision between two such balls. At large velocities, the collision is found to be elastic, which turns out to be inelastic as the velocity is lowered. We consider the numerical solution of a beyond-mean-field model for the binary ball as well as a single-mode variational approximation to it in this study.

Effective equation for quasi-one dimensional tube-shaped Bose–Einstein condensates

In this letter we derive an effective 1D equation that describes the axial dynamics of a tube-shaped Bose–Einstein condensate. The dimensional reduction is achieved by using a variational approach starting from the 3D Gross–Pitaevskii equation (GPE) in presence of a nonharmonic external potential in the transverse direction and generic in the axial one. The resulting equation is a time-dependent 1D nonpolynomial Schrödinger equation (NPSE). In view to check the accuracy of our 1D-NPSE, we numerically investigated the ground state properties of such a system that are in perfect agreement with the results produced by the 3D-GPE. We also compare the results with those from an 1D cubic nonlinear Schrödinger equation and the Thomas–Fermi approximation. Finally, the dynamics of ground states obtained from our 1D-NPSE is verified numerically by considering a small change in the strength of the axial confining potential.

Bifurcation of nonlinear bound states in the periodic Gross-Pitaevskii equation with 𝒫𝒯-symmetry

AbstractThe stationary Gross–Pitaevskii equation in one dimension is considered with a complex periodic potential satisfying the conditions of the 𝒫𝒯 (parity-time reversal) symmetry. Under rather general assumptions on the potentials, we prove bifurcations of 𝒫𝒯-symmetric nonlinear bound states from the end points of a real interval in the spectrum of the non-selfadjoint linear Schrödinger operator with a complex 𝒫𝒯-symmetric periodic potential. The nonlinear bound states are approximated by the effective amplitude equation, which bears the form of the cubic nonlinear Schrödinger equation. In addition, we provide sufficient conditions for the appearance of complex spectral bands when the complex 𝒫𝒯-symmetric potential has an asymptotically small imaginary part.