Lobatto deferred correction for stiff two-point boundary value problems

An iterated deferred correction algorithm based on Lobatto Runge-Kutta formulae is developed for the efficient numerical solution of nonlinear stiff two-point boundary value problems. An analysis of the stability properties of general deferred correction schemes which are based on implicit Runge-Kutta methods is given and results which are analogous to those obtained for initial value problems are derived. A revised definition of symmetry is presented and this ensures that each deferred correction produces an optimal increase in order. Finally, some numerical results are given to demonstrate the superior performance of Lobatto formulae compared with mono-implicit formulae on stiff two-point boundary value problems. (C) 1998 Elsevier B.V. Ltd. All rights reserved.

Stabilization of a (3+1)-dimensional soliton in a Kerr medium by a rapidly oscillating dispersion coefficient

Using the numerical solution of the nonlinear Schrodinger equation and a variational method it is shown that (3 + 1)-dimensional spatiotemporal optical solitons can be stabilized by a rapidly oscillating dispersion coefficient in a Kerr medium with cubic nonlinearity. This has immediate consequence in generating dispersion-managed robust optical soliton in communication as well as possible stabilized Bose-Einstein condensates in periodic optical-lattice potential via an effective-mass formulation. We also critically compare the present stabilization with that obtained by a rapid sinusoidal oscillation of the Kerr nonlinearity parameter.

Resonance in Bose-Einstein condensate oscillation from a periodic variation in scattering length

Using the explicit numerical solution of the axially symmetric Gross-Pitaevskii equation, we study the oscillation of the Bose-Einstein condensate (BEC) induced by a periodic variation in the atomic scattering length a. When the frequency of oscillation of a is an even multiple of the radial or axial trap frequency, respectively, the radial or axial oscillation of the condensate exhibits resonance with a novel feature. In this nonlinear problem without damping, at resonance in the steady state the amplitude of oscillation passes through a maximum and minimum. Such a growth and decay cycle of the amplitude may keep on repeating. Similar behaviour is also observed in a rotating BEC.

The critical number of atoms in an attractive Bose-Einstein condensate on optical plus harmonic traps

The stability of an attractive Bose-Einstein condensate on a joint one-dimensional optical lattice and an axially symmetrical harmonic trap is studied using the numerical solution of the time-dependent mean-field Gross-Pitaevskii equation and the critical number of atoms for a stable condensate is calculated. We also calculate this critical number of atoms in a double-well potential which is always greater than that in an axially symmetrical harmonic trap. The critical number of atoms in an optical trap can be made smaller or larger than the corresponding number in the absence of the optical trap by moving a node of the optical lattice potential in the axial direction of the harmonic trap. This variation of the critical number of atoms can be observed experimentally and compared with the present calculations.

Positronium atom scattering by H-2 in a coupled-channel framework

The scattering of ortho-positronium (Ps) by H-2 has been investigated using a three-Ps-state (Ps(1s,2s, 2p)H-2(X (1)Sigma(g)(+))) coupled-channel model and using the Born approximation for higher excitations and ionization of Ps and B (1)Sigma(u)(+) and b (3)Sigma(u)(+) excitations of H-2. We employ a recently proposed time-reversal-symmetric non-local electron-exchange model potential. We present a calculational scheme for solving the body-frame fixed-nuclei coupled-channel scattering equations for Ps-H-2, which simplifies the numerical solution technique considerably. Ps ionization is found to have the leading contribution to target-elastic and all target-inelastic processes. The total cross sections at low and medium energies are in good agreement with experiment.

Dynamics of collapsing and exploding Bose-Einstein condensate

Recently, Donley et al. performed an experiment on the dynamics of collapsing and exploding Bose-Einstein condensates by suddenly changing the scattering length of atomic interaction to a large negative value on a preformed repulsive condensate of Rb-85 atoms in an axially symmetric trap. Consequently, the condensate collapses and ejects atoms via explosions, We show that the accurate numerical solution of the time-dependent Gross-Pitaevskii equation with axial symmetry can explain some aspects of the dynamics of the collapsing condensate. (C) 2002 Published by Elsevier B.V. B.V.

Dynamical collapse in a degenerate binary fermion mixture using a hydrodynamic model

We use a time-dependent dynamical hydrodynamic model to study a collapse in a degenerate fermion-fermion mixture ( DFFM) of different atoms. Due to a strong Pauli-blocking repulsion among identical spin-polarized fermions at short distances, there cannot be a collapse for repulsive interspecies fermion fermion interaction. However, there can be a collapse for a sufficiently attractive interspecies fermion-fermion interaction in a DFFM of different atoms. Using a variational analysis and numerical solution of the hydrodynamic model, we study different aspects of collapse in such a DFFM initiated by a jump in the interspecies fermion-fermion interaction ( scattering length) to a large negative ( attractive) value using a Feshbach resonance. Suggestion for experiments of collapse in a DFFM of distinct atoms is made.

