The central aim of this thesis is the study of the low-energy and low-temperature properties of strongly correlated systems of bosonic particles interacting via finite- and long-range potentials, and relevant to experimental realization with cold atomic gases. This study is carried out with a combination of state-of-the-art numerical techniques such as Path Integral Monte Carlo and analytical techniques. The main result of my work is the demonstration of the existence of a stripe supersolid phase and of a rare transition between isotropic and anisotropic supersolids in a finite-range interacting model of hard-bosons on a square lattice. I also investigate the out-of-equilibrium scenarios of such models via simulated temperature quenches. Finally, I investigate how restoring energy extensivity in long-range interacting systems can have a profound incidence on the low-energy properties in the thermodynamic limit.
The central aim of this thesis is the study of the low-energy and low-temperature properties of strongly correlated systems of bosonic particles interacting via finite- and long-range potentials, and relevant to experimental realization with cold atomic gases. This study is carried out with a combination of state-of-the-art numerical techniques such as Path Integral Monte Carlo and analytical techniques. The main result of my work is the demonstration of the existence of a stripe supersolid phase and of a rare transition between isotropic and anisotropic supersolids in a finite-range interacting model of hard-bosons on a square lattice. I also investigate the out-of-equilibrium scenarios of such models via simulated temperature quenches. Finally, I investigate how restoring energy extensivity in long-range interacting systems can have a profound incidence on the low-energy properties in the thermodynamic limit.
Disorder is ubiquitous in condensed matter systems, but a thorough understanding of these systems is hindered by poor control over the nature of disorder and competing interactions. Ultracold gases trapped in optical lattices, on the other hand, offer an unprecedented level of control over interactions and disorder, and are ideal candidates to test theoretical predictions and study fascinating phenomena resulting from the interplay between disorder, interaction, and quantum degeneracy. In this dissertation, I present my research on strongly correlated many-body lattice bosons in the presence of disordered potentials. My results are based on quantum Monte Carlo simulations by the Worm algorithm. Experimentally, lattice bosons can be realized by ultracold gases in optical lattices. I studied the ground state phase diagram of the above systems in the presence of diagonal (quasi-periodic disordered potential or random disordered potential) and off-diagonal disorder (in the form of site dilution). First, I will discuss the ground state phase diagram of two-dimensional bosons in quasi-periodic lattices. Then, I will discuss the phase diagram as well as finite temperature thermodynamic properties of two-dimensional dipolar lattice bosons in the presence of random diagonal disordered potentials. Finally, I will discuss the phase diagram of two-dimensional hard-core lattice bosons system in the presence of long-range hopping and long-range two-body interactions both in the absence of disorder and in the presence of site dilution (equivalent to off-diagonal disorder). This system is equivalent to the long-range XXZ model as realized by polar molecules trapped in a deep optical lattice. These results provide theoretical guidance for experimentalists to find quantum phase transitions as well as exotic quantum phases in strongly correlated bosonic lattice systems.
The study of out-of-equilibrium dynamics in ultra-cold gases is a new and exciting field, driven largely by the recent experimental advances in controlling and imaging cold clouds. The experimental and theoretical work thus far has been somewhat exploratory and largely numerical in nature, as the very paradigms for thinking about these systems are not well established. In this thesis I consider several different scenarios of ultra-cold bosonic and fermionic gases driven out of equilibrium and study their properties. In Chapter 1, I provide an overview of the phenomenology of ultra-cold gases, highlighting the timescales governing these systems and how the experimentalist can tune them. I discuss how cold gases can be cooled and trapped and discuss the basic physics behind optical lattices. I also discuss experimental probes of these gases, in particular the new high resolution imaging techniques developed recently at Chicago, Munich and Harvard. In Chapter 2, I discuss an early experiment (circa. 2008) which observed long lived spin dynamics in a thermal spin-1/2 Fermi gas. This experiment is an nice illustration of interesting physics resulting from the separation of timescales between spin and collisional dynamics. In my opinion, it is an excellent example of why cold gases are naturally suited to studying non-equilibrium dynamics. I simulate the experiment numerically using a collisionless Boltzmann equation and explain the observed spin dynamics both qualitatively and quantitively. In Chapter 3, I continue the discussion of spin waves in thermal gases by extending previous works on spin-1/2 gases to spin-1 Bose gases. In contrast to Chapter 2, the bulk of the work in this Chapter is analytic in nature. In particular, I find a spin wave instability in the thermal spin-1 Bose gas, which is the high temperature analog of the polar to ferromagnetic transition in a spin-1 Bose Einstein condensate. In Chapter 3, I turn my attention to bosonic systems and briefly review the the Bogoliubov mean-field theory. I calculate the momentum distribution and density-density correlation function of an interacting Bose gas within the Bogoliubov framework. Then I consider bosons in an optical lattice, and introduce the Bose Hubbard model. I calculate the mean-field phase diagram of the Bose Hubbard model and then consider fluctuations about the mean field, and derive the excitation spectrum of the lattice gas in the superfluid and insulating regimes. In Chapter 4, I ask what we learn by studying the dynamics of correlation functions following a sudden change in the interactions in a superfluid. Using the Bogoliubov theory developed in Chapter 3, I will show how the underlying excitation spectrum influences the long and short time behavior of the correlation functions. By considering a lattice dispersion, I study the analogous problem in a weak optical lattice and discuss how the lattice dispersion leads to additional features in the correlation functions. I will also discuss the timescale governing the revival of the condensate fraction in a quantum depleted gas. In Chapter 5, I derive equations of motion governing the dynamics of one and two body correlation functions in the single-band Bose Hubbard model, applicable to bosons in deep lattices. I then consider a simple quench from a Mott insulating initial state to a weakly interacting final state and produce analytic expressions describing the dynamics of correlations following such a quench. I discuss the timescale for the development of long range order following such a quench. I study the problem of chapter 4 using an equations of motion approach. This approach complements the Bogoliubov approach of Chapter 4. First, I derive exact expressions for a quench to a non-interacting state. I then consider how interactions redistribute quasi-momentum to first order in perturbation theory in different dimensions. In Chapter 6, I calculate the relevant timescales for local and global dynamics in trapped lattice Bose gases, a work done in collaboration with Dr. Kaden R. A Hazzard. Using a time-dependent Gutzwiller mean-field theory, I show that the timescale for local equilibration in these systems is fast in experimental terms. I then show that due to the spatial inhomogeneities inherent to cold gases, achieving global equilibrium can be quite complicated, sometimes taking longer than the lifetime of the experiment, an issue of practical importance to current day experiments. I continue this discussion in Chapter 7 which is a collaborative work with experimentalists David McKay and Prof. Brian DeMarco from the University of Toronto and the University of Illinois, Urbana Champaign. Using experimental and numerical methods, we show that the rapid timescales for local dynamics in interacting systems invalidates a frequently used cold atom technique for mapping out the momentum distribution of atoms in an optical lattice.
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.
Pairing and topological phases in cold atoms with long range interactions Chapter 1 Introduction The systematic study and classification of phase transitions has become a popular research area in physics over the past decades. The most familiar phenomena associated to the concept of phase transition are the ones giving rise to macroscopic changes in the properties of strongly correlated many-body systems due to thermal fluctuations. The most prominent examples of this kind of behaviour are the ferromagnetic-paramagnetic transition in lattice spin models and the liquid-gas transition. On the other hand, when dealing with quantum-mechanical systems changes of matter from one state to another can arise as the result of the variation of parameters other than temperature. In particular, when a given many-body quantum problem is studied at zero temperature, thereby suppressing thermal fluctuations, it is the interplay between kinetic energy and interactions that can cause the system to jump into different phases as the model parameters are varied. Phase transitions occurring at zero temperature are widely known as quantum phase transitions. In the framework of the physics of phase transitions, the role of dimensionality is a noticeable one. The major challenge in this perspective is the treatment of low dimensional systems, where standard mean-field treatments and ordinary perturbation theory are known to fail due to the enhancement of fluctuations. The goal of describing collective behaviours in reduced dimensionality is usually pursued in one-dimensional (1D) systems, which turn out to be simple enough to be studied in an effective way through both numerical and analytical techniques, while exhibiting a rich and deeply fascinating phenomenology. As a matter of fact, since particles confined to one dimension cannot avoid interactions with each other, the typical low energy excitations of the system are represented by collective density waves, which substitute nearly free quasiparticle excitations characterizing the behaviour of three-dimensional Fermi liquids in the description of the effects of the interplay between kinetic fluctuations and interactions. On the other hand, the tools employed in the activity of tackling 1D manybody problems are efficient and well developed. The analytical techniques aremainly based on the bosonization technique, which allows for the reformulation of strongly correlated 1D systems onto free bosonic theories, which can then be analyzed by means of standard path integral techniques. Meanwhile, the density matrix renormalization group (DMRG) algorithm represents the state-of-theart numerical machinery in the study of 1D physics, since it allows to efficiently extract the low energy properties of the model Hamiltonian and consequently characterize the corresponding phase diagram by means of the behaviour of properly chosen order parameters. As a final remark pointing out the relevance of the research effort in the direction of a better understanding of 1D many-body quantum systems, it is worth noticing that 1D systems do not represent only a playground for theorists, but have been experimentally realized by means of setups with cold atoms trapped in optical lattices, where the experimentalist is able to tune the model parameters by varying the characteristic features of the trapping laser beams and meanwhile has a comfortable access to the most relevant observables.
This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems.
We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the 1/k4 momentum tail.
An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling.
In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.
During this Ph.D, we studied one-dimensional systems with long-range couplings. In the first part, we demonstrate that power-law couplings lead to an algebraic decay of correlations at long distances in disordered quantum wires. In the second chapter, we analyzed an extended Hubbard model where particles interact via a finite-range potential that induces frustration and new exotic phases. In the third chapter, we demonstrated that restoring energy extensivity has an influence on the low-energy properties of quantum model in the thermodynamic limit. Finally, we provide preliminary results on the modification of Anderson localization due to the coupling to a cavity mode.