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.

Effect of anharmonicities in the critical number of trapped condensed atoms with an attractive two-body interaction

The quantitative effect in the maximum number of particles and other static observables was determined. A deviation in the harmonic trap potential that is effective only outside the central part of the potential, with the addition of a term that is proportional to a cubic or quartic power of the distance was considered. Results showed that this study could be easily transferred to other trap geometries to estimate anharmonic effects.

Liquid-gas phase transition in Bose-Einstein condensates with time evolution

We study the effects of a repulsive three-body interaction on a system of trapped ultracold atoms in a Bose-Einstein condensed state. The stationary solutions of the corresponding s-wave nonlinear Schrödinger equation suggest a scenario of first-order liquid-gas phase transition in the condensed state up to a critical strength of the effective three-body force. The time evolution of the condensate with feeding process and three-body recombination losses has a different characteristic pattern. Also, the decay time of the dense (liquid) phase is longer than expected due to strong oscillations of the mean-squared radius.

Chaos in collapsing Bose-condensed gas

We reinvestigate the dynamics of the grow and collapse of Bose-Einstein condensates in a system of trapped ultracold atoms with negative scattering lengths, and found a new behavior in the long time scale evolution: the number of atoms can go far beyond the static stability limit. The condensed state is described by the solution of the time-dependent nonlinear Schrödinger equation, in a model that includes atomic feeding and three-body dissipation. Our results for the model show that, by changing the feeding parameter and when a substantial depletion of the ground-state exists, a chaotic behavior is found. We consider a criterion proposed by Deissler and Kaneko [Phys. Lett. A 119, 397 (1987)] to diagnose spatiotemporal chaos. ©2000 The American Physical Society.

Mean-field description of a dynamical collapse of a fermionic condensate in a trapped boson-fermion mixture

We suggest a time-dependent dynamical mean-field-hydrodynamic model for the collapse of a trapped boson-fermion condensate and perform numerical simulation based on it to understand some aspects of the experiment by Modugno et al. [Science 297, 2240 (2002)] on the collapse of the fermionic condensate in the K-40-Rb-87 mixture. We show that the mean-field model explains the formation of a stationary boson-fermion condensate at zero temperature with relative sizes compatible with experiment. This model is also found to yield a faithful representation of the collapse dynamics in qualitative agreement with experiment. In particular we consider the collapse of the fermionic condensate associated with (a) an increase of the number of bosonic atoms as in the experiment and (b) an increase of the attractive boson-fermion interaction using a Feshbach resonance. Suggestion for experiments of fermionic collapse using a Feshbach resonance is made.

Stability and collapse of a hybrid Bose-Einstein condensate of atoms and molecules

The dynamics of stability and collapse of a trapped atomic Bose-Einstein condensate (BEC) coupled to a molecular one is studied using the time-dependent Gross-Pitaevskii (GP) equation including a nonlinear interaction term which can transform two atoms into a molecule and vice versa. We find an interesting oscillation of the number of atoms and molecules for a BEC of fixed mass. This oscillation is a consequence of continuous transformation in the condensate of two atoms into a molecule and vice versa. For the study of collapse an absorptive contact interaction and an imaginary quartic three-body recombination term are included in the GP equation. It is possible to have a collapse of one or both components when one or more of the nonlinear terms in the GP equation are attractive in nature, respectively.

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.

BCS-BEC crossover in a trapped Fermi super-fluid using a density-functional equation

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

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

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

Solitons in Bose-Einstein condensates trapped in a double-well potential

We investigate, analytically and numerically, families of bright solitons in a system of two linearly coupled nonlinear Schrodinger/Gross-Pitaevskii equations, describing two Bose-Einstein condensates trapped in an asymmetric double-well potential, in particular, when the scattering lengths in the condensates have arbitrary magnitudes and opposite signs. The solitons are found to exist everywhere where they are permitted by the dispersion law. Using the Vakhitov-Kolokolov criterion and numerical methods, we show that, except for small regions in the parameter space, the solitons are stable to small perturbations. Some of them feature self-trapping of almost all the atoms in the condensate with no atomic interaction or weak repulsion is coupled to the self-attractive condensate. An unusual bifurcation is found, when the soliton bifurcates from the zero solution with vanishing amplitude and width simultaneously diverging but at a finite number of atoms in the soliton. By means of numerical simulations, it is found that, depending on values of the parameters and the initial perturbation, unstable solitons either give rise to breathers or completely break down into incoherent waves (radiation). A version of the model with the self-attraction in both components, which applies to the description of dual-core fibers in nonlinear optics, is considered too, and new results are obtained for this much studied system. (C) 2003 Elsevier B.V. All rights reserved.