Expansion of a Bose-Einstein condensate formed on a joint harmonic and one-dimensional optical-lattice potential

We study the expansion of a Bose-Einstein condensate trapped in a combined optical-lattice and axially-symmetric harmonic potential using the numerical solution of the mean-field Gross-Pitaevskii equation. First, we consider the expansion of such a condensate under the action of the optical-lattice potential alone. In this case the result of numerical simulation for the axial and radial sizes during expansion is in agreement with two experiments by Morsch et al (2002 Phys. Rev. A 66 021601(R) and 2003 Laser Phys. 13 594). Finally, we consider the expansion under the action of the harmonic potential alone. In this case the oscillation, and the disappearance and revival of the resultant interference pattern is in agreement with the experiment by Muller et al (2003 J. Opt. B: Quantum Semiclass. Opt. 5 S38).

Stabilization of a light bullet in a layered Kerr medium with sign-changing nonlinearity

Using the numerical solution of the nonlinear Schrodinger equation and a variational method, it is shown that (3+1)-dimensional spatiotemporal optical solitons, known as light bullets, can be stabilized in a layered Kerr medium with sign-changing nonlinearity along the propagation direction.

Bose-Einstein condensation thermodynamics of a trapped gas with attractive interaction

We study the Bose-Einstein condensation of an interacting gas with attractive interaction confined in a harmonic trap using a semiclassical two-fluid mean-field model. The condensed state is described by the converged numerical solution of the Gross-Pitaevskii equation. By solving the system of coupled equations of this model iteratively we obtain the converged results for the temperature dependencies of the condensate fraction, chemical potential, and internal energy for the Bose-Einstein condensate of Li-7 atoms. (C) 2000 Elsevier B.V. B.V. All rights reserved.

Collapse of attractive Bose-Einstein condensed vortex states in a cylindrical trap

The quantized vortex states of a weakly interacting Bose-Einstein condensate of atoms with attractive interatomic interaction in an axially symmetric harmonic oscillator trap are investigated using the numerical solution of the time-dependent Gross-Pitaevskii equation obtained by the semi-implicit Crank-Nicholson method. The collapse of the condensate is studied in the presence of deformed traps with the larger frequency along either the radial or the axial direction. The critical number of atoms for collapse is calculated as a function of the vortex quantum number L. The critical number increases with increasing angular momentum L of the cortex state but tends to saturate for large L.

Asymptotic soliton train solutions of the defocusing nonlinear Schrodinger equation

Asymptotic behavior of initially large and smooth pulses is investigated at two typical stages of their evolution governed by the defocusing nonlinear Schrodinger equation. At first, wave breaking phenomenon is studied in the limit of small dispersion. A solution of the Whitham modulational equations is found for the case of dissipationless shock wave arising after the wave breaking point. Then, asymptotic soliton trains arising eventually from a large and smooth initial pulse are studied by means of a semiclassical method. The parameter varying along the soliton train is calculated from the generalized Bohr-Sommerfeld quantization rule, so that the distribution of eigenvalues depends on two functions-intensity rho(0)(x) of the initial pulse and its initial chirp v(0)(x). The influence of the initial chirp on the asymptotic state is investigated. Excellent agreement of the numerical solution of the defocusing NLS equation with predictions of the asymptotic theory is found.

Dynamical classical superfluid-insulator transition in a Bose-Einstein condensate on an optical lattice

We predict a dynamical: classical superfluid-insulator transition in a Bose-Einstein condensate (BEC) trapped in combined optical and axially symmetrical harmonic potentials initiated by the periodic modulation of the radial trapping potential. The transition is marked by a loss of phase coherence in the BEC and a subsequent destruction of the interference pattern upon free:expansion. For a weak modulation of the radial potential the phase coherence is maintained. For a stronger modulation and a longer holding time in the modulated trap, the phase coherence is destroyed thus signalling a classical superfluid-insulator transition. The results are illustrated by a complete numerical solution of the axially symmetrical mean-field Gross-Pitaevskii equation for a repulsive BEC. Suggestions for future experimentation are-made.

Josephson oscillation of a superfluid Fermi gas

Using the complete numerical solution of a time-dependent three-dimensional rnean-field model we study the Josephson oscillation of a superfluid Fermi gas (SFG) at zero temperature formed in a combined axially-symmetric harmonic plus one-dimensional periodic optical-lattice (OL) potentials after displacing the harmonic trap along the axial OL axis. We study the dependence of Josephson frequency on the strength of the OL potential. The Josephson frequency decreases with increasing strength as found in the experiment of Cataliotti et al. [Science 293, 843 (2001)] for a Bose-Einstein condensate and of the experiment of Pezze et al. [Phys. Rev. Lett. 93, 120401 (2004)] for an ideal Fermi gas. We demonstrate a breakdown of Josephson oscillation in the SFG for a large displacement of the harmonic trap. These features of Josephson oscillation of a SFG can be tested experimentally.

Study of a degenerate dipolar Fermi gas of Dy-161 atoms

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)