Numerical study of the spherically symmetric gross-pitaevskii equation in two space dimensions

Hubbard Model for Atomic Impurities Bound by the Vortex Lattice of a Rotating Bose-Einstein Condensate

We investigate cold bosonic impurity atoms trapped in a vortex lattice formed by condensed bosons of another species. We describe the dynamics of the impurities by a bosonic Hubbard model containing occupation-dependent parameters to capture the effects of strong impurity-impurity interactions. These include both a repulsive direct interaction and an attractive effective interaction mediated by the Bose-Einstein condensate. The occupation dependence of these two competing interactions drastically affects the Hubbard model phase diagram, including causing the disappearance of some Mott lobes.

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 quasielastic collision at medium velocities. Possibilities of forming such a 2D spatial soliton in the core of a vortical beam are discussed.

Stationary solutions for the nonlinear Schrödinger equation modeling three-dimensional spherical Bose-Einstein condensates in general potentials

Stationary solutions for the cubic nonlinear Schrödinger equation modeling Bose-Einstein condensates (BECs) confined in three spatial dimensions by general forms of a potential are studied through a perturbation method and also numerically. Note that we study both repulsive and attractive BECs under similar frameworks in order to deduce the effects of the potentials in each case. After outlining the general framework, solutions for a collection of specific confining potentials of physical relevance to experiments on BECs are provided in order to demonstrate the approach. We make several observations regarding the influence of the particular potentials on the behavior of the BECs in these cases, comparing and contrasting the qualitative behavior of the attractive and repulsive BECs for potentials of various strengths and forms. Finally, we consider the nonperturbative where the potential or the amplitude of the solutions is large, obtaining various qualitative results. When the kinetic energy term is small (relative to the nonlinearity and the confining potential), we recover the expected Thomas-Fermi approximation for the stationary solutions. Naturally, this also occurs in the large mass limit. Through all of these results, we are able to understand the qualitative behavior of spherical three-dimensional BECs in weak, intermediate, or strong confining potentials.

Spontaneous formation and nonequilibrium dynamics of a soliton-shaped Bose-Einstein condensate in a trap

The Bose-stimulated self-organization of a quasi-two-dimensional nonequilibrium Bose-Einstein condensate in an in-plane potential is proposed. We obtained the solution of the nonlinear, driven-dissipative Gross-Pitaevskii equation for a Bose-Einstein condensate trapped in an external asymmetric parabolic potential within the method of the spectral expansion. We found that, in sharp contrast to previous observations, the condensate can spontaneously acquire a solitonlike shape for spatially homogeneous pumping. This condensate soliton performs oscillatory motion in a parabolic trap and, also, can spontaneously rotate. Stability of the condensate soliton in the spatially asymmetric trap is analyzed. In addition to the nonlinear dynamics of nonequilibrium Bose-Einstein condensates of ultracold atoms, our findings can be applied to the condensates of quantum well excitons and cavity polaritons in semiconductor heterostructure, and to the condensates of photons.

Stationary solutions for the 2+1 nonlinear Schrödinger equation modeling Bose-Einstein condensates in radial potentials

Stationary solutions for the 2+1 cubic nonlinear Schrödinger equation modeling Bose-Einstein condensates (BEC) in a small potential are obtained via a form of perturbation. In particular, perturbations due to small potentials which either confine or repel the BECs are studied, and under arbitrary piecewise continuous potentials, we obtain the general representation for the perturbation theory of radial BEC solutions. Numerical results are also provided for regimes where perturbative results break down (i.e., the large-potential regime). Both repulsive and attractive BECs are considered under this framework. Solutions for many specific confining potentials of physical relevance to experiments on BECs are provided in order to demonstrate the approach. We make several observations regarding the influence of the particular small potentials on the behavior of the BECs.

Stationary solutions for the 1+1 nonlinear Schrödinger equation modeling attractive Bose-Einstein condensates in small potentials

Stationary solutions for the 1+1 cubic nonlinear Schrödinger equation (NLS) modeling attractive Bose-Einstein condensates (BECs) in a small potential are obtained via a form of nonlinear perturbation. The focus here is on perturbations to the bright soliton solutions due to small potentials which either confine or repel the BECs: under arbitrary piecewise continuous potentials, we obtain the general representation for the perturbation theory of the bright solitons. Importantly, we do not need to assume that the nonlinearity is small, as we perform a sort of nonlinear perturbation by allowing the zeroth-order perturbation term to be governed by a nonlinear equation. This is useful, in that it allows us to consider perturbations of bright solitons of arbitrary size. In some cases, exact solutions can be recovered, and these agree with known results from the literature. Several special cases are considered which involve confining potentials of specific relevance to BECs. We make several observations on the influence of the small potentials on the behavior of the perturbed bright solitons. The results demonstrate the difference between perturbed bright solitons in the attractive NLS and those results found in the repulsive NLS for dark solitons, as discussed by Mallory and Van Gorder, [Phys. Rev. E 88 013205 (2013)]. Extension of these results to more spatial dimensions is mentioned.

Stationary solutions for the 1+1 nonlinear Schrödinger equation modeling repulsive Bose-Einstein condensates in small potentials

