Ultracold atoms have revolutionized atomic and condensed matter physics. In addition to having clean, controllable Hamiltonians, ultracold atoms are near-perfect realizations of isolated quantum systems, in which weak environmental coupling can be neglected on experimental time scales. This opens new opportunities to explore these systems not just in thermal equilibrium, but out of equilibrium as well. In this dissertation, we investigate some properties of closed quantum systems, utilizing a combination of numerical and analytical techniques. We begin by applying full configuration-interaction quantum Monte Carlo (FCIQMC) to the Fermi polaron [1], which we use as a test bed to improve the algorithm. In addition to adapting standard QMC techniques, we introduce novel controlled approximations that allow mitigation of the sign problem and simulation directly in the thermodynamic limit. We also contrast the sign problem of FCIQMC with that of more standard techniques, focusing on FCIQMC’s capacity to work in a second quantized determinant space. Next, we discuss nonequilibrium dynamics near a quantum critical point, focusing on the one-dimensional transverse-field Ising (TFI) chain. We show that the TFI dynamics exhibit critical scaling, within which the spin correlations exhibit qualitatively athermal behavior. We provide strong numerical evidence for the universality of dynamic scaling by utilizing time-dependent matrix product states to simulate a nonintegrable model in the same equilibrium universality class. As this non-integrable model has been realized experimentally [2], we investigate the robustness of our predictions against the presence of open boundary conditions and disorder. We find that the qualitatively athermal correlations remain visible, although other phenomena such as even/odd effects become relevant within the finite size scaling theory. Finally, we investigate the properties of the integrable TFI model upon varying the strength of a non-integrable perturbation, attempting to understand the finite-size scaling behavior of the eigenstate thermalization crossover. We numerically estimate iii the crossover scale of this perturbation strength using exact diagonalization, finding that the crossover scale decreases strongly with system size, and analytically provide a lower bound on the finite-size scaling of this crossover. However, we are unable to solve for the thermodynamic limit of this crossover; doing so will require larger systems and more comprehensive theory.
Ultracold atoms have revolutionized atomic and condensed matter physics. In addition to having clean, controllable Hamiltonians, ultracold atoms are near-perfect realizations of isolated quantum systems, in which weak environmental coupling can be neglected on experimental time scales. This opens new opportunities to explore these systems not just in thermal equilibrium, but out of equilibrium as well. In this dissertation, we investigate some properties of closed quantum systems, utilizing a combination of numerical and analytical techniques. We begin by applying full configuration-interaction quantum Monte Carlo (FCIQMC) to the Fermi polaron [1], which we use as a test bed to improve the algorithm. In addition to adapting standard QMC techniques, we introduce novel controlled approximations that allow mitigation of the sign problem and simulation directly in the thermodynamic limit. We also contrast the sign problem of FCIQMC with that of more standard techniques, focusing on FCIQMC’s capacity to work in a second quantized determinant space. Next, we discuss nonequilibrium dynamics near a quantum critical point, focusing on the one-dimensional transverse-field Ising (TFI) chain. We show that the TFI dynamics exhibit critical scaling, within which the spin correlations exhibit qualitatively athermal behavior. We provide strong numerical evidence for the universality of dynamic scaling by utilizing time-dependent matrix product states to simulate a nonintegrable model in the same equilibrium universality class. As this non-integrable model has been realized experimentally [2], we investigate the robustness of our predictions against the presence of open boundary conditions and disorder. We find that the qualitatively athermal correlations remain visible, although other phenomena such as even/odd effects become relevant within the finite size scaling theory. Finally, we investigate the properties of the integrable TFI model upon varying the strength of a non-integrable perturbation, attempting to understand the finite-size scaling behavior of the eigenstate thermalization crossover. We numerically estimate iii the crossover scale of this perturbation strength using exact diagonalization, finding that the crossover scale decreases strongly with system size, and analytically provide a lower bound on the finite-size scaling of this crossover. However, we are unable to solve for the thermodynamic limit of this crossover; doing so will require larger systems and more comprehensive theory.
Many strongly correlated quantum systems are difficult to study numerically because of the exponential growth of the Hilbert space combined with the failure of numerical methods in important parameter regimes. This includes the long-standing problems of high-temperature superconductivity and frustrated magnetism. A promising solution to this dilemma are quantum simulators, well-controlled (synthetic) quantum systems, emulating model Hamiltonians not only of condensed matter systems, where parameters can be tuned at will and model complexity can be added step by step. This dissertation presents numerical simulation studies that are motivated by quantum simulators. Using large-scale quantum Monte Carlo computer simulations and mean-field theory the semiquantitative ground state phase diagram of a Penning ion trap quantum simulator, realizing a two-dimensional model of quantum magnetism with long-range interactions, is obtained. We study also thermal and quantum phase transitions in one-dimensional transverse-field Ising models with long-range interactions, which exhibit a Kosterlitz-Thouless transition and universality classes with continuously varying critical exponents. The second part of this thesis is concerned with the Fermi-Hubbard model, which has been realized recently with ultracold atoms in optical lattices, both for repulsive and attractive interactions. Equipped with a quantum gas microscope, which allows single-site single-atom detection, these experimental setups can collect histograms of the particle number from repeated projective measurements. This dissertation presents a method to compute such full counting statistics for quadratic operators such as particle number or magnetization in the framework of determinantal quantum Monte Carlo. From an analysis of the full counting statistics of the particle number for different subsystem sizes the size of preformed pairs or Cooper pairs in the attractive Hubbard model is inferred. Furthermore, excellent agreement of the numerically computed probability distribution of the staggered magnetization with a recent quantum gas experiment demonstrates that current experiments are capable of resolving the differences between the Hubbard model and its limiting Heisenberg model. An extensive method section reviews the stochastic series expansion quantum Monte Carlo method for Ising models in a transverse field, which is free of the sign problem for arbitrary interactions, with special emphasis on the efficient simulation of long-range interactions. Also a self-contained introduction to the determinantal quantum Monte Carlo method is provided.
In this Thesis we discussed several issues related to the dynamics of isolated quantum many-body systems. The brief overview of Chapter 1 highlighted a series of open fundamental theoretical questions which are crucial for the investigation of nonequilibrium phenomena. Among them, we mainly focused on the emergence of thermal behaviors and on the characterization of dynamical phase transitions and of the associated universal dynamics. To this end, we analyzed some paradigmatic models, such as the Fermi-Hubbard model, the quantum Ising chain, and the O(N) model. -Chapter 2: Absence of thermalization in a Fermi liquid. In this Chapter we report two results that question the common expectation of thermalization in a Landau-Fermi liquid after an interaction quench. We first consider the perturbative expansion in the interaction strength of the long-wavelength structure factor S(q) and show that it does not satisfy the hypothesis that steady-state averages correspond to thermal ones. In particular, S(q) has an analytical component ∼ const. + O(q^2) for q → 0, compatible with thermalization, but it retains also a nonanalytic term ∼ |q| characteristic of a Fermi liquid at zero-temperature. In real space, this nonanalyticity turns in a power-law decay of the density-density correlations, in contrast with the exponential decay associated to thermalization. We next analyze the low-density case, where one can obtain results valid at any order in interaction but at leading in the density, and find that in the steady state the momentum distribution jump at the Fermi surface is strictly finite, though smaller than at equilibrium. -Chapter 3: Nonadiabatic stationary behavior in a driven Ising chain. In this Chapter we discuss the emergence of a nonadiabatic behavior in the dynamics of the order parameter in a low-dimensional system subject to a linear ramp of one of its parameters within a gapped phase, which should be the most favorable situation for an adiabatic evolution. We study in details the dynamics of the spontaneous magnetization m_x(t) in a quantum Ising chain after a linear variation in time of the transverse field within the ordered phase. In particular, focusing on the value of the order parameter at the end of the ramp m_x(τ), we show that the smaller the switching rate of the transverse field is the closer m_x(τ) gets to its ground state value. Nonetheless, this apparently small correction to adiabaticity eventually leads to the disruption of the order exponentially fast in the subsequent time evolution, no matter how slowly the ramp is performed. -Chapter 4: Aging and coarsening in isolated quantum systems after a quench. This Chapter analyzes the nonequilibrium dynamics of an isolated quantum system after a sudden quench to the dynamical critical point, where the emergence of scaling and universal exponents is expected, due to the absence of time and length scales. We explore these features for a bosonic interacting scalar field theory with O(N) symmetry in the large-N limit, where the model is exactly solvable and the exponents and scaling forms of the relevant two-time correlation functions can be analytically calculated. Moreover, we provide evidence of the emergence of a dynamical scaling behavior also for quenches below the dynamical critical point, associated with coarsening. We find that the latter case is characterized by the same scaling functions as those describing the critical case, yet with different exponents. -Chapter 5: Dynamical transitions and statistics of excitations. In this Chapter we study the dynamics of the O(N) model (introduced in the previous Chapter) resulting from a different protocol: a linear ramp of one of its parameters. We find that the presence of a dynamical phase transition, as well as its critical properties, are robust against the change of the protocol. We show that a characterization based on the critical dimensions and exponents would suggest that the dynamical phase transition is analogous to the equilibrium thermal one. However, its nonequilibrium nature becomes evident analyzing the statistics of excitations produced in the ramp process. In particular, the critical properties are encoded in the fluctuations in the number of excitations, which display qualitatively different behaviors depending on the ramp being performed above, at, or below the dynamical critical point. These behaviors bear no dependence on the duration of the protocol.
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.
Understanding how quantum matter behaves when driven out of equilibrium is one of the key focuses in quantum physics. Thanks to impressive progress in the control and precision achieved in quantum synthetic matter over the past decades, the nonequilibrium quantum many-body physics has become one of the most active research areas today, especially after the experimental realization of Bose-Einstein condensates and optical lattices, which allows us to directly observe and study nonequilibrium quantum matter with great accuracy and controllability. In this dissertation, I explore the rich landscape of nonequilibrium quantum many-body physics and how quantum phase transitions, both symmetry-breaking and topological, can be extended to the nonequilibrium setting. In the first part of the dissertation, I focus on spinor Bose-Einstein condensates as an isolated quantum many-body system, and reveal their various dynamical behaviors, including quantum collapse and revivals, thermalization and nonthermal equilibration with no revival even though the system has finite degrees of freedom. In contrast to typical integrable systems, which usually do not thermalize, we find that spinor condensates have a parameter range in which the system thermalizes via the Eigenstate Thermalization Hypothesis (ETH). We show that this observation is linked to the presence of rare nonthermal states whose fraction vanishes with system size, and contributes to the notion of thermalization via weak ETH. Next, I explore a dynamical process that is complementary to thermalization in isolated quantum systems: information scrambling, which could be probed via out-of-time-order correlators (OTOC). I propose a nonintegrable, disordered and quasi-1D spin model, the ladder-XX model, for a feasible detection of information scrambling in a cold atom simulator. This chapter poses a fundamental question: "What are the signatures of quantum phases and phase transitions in isolated interacting systems driven out-of-equilibrium?" I study the ladder-XX model in both clean and disordered potentials, and characterize different nonequilibrium phases, i.e., ergodic and many-body localized, of the model based on the decay properties of OTOCs. Emergent light cone shows sublinear behavior, while the butterfly cones drastically differ from the light cone by demonstrating superlinear spread of information with a velocity that is bounded by the light cone velocity. In the second part of the dissertation, I continue to search for answers to the question posed above, however this time with a particular focus on symmetry-breaking and topological quantum phase transitions. I pin down a universal mechanism underlying the relation between information scrambling at any temperature and quantum phases at low temperatures. Our method points to key ingredients to dynamically detect long-range order in gapped phases through OTOCs for symmetry-breaking quantum phase transitions and Z2 topological order associated with Majorana zero modes localized at the edges. Our results pave the way to an intriguing observation that phases of quantum matter could protect the information from scrambling and thermalization, even when the system is interacting and nonintegrable. Finally, I explore and propose utilizing short-time transient temporal regimes and single-site probes to detect the phases and phase transitions in quantum matter. These studies reveal a dynamical crossover and a dynamical phase transition, respectively for periodic and open-boundary chains. In both cases, a nonequilibrium scaling law appears in the vicinity of the crossover/transition with associated exponents that differ from the analytical predictions for long times. Feasibility of detecting such dynamical criticality in experimental systems are discussed.
Strongly correlated quantum systems display a myriad of complex, fascinating behaviour. The competition between charge, spin, orbital, or lattice degrees of freedom and their respective energy scales gives rise to rich and complex emergent phases. More recently, the non-equilibrium properties of correlated quantum many-body systems have come into the focus of active research, with remarkable advances. While the advent of tensor network based simulations on the theoretical side has allowed access to the full time-evolution of the quantum state as it traverses the Hilbert space, the experimental progress in the ultracold atom community has allowed physicists to enter the paradigm of analog quantum simulation. Their tunability and isolation from the environment makes them ideal platforms to emulate and study open problems in condensed matter physics and go beyond the simulational capabilities of state of the art classical simulations. In this thesis, we explore the dynamics of correlated Fermi gases. We develop and present different facets of applicability of radiofrequency modulation techniques to drive the system out of equilibrium. We begin by considering a homogeneous, three-dimensional Fermi gas in the BCS-BEC crossover and develop an off-resonant radiofrequency transfer scheme to excite the Higgs mode of the superfluid. In a subsequent study we investigate the tunability of the coherence between Cooper pairs by tuning the duration of interaction ramps of the internal interaction strength. These works highlight the possibility to externally tune and stabilise complex quantum many-body states away from equilibrium. Using quasi-exact matrix product state simulations, we study the response of an attractively interacting, one-dimensional Fermi-Hubbard model to weak radiofrequency perturbations. We reveal the emergence of two distinct dynamical regimes in the time evolution and are able to relate it back to the underlying excitation spectrum of the system. Finally, we explore the possibility to perform quantum quenches through radiofrequency π-pulses, achieving nearly complete population transfer between different interacting fermionic states. We reveal non-trivial dynamical effects in the pair correlation of the final state, as well as the excitation of a collective trap mode of the system.
In this thesis we examine universal scaling properties of strongly-correlated systems near and far from equilibrium. We discuss quantum phase transitions at vanishing temperature, multicritical and dynamic critical behavior near thermal equilibrium, and scaling properties of nonequilibrium steady states. We employ nonperturbative methods including the functional renormalization group as well as Monte Carlo simulations. A general outline of the functional renormalization group is given in the introductory chapters. \n \nIn the first part of this thesis, we investigate spinless fermions on the honeycomb lattice interacting via short-range repulsive interactions. Such a system can be seen as a simple model for suspended graphene. The short-range interactions control the ground state properties of the system that may lead to a chiral phase transition from the semimetal to the charge density wave (CDW)/Kekulé ordered state. We determine the universal scaling properties at the chiral transition, and establish the presence of large anomalous dimensions indicating the importance of strong fluctuations. \n \nThe competition of two nonvanishing order parameters and their corresponding multicritical behavior are investigated in the subsequent chapter. We characterize the bicritical and tetracritical behavior in the purely bosonic O(N_1) + O(N_2) symmetric model and comment on possible applications to condensed-matter and high-energy physics. \n \nIn the following chapter we discuss the long-time relaxational behavior at criticality of an order parameter with O(N) symmetry coupled to an additional conserved density. We find an anomalous diffusion phase with new dynamic scaling properties. Using the functional renormalization group we determine the complete dynamic critical behavior of the model in 2 < d < 4 dimensions and compare our results to experiments. \n \nFinally, we investigate the scaling properties of stationary states far from equilibrium. At the example of the one-dimensional Burgers’ equation we develop a novel approach to hydrodynamic turbulence using lattice Monte Carlo methods. We apply these techniques to determine the statistical properties of small-scale fluctuations in this model and identify the anomalous scaling behavior. \n