This thesis presents a series of theoretical studies of ultra cold atomic systems which model and propose experiments, and develop new computational techniques in order to elucidate aspects of many-body physics and non-equilibrium dynamics. In the first two studies I model the dynamics of nonlinear solitonic excitations in ultracold fermionic superfluids: the first simulates recent experiments and supports the hypothesis that the solitons generated in those experiments are unstable to the formation of vortex rings; the second demonstrates how population imbalance between up and down spin fermions can be used to prevent this instability. In the next study I discuss a method for generating and probing topologically protected edge states using periodically driven optical lattices potentials. Next I use a perturbative approach to study the spectral density of fermions with strong attractive interactions in the normal phase. After that I develop a novel cluster expansion technique to model the dynamics of interacting fermions in a disordered optical lattice. Finally I apply a Ginzurg-Landau theory to model experimental studies of superfluid 3He embedded in nematically ordered aerogel, finding evidence for a new phase of matter --the ``polar phase"-- which is not seen in bulk 3He.
This thesis presents a series of theoretical studies of ultra cold atomic systems which model and propose experiments, and develop new computational techniques in order to elucidate aspects of many-body physics and non-equilibrium dynamics. In the first two studies I model the dynamics of nonlinear solitonic excitations in ultracold fermionic superfluids: the first simulates recent experiments and supports the hypothesis that the solitons generated in those experiments are unstable to the formation of vortex rings; the second demonstrates how population imbalance between up and down spin fermions can be used to prevent this instability. In the next study I discuss a method for generating and probing topologically protected edge states using periodically driven optical lattices potentials. Next I use a perturbative approach to study the spectral density of fermions with strong attractive interactions in the normal phase. After that I develop a novel cluster expansion technique to model the dynamics of interacting fermions in a disordered optical lattice. Finally I apply a Ginzurg-Landau theory to model experimental studies of superfluid 3He embedded in nematically ordered aerogel, finding evidence for a new phase of matter --the ``polar phase"-- which is not seen in bulk 3He.
Advances in experiments on dilute ultracold atomic gases have given us access to highly tunable quantum systems. In particular, there have been substantial improvements in achieving different kinds of interaction between atoms. As a result, utracold atomic gases oer an ideal platform to simulate many-body phenomena in condensed matter physics, and engineer other novel phenomena that are a result of the exotic interactions produced between atoms. In this dissertation, I present a series of studies that explore the physics of dilute ultracold atomic gases in different settings. In each setting, I explore a different form of the inter-particle interaction. Motivated by experiments which induce artificial spin-orbit coupling for cold fermions, I explore this system in my first project. In this project, I propose a method to perform universal quantum computation using the excitations of interacting spin-orbit coupled fermions, in which effective p-wave interactions lead to the formation of a topological superfluid. Motivated by experiments which explore the physics of exotic interactions between atoms trapped inside optical cavities, I explore this system in a second project. I calculate the phase diagram of lattice bosons trapped in an optical cavity, where the cavity modes mediates effective global range checkerboard interactions between the atoms. I compare this phase diagram with one that was recently measured experimentally. In two other projects, I explore quantum simulation of condensed matter phenomena due to spin-dependent interactions between particles. I propose a method to produce tunable spin-dependent interactions between atoms, using an optical Feshbach resonance. In one project, I use these spin-dependent interactions in an ultracold Bose-Fermi system, and propose a method to produce the Kondo model. I propose an experiment to directly observe the Kondo effect in this system. In another project, I propose using lattice bosons with a large hyperfine spin, which have Feshbach-induced spin-dependent interactions, to produce a quantum dimer model. I propose an experiment to detect the ground state in this system. In a final project, I develop tools to simulate the dynamics of fermionic superfluids in which fermions interact via a short-range interaction.
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.
This thesis covers most of my work in the field of ultracold atoms loaded in optical lattices. It makes a route through the physics of cold atoms in periodic potentials starting from the simple noninteracting system and going into the many-body physics that describes the strongly correlated Mott insulator regime. Even though this thesis is a theoretical work all the chapters are linked either with experiments already done or with ongoing experimental efforts. This thesis can be divided into four different parts. The first part comprises chapters 1 to 3. In these chapters, after a brief introduction to the field of optical lattices I review the fundamental aspects pertaining to the physics of systems in periodic potentials. The second part deals with the superfluid weakly interacting regime where standard mean field techniques can be applied. This is covered in chapters 4 and 5. Specifically, chapter 4 introduces the discrete nonlinear Schroedinger equation (DNLSE) and uses it to model some experiments. In chapter 5 I go one step further and include the small quantum fluctuations neglected in the DNLSE by studying quadratic approximations of the Bose-Hubbard Hamiltonian. Chapters 6 to 8 can be grouped as the third part of the thesis. In them I adopt an effective action formalism, the so called two particle irreducible effective action (2PI) together with the closed time path (CTP) formalism to study far-from-equilibrium dynamics. The many-body techniques discussed in these chapters systematically include higher order quantum corrections, not included in the quadratic approximations of the Hamiltonian, which we show are crucial for a correct description of the quantum dynamics outside the very weakly interacting regime. Finally, chapter 9 to 11 are devoted to study the Mott insulator phase. In these chapters using perturbation theory I study the Mott insulator ground state and its excitation spectrum, the response of the system to Bragg spectroscopy, and propose a mechanism to correct for the residual quantum coherences inherent to the Mott insulator ground state. Even though small these are not ideal for the use of neutral atoms in optical lattice as a tool for quantum computation.
A new millennia ushered in novel experimental capabilities to investigate engineered quantum systems. Gases consisting of charge neutral atoms are cooled to the ground state and by exploiting the dipole interaction can be trapped by laser light. In ultracold atomic systems the tunneling and interactions are controllable, and the system is absent of disorder. Disorder and interactions in solid-state systems inhibits a complete understanding of transport phenomena. Theoretical investigations into non-equilibrium and equilibrium dynamics of fermions in discretized potentials and in continuum are presented. The affect of controllable, external system parameters is evaluated and exploited to reveal new properties of fermion dynamics. My thesis will begin with an overview of the physics associated with many-body fermion and the models used to obtain results, followed by my thesis research on how matter-waves can be manipulated by lattice geometry and tunneling strengths. Next, I present how dissipation affects current in a continuous ring, demonstrating rate-dependent hysteresis as the particles are driven by an artificial gauge field. After presenting the continuous ring, I will then focus on persistent current and corresponding quantum current fluctuations in a discretized ring. Lissajous curves are found in the system whose area strongly depends on the presence and absence of interactions. Lastly, I will discuss dynamical detection methods for one-dimensional topological insulators using ultracold fermions. All of the results are supported with discussion on how they can be realized using current experimental technology.
Predicting the quantum dynamics of charged matter interacting with dynamical gauge fields poses an outstanding challenge in theoretical physics. Lacking a generally applicable computational method, quantum simulators offer a promising alternative. \n \nIn this thesis, we contribute to the quantum simulation of high-energy physics, focusing on the platform of ultracold atoms. Using Wilson fermions, we propose to improve implementations of lattice gauge theories based on mixtures of cold atoms in optical lattices. Numerical benchmarks indicate that this makes the realization of Schwinger pair production feasible with current technology. Our proposal is modular and an elementary building is demonstrated experimentally. We further identify dynamical topological transitions, which we discovered in the massive Schwinger model, as a suitable target for quantum simulators. Defining a gauge-invariant order parameter, these transitions are shown to persist beyond weak coupling. \n \nIn the second part of this thesis, we develop a framework for analyzing quantum simulators in terms of experimentally accessible irreducible correlation functions at equal times. We verify this approach numerically for the sine-Gordon model in thermal equilibrium, quantum simulated by two tunnel-coupled superfluids. Finally, we apply our analysis to the non-equilibrium dynamics of a spinor Bose gas, revealing suppressed effective interactions in a strongly-correlated infrared regime.
The study of quantum dynamics is an increasingly important aspect of atomic, molecular and optical physics, as well as condensed-matter physics and other disciplines. It is essential in order to understand non-equilibrium phenomena which do not fit within the framework of equilibrium thermodynamics, such as many-body localization, dynamical quantum phase transitions, and non-equilibrium steady states in driven-dissipative systems. However, the exact calculation of dynamical observables is extremely difficult in general for several reasons, which makes it challenging for theoretical predictions and analyses to keep pace with rapidly advancing experiments. In this thesis I use and develop several computational methods to investigate the dynamics of quantum many-body systems, with the two goals of better understanding the evolution of local observables and correlation functions in spin and Fermion systems as well as creating and testing more effective numerical methods for solving quantum dynamics problems. I explore quantum dynamics by using Hamiltonian parameter quenches in several commonly-used models in AMO physics and elsewhere: the Fermi-Hubbard model, the XXZ model, and the transverse-Ising model. The first represents a system of short-range interacting Fermions, while the latter two represent spin systems with different spin-spin couplings and magnetic field terms. All of these models can be realized in current experiments. First, I derive an exact solution for an extreme interaction quench in the Fermi-Hubbard model, from a strongly-interacting initial state to a noninteracting final Hamiltonian. From this result I show that there are nontrivial transient connected correlations, despite the fact that the system has very high temperature and a noninteracting Hamiltonian. Those features would suppress such correlations in equilibrium, but the dynamical correlations are still significant. Then, I develop and examine two numerical methods for solving dynamics of the spin expectation values and correlation functions in the transverse-Ising and XXZ models. Both methods rely on approximating a thermodynamically large system with a much smaller and more computationally tractable cluster of lattice sites, or a series of clusters. The first method is a dynamical numerical linked cluster expansion (NLCE). NLCEs have already been used successfully on systems in thermodynamic equilibrium, and I extended the technique to directly compute the dynamics of an observable after a quantum quench. I review the general NLCE method, describe my dynamical variant of it, and compute the quench dynamics of spin and correlation observables in the transverse-Ising and XXZ models. The second method, which I term the adaptive-boundary cluster method, relies on recently refined bounds on information propagation in quantum systems to choose a cluster of lattice sites that minimizes the finite size error bound of a dynamical computation. I briefly review these bounds, and describe a new algorithm which uses them to construct a cluster whose boundary is chosen according to a related error estimate. Such an adaptive-boundary cluster has the property that there exists a computable upper bound on the finite size error of observables measured at the cluster center, and this error bound is nearly minimal when compared to that of other clusters of the same size. I then compute the quench dynamics of spin and correlation observables in the transverse-Ising model using the adaptive boundary clusters. I apply this method to both uniform and strongly disordered systems. To evaluate the effectiveness of both dynamical NLCE and adaptive-boundary clusters, I compare the results obtained using them to results from a standard and general purpose numerical method: exact solvers applied a rectangular cluster with periodic boundary conditions. Both of the new methods usually converge faster and are more accurate than the standard method for comparable cluster sizes. This thesis illustrates that choosing numerical methods which take advantage of the geometric and entanglement structure of the initial state and Hamiltonian can substantially enhance our ability to study quantum dynamics, by reducing the computational cost of numerically exact solutions. I conclude with a discussion of various avenues for future research to build on this work.