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.
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 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 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.
The physics of strongly interacting particles is often extraordinarily rich, leading to the \textit{emergence} of fundamentally new phases and phenomena that cannot exist in the absence of correlations. This thesis is devoted to the study of different classes of quantum systems with such strong correlations. Programmable quantum simulators based on Rydberg atom arrays have recently emerged as versatile platforms to probe quantum many-body phases and their dynamics. In Chapter 2, we start with theories of quantum phase transitions (QPTs) in one-dimensional chains of Rydberg atoms. We demonstrate the existence of a novel QPT in the universality class of the Z3 chiral clock model, calculate its critical exponents, and discuss their measurement with the quantum Kibble-Zurek mechanism. Proceeding to two dimensions, Chapter 3 describes the zero-temperature phase diagram of a square-lattice Rydberg atom array, revealing several experimentally realizable density-wave-ordered phases and QPTs. For a different geometry, namely, the kagome lattice studied in Chapter 4, we identify a regime that constitutes a likely candidate for hosting a phase with long-range quantum entanglement and topological order. Based on mappings to dimer models and gauge theories, we argue for the existence of a quantum spin liquid phase that has proved elusive in solid-state materials. In the next part of this thesis, we focus on similarly strongly correlated phases in the realm of quantum materials. Motivated by recent observations of an enhanced thermal Hall response in the high-temperature cuprate superconductors, in Chapter 5, we study the thermal Hall conductivity, κxy, of chiral spin-liquid ans\"{a}tze on the square lattice, focusing on the contribution from neutral (bosonic or fermionic) spinons. Thereafter, going beyond mean-field theory, we consider the corrections due to emergent gauge fields, which yield a response of opposite sign to that of the fermionic matter. We additionally consider the contribution of phonons to κxy and establish its relation to the spinon Hall viscosity. Thermal Hall measurements also leave intriguing fingerprints in Kitaev materials like α-RuCl3, and in this context, Chapter 6 investigates the unquantized κxy obtained in a gapless spinon Fermi surface state.
Materials with strongly correlated electrons show a multitude of unconventional phenomena, such as superconductivity or magnetism, that cannot be explained in a single-particle picture and necessitate to take many-body effects into account. A route to explain the origin of these phenomena consists in constructing effective models for these materials that are then solved, often using numerical techniques. Paradigmatic models for two different classes of strongly correlated electron systems are investigated in this thesis: The Hubbard model is a paradigmatic model to describe the Mott metal-insulator transition (MIT) and the Kondo lattice model describes the low-energy physics in heavy fermion compounds. Both models are analyzed numerically using different variations of a quantum cluster technique, the variational cluster approximation (VCA). The first part focuses on the MIT, as it is investigated experimentally in quasi-two-dimensional charge-transfer salts. Its nature and associated universality class are still heavily discussed and both experiments and theoretical approaches come to different conclusions. In this thesis, the Mott transition is investigated on an anisotropic two-dimensional Hubbard model at half-filling. By using a control parameter, which induces this anisotropy, strong evidence for Mott quantum criticality is found in weakly coupled Hubbard chains. The results at zero temperature show that the second-order critical end point Tc of the interaction-driven metal-insulator transition can be tuned down to zero at strong anisotropy. Further results for the antiferromagnetic phase suggest a similar picture and motivate adding the anisotropy as a new axis in the phase diagram to account for a low-temperature critical end point of the Mott transition. The second part of this thesis focuses on unconventional phases that emerge in a paradigmatic model for heavy fermion systems, the Kondo lattice model. Using the VCA, the two-dimensional Kondo lattice model is investigated in the paramagnetic phase as well as in phases with broken symmetry. An antiferromagnetic phase at weak coupling and a phase with Kondo-singlet formation at strong coupling are analyzed. Within the antiferromagnetically ordered region, two phases with different Fermi-surface topology, which are separated by a discontinuous transition, are identified. The model is also tested for s-wave superconductivity as found recently using dynamical mean-field theory, but no indications for robust local pairing are seen within VCA. Instead, nodal d-wave superconductivity is found and analyzed in a large range of couplings and electron fillings, and its interplay with antiferromagnetic order at weak coupling is discussed. This motivates further studies on extended models that might allow for a close comparison with experiments on heavy fermion systems. Understanding the formation of unconventional phases in strongly correlated electron systems is not only of fundamental interest, but it is also essential for designing new functional materials. Although the results of approximate cluster techniques have to be handled with care, the new aspects presented in this thesis pave the way for further investigations in this direction.
The thesis contains theoretical study of number of systems in low dimension, which are related to quantum Ising model. The main results are emergence of Majorana fermions in condensed matter system. We discuss an array of cavity-QED system, describing atom-photon (spin-boson) interacting problem, namely Rabi lattice model, and calculated its phase diagram. More importantly we showed that in any dimension, the system has Ising dynamics in strong atom-photon coupling limit. For 1D, a special relation between end-to-end spin-spin correlation is presented, which is very much different form the bulk spin-spin correlation. The entanglement between cavities is calculated and shown to have maximum value at the phase transition point. Other problems we studied were the frustrated quantum Ising model and its phase diagram as well as existence of Majorana like edge modes in different phases of the model. The fermionic version of the frustrated quantum Ising problem is also presented. On a different but related note, we present some beautiful results on Falicov Kimball (FK) model, a fermionic model similar to Hubbard model, on a square lattice. Using a canonical representation for electrons (invented by Brijesh Kumar), the FK model is mapped to classical Ising model (at infinite Hubbard interaction limit) on 2D lattice. We performed classical Monte Carlo simulation to study magnetization as a function of Hubbard interaction and also as a function of temperature. The motivation was to study (quantum) charge fluctuation (thus, quantum fluctuation of top of classical Ising thermodynamics) due to finite Hubbard interaction term. We also devote a chapter for basic numerical techniques, such as Density matrix Renormalization Group method, Exact diagonalization for spin system and classical Monte Carlo simulation (for classical 2D Ising model), which we have used in this thesis.
The study of collective behaviour in many-body systems often explores fundamentally new ideas absent from the mere constituents of such a system. A paradigmatic model for these studies is the spin-1/2 XXZ chain and its fermionic equivalent. This thesis can be broadly divided into the study of two fundamental aspects of this model. Firstly, we discuss localisation phenomena in one dimensional lattices as often experimentally realised in cold atom systems. Secondly, we investigate how disorder and symmetry influence heat transport in spin chains. More specifically, in the first part we consider a system of non-interacting fermions in one dimension subject to a single-particle potential consisting of a strong optical lattice, a harmonic trap, and uncorrelated on-site disorder. We investigate a global inhomogeneous quantum quench and present numerical and analytical results for static and dynamical properties. We show that the approach to the non-thermal equilibrium state is extremely slow and that it implies a sensitivity to disorder parametrically stronger than that expected from Anderson localisation. We also consider the above system in a strong non-uniform electric field. In the non-interacting case, due to Wannier-Stark localisation, the single-particle wave functions are exponentially localised without quenched disorder. We show that this system remains localised in the presence of nearest-neighbour interactions and exhibits physics analogous to models of conventional many-body localisation. The second part explores the hydrodynamics of the disordered XYZ spin chain. Using time-evolving block decimation on open chains of up to 400 spins attached to thermal baths, we probe the energy transport of this system. Our principal findings are as follows. For weak disorder there is a stable diffusive region that persists up to a critical disorder strength that depends on the XY anisotropy. Then, for disorder strengths above this critical value energy transport becomes increasingly subdiffusive.