In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.
In the first part of the Thesis we mostly concentrate on spectral properties of strongly correlated systems and on their equilibrium properties. This is accomplished by the general concept of imaginary-time dynamics which we apply to a number of different problems in which different strengths of this approach emerge. In Chapter 1 we introduce the formalism that allows for a connection between the quantum and the classical worlds. The connection is established by means of the imaginary-time quantum evolution which, under certain circumstances, is shown to be equivalent to a classical stochastic process. It is further shown that exact static and spectral properties of correlated systems can be obtained when this mapping is feasible. The relationship between the imaginary-time dynamics in different frameworks such as the path-integral and the perturbative one is also underlined. In Chapter 2 we present a specific implementation of the general ideas previously presented. In particular we introduced an extension to lattice systems of the Reptation Monte Carlo algorithm [30] which benefits of a sampling scheme based on directed updates. Specific improvements over the existing methodologies consist in the unbiased evaluation of the imaginary-time path integrals for bosons and a systematic scheme to improve over the Fixed-node approximation for fermions. Applications to the Hubbard and the Heisenberg models are presented. In Chapter 3 we demonstrate the application of the imaginary-time dynamics to the exact study of spectral properties. Subject of our attention is a highly anharmonic and correlated quantum crystal such as Helium 4 at zero temperature.[33] Concerning this system, we have obtained the first ab-initio complete phonon dispersion in good agreement with neutron spectroscopy experiments. Moreover, we have also studied the density excitations of solid helium in a region of wave-vectors in between the collective (phonon) and the single-particle regimes, where the presence of residual coherence in the dynamics shows analogies between the highly anharmonic crystal and the superfluid phase. In Chapter 4 we introduce a novel method, based on the imaginary-time dynamics, to obtain unbiased estimates of fermionic properties.[34] By means of this method and of a very accurate variational state, we provide strong evidence for the stability of a saturated ferromagnetic phase in the high-density regime of the two-dimensional infinite-U Hubbard model. By decreasing the electron density, we observe a discontinuous transition to a paramagnetic phase, accompanied by a divergence of the susceptibility on the paramagnetic side. This behavior, resulting from a high degeneracy among different spin sectors, is consistent with an infinite-order phase transition scenario. In Chapter 5 the use of imaginary-time dynamics in the context of finite-temperature response functions is highlighted. As an application, we study an intriguing quantum phase featuring both glassy order and Bose-Einstein condensation. [35] We introduce and validate a model for the role of geometrical frustration in the coexistence of off-diagonal long range order with an amorphous density profile. The exact characterization of the response of the system to an external density perturbation is what allows here to establish the existence of a spin-glass phase. The differences between such a phase and the otherwise insulating "Bose glasses" are further elucidated in the Chapter. In the second part of the Thesis we focus our attention on the dynamics of closed systems out of equilibrium. This is accomplished by both non-stochastic exact methods for the dynamics and the introduction of a novel time-dependent Variational Monte Carlo scheme. In Chapter 6 exact diagonalization schemes and renormalization-based methods for one-dimensional systems are introduced. We identify key phenomenological traits resulting from the many-body correlation in closed systems driven sufficiently away from equilibrium.[31] We provide evidences that the dynamics of interacting lattice bosons away from equilibrium can be trapped into extremely long-lived inhomogeneous metastable states. The slowing down of incoherent density excitations above a threshold energy, much reminiscent of a dynamical arrest on the verge of a glass transition, is identified as the key feature of this phenomenon. In Chapter 7 we present an extension to dynamical properties of the Variational Quantum Monte Carlo method.[32] This is accomplished by introducing a general class of time-dependent variational states which is based on the mapping of the many-body dynamics onto an instantaneous ground-state problem. The application of the method to the experimentally relevant quantum quenches of interacting bosons reveals the accuracy and the reliability of the introduced numerical scheme. We indeed obtain for the first time a consistent variational description of the approach to the equilibrium of local observables and underline the origin of the metastability and glassy behavior previously identified. In the very last part we draw our conclusions and show some possible paths for stimulating future research.
The investigation of quantum many-particle systems is a central goal of modern condensed-matter physics. Complex physical phenomena such as high-temperature superconductivity, heavy fermion behavior, and correlation-driven metal-insulator transitions are central aspects of this particular field of science. The Hubbard model is a paradigm of correlated-electron physics. It is one of the simplest quantum-mechanical lattice models capable of capturing relevant physical aspects of strongly correlated electron systems. The Hubbard model also plays an important role in the explanation of many-body phenomena observed in ultra-cold atoms trapped in optical lattices. In this thesis, we study the three-dimensional Hubbard model using quantum Monte Carlo simulations (QMC). Our aim is a precise numerical study of its properties in the thermodynamic limit. However, a direct investigation of lattice models in the thermodynamic limit is generally impossible because QMC simulations cannot treat infinite lattice sizes. Therefore, approximative schemes that allow calculations directly in the thermodynamic limit are advantageous. A widely applied approximation of this kind is the dynamical mean-field theory (DMFT) which maps the infinite-lattice problem onto a single-site impurity model embedded in a mean field. We apply cluster mean-field theories which represent a systematic extension of the DMFT by expanding the single impurity to a finite cluster. The second focus of this thesis is the calculation of spectral properties of the Hubbard model. QMC methods map the quantum-mechanical lattice model on a classical one at the expense of an additional imaginary time dimension. For this reason, QMC algorithms can only provide dynamical data on the imaginary time or frequency axis. The necessary analytic continuation to the physically relevant real axis has proven to be a difficult problem that has to be approached by advanced data analysis tools. The maximum-entropy method (MEM) is the standard technique to handle this problem. We apply the MEM to the problem of extracting momentum-resolved single particle spectra from QMC simulations of the 3D Hubbard model. Furthermore, we discuss and expand algorithms that were recently proposed as alternatives to the standard MEM approach.
Strongly correlated systems, i.e., quantum materials for which the interactions between its constituents are strong, are good candidates for the development of applications based on quantum-mechanical principles, such as quantum computers. Two paradigmatic models of strongly correlated systems are heavy-fermionic systems and one-dimensional spin-12 systems, with and without quenched disorder. In the past decade, improvement in computational methods and a vast enhancement in computational power has made it possible to study these systems in a a non-perturbative manner. In this thesis we present state-of-the-art numerical methods to investigate the properties of strongly correlated systems, and we apply these methods to solve a couple of selected problems in quantum condensed matter theory. We start by revisiting the phase diagram of the Falicov-Kimball model on the square lattice which can be considered as a heavy-fermionic systems. This model describes an interplay between conduction electrons and heavy electrons and reveals several distinct metal-insulator phase transitions. Using a lattice Monte-Carlo method, we study the transport properties of the model. Our analysis describes the role of temperature and interaction strength on the metal-insulator phase transitions in the Falicov-Kimball model. The second part of the thesis investigate the spatial structure of the entanglement in ground and thermal statesof the transverse-field Ising chain. We use the logarithmic negativity as a measure for the entanglement between two disjoint blocks. We investigate how logarithmic negativity depends on the spatial separation between two blocks, which can be viewed as the entanglement analog of a spatial correlation function. We find sharp entanglement thresholds at a critical distance beyond which the logarithmic negativity vanishes exactly and thus the two blocks become unentangled. Our results hold even in the presence of long-ranged quantum correlations, i.e., at the system’s quantum critical point. Using Time-Evolving Block Decimation (TEBD), we explore this feature as a function of temperature and size of the two blocks. We present a simple model to describe our numerical observations. In the last part of this thesis, we introduce an order parameter for a many-body localized spin-glass (MBL-SG) phase. We show that many-body localized spin-glass order can also be detected from two-site reduced density matrices, which we use to construct an eigenstate spin-glass order parameter. We find that this eigenstate spin-glass order parameter captures spin-glass phases in random Ising chains, both in many-body eigenstates as well as in the nonequilibrium dynamics, from a local in time measurement. We discuss how our results can be used to observe MBL-SG order within current experiments in Rydberg atoms and trapped ion systems.
In this thesis, strongly correlated quantum many-body systems in equilibrium and in out-of-equilibrium situations are investigated. This is done by applying and developing well established numerical methods. The focus of the thesis lies in the development and application of the density matrix renormalization group method (DMRG) to quantum many-body systems out of equilbrium. In this thesis, in addition to the DMRG, methods for the exact diagonalization of the Hamiltonian of the system, like the Lanczos- or the Jacobi-Davidson method, are treated. An extension of the Lanczos method makes it possible to treat the time evolution of strongly correlated quantum systems with an accuracy comparable to machine precision. This method is the basis for a possible extension of the DMRG for the treatment of systems out of equilibrium, the so-called adaptive time-dependent DMRG ("adaptive t-DMRG"). A second variant of the adaptive t-DMRG uses the Suzuki-Trotter decomposition of the time-evolution operator. An error analysis demonstrates that both methods, for a suitable choice of control parameters, have errors < 1% at the end of the time evolution in complicated quantities like, e.g., the momentum distribution function. We apply these numerical methods to investigate the quantum critical behavior of a variant of the Hubbard model and to treat two non-equilibrium situations of current interest. Extensive use of the DMRG makes it possible to clarify the quantum-critical behavior of the ionic Hubbard model; in particular, the numerical results demonstrate quite clearly that the correct scenario has two critical points; at the first critical point, only the charge degrees of freedom are critical, while at the second one only the spin degrees of freedom are critical. Quite surprisingly, the Mott-insulator phase in the strong-coupling regime shows a divergent electrical susceptibility. Next, the dynamics of a system of so-called soft-core bosons is treated. The particles initially are trapped in a deep box-potential and the dynamics is investigated after releasing them from this trapping potential. Similar to exact results obtained for so-called hard-core bosons, (quasi-)coherent matter waves emerge, demonstrating the possibility of realizing an atom laser in experiments on optical lattices. We find that the wave vector of the emerging matter wave can be tuned by changing the strength of the interaction between the particles. The second non-equilibrium situation is a so-called "quantum-quench" of a system of strongly correlated fermions, i.e., the evolution of the system after suddenly changing an intrinsic parameter like, e.g., the interaction strength between the particles, is investigated. In this case, we are mainly interested in the long-time behavior; in particular, the question arises whether, relying on Boltzmann's ergodic hypothesis, relaxation to a thermal state is obtained. We find that, in general, non-thermal final states are reached, which can be described by generalized Gibbs-Boltzmann ensembles.
We study non-equilibrium and strongly-correlated dynamics in two contexts. We begin by analysing quantum many-body systems out of equilibrium through the lens of cold atomic impurities in Bose gases. Such highly-imbalanced mixtures provide a controlled arena for the study of interactions, dissipation, decoherence and transport in a many-body quantum environment. Specifically we investigate the oscillatory dynamics of a trapped and initially highly-localised impurity interacting with a weakly-interacting trapped quasi low-dimensional Bose gas. This relates to and goes beyond a recent experiment by the Inguscio group in Florence. We witness a delicate interplay between the self-trapping of the impurity and the inhomogeneity of the Bose gas, and describe the dissipation of the energy of the impurity through phononic excitations of the Bose gas. We then study the transport of a driven, periodically-trapped impurity through a quasi one-dimensional Bose gas. We show that placing the weakly-interacting Bose gas in a separate periodic potential leads to a phononic excitation spectrum that closely mimics those in solid state systems. As a result we show that the impurity-Bose gas system exhibits phonon-induced resonances in the impurity current that were predicted to occur in solids decades ago but never clearly observed. Following this, allowing the bosons to interact strongly, we predict the effect of different strongly-correlated phases of the Bose gas on the motion of the impurity. We show that, by observing the impurity, properties of the excitation spectrum of the Bose gas, e.g., gap and bandwidth, may be inferred along with the filling of the bosonic lattice. In other words the impurity acts as a probe of its environment. To describe the dynamics of such a strongly-correlated system we use the powerful and near-exact time-evolving block decimation (TEBD) method, which we describe in detail. The second part of this thesis then analyses, for the first time, the performance of this method when applied to simulate non-equilibrium classical stochastic processes. We study its efficacy for a well-understood model of transport, the totally-asymmetric exclusion process, and find it to be accurate. Next, motivated by the inefficiency of sampling-based numerical methods for high variance observables we adapt and apply TEBD to simulate a path-dependent observable whose variance increases exponentially with system size. Specifically we calculate the expected value of the exponential of the work done by a varying magnetic field on a one-dimensional Ising model undergoing Glauber dynamics. We confirm using Jarzynski's equality that the TEBD method remains accurate and efficient. Therefore TEBD and related methods complement and challenge the usual Monte Carlo-based simulators of non-equilibrium stochastic processes.
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
The theory of statistical mechanics provides a powerful conceptual framework within which the relevant (macroscopic) features of systems at equilibrium can be described. As there is currently no equivalent capable of encompassing the much richer class of non-equilibrium phenomena, research in this direction proceeds mainly on an instance-by-instance basis. The aim of this Thesis is to describe in some detail three such attempts, which involve different dynamical aspects of classical and quantum systems. As summarised below, each of the last three Chapters of this document delves into one of these different topics, while Chapter 2 provides a brief introduction on the study of non-equilibrium dynamics. In Chapter 3 we investigate the purely relaxational dynamics of classical critical ferromagnetic systems in the proximity of surfaces, paying particular attention to the effects that the latter induce on the early stages of the evolution following an abrupt change in the temperature of the sample. When the latter ends close enough to the critical value which separates the paramagnetic from the ferromagnetic phase, it effectively introduces a temporal boundary which can be treated as if it were a surface. Within this picture, we highlight the emergence of novel effects near the effective edge formed by the intersection of the two spatial and temporal boundaries. Our findings are apparently in disagreement with previous predictions which were based on the assumption that the presence of such an edge would not affect the scaling behaviour of observables; in order to explain this discrepancy, we propose an alternative for the original power-counting argument which, at least, correctly predicts the emergence of novel field-theoretical divergences in our one-loop calculations. We show that said singularities are associated with the scaling at the edge. Moreover, by encoding our findings in a boundary renormalisation group framework, we argue that the new predicted behaviour represents a universal feature associated to the short-distance expansion of the order parameter of the transition near the edge; we also calculate explicitly its anomalous dimension at the first-order in a dimensional expansion. As a qualitative feature, this anomalous dimension depends on the type of phase transition occurring at the surface. We exploit this fact in order to provide numerical support to our predictions via Monte Carlo simulations of the dynamical behaviour of a three-dimensional Ising model. The main results reported in Chap. 3 have appeared in Ref. [1]. In Chapter 4 we revisit the Euclidean mapping to imaginary times which has been recently proposed [2, 3] as an alternative for approaching the problem of quantum dynamics following a quench. This is expected to allow one to reformulate the original problem as a static one confined in a film geometry. We show that this interpretation actually holds only if the initial state of the dynamics is pure. Statistical mixtures, instead, intertwine the effects due to the two boundaries, which therefore cannot be regarded as being independent. We emphasize that, although the aforementioned reinterpretation as a confined static problem fails, one is still able, in principle, to write down and solve the corresponding equations. We also discuss in some detail the relation between this approach and the real-time field-theoretical one which makes use of the two-time Keldysh contour. For this purpose, we study the analytical structure of relevant observables — such as correlation functions — in the complex plane of times, identifying a subdivision of this domain into several sectors which depend on the ordering of the imaginary parts of the involved time coordinates. Within each of these subdomains, the analytic continuation to the real axis provides in principle a different result. This feature allows one to reconstruct from the Euclidean formalism all possible non-time-ordered functions, which in particular include all those which can be calculated via the Keldysh two-time formalism. Moreover, we give a prescription on how to retrieve response functions, discussing some simple examples and rationalising some recent numerical data obtained for one of these observables in a one-dimensional quantum Ising chain [4]. We also highlight the emergence of a light-cone effect fairly similar to the one previously found for correlation functions [2], which therefore provides further confirmation to the fact that information travels across the system in the form of the entanglement of quasi-particles produced by the quenching procedure. We have reported part of this analysis in Ref. [5]. Chapter 5 presents part of our recent work on effective relaxation in quantum systems following a quench and on the observed prethermalisation. We analyse the effects caused by the introduction of a long-range integrability-breaking interaction in the early stages of the dynamics of an otherwise integrable quantum spin chain following a quench in the magnetic field. By employing a suitable transformation, we redefine the theory in terms of a fully-connected model of hard-core bosons, which allows us to exploit the (generically) low density of excitations for rendering our model exactly solvable (in a numerical sense, i.e., by numerically diagonalising an exact matrix). We verify that, indeed, as long as the parameters of the quench are not too close to the critical point, the low-density approximation captures the dynamical features of the elementary operators, highlighting the appearance of marked plateaux in their dynamics, which we reinterpret as the emergence of a prethermal regime in the original model. As expected, the latter behaviour is reflected also on extensive observables which can be constructed as appropriate combinations of the mode populations. For these quantities, the typical approach to the quasi-stationary value is algebraic with exponent a ≈ 3, independently of the size of the system, the strength of the interaction and the amplitude of the magnetic field (as long as it is kept far from the critical point). The plateaux mentioned above last until a recurrence time — which can be approximately identified with tR ≈ N/2 for single modes and t′R ≈ N/4 for extensive quantities — after which quantum oscillations due to the finite size of the chain reappear. Our procedure allows us to shed some light over prethermal features without having to considerably limit the size of the system, which we can choose to be quite large, as we discuss in Ref. [6].