Numerical solution of the two-dimensional Gross-Pitaevskii equation for trapped interacting atoms

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

The complex sine-Gordon equation as a symmetry flow of the AKNS hierarchy

It is shown how the complex sine-Gordon equation arises as a symmetry flow of the AKNS hierarchy. The AKNS hierarchy is extended by the 'negative' symmetry flows forming the Borel loop algebra. The complex sine-Gordon and the vector nonlinear Schrodinger equations appear as lowest-negative and second-positive flows within the extended hierarchy. This is fully analogous to the well known connection between the sine-Gordon and mKdV equations within the extended mKdV hierarchy. A general formalism for a Toda-like symmetry occupying the 'negative' sector of the sl(N) constrained KP hierarchy and giving rise to the negative Borel sl(N) loop algebra is indicated.

Lie symmetry analysis and reductions of a two-dimensional integrable generalization of the Camassa-Holm equation

In this Letter we investigate Lie symmetries of a (2 + 1)-dimensional integrable generalization of the Camassa-Holm (CH) equation. Through the similarity reductions we obtain four different (1 + 1)-dimensional systems of partial differential equations in which one of them turns out to be a (1 + 1)-dimensional CH equation. We establish their integrability by providing the Lax pair for all of them. Further, we present a brief analysis for some types of particular solutions which include the cuspon, peakon and soliton solutions for the two-dimensional generalization of the CH equation. (C) 2000 Published by Elsevier B.V. B.V.

Stability of the trapped nonconservative Gross-Pitaevskii equation with attractive two-body interaction

The dynamics of a nonconservative Gross-Pitaevskii equation for trapped atomic systems with attractive two-body interaction is numerically investigated, considering wide variations of the nonconservative parameters, related to atomic feeding and dissipation. We study the possible limitations of the mean-field description for an atomic condensate with attractive two-body interaction, by defining the parameter regions, where stable or unstable formation can be found. The present study is useful and timely considering the possibility of large variations of attractive two-body scattering lengths, which may be feasible in recent experiments.

Classical integrable super sinh-Gordon equation with defects

The introduction of defects is discussed under the Lagrangian formalism and Backlund transformations for the N = 1 super sinh-Gordon model. Modified conserved momentum and energy are constructed for this case. Some explicit examples of different Backlund soliton solutions are discussed. The Lax formulation within the space split by the defect leads to the integrability of the model and henceforth to the existence of an infinite number of constants of motion.

Symmetry analysis of an integrable reaction-diffusion equation

In this paper, we investigate the invariance and integrability properties of an integrable two-component reaction-diffusion equation. We perform Painleve analysis for both the reaction-diffusion equation modelled by a coupled nonlinear partial differential equations and its general similarity reduced ordinary differential equation and confirm its integrability. Further, we perform Lie symmetry analysis for this model. Interestingly our investigations reveals a rich variety of particular solutions, which have not been reported in the literature, for this model. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

Bright solitons in coupled defocusing NLS equation supported by coupling: Application to Bose-Einstein condensation

We demonstrate the formation of bright solitons in coupled self-defocusing nonlinear Schrodinger (NLS) equation supported by attractive coupling. As an application we use a time-dependent dynamical mean-field model to study the formation of stable bright solitons in two-component repulsive Bose-Einstein condensates (BECs) supported by interspecies attraction in a quasi one-dimensional geometry. When all interactions are repulsive, there cannot be bright solitons. However, bright solitons can be formed in two-component repulsive BECs for a sufficiently attractive interspecies interaction, which induces an attractive effective interaction among bosons of same type. (c) 2005 Elsevier B.V. All rights reserved.

An exact equation for the free surface of a fluid in a porous medium

We study the problem of the evolution of the free surface of a fluid in a saturated porous medium, bounded from below by a. at impermeable bottom, and described by the Laplace equation with moving-boundary conditions. By making use of a convenient conformal transformation, we show that the solution to this problem is equivalent to the solution of the Laplace equation on a fixed domain, with new variable coefficients, the boundary conditions. We use a kernel of the Laplace equation which allows us to write the Dirichlet-to-Neumann operator, and in this way we are able to find an exact differential-integral equation for the evolution of the free surface in one space dimension. Although not amenable to direct analytical solutions, this equation turns out to allow an easy numerical implementation. We give an explicit illustrative case at the end of the article.

Soliton-cuspon interaction for the Camassa-Holm equation

The interaction of different kinds of solitary waves of the Camassa-Holm equation is investigated. We consider soliton-soliton, soliton-cuspon and cuspon-cuspon interactions. The description of these solutions had previously been shown to be reducible to the solution of an algebraic equation. Here we give explicit examples, numerically solving these algebraic equations and plotting the corresponding solutions. Further, we show that the interaction is elastic and leads to a shift in the position of the solitons or cuspons. We give the analytical expressions for this shift and represent graphically the coupled soliton-cuspon, soliton-soliton and cuspon-cuspon interactions.

Short-wave instabilities in the Benjamin-Bona-Mahoney-Peregrine equation: theory and numerics

In this paper we discuss the nonlinear propagation of waves of short wavelength in dispersive systems. We propose a family of equations that is likely to describe the asymptotic behaviour of a large class of systems. We then restrict our attention to the analysis of the simplest nonlinear short-wave dynamics given by U-0 xi tau, = U-0 - 3(U-0)(2). We integrate numerically this equation for periodic and non-periodic boundary conditions, and we find that short waves may exist only if the amplitude of the initial profile is not too large.

Nonlinear Schrodinger equation with power-law confining potential

Critical limits of a stationary nonlinear three-dimensional Schrodinger equation with confining power-law potentials (similar to r(alpha)) are obtained using spherical symmetry. When the nonlinearity is given by an attractive two-body interaction (negative cubic term), it is shown how the maximum number of particles N-c in the trap increases as alpha decreases. With a negative cubic and positive quintic terms we study a first order phase transition, that occurs if the strength g(3) of the quintic term is less than a critical value g(3c). At the phase transition, the behavior of g(3c) with respect to alpha is given by g(3c)similar to 0.0036+0.0251/alpha+0.0088/alpha(2).

Stability analysis of the D-dimensional nonlinear Schrodinger equation with trap and two- and three-body interactions

Considering the static solutions of the D-dimensional nonlinear Schrodinger equation with trap and attractive two-body interactions, the existence of stable solutions is limited to a maximum critical number of particles, when D greater than or equal to 2. In case D = 2, we compare the variational approach with the exact numerical calculations. We show that, the addition of a positive three-body interaction allows stable solutions beyond the critical number. In this case, we also introduce a dynamical analysis of the conditions for the collapse. (C) 2000 Published by Elsevier B.V. B.V. All rights reserved.

Darboux-Egoroff metrics, rational Landau-Ginzburg potentials and the Painleve VI equation

We present a class of three-dimensional integrable structures associated with the Darboux-Egoroff metric and classical Euler equations of free rotations of a rigid body. They are obtained as canonical structures of rational Landau-Ginzburg potentials and provide solutions to the Painleve VI equation.

Bose-Einstein condensation dynamics from the numerical solution of the Gross-Pitaevskii equation

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.