I present simple analytical methods for computing the properties of ground and excited states of Bose-Einstein condensates, and compare their results to extensive numerical simulations.I consider the effect of vortices in the condensate for both positive and negative scattering lengths, a , and find an analytical expression for the large-N0 limit of the vortex critical frequency for a > 0, by approximate solution of the timeindependent nonlinear Schro ¨dinger equation.
I present simple analytical methods for computing the properties of ground and excited states of Bose-Einstein condensates, and compare their results to extensive numerical simulations.I consider the effect of vortices in the condensate for both positive and negative scattering lengths, a , and find an analytical expression for the large-N0 limit of the vortex critical frequency for a > 0, by approximate solution of the timeindependent nonlinear Schro ¨dinger equation.
For Bose–Einstein condensation of neutral atoms in anisotropic traps at zero temperature, we present simple analytical methods for computing the properties of ground state and single vortex of Bose–Einstein condensates, and compare those results to extensive numerical simulations. The critical angular velocity for production of vortices is calculated for both positive and negative scattering lengths a, and find an analytical expression for the large- N limit of the vortex critical angular velocity for a >0, and the critical number for condensate population approaches the point of collapse for a <0, by using approximate variational method.
Numerical solution of nonlinear Schrodinger equation (NLSE) can be used to describe an inhomogeneous,weakly interact Bose-Einstein condensate in a small harmonic trap potential at zero temperature. In one dimension situation, we find solutions for the NLSE for ground-state condensate wave function, and obtain the ground-state condensate energy per particle. Furthermore, with the time-dependent NLSE, we can examine the time evolution of the macroscopic wave function even when the trap potential is changed on a time scale comparable to that of the condensate dynamics, a situation that can be easily achieved in magneto-optical traps. We show that there are stable solutions for atomic species with both positive and negative s-wave scattering length in one-dimensional (1D) system for a fixed number of atoms, and the negative scattering length solutions have solitonlike properties.
In this paper, we study dynamics of the ground state and central vortex states in Bose–Einstein condensation (BEC) analytically and numerically. We show how to define the energy of the Thomas–Fermi (TF) approximation, prove that the ground state is a global minimizer of the energy functional over the unit sphere and all excited states are saddle points in linear case, derive a second-order ordinary differential equation (ODE) which shows that time-evolution of the condensate width is a periodic function with/without a perturbation by using the variance identity, prove that the angular momentum expectation is conserved in two dimensions (2D) with a radial symmetric trap and 3D with a cylindrical symmetric trap for any initial data, and study numerically stability of central vortex states as well as interaction between a few central vortices with winding numbers ±1 by a fourth-order time-splitting sine-pseudospectral (TSSP) method. The merit of the numerical method is that it is explicit, unconditionally stable, time reversible and time transverse invariant. Moreover, it conserves the position density, performs spectral accuracy for spatial derivatives and fourth-order accuracy for time derivative, and possesses "optimal" spatial/temporal resolution in the semiclassical regime. Finally we find numerically the critical angular frequency for single vortex cycling from the ground state under a far-blue detuned Gaussian laser stirrer in strong repulsive interaction regime and compare our numerical results with those in the literatures.
Amplifying on a proposal by O'Dell et al. for the realization of Bose-Einstein condensates of neutral atoms with attractive 1∕r interaction, we point out that the instance of self-trapping of the condensate, without an external trap potential, is physically best understood by introducing appropriate ``atomic'' units. This reveals a remarkable scaling property: the physics of the condensate depends only on the two parameters N2a∕au and $\ensuremath{\gamma}∕{N}^{2}$, where N is the particle number, a the scattering length, au the ``Bohr'' radius, and $\ensuremath{\gamma}$ the trap frequency in atomic units. We calculate accurate numerical results for self-trapping wave functions and potentials, and for energies, sizes, and peak densities, and compare with previous variational results. We point out the existence of a second solution of the extended Gross-Pitaevskii equation for negative scattering lengths, with and without trapping potential, which is born together with the ground state in a tangent bifurcation. This indicates the existence of an unstable collectively excited state of the condensate for negative scattering lengths.
We numerically study properties of the dynamics of vortices in nonrotating Bose-Einstein condensates in the Thomas-Fermi regime. On the one hand, we compute the vortex energy as a function of its position and we predict, using the expression of the Magnus force, the vortex precession velocity. On the other hand, we calculate the temporal evolution of the vortex-state and test the accuracy of the previous prediction. We also investigate the validity of analytical formulae of this velocity involving the healing length. In addition, we analyze the velocity field and the angular momentum and we compare them to available analytical expressions.
This dissertation focuses on the properties of nonlinear waves in Bose-Einstein condensates (BECs). The fundamental model here is the nonlinear Schrodinger equation, the so-called Gross-Pitaevskii (GP) equation, which is a mean-field description of BECs. The systematic analysis begins by considering the dark-bright (DB)-soliton interactions and multiple-dark-bright-soliton complexes in atomic two-component BECs. The interaction between two DB solitons in a homogeneous condensate and at the presence of the trap are both considered. Our analytical approximation relies in a Hamiltonian perturbation theory, which leads to an equation of motion of the centers of DB-soliton interacting pairs. Employing this equation, we demonstrate the existence of robust DB-soliton molecules, in the form of stationary two- and three-DB-soliton states. Also the equilibrium distance of the constituent solitons and the corresponding oscillation frequencies are found semianalytically, where the latter corresponds to the characteristic anomalous modes' eigenfrequencies that we numerically computed via a so called Bogoliubov-de Gennes (BdG) linearization analysis. Those studies are discussed in Chapter 2. Then, we extend our studies to the dynamics of dark-bright (DB) solitons in binary BECs at finite temperature using a system of two-coupled dissipative GPs. We show that the effect of the bright soliton is to partially stabilize dark solitons against temperature-induced dissipation, thus providing longer lifetimes in Chapter 3. Furthermore, the dark-dark (DD) solitons as a prototypical coherent structure that emerges in two-component BECs are studied and are connected to dark-bright (DB) solitons via SO(2) rotation. We obtained their beating frequency and their frequency of oscillation inside a parabolic trap. They are identified as exact periodic orbits in the Manakov limit of equal inter- and intra- species nonlinearity strengths with and without the trap and we showcase the persistence of such states upon weak deviations from this limit. Also we investigated in detail the effect of the deviation from the Manakov case by considering different from unity scattering length ratios in Chapter 4. Next, we revisited Hamiltonian eigenvalue problems that typically arise in the linearization around a stationary state of a Hamiltonian nonlinear PDE. Also we presented a overview of the known facts for the eigenvalue counts of the corresponding unstable spectra. In particular, we focused on a straightforward plan to implement finite-dimensional techniques for locating this spectrum via the singular points of the meromorphic Krein Matrix and illustrated the value of the approach by considering realistic problems for recently observed experimentally multivortex and multisoliton solutions in atomic Bose-Einstein condensates in Chapter 5. In the two dimensional scenario, we also examine the stability and dynamics of vortices under the effect of dissipation used as a simplified model for the inclusion of the effect of finite temperatures in atomic BECs, which enables an analytical prediction that can be compared directly to numerical results in Chapter 6. In all the above studies, our analytical prediction from the equation of motion are in good agreement with the numerical results from the BdG analysis.