Numerical study of the spherically symmetric Gross-Pitaevskii equation in two space dimensions

The Gross-Pitaevskii equation for Bose-Einstein condensation (BEC) in two space dimensions under the action of a harmonic oscillator trap potential for bosonic atoms with attractive and repulsive interparticle interactions was numerically studied by using time-dependent and time-independent approaches. In both cases, numerical difficulty appeared for large nonlinearity. Nonetheless, the solution of the time-dependent approach exhibited intrinsic oscillation with time iteration which is independent of space and time steps used in discretization.

Two-dimensional supersonic nonlinear Schrodinger flow past an extended obstacle

Supersonic flow of a superfluid past a slender impenetrable macroscopic obstacle is studied in the framework of the two-dimensional (2D) defocusing nonlinear Schroumldinger (NLS) equation. This problem is of fundamental importance as a dispersive analog of the corresponding classical gas-dynamics problem. Assuming the oncoming flow speed is sufficiently high, we asymptotically reduce the original boundary-value problem for a steady flow past a slender body to the one-dimensional dispersive piston problem described by the nonstationary NLS equation, in which the role of time is played by the stretched x coordinate and the piston motion curve is defined by the spatial body profile. Two steady oblique spatial dispersive shock waves (DSWs) spreading from the pointed ends of the body are generated in both half planes. These are described analytically by constructing appropriate exact solutions of the Whitham modulation equations for the front DSW and by using a generalized Bohr-Sommerfeld quantization rule for the oblique dark soliton fan in the rear DSW. We propose an extension of the traditional modulation description of DSWs to include the linear ""ship-wave"" pattern forming outside the nonlinear modulation region of the front DSW. Our analytic results are supported by direct 2D unsteady numerical simulations and are relevant to recent experiments on Bose-Einstein condensates freely expanding past obstacles.

Local approximations for polarization potentials

We discuss the derivation of an equivalent polarization potential independent of angular momentum l for use in the optical Schrodinger equation that describes the elastic scattering of heavy ions. Three different methods are used for this purpose. Application of our theory to the low energy scattering of light heavy-ion systems at near-barrier energies is made. It is found that the notion of an l-independent polarization potential has some validity but cannot be a good substitute for the l-dependent local equivalent Feshbach polarization potential.

Decoherence-free evolution of time-dependent superposition states of two-level systems and thermal effects

In this paper we detail some results advanced in a recent letter [Prado et al., Phys. Rev. Lett. 102, 073008 (2009).] showing how to engineer reservoirs for two-level systems at absolute zero by means of a time-dependent master equation leading to a nonstationary superposition equilibrium state. We also present a general recipe showing how to build nonadiabatic coherent evolutions of a fermionic system interacting with a bosonic mode and investigate the influence of thermal reservoirs at finite temperature on the fidelity of the protected superposition state. Our analytical results are supported by numerical analysis of the full Hamiltonian model.

Simulations of Bose fields at finite temperature

We introduce a time-dependent projected Gross-Pitaevskii equation to describe a partially condensed homogeneous Bose gas, and find that this equation will evolve randomized initial wave functions to equilibrium. We compare our numerical data to the predictions of a gapless, second order theory of Bose-Einstein condensation [S. A. Morgan, J. Phys. B 33, 3847 (2000)], and find that we can determine a temperature when the theory is valid. As the Gross-Pitaevskii equation is nonperturbative, we expect that it can describe the correct thermal behavior of a Bose gas as long as all relevant modes are highly occupied. Our method could be applied to other boson fields.

Chebyshev real wave packet propagation: H+O-2 (J=0) state-to-state reactive scattering calculations

In this paper we explore the relative performance of two recently developed wave packet methodologies for reactive scattering, namely the real wave packet Chebyshev domain propagation of Gray and Balint-Kurti [J. Chem. Phys. 108, 950 (1998)] and the Lanczos subspace wave packet approach of Smith [J. Chem. Phys. 116, 2354 (2002); Chem. Phys. Lett. 336, 149 (2001)]. In the former method, a modified Schrodinger equation is employed to propagate the real part of the wave packet via the well-known Chebyshev iteration. While the time-dependent wave packet from the modified Schrodinger equation is different from that obtained using the standard Schrodinger equation, time-to-energy Fourier transformation yields wave functions which differ only trivially by normalization. In the Lanczos subspace approach the linear system of equations defining the action of the Green operator may be solved via either time-dependent or time-independent methods, both of which are extremely efficient due to the simple tridiagonal structure of the Hamiltonian in the Lanczos representation. The two different wave packet methods are applied to three dimensional reactive scattering of H+O-2 (total J=0). State-to-state reaction probabilities, product state distributions, as well as initial-state-resolved cumulative reaction probabilities are examined. (C) 2002 American Institute of Physics.

Real spectra for the non-Hermitian Dirac equation in 1+1 dimensions with the most general coupling

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Solutions of the Fokker-Planck equation for a Morse isospectral potential

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Variational supersymmetric approach to evaluate Fokker-Planck probability

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A General Solution of the Fokker-Planck Equation

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Atomic Bose-Einstein condensation with three-body interactions and collective excitations

The stability of a Bose-Einstein condensed state of trapped ultra-cold atoms is investigated under the assumption of an attractive two-body and a repulsive three-body interaction. The Ginzburg-Pitaevskii-Gross (GPG) nonlinear Schrodinger equation is extended to include an effective potential dependent on the square of the density and solved numerically for the s-wave. The lowest collective mode excitations are determined and their dependences on the number of atoms and on the strength of the three-body force are studied. The addition of three-body dynamics can allow the number of condensed atoms to increase considerably, even when the strength of the three-body force is very small compared with the strength of the two-body force. We study in detail the first-order liquid-gas phase transition for the condensed state, which can happen in a critical range of the effective three-body force parameter.

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.

Effective nonlinear Schrodinger equations for cigar-shaped and disc-shaped Fermi superfluids at unitarity

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

Solving Schrödinger equation for two dimensional potential using supersymmetry

The formalism of supersymmetric Quantum Mechanics can be extended to arbitrary dimensions. We introduce this formalism and explore its utility to solve the Schodinger equation for a bidimensional potential. This potential can be applied in several systens in physical and chemistry context, for instance, it can be used to study benzene molecule.

Wide localized solitons in systems with time- and space-modulated nonlinearities

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)