Advances in experimental atomic systems have given us access to highly tunable quantum systems, and to an unprecedented range of observables of these systems. One fundamental system that has been made accessible in this way is a gas of bosons trapped in a periodic potential. We present here a several studies of the many-body physics of bosons in optical lattices. We discuss the ferromagnetic effects of a single vacancy in a two-species gas on a lattice. We present a derivation of the superfluid density, the order parameter for the superfluid state, and discuss the mathematical subtleties of calculating it. Finally, we calculate the dynamics of bosonic lattice systems, in the presence of inelastic light scattering used as a density measurement, and during a ramp of the interaction strength from the Mott to the superfluid phase.
Advances in experimental atomic systems have given us access to highly tunable quantum systems, and to an unprecedented range of observables of these systems. One fundamental system that has been made accessible in this way is a gas of bosons trapped in a periodic potential. We present here a several studies of the many-body physics of bosons in optical lattices. We discuss the ferromagnetic effects of a single vacancy in a two-species gas on a lattice. We present a derivation of the superfluid density, the order parameter for the superfluid state, and discuss the mathematical subtleties of calculating it. Finally, we calculate the dynamics of bosonic lattice systems, in the presence of inelastic light scattering used as a density measurement, and during a ramp of the interaction strength from the Mott to the superfluid phase.
We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.
Ultracold atomic gases in optical lattices provide an unique framework to study quantum phenomena in strongly correlated systems. In addition to the precise control over all relevant parameters in the experiment, these experiments can be mapped to a fundamental model from solid-state physics. For the bosonic version of this model, the so-called Bose-Hubbard model, a phase transition from a superfluid to a Mott insulator was theoretically predicted and later experimentally observed in an ultracold gas of Rb-87 atoms in three-dimensional as well as in one-dimensional optical lattices. Apart from homogeneous optical lattices one can introduce more complex lattice topologies such as two-color superlattices which give rise to a rich phase diagram including more exotic phases like the Bose-glass. We employ various powerful many-body techniques like exact diagonalization in complete and truncated Hilbert spaces and the Density-Matrix Renormalization Group (DMRG) algorithm to study the phase diagrams of the one-dimensional Bose-Hubbard and the one-dimensional Bose-Fermi-Hubbard Hamiltonian. Most theoretical studies of these systems discuss the phase diagrams with respect to the generic parameters of the Hubbard model. These Hubbard parameters, however, depend non-trivially on the control parameters used in experiments. The focus of this work is on the ab-initio calculation of the phase diagram of ultracold Rb-87 in one-dimensional optical superlattices starting directly from the experimental setup. To this end, we first employ band-structure calculations to extract the Hubbard parameters from the experimental parameters. Then, we use state-of-the-art DMRG techniques to solve the many-body problem for realistic particle numbers and lattice sizes that occur in experiments. Our results show that by using the intensities of the two laser fields forming the two-color superlattice as control parameters while keeping all other experimental parameters fixed, it is possible to access all relevant quantum phases of the system. Furthermore, we found out that the critical values of the laser intensities for the different phase transitions depend strongly on a third parameter that has to be included for a realistic description of the experiment. This third parameter is the strength of a harmonic trapping potential which accounts for the Gaussian shape of the laser fields and an additional magnetic potential used to confine the atoms in the center of the trap.
Ultracold atoms have emerged as an indispensable setting to study quantum many-body systems. Recent experimental and theoretical work has explored the curious phases and novel properties of Bose-Einstein Condensate with optical lattices, and Bose-Einstein Condensate with light-induced artificial spin-orbit coupling. In this thesis, we report our research on these two types of boson systems. In the first topic, in contrast with calculations of bosons in optical lattices that focus on the tight-binding regime, we note that the single-particle states of bosons in a periodic potential generally satisfy the Mathieu equation, and have developed a formalism for studying bosons in an optical lattice using the Mathieu equation. Moreover, based on this formalism, we have proposed a self-consistent scheme for describing interacting bosons in an optical lattice using Hartree Fock approximation. We apply this scheme to quantify the effects of inter-atomic interactions on the properties of bosons in an optical lattice, as exhibited in the comparison between observables of non-interacting and interacting systems, such as the superfluid transition temperature and momentum distribution as probed in time-of-flight expansion. In the second topic, the phases of a Bose-Einstein condensate with light-induced spin-orbit coupling are studied within the mean-field approximation. We obtain the phase diagram at fixed chemical potential and at fixed density for bosons with spin-orbit coupling, finding a regime of phase separation and a regime in which the bosons condensed into a mixed phase. We determine how this phase diagram evolves as a function of the atom interaction parameters and as a function of the strength of light-atom coupling. The mixed phase is found to be stable for sufficiently small light-atom coupling. Specifically, we show that the structure of the phase diagram at fixed chemical potential suggests an unusual density dependence for the mixed phase in a harmonic trapping potential, in which the density of one spin increases with increasing radius, suggesting a unique experimental signature of this state. The collective Bogoliubov sound mode is shown to also provide a signature of the mixed phase, vanishing as the boundary to the regime of phase separation is approached. Together, in these two topics we address the need to enhance the understanding of the unconventional physical properties of Bose-Einstein Condensate in a controlled electromagnetic environment (optical lattices, Raman lasers, etc.), and provide predictions for possible experimental findings.
Motivated by recent experiments on trapped ultracold bosonic atoms in an optical lattice potential, we consider the nonequilibrium dynamic properties of such bosonic systems for a number of experimentally relevant situations. When the number of bosons per lattice site is large, there is a wide parameter regime where the effective boson interactions are strong, but the ground state remains a superfluid (and not a Mott insulator): we describe the conditions under which the dynamics in this regime can be described by a discrete Gross-Pitaevskii equation. We describe the evolution of the phase coherence after the system is initially prepared in a Mott insulating state, and then allowed to evolve after a sudden change in parameters places it in a regime with a superfluid ground state. We also consider initial conditions with a ``\ensuremath{\pi} phase'' imprint on a superfluid ground state (i.e., the initial phases of neighboring wells differ by \ensuremath{\pi}), and discuss the subsequent appearance of the density wave order and ``Schr\"odinger cat,'' i.e., macroscopic quantum interference, states.
Two models are discussed here to illustrate the quantum many-body phenomena in mixtures of ultra-cold atoms in optical lattices. The first model describes a mixture of two species of bosonic atoms of equal masses in optical lattices and the second describes a mixture of heavy bosonic atoms and light fermionic atoms in optical lattices. For both models, we assume the trap is present and use parameters typical in experiment. For the first model, the discussion is aimed at providing a thorough description of the collective behavior of the binary mixture in various interaction regions, with emphasis on two many-body phenomena, pairing and anti-pairing, as a result of the inter-species interaction. The pairing leads to a new type of superfluid order, called the paired superfluid (PSF) and the anti-pairing leads to another type of superfluid order, called the counter-flow superfluid (CFSF). In addition, we discuss the coexistence of charge density wave order with the three superfluid orders in the strong interaction region. We use both Luttinger liquid theory and the time evolving block decimation (TEBD) method to study this model in one dimension. The discussion is organized in three parts: the phase diagram and the correlation functions; the noise correlation functions; and the transport properties. Two phase diagrams are constructed to map the different orders in the parameter space. The correlation functions, include noise correlations, are carefully examined for the determination of the orders and for possible detection methods. In the end, the transport properties of the PSF and CFSF orders are studied through the dipole oscillation induced by trap displacement. For the second model, examining a mixture of heavy bosons and light fermions, the discussion is oriented toward determining the thermal properties of the mixture for attractive inter-species interactions. This work is motivated by experiments creating artificial molecules through optical and magnetic control of ultra-cold atoms. We use the strong coupling (SC) expansion method to evaluate the density profile, the onsite inter-species correlations, the density fluctuations and the entropy per particle. Analytical expressions are derived for all the quantities above as well as the partition function. To benchmark the accuracy, the SC calculations are compared with inhomogeneous dynamical mean field theory (IDMFT) and Monte Carlo (MC) simulation. From the calculations, we find that 1) the efficiency of creating pre-formed molecules is significantly increased by confining the mixtures onto optical lattices; 2) the temperature of the mixtures in optical lattices can be reliably estimated through the density gradient and the density fluctuations.
Bosonic atoms confined in optical lattices are described by the Bose-Hubbard model and can exist in two different phases, Mott insulator or superfluid, depending on the strength of the system parameters. In the vicinity of the phase boundary, there are degeneracies that occur between every two adjacent Mott lobes. Because of this, nondegenerate perturbation theory fails to give meaningful results for the condensate density: it predicts a phase transition in a point of the phase diagram where no transition occurs. Motivated by this, we develop two different degenerate perturbative methods to solve the degeneracy-related problems. Moreover, we study the one-dimensional repulsively interacting Bose gas under harmonic confinement, with special attention to the asymptotic behavior of the momentum distribution, which is a universal k−4 decay characterized by the Tan's contact. The latter constitutes a direct signature of the short-range correlations in such an interacting system and provides valuable insights about the role of the interparticle interactions. We investigate the system constituted of N interacting particles in the strongly interacting limit. In such a regime, the strong interparticle interaction makes the bosons behave similarly to the ideal Fermi gas. Because of the difficulty in analytically solving the system for N particles at finite interaction, the Tonks-Girardeau regime provides a favorable scenario to probe the contact. Therefore, we are able to provide an analytical formula for the Tan's contact. Furthermore, we analyze the scaling properties of the Tan's contact in terms of N in the high-temperature regime as well as in the strongly interacting regime. Finally, we compare our analytical calculations of the Tan's contact to quantum Monte Carlo simulations and discuss some fundamental differences between the canonical and the grand-canonical ensembles.