Stationary solutions for the 1+1 cubic nonlinear Schrödinger equation modeling repulsive Bose-Einstein condensates (BEC) in a small potential are obtained through a form of nonlinear perturbation. In particular, for sufficiently small potentials, we determine the perturbation theory of stationary solutions, by use of an expansion in Jacobi elliptic functions. This idea was explored before in order to obtain exact solutions [Bronski, Carr, Deconinck, and Kutz, Phys. Rev. Lett. 86, 1402 (2001)], where the potential itself was fixed to be a Jacobi elliptic function, thereby reducing the nonlinear ODE into an algebraic equation, (which could be easily solved). However, in the present paper, we outline the perturbation method for completely general potentials, assuming only that such potentials are locally small. We do not need to assume that the nonlinearity is small, as we perform a sort of nonlinear perturbation by allowing the zeroth-order perturbation term to be governed by a nonlinear equation. This allows us to consider even poorly behaved potentials, so long as they are bounded locally. We demonstrate the effectiveness of this approach by considering a number of specific potentials: for the simplest potentials, and we recover results from the literature, while for more complicated potentials, our results are new. Dark soliton solutions are constructed explicitly for some cases, and we obtain the known one-soliton tanh-type solution in the simplest setting for the repulsive BEC. Note that we limit our results to the repulsive case; similar results can be obtained for the attractive BEC case.

Unified approach to split absorbing boundary conditions for nonlinear Schrödinger equations: Two-dimensional case

This paper aims to design local absorbing boundary conditions (LABCs) for the two-dimensional nonlinear Schrödinger equations on a rectangle by extending the unified approach. Based on the time-splitting idea, the main process of the unified approach is to approximate the kinetic energy part by a one-way equation, unite it with the potential energy equation, and then obtain the well-posed and accurate LABCs on the artificial boundaries. In the corners, we use the (1,1)-Padé approximation to the kinetic term and also unite it with the nonlinear term to give some local corner boundary conditions. Numerical tests are given to verify the stable and tractable advantages of the method.

Unified approach to split absorbing boundary conditions for nonlinear Schrödinger equations

An efficient method is proposed for numerical solutions of nonlinear Schrödinger equations on an unbounded domain. Through approximating the kinetic energy term by a one-way equation and uniting it with the potential energy equation, absorbing boundary conditions are designed to truncate the unbounded domain, which are in nonlinear form and can perfectly absorb waves outgoing from the boundaries of the truncated computational domain. The stability of the induced initial boundary value problem defined on the computational domain is examined by a normal mode analysis. Numerical examples are given to illustrate the stable and tractable advantages of the method.

Quantum lattice Boltzmann simulation of expanding Bose-Einstein condensates in random potentials

The phenomenon of Anderson localization in expanding one-dimensional Bose-Einstein condensates is investigated by numerically solving the Gross-Pitaevskii equation with a random speckle potential. To this purpose, a quantum lattice Boltzmann (QLB) method is used, and compared with a standard Crank-Nicolson scheme. The QLB simulations show evidence of Anderson localization even for relatively low-energy condensates, with a healing length as large as one-tenth of the Thomas-Fermi length. Moreover, very long-time simulations, lasting up to 15 000 optical confinement periods, indicate that the Anderson localization degrades in time, although at a very slow pace. In particular, the inverse localization length is found to decay according to a t;{-1/3} law. This lends support to the idea that localized wave functions, although not strictly ground states, represent extremely long-lived metastable states of the expanding condensate.

Ground-state computation of Bose-Einstein condensates by an imaginary-time quantum lattice Boltzmann scheme

The multidimensional formulation of the quantum lattice Boltzmann (qLB) scheme is extended to the case of nonlinear quantum wave equations. More specifically, imaginary-time formulations of the qLB scheme are developed and applied to the numerical computation of the ground state of the Gross-Pitaevskii equation in one and two spatial dimensions. The calculation is validated through detailed comparison with other numerical methods, as well as with analytical results based on the Thomas-Fermi approximation. The linear scaling of the time-step size with the spatial mesh spacing, a distinctive feature of the present quantum kinetic approach, is also numerically demonstrated.

Absorbing boundary conditions for nonlinear Schrödinger equations

A local time-splitting method (LTSM) is developed to design absorbing boundary conditions for numerical solutions of time-dependent nonlinear Schrödinger equations associated with open boundaries. These boundary conditions are significant for numerical simulations of propagations of nonlinear waves in physical applications, such as nonlinear fiber optics and Bose-Einstein condensations. Numerical examples are implemented to demonstrate the attractive features of using the LTSM.

Spectral method for the time-dependent Gross-Pitaevskii equation with a harmonic trap

We study the numerical resolution of the time-dependent Gross-Pitaevskii equation, a nonlinear Schrödinger equation used to simulate the dynamics of Bose-Einstein condensates. Considering condensates trapped in harmonic potentials, we present an efficient algorithm by making use of a spectral-Galerkin method, using a basis set of harmonic-oscillator functions, and the Gauss-Hermite quadrature. We apply this algorithm to the simulation of condensate breathing and scissor modes.

Numerical study of the coupled time-dependent Gross-Pitaevskii equation: application to Bose-Einstein condensation

We present a numerical study of the coupled time-dependent Gross-Pitaevskii equation, which describes the Bose-Einstein condensate of several types of trapped bosons at ultralow temperature with both attractive and repulsive interatomic interactions. The same approach is used to study both stationary and time-evolution problems. We consider up to four types of atoms in the study of stationary problems. We consider the time-evolution problems where the frequencies of the traps or the atomic scattering lengths are suddenly changed in a stable preformed condensate. We also study the effect of periodically varying these frequencies or scattering lengths on a preformed condensate. These changes introduce oscillations in the condensate, which are studied in detail. Good convergence is obtained in all cases studied.