Stability of trapped Bose-Einstein condensates

In three-dimensional trapped Bose-Einstein condensate (BEC), described by the time-dependent Gross-Pitaevskii-Ginzburg equation, we study the effect of initial conditions on stability using a Gaussian variational approach and exact numerical simulations. We also discuss the validity of the criterion for stability suggested by Vakhitov and Kolokolov. The maximum initial chirp (initial focusing defocusing of cloud) that can lead a stable condensate to collapse even before the number of atoms reaches its critical limit is obtained for several specific cases. When we consider two- and three-body nonlinear terms, with negative cubic and positive quintic terms, we have the conditions for the existence of two phases in the condensate. In this case, the magnitude of the oscillations between the two phases are studied considering sufficient large initial chirps. The occurrence of collapse in a BEC with repulsive two-body interaction is also shown to be possible.

Dipolar droplet bound in a trapped Bose-Einstein condensate

We study the statics and dynamics of a dipolar Bose-Einstein condensate (BEC) droplet bound by interspecies contact interaction in a trapped nondipolar BEC. Our findings are demonstrated in terms of stability plots of a dipolar 164Dy droplet bound in a trapped nondipolar 87Rb BEC with a variable number of 164Dy atoms and interspecies scattering length. A trapped nondipolar BEC of a fixed number of atoms can bind only a dipolar droplet containing fewer atoms than a critical number for the interspecies scattering length between two critical values. The shape and size (statics) as well as the small breathing oscillation (dynamics) of the dipolar BEC droplet are studied using numerical and variational solutions of a mean-field model. We also suggest an experimental procedure for achieving such a 164Dy droplet by relaxing the trap on the 164Dy BEC in a trapped binary 87Rb-164Dy mixture. © 2013 American Physical Society.

Stability of trapped degenerate dipolar Bose and Fermi gases

Trapped degenerate dipolar Bose and Fermi gases of the cylindrical symmetry with the polarization vector along the symmetry axis are only stable for the strength of dipolar interaction below a critical value. In the case of bosons, the stability of such a dipolar Bose-Einstein condensate (BEC) is investigated for different strengths of contact and dipolar interactions using a variational approximation and a numerical solution of a mean-field model. In the disc shape, with the polarization vector perpendicular to the plane of the disc, the atoms experience an overall dipolar repulsion and this fact should contribute to the stability. However, a complete numerical solution of the dynamics leads to the collapse of a strongly disc-shaped dipolar BEC due to the long-range anisotropic dipolar interaction. In the case of fermions, the stability of a trapped single-component degenerate dipolar Fermi gas is studied including the Hartree-Fock exchange and Brueckner-Goldstone correlation energies in the local-density approximation valid for a large number of atoms. Estimates for the maximum allowed number of polar Bose and Fermi molecules in the BEC and degenerate Fermi gas are given. © 2013 IOP Publishing Ltd.

Neutral Higgs bosons in the SU(3)(L)circle times U(1)(N) model

The SU(3)(L)circle times U(1)(N) electroweak model predicts new Higgs bosons beyond the one of the standard model. In this work we investigate the signature and production of neutral SU(3)(L)circle times U(1)(N) Higgs bosons in the e(-)e(+) Next Linear Collider and in the CERN Linear Collider . We compute the branching ratios of two of the SU(3)(L)circle times U(1)(N) neutral Higgs bosons and study the possibility to detect them and the Z' extra neutral boson of the model.

Distribution of neutral prilocaine in a phospholipid bilayer: Insights from molecular dynamics simulations

In this work, we report a 20-ns constant pressure molecular dynamics simulation of prilocaine (PLC), in amine-amide local anesthetic, in a hydrated liquid crystal bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylcholine. The partition of PLC induces the lateral expansion of the bilayer and a concomitant contraction in its thickness. PLC molecules are preferentially found in the hydrophobic acyl chains region, with a maximum probability at similar to 12 angstrom from the center of the bilayer (between the C(4) and C(5) methylene groups). A decrease in the acyl chain segmental order parameter, vertical bar S-CD vertical bar, compared to neat bilayers, is found, in good agreement with experimental H-2-NMR studies. The decrease in vertical bar S-CD vertical bar induced by PLC is attributed to a larger accessible volume per lipid in the acyl chain region. (C) 2008 Wiley Periodicals, Inc.