Title: Modelling of Ultracold Gases in Multidimensional Optical Lattices Author: Miroslav Urbanek Department: Department of Chemical Physics and Optics Supervisor: doc. Ing. Pavel Soldán, Dr. Abstract: Optical lattices are experimental devices that use laser light to confine ultracold neutral atoms to periodic spatial structures. A system of bosonic atoms in an optical lattice can be described by the Bose-Hubbard model. Although there exist powerful analytic and numerical methods to study this model in one dimension, their extensions to multiple dimensions have not been as successful yet. I present an original numerical method based on tree tensor networks to simulate time evolution in multidimensional lattice systems with a focus on the two-dimensional Bose-Hubbard model. The method is used to investigate phenomena accessible in current experiments. In particular, I have studied phase collapse and revivals, boson expansion, and many-body localization in two-dimensional optical lattices. The outcome of this work is TEBDOL - a program for modelling one-dimensional and two-dimensional lattice systems. Keywords: Bose-Hubbard model, multidimensional system, optical lattice, tensor network
Title: Modelling of Ultracold Gases in Multidimensional Optical Lattices Author: Miroslav Urbanek Department: Department of Chemical Physics and Optics Supervisor: doc. Ing. Pavel Soldán, Dr. Abstract: Optical lattices are experimental devices that use laser light to confine ultracold neutral atoms to periodic spatial structures. A system of bosonic atoms in an optical lattice can be described by the Bose-Hubbard model. Although there exist powerful analytic and numerical methods to study this model in one dimension, their extensions to multiple dimensions have not been as successful yet. I present an original numerical method based on tree tensor networks to simulate time evolution in multidimensional lattice systems with a focus on the two-dimensional Bose-Hubbard model. The method is used to investigate phenomena accessible in current experiments. In particular, I have studied phase collapse and revivals, boson expansion, and many-body localization in two-dimensional optical lattices. The outcome of this work is TEBDOL - a program for modelling one-dimensional and two-dimensional lattice systems. Keywords: Bose-Hubbard model, multidimensional system, optical lattice, tensor network
Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.
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 in optical lattices are one of the most promising experimental setups to simulate strongly correlated systems. However, efficient numerical algorithms able to benchmark experiments at low-temperatures in interesting 3d lattices are lacking. To this aim, here we introduce an efficient tensor network algorithm to accurately simulate thermal states of local Hamiltonians in any infinite lattice, and in any dimension. We apply the method to simulate thermal bosons in optical lattices. In particular, we study the physics of the (soft-core and hard-core) Bose-Hubbard model on the infinite pyrochlore and cubic lattices with unprecedented accuracy. Our technique is therefore an ideal tool to benchmark realistic and interesting optical-lattice experiments.
Experiments with ultracold atoms in optical lattice have become a versatile testing ground to study diverse quantum many-body Hamiltonians. A single-band Bose-Hubbard (BH) Hamiltonian was first proposed to describe these systems in 1998 and its associated quantum phase-transition was subsequently observed in 2002. Over the years, there has been a rapid progress in experimental realizations of more complex lattice geometries, leading to more exotic BH Hamiltonians with contributions from excited bands, and modified tunneling and interaction energies. There has also been interesting theoretical insights and experimental studies on “un- conventional” Bose-Einstein condensates in optical lattices and predictions of rich orbital physics in higher bands. In this thesis, I present our results on several multi- band BH models and emergent quantum phenomena. In particular, I study optical lattices with two local minima per unit cell and show that the low energy states of a multi-band BH Hamiltonian with only pairwise interactions is equivalent to an effec- tive single-band Hamiltonian with strong three-body interactions. I also propose a second method to create three-body interactions in ultracold gases of bosonic atoms in a optical lattice. In this case, this is achieved by a careful cancellation of two contributions in the pair-wise interaction between the atoms, one proportional to the zero-energy scattering length and a second proportional to the effective range. I subsequently study the physics of Bose-Einstein condensation in the second band of a double-well 2D lattice and show that the collision aided decay rate of the con- densate to the ground band is smaller than the tunneling rate between neighboring unit cells. Finally, I propose a numerical method using the discrete variable repre- sentation for constructing real-valued Wannier functions localized in a unit cell for optical lattices. The developed numerical method is general and can be applied to a wide array of optical lattice geometries in one, two or three dimensions.
In this dissertation I present results for lattice boson systems based on quantum Monte Carlo simulations by the Worm algorithm. Lattice bosons can be realized by ultracold gases in optical lattices and, under certain conditions, they are described by the Bose-Hubbard model. I will discuss the ground state, i.e. the phase diagram and properties of elementary excitations, as well as finite temperature thermodynamic properties of the two- and three-dimensional systems. These results provide bench-marks and guidance for experimentalists. I will also discuss the phase diagram of the one-dimensional hard-core Bose-Hubbard model with three-body interactions. The model is realized by cold polar molecules in optical lattices. The ground state presents a novel solid phase at unconventional filling, characterized by coexisting charge density wave and bond orders. Such phase has not been predicted by the preexisting mean field theory.
A major focus in condensed matter physics is to study the origin of exotic quantum phases such as coexistent and inhomogeneous phases, quantum criticality, and secondary ordered phases close to quantum critical points. Exotic phenomena in strongly correlated systems occur due to competing complex interactions of spin, charge, lattice, and orbital degrees of freedom. Complex quantum phases in strongly correlated systems are challenging, but they might be very useful due to their possible functionality to make advance devices. In order to understand, utilize, and optimize such behaviors, we need to improve our understanding of these systems. Studies in cold atom systems are of interest since cold atom experiments provide a control on model parameters. In this thesis, we use novel analytical and computational techniques to treat strongly interacting bosonic systems. Taking advantage of the very versatile quantum Monte Carlo Stochastic Green Function algorithm, we studied several interesting problems. First, we explore a recently developed new confining method for cold atoms on optical lattices. Atoms are confined via a hopping integral that decreases as a function of the distance from the center of the lattice. This method might lead to lower temperatures than existing diagonal confinement methods. Next, we study the ground state phase diagram of interacting bosons on a ring-shape lattice with a region of weak hopping integrals. The model, an extension of the well known Bose-Hubbard model, develops a novel local Mott phase in addition to the usual Mott and superfluid phases in the homogeneous system. This might provide a new insight to the description of atomtronics applications. Finally, we study the two species Bose Hubbard model in a two-dimensional lattice. This model presents novel phases due to the complexity associated with multiple species. Its phase diagram shows ordered and coexistence phases including a ferromagnetic phase separated phase with high entropy. This phase might be accessible experimentally. The novel phases found from our studies are linked to experiments on ultra-cold atoms trapped by laser beams.