Presented is a study of a spin-2 Bose-Einstein condensate (BEC) by unitary quantum simulations of time-dependent soliton-soliton scattering. The quantum simulation method is based on a quantum lattice algorithm which is designed for implementation on a future digital quantum computer but is tested today using a parallel computing architecture based on graphical processing units (GPUs). We analytically solve the spin-2 BEC equations of motion, a nonlinear system of 5 coupled Gross-Pitiaevskii (GP) equations, in one- and two-spatial dimensions. In 1D there are 16 bright soliton and 16 dark soliton solutions. In 2D there are 3 dark soliton solutions Pade approximation solutions, for m_f=+/-2, m_f=+/-1 and m_f=0, corresponding to quantum vortices. We report on the implementation the unitary quantum lattice gas algorithm for spinor superfluid and establish its efficacy by validating the stability of the 1D and 2D energy eigenstate solutions of the spin-2 BEC Hamiltonian. Using the calibrated quantum lattice gas algorithm, the highly nonlinear physics in the nonintegrable regime of the spin-2 BEC is studied by performing soliton-soliton scattering experiments. The scattering of topological solitons produces breathers and complex quantum vortices characterized by local entanglement across multiple m_f-hyperfine states of the Zeeman manifold.
Presented is a study of a spin-2 Bose-Einstein condensate (BEC) by unitary quantum simulations of time-dependent soliton-soliton scattering. The quantum simulation method is based on a quantum lattice algorithm which is designed for implementation on a future digital quantum computer but is tested today using a parallel computing architecture based on graphical processing units (GPUs). We analytically solve the spin-2 BEC equations of motion, a nonlinear system of 5 coupled Gross-Pitiaevskii (GP) equations, in one- and two-spatial dimensions. In 1D there are 16 bright soliton and 16 dark soliton solutions. In 2D there are 3 dark soliton solutions Pade approximation solutions, for m_f=+/-2, m_f=+/-1 and m_f=0, corresponding to quantum vortices. We report on the implementation the unitary quantum lattice gas algorithm for spinor superfluid and establish its efficacy by validating the stability of the 1D and 2D energy eigenstate solutions of the spin-2 BEC Hamiltonian. Using the calibrated quantum lattice gas algorithm, the highly nonlinear physics in the nonintegrable regime of the spin-2 BEC is studied by performing soliton-soliton scattering experiments. The scattering of topological solitons produces breathers and complex quantum vortices characterized by local entanglement across multiple m_f-hyperfine states of the Zeeman manifold.
The nonlinear Schrodinger (NLS) equation in a self-defocusing Kerr medium supports dark solitons. Moreover the mean field description of a dilute Bose-Einstein condensate (BEC) is described by the Gross-Pitaevskii equation, which for a highly anisotropic (cigar-shaped) magnetic trap reduces to a one-dimensional (1D) cubic NLS in an external potential. A quantum lattice algorithm is developed for the dark solitons. Simulations are presented for both black (stationary) solitons as well as (moving) dark solitons. Collisions of dark solitons are compared with the exact analytic solutions and coupled dark-bright vector solitons are examined. The quantum algorithm requires 2 qubits per scalar field at each spatial node. The unitary collision operator quantum mechanically entangles the on-site qubits, and this transitory entanglement is spread throughout the lattice by the streaming operators. These algorithms are suitable for a Type-II quantum computers, with wave function collapse induced by quantum measurements required to determine the coupling potentials.
Bose–Einstein Condensation (BEC) of atomic gases has attracted a renewed theoretical and experimental interest in quantum many-body systems at extremely low temperatures (Pethick & Smith; 2002). This excitement stems from two favorable features: (1) by applying magnetic fields and lasers, most of the system parameters, such as the shape, dimensionality, internal states of the condensates, and even the strength of the interatomic interactions, are controllable; (2) due to the diluteness, the mean-field theory explains experiments quite well. In particular, the Gross–Pitaevskii (GP) equation demonstrates its validity as a basic equation for the condensate dynamics. The GP equation is the counterpart of the nonlinear Schrodinger (NLS) equation in nonlinear optics. Thus, a study based on nonlinear analysis is possible and important. In nonlinear physics, a soliton is remarkable object not only for the fact that exact solutions can be obtained but also for its usefulness as a communications tool due to its robustness. In general, solitons are formed under the balance between nonlinearity and dispersion. For atomic condensates, the former is attributed to the interatomic interactions, while the latter comes from the kinetic energy. Either dark or bright solitons are allowable depending on the positive or negative sign of the interatomic coupling constants g, respectively, and indeed have been observed in a quasione-dimensional (q1D) optically constructed waveguide (Strecker et al.; 2002) (Khaykovich et al.; 2002). Such matter-wave solitons are expected in atom optics for applications in atom laser, atom interferometry, and coherent atom transport (Meystre; 2002). In this chapter, we extend the analysis of the matter-wave solitons to a multicomponent case by considering the so-called spinor condensate (Stenger et al.; 1998) whose spin degrees of freedom are liberated under optical traps. Based on theoretical and experimental results, we introduce a new integrable model which describes the dynamical properties of the matter-wave soliton of spinor condensates (Ieda et al.; 2004a). We employ the inverse scattering method to solve this model exactly. As a result, we predict the occurrence of undiscovered physical phenomena such as macroscopic spin precession and spin switching. The chapter is organized as follows. In Sec. 2, the mean field theory of condensate is briefly reviewed. Section 3 introduces an effective interatomic coupling in a q1D condensate. Using these results, we consider a spinor condensate in q1D regime in Sec. 4. Then, in Sec. 5, we
Presented is an unitary operator splitting method for handling the spin-density interaction in spinor Bose-Einstein condensates. The zero temperature behavior of a spinor BEC is given by mean field theory, where the Hamiltonian includes a nonlinear hyperfine spin interaction. This hyperfine interaction has a diagonal probability-density term (leading to the usual Gross-Pitaevskii type equation of motion) but also has a nondiagonal spin-density term. Since the F=2 spinor BEC (spin-2 BEC) has a non-Abelian superfluid phase (nonperturbative cyclic phase in the strong spin-density coupling regime), an infinite-order expansion of the quantum evolution operator is needed for quantum simulation applications. An infinite-order expansion, obtained by analytical continuation and expressed in analytically closed form, for the spin-2 BEC is presented.
We investigate the different types of matter-wave solitons in spin-orbit-coupled spin-2 spinor Bose-Einstein condensates. Using mean-field theory and adopting the multiscale perturbation method, the original five-component Gross-Pitaevskii spin-orbit-coupled spin-2 spinor Bose-Einstein condensate model can be reduced to a single effective nonlinear Schrödinger equation, which allows us to find analytical soliton solutions of this system. In this way, for different regimes of the spin-orbit coupling, Raman coupling, and interatomic interactions, we find approximate bright and dark soliton solutions. Particularly, the type of solitons depends on the dispersion properties of the system. When the lowest-energy band has a single-well structure, we find there only exist positive mass bright or dark solitons due to the dispersion coefficient of effective nonlinear Shrödinger equation always positive. However, when the lowest-energy band has a double-well structure, there will appear positive (negative) mass bright or dark solitons because the sign of the dispersion coefficient can be positive (negative) under different momentum. We employ direct numerical simulation of the original five-component Gross-Pitaevskii equations to confirm the analytical results.
Spinor Bose Einstein Condensates are intriguing because of their vast range of different topological vortices. These states occur when a BEC gas is trapped in an optical lattice rather than in a magnetic well (which would result in scalar BEC vortices). A spinor BEC states also occur in a quantum gas when several hyperfine states of the atom co-exist in the same trap. A unitary quantum lattice algorithm that is ideally parallelized to all available processors is used to solve the evolution of non-eigenstate Skyrmions in a coupled BEC system. The incompressible kinetic energy spectrum of the inner quantum vortex ring core rapidly deviates from the k-3 spectrum found in the evolution of scalar BECs.
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.