Understanding the collective behavior of a quantum many-body system, a system composed of a large number of interacting microscopic degrees of freedom, is a key aspect in many areas of contemporary physics. However, as a direct consequence of the difficultly of the so-called many-body problem, many exotic quantum phenomena involving extended systems, such as high temperature superconductivity, remain not well understood on a theoretical level. Entanglement renormalization is a recently proposed numerical method for the simulation of many-body systems which draws together ideas from the renormalization group and from the field of quantum information. By taking due care of the quantum entanglement of a system, entanglement renormalization has the potential to go beyond the limitations of previous numerical methods and to provide new insight to quantum collective phenomena. This thesis comprises a significant portion of the research development of ER following its initial proposal. This includes exploratory studies with ER in simple systems of free particles, the development of the optimisation algorithms associated to ER, and the early applications of ER in the study of quantum critical phenomena and frustrated spin systems.
Understanding the collective behavior of a quantum many-body system, a system composed of a large number of interacting microscopic degrees of freedom, is a key aspect in many areas of contemporary physics. However, as a direct consequence of the difficultly of the so-called many-body problem, many exotic quantum phenomena involving extended systems, such as high temperature superconductivity, remain not well understood on a theoretical level. Entanglement renormalization is a recently proposed numerical method for the simulation of many-body systems which draws together ideas from the renormalization group and from the field of quantum information. By taking due care of the quantum entanglement of a system, entanglement renormalization has the potential to go beyond the limitations of previous numerical methods and to provide new insight to quantum collective phenomena. This thesis comprises a significant portion of the research development of ER following its initial proposal. This includes exploratory studies with ER in simple systems of free particles, the development of the optimisation algorithms associated to ER, and the early applications of ER in the study of quantum critical phenomena and frustrated spin systems.
In contrast to classical many-particle systems for which the number of degrees of freedom grows linearly with the system size, it grows exponentially for quantum many-particle systems. On the one hand, the resulting complexity of the quantum systems has a high potential for technological advances if it can be exploited directly. For example, in recent years, the fundamentals for the realization of quantum computers and their possible applications have been investigated intensely. On the other hand, the high complexity, especially in strongly correlated regimes, often makes the analytical or numerical analysis of underlying physical phenomena very difficult. This is why some systems of high technological relevance, like the high-temperature superconductors, are still unsatisfactorily understood. A main topic of this thesis is the development of efficient numerical methods for the simulation of strongly correlated quantum lattice models. For one-dimensional systems, the density-matrix renormalization-group (DMRG) is such a very successful method. The physical states of interest are approximated within a certain class of ansatz states. These ansatz states are designed in a way that the number of degrees of freedom are prevented from growing exponentially. With the DMRG method, ground states of interacting one-dimensional systems can often be determined with very high precision. This observation indicates already that ground states actually do not exhaust the available exponential number of degrees of freedom but can be described efficiently with much fewer parameters. The ansatz states used in DMRG are the so-called matrix product states. A very useful criterion for the construction of corresponding appropriate ansatz states for higher-dimensional systems is the quantum nonlocality as for example quantified by the entanglement entropy. The ansatz states should be able to reproduce the scaling properties of the quantum nonlocality in the physical states of interest with as few parameters as possible. The first part of the thesis, therefore, provides analytical and numerical analysis of the scaling of quantum nonlocality with the system size or time in different, physically relevant scenarios. For example, the scaling of Renyi entropies and their dependence on boundary conditions is derived within the 1+1-dimensional conformal field theory. Conjectures and analytical indications concerning the properties of entanglement entropy in critical fermionic and bosonic systems are confirmed numerically with high precision. For integrable models in the thermodynamic limit, general preconditions are derived under which subsystems converge to steady states. These steady states are non-thermal and retain information about the initial state. It is shown that the entanglement entropy in such steady states is extensive. For short times, the entanglement entropy grows typically linearly with time, causing an exponential increase in computation costs for the DMRG method. Consequently, time-dependent DMRG simulations usually have no direct access to the long time regime. The second part of the thesis focuses on the development and improvement of the above-mentioned numerical techniques. The time-dependent DMRG is complemented with an extrapolation technique for the evaluated observables. In this way, the problem of the entropy increase can be circumvented, allowing for a precise determination of spectral functions. The method is demonstrated using the example of the Heisenberg antiferromagnet and results are compared to Bethe-Ansatz data for T=0 and quantum Monte Carlo data for T>0. For the simulation of higher-dimensional systems, projected-entangled pair-states (PEPS) and the multiscale entanglement renormalization ansatz (MERA) turn out to be appropriate wavefunction classes. Initially, they could only be used for spin systems. Within this thesis, variants of those state classes are being suggested for the fermionic case. It is shown how the known algorithms for spin systems can be translated into corresponding algorithms for the fermionic systems. The occurring computation costs increase only by a marginal overhead. First numerical benchmarking tests proceeded successfully. Thus, new techniques are available which can for example be used to examine the two-dimensional Hubbard model. This model is a candidate for the description of high-temperature superconductivity.
Quantum many-body systems are comprised of complex networks of microscopic interactions that work together to produce novel collective phases and phenomena. \nFor strongly-interacting systems, the connection between the underlying atomic-scale behaviour and the resulting emergent observable phenomena is exponentially (and, in practice, infinitely) complicated such that it becomes impossible to fully understand the connection between the physics at the microscopic and macroscopic levels. \nThe most promising theoretical approaches to addressing this infinite complexity utilize microscopic coarse-graining along with powerful computational simulations and strategies. \nIn this thesis, we apply numerical methods to effective lattice models and field theories with the goal of shedding light onto universal critical behaviour and exotic low-temperature phases in condensed matter physics. \n \nWe start by exploring universal features of critical non-interacting systems, for which several analytical strategies are readily available for calculating certain observables. \nWe study the behaviour of the system's entanglement for various entangling geometries, and utilize numerical techniques in order to isolate new universal features of Gaussian fixed points, which provide insight into the underlying critical theory's renormalization group flow. \nWe then proceed to examine computational strategies and, in particular, Monte Carlo simulations for studying more general interacting models. \nWe focus on the application of such strategies to the field of high-temperature superconductivity, for which we develop a coarse-grained model, simulate macroscopic observables and compare our results with those of recent experiments. \nWe conclude that lattice field theories together with innovative computational methods offer new perspectives on both universality and emergent phenomena in quantum matter.
Emergent phenomena of interacting quantum many-body systems arguably pose some of the greatest challenges to modern theoretical physics. Many of these systems are characterized by drastic changes in their properties as they are driven into a regime where quantum effects and interactions become relevant. Examples are the fractional quantum Hall effect and high-temperature superconductors, to name just two. New theoretical tools are necessary to study such systems, as even simple model Hamiltonians are not fully understood after decades of intense research. Numerical simulations are emerging as an indispensable tool alongside more traditional analytic approaches and have seen immense progress in recent years. In this thesis, we study tensor network states, a new class of numerical methods that have been developed in recent years drawing on insights from quantum information theory. Combining renormalization group methods with knowledge about the entanglement structure of ground states of quantum systems, tensor network states aim to efficiently describe ground states of strongly correlated quantum systems. Center stage is taken by the projected entangled-pair states. This ansatz for two-dimensional systems is discussed in detail and its capabilities are assessed by comparison to established methods and by exploring a frustrated system. We then develop a general formalism to exploit Abelian symmetries in a tensor network algorithm and demonstrate its application to projected entangled-pair states. A large part of this thesis is devoted to the discussion of a class of supersymmetric models for interacting lattice fermions. After introducing the concept of supersymmetry in quantum mechanics and an extensive overview of conformal field theory and its relation to numerical simulations, we study the critical theory of the supersymmetric model on the square ladder. We use two tensor network state algorithms, the density-matrix renormalization group and multi-scale entanglement renormalization along with exact diagonalization to explore this challenging system. In the final two chapters, we discuss two applications of tensor network states, namely simulations of the SU(3) Heisenberg model in two dimensions using two different tensor network state algorithms with the goal of illuminating the nature of its ground state, and the indistinguishability, a measure to detect phase transitions and classify wave functions which is particularly suitable for calculations based on tensor network states.
Quantum many-body calculation methods, including exact diagonalization, quantum Monte Carlo, density matrix and tensor renormalization, and dynamic mean field theory, as well as emerging methods harnessing the power of artificial intelligence and quantum computing, can be used to accurately and efficiently calculate the physical properties of correlated quantum many-body systems. The quantum particles in the interacting, non-perturbative many-body system are highly entangled, and the mean field theory often lacks sufficient accuracy or may even break down for certain problems. It is thus very necessary to develop new methods for studying the exotic states and emergent phenomena in correlated systems, such as high-temperature superconductivity and frustrated quantum magnetism. Recently, the rapid progress in the field of many-body calculation has shown the characteristics of interdisciplinary cross, and there is increasingly significant interplay and integration with machine learning, materials science, quantum simulation and computation, among other fields. This paper outlines the brief history of many-body calculations, major approaches and their current status, as well as main challenges in this field.
At the interface of quantum information and condensed matter physics, the study of entanglement in quantum many-body systems requires a new toolset which combines concepts from each. This thesis introduces a set of computational methods to study phases and phase transitions in lattice models of quantum systems, using the Renyi entropies as a means of quantifying entanglement. The scaling of entanglement entropy can give valuable insight into the phase of a condensed matter system. It can be used to detect exotic types of phases, to pinpoint transitions between phases, and can give us universal information about a system. \nThe first approach in this thesis is a technique to measure entanglement in finite size \nlattice systems using zero-temperature quantum Monte Carlo simulations. The algorithm \nis developed, implemented, and used to explore anomalous entanglement scaling terms in \nthe spin-1/2 Heisenberg antiferromagnet. \nIn the second part of this thesis, a new and complementary numerical technique is introduced to study entanglement not just in finite size systems, but as we approach the thermodynamic limit. This “numerical linked-cluster expansion” is used to study two different systems at their quantum critical points — continuous phase transitions occurring at zero temperature, at which these systems exhibit universal properties. Remarkably, these universal properties can be reflected in the scaling of entanglement. \nEntanglement offers a new perspective on condensed matter systems, one which takes us closer to genuinely understanding what goes on in these materials at the quantum mechanical level. This thesis demonstrates the first steps in developing an extensive list of computational tools that can be used to study entanglement over a wide range of interacting quantum many-body systems. With the ever increasing computational power available, it may be only a matter of time before these tools are used to create a comprehensive framework for the characterization of condensed matter phases and phase transitions.
The collective behavior exhibited by a large number of microscopic quantum particles is at the heart of some of the most striking phenomena in condensed matter physics such as Bose-Einstein condensation and superconductivity. Physicists and mathematicians have made great progress in understanding when and how these collective phenomena emerge through the interplay of particle statistics, particle interaction and the value of thermodynamic parameters like the temperature or the chemical potential. Due to the extreme complexity of realistic many-body systems, it is natural to introduce appropriate simplifications to render their analysis feasible. Three examples of such simplifications which have proven themselves as viable starting points for a fruitful and mathematically rigorous analysis of many-body systems are the following: (a) the study of integrable models; (b) the derivation of effective theories, valid on a macroscopic scale, from more fundamental microscopic theories under appropriate coarse-graining; and (c) the use of quantum information theory to understand general connections between correlation, entanglement and particle statistics. In this thesis, we present mathematically rigorous results that were obtained in these three directions. (1) We prove anomalous quantum many-body transport in XY quantum spin chains for certain choices of the external magnetic field. The anomalous transport is described via new kinds of anomalous Lieb-Robinson bounds, including one of power-law type. We note that the XY spin chain is integrable as it can be mapped to free fermions via the non-local Jordan-Wigner transformation. (2) We derive effective macroscopic theories of Ginzburg-Landau type from the microscopic BCS theory of superconductivity in certain circumstances. We study the case of a multi-component order parameter for translation-invariant systems and the condensation of fermion pairs at zero temperature in a domain with a hard boundary. (3) We use techniques from quantum information-theory to derive bounds on the entropy of fermionic reduced density matrices, a measure of the entanglement inherent to a fermionic quantum state.