I first give an overview of the thesis and Matrix Product States (MPS) representation of quantum spin chains with an improvement on the conventional notation. The rest of this thesis is divided into two parts. The first part is devoted to eigenvalues of quantum many-body systems (QMBS). I introduce Isotropic Entanglement, which draws from various tools in random matrix theory and free probability theory (FPT) to accurately approximate the eigenvalue distribution of QMBS on a line with generic interactions. Next, I discuss the energy distribution of one particle hopping random Schrödinger operator in 1D from FPT in context of the Anderson model. The second part is devoted to ground states and gap of QMBS. I first give the necessary background on frustration free (FF) Hamiltonians, real and imaginary time evolution within MPS representation and a numerical implementation. I then prove the degeneracy and FF condition for quantum spin chains with generic local interactions, including corrections to our earlier assertions. I then summarize my efforts in proving lower bounds for the entanglement of the ground states, which includes some new results, with the hope that they inspire future work resulting in solving the conjecture given therein. Next I discuss two interesting measure zero examples where FF Hamiltonians are carefully constructed to give unique ground states with high entanglement. One of the examples (i.e., $d=4$) has not appeared elsewhere. In particular, we calculate the Schmidt numbers exactly, entanglement entropies and introduce a novel technique for calculating the gap which may be of independent interest. The last chapter elaborates on one of the measure zero examples (i.e., $d=3$) which is the first example of a FF translation-invariant spin-1 chain that has a unique highly entangled ground state and exhibits signatures of a critical behavior.
I first give an overview of the thesis and Matrix Product States (MPS) representation of quantum spin chains with an improvement on the conventional notation. The rest of this thesis is divided into two parts. The first part is devoted to eigenvalues of quantum many-body systems (QMBS). I introduce Isotropic Entanglement, which draws from various tools in random matrix theory and free probability theory (FPT) to accurately approximate the eigenvalue distribution of QMBS on a line with generic interactions. Next, I discuss the energy distribution of one particle hopping random Schrödinger operator in 1D from FPT in context of the Anderson model. The second part is devoted to ground states and gap of QMBS. I first give the necessary background on frustration free (FF) Hamiltonians, real and imaginary time evolution within MPS representation and a numerical implementation. I then prove the degeneracy and FF condition for quantum spin chains with generic local interactions, including corrections to our earlier assertions. I then summarize my efforts in proving lower bounds for the entanglement of the ground states, which includes some new results, with the hope that they inspire future work resulting in solving the conjecture given therein. Next I discuss two interesting measure zero examples where FF Hamiltonians are carefully constructed to give unique ground states with high entanglement. One of the examples (i.e., d=4) has not appeared elsewhere. In particular, we calculate the Schmidt numbers exactly, entanglement entropies and introduce a novel technique for calculating the gap which may be of independent interest. The last chapter elaborates on one of the measure zero examples (i.e., d=3) which is the first example of a FF translation-invariant spin-1 chain that has a unique highly entangled ground state and exhibits signatures of a critical behavior.
Classical simulation of quantum many-body systems is in general a challenging problem for the simple reason that the dimension of the Hilbert space grows exponentially with the system size. In particular, merely encoding a generic quantum many-body state requires an exponential number of bits. However, condensed matter physicists are mostly interested in local Hamiltonians and especially their ground states, which are highly non-generic. Thus, we might hope that at least some physical systems allow efficient classical simulation.Starting with one-dimensional (1D) quantum systems (i.e., the simplest nontrivial case), the first basic question is: Which classes of states have efficient classical representations? It turns out that this question is quantitatively related to the amount of entanglement in the state, for states with ``little entanglement'' are well approximated by matrix product states (a data structure that can be manipulated efficiently on a classical computer). At a technical level, the mathematical notion for ``little entanglement'' is area law, which has been proved for unique ground states in 1D gapped systems. We establish an area law for constant-fold degenerate ground states in 1D gapped systems and thus explain the effectiveness of matrix-product-state methods in (e.g.) symmetry breaking phases. This result might not be intuitively trivial as degenerate ground states in gapped systems can be long-range correlated.Suppose an efficient classical representation exists. How can one find it efficiently? The density matrix renormalization group is the leading numerical method for computing ground states in 1D quantum systems. However, it is a heuristic algorithm and the possibility that it may fail in some cases cannot be completely ruled out. Recently, a provably efficient variant of the density matrix renormalization group has been developed for frustration-free 1D gapped systems. We generalize this algorithm to all (i.e., possibly frustrated) 1D gapped systems. Note that the ground-state energy of 1D gapless Hamiltonians is computationally intractable even in the presence of translational invariance.It is tempting to extend methods and tools in 1D to two and higher dimensions (2+D), e.g., matrix product states are generalized to tensor network states. Since an area law for entanglement (if formulated properly) implies efficient matrix product state representations in 1D, an interesting question is whether a similar implication holds in 2+D. Roughly speaking, we show that an area law for entanglement (in any reasonable formulation) does not always imply efficient tensor network representations of the ground states of 2+D local Hamiltonians even in the presence of translational invariance. It should be emphasized that this result does not contradict with the common sense that in practice quantum states with more entanglement usually require more space to be stored classically; rather, it demonstrates that the relationship between entanglement and efficient classical representations is still far from being well understood.Excited eigenstates participate in the dynamics of quantum systems and are particularly relevant to the phenomenon of many-body localization (absence of transport at finite temperature in strongly correlated systems). We study the entanglement of excited eigenstates in random spin chains and expect that its singularities coincide with dynamical quantum phase transitions. This expectation is confirmed in the disordered quantum Ising chain using both analytical and numerical methods.Finally, we study the problem of generating ground states (possibly with topological order) in 1D gapped systems using quantum circuits. This is an interesting problem both in theory and in practice. It not only characterizes the essential difference between the entanglement patterns that give rise to trivial and nontrivial topological order, but also quantifies the difficulty of preparing quantum states with a quantum computer (in experiments).
We prove that generic quantum local Hamiltonians are gapless. In fact, we prove that there is a continuous density of states above the ground state. The Hamiltonian can be on a lattice in any spatial dimension or on a graph with a bounded maximum vertex degree. The type of interactions allowed for include translational invariance in a disorder (i.e., probabilistic) sense with some assumptions on the local distributions. Examples include many-body localization and random spin models. We calculate the scaling of the gap with the system's size when the local terms are distributed according to a Gaussian β orthogonal random matrix ensemble. As a corollary, there exist finite size partitions with respect to which the ground state is arbitrarily close to a product state. When the local eigenvalue distribution is discrete, in addition to the lack of an energy gap in the limit, we prove that the ground state has finite size degeneracies. The proofs are simple and constructive. This work excludes the important class of truly translationally invariant Hamiltonians where the local terms are all equal.
This thesis is focused on many-body localization (MBL) and the development of algorithms using the tensor networks representation for many-body localized systems. Chapter 1 is a detailed discussion on the MBL phenomenon. Chapter 2 gives an introduction to the numerical techniques that are frequently used in studying many-body quantum systems. Chapter 3 studies the shift-and-invert matrix product state method for obtaining highly excited states of MBL systems. Chapter 4 studies an efficient full diagonalization method using Wegner unitary gates. Chapter 5 studies the distribution of the entanglement entropy in the MBL-ergodic transition region. Chapter 6 studies the the influence of pseudo-spin SU (2) symmetry on the MBL phase. Chapter 7 is a side project that studies the influence of the duality-twisted boundary conditions on the critical twist defect chains. The MBL-ergodic transition is a dynamical phase transition which happens at finite energy density for a disordered and isolated many-body interacting system. It is a natural extension of Anderson localization where the localized phase survives the inter-particle interaction. In Chapter 1, we describe its relation to quantum thermalization, typical lattice models that harbor a MBL phase, MBL system’s phenomenological description, the characteristics of the MBL phase, experimental studies of MBL phenomenon in 1D and higher dimensions, the latest developments and the open questions on the nature of the MBL-ergodic transition. In Chapter 2, we provide a short but self-contained discussion of basic numerical methods that are frequently used in studying MBL physics. We mainly focus on the exact diagonalization methods and the matrix product state (MPS) representation, and explain their advantages and disadvantages. In Chapter 3, because obtaining the interior eigenstates of MBL Hamiltonians is a key step to studying the MBL phases and the transitions, we develop algorithms that can capture individual excited states to high fidelity. Using the fact that eigenstates of many-body localized systems have area-law entanglement and can therefore be efficiently represented as a matrix product state (MPS), we designed two algorithms that can generate excited states in the MPS representation, and use them to test the basic properties of MBL in the regime of large one-dimensional random field Heisenberg chains that were previously inaccessible due to the limitations of exact diagonalization (ED). In Chapter 4, we investigate an efficient way to diagonlize a fully many-body localized (FMBL) Hamiltonian. Using the Wegner flow technique for subsystems, we constructively build unitary gates that can greedily transform a local part of the FMBL Hamiltonian with the aim of lowering its average energy variance. Applying the Wegner unitary successively can make the system “flow” towards the diagonalized form. We find that the performance of this method is mainly controlled by the length of each individual Wegner unitary gate, and gets better with increasing length of the gate. We compare the performance of this constructive method with previous methods based on quasi-Newtonian optimization methods, and find that they have similar performance, while our constructive method is far more efficient. We also comment the possibility of constructing a tensor network state that resembles a multi-scale entanglement renormalization ansatz (MERA) using this constructive method. In Chapter 5, using strong subadditivity (SSA) theorem, we develop two order parameters – cut-averaged entanglement entropy (CAEE) and its slope (SCAEE), and show that they can be used to directly identify the volume or area law scaling in single eigenstates. We study the distribution of the SCAEE over disorder realizations. This distribution appears Gaussian at weak disorder, while on the MBL side, the distribution of the SCAEE is peaked at zero slope and has an exponential tail. In the critical regime we find that the distribution of the SCAEE is bimodal both over multiple disorder realizations as well as for single disorder realizations. The variance of the SCAEE distribution in the transition region seems to grow with system size when considered over disorder realizations. Our system sizes are too small to pin down its maximal value, but they are consistent with (among other possibilities) the maximal variance possible which would lead to half the states having zero SCAEE and half having maximal SCAEE. This scenario would lead to an entanglement entropy at the transition which scales as a volume law with half its thermal value. In Chapter 6, we consider the strongly disordered one-dimensional Hubbard model with spin disorder, which, as we argue, cannot be categorized as a conventional MBL phase. Though the spin rotation SU (2) symmetry is broken by the spin disorder, this model still preserves the pseudo-spin SU (2) symmetry under periodic boundary conditions. On the theoretical side, using the pseudo-spin algebra we show that a significant number of the excited states at any disorder strength have logarithmic correction to their von Neumann entanglement entropy, which violates the area-law entanglement for a typical MBL phase. A common feature of this group of excited states is that they usually have far more double occupancies than single occupancies. On the numerical side, we studied the time evolution after quantum quench from two extreme cases of product states – all single occupancies at quarter filling, and all double occupancies at half filling. We find in the former case entanglement entropy behaves as in a typical MBL system, while the latter is clearly delocalized, which suggests the existence of incomplete set of local integrals of motion. With the above evidence, it is convincing that this model provides a playground for studying a non-ergodic, non-MBL phase, and its relation to continuous SU (2) symmetry. Chapter 7 deviates from the common scene of MBL physics. In this chapter, we consider a series of one-dimensional critical twist defect models with various boundary conditions. These twist defects are embedded in the background of a generalization of Kitaev Z 2 toric code model, and can also be understood as a discrete Z k gauge theory in its deconfined phase. We analytically construct the duality transformation of these critical models, and use numerical methods to extract information on the central charges, sound velocities, and conformal dimensions. With the numerical data, we show the relation between the CFT contents of the critical chains with periodic boundary conditions and the critical chains with duality-twisted boundary conditions.
This thesis treats the classical simulation of strongly-interacting many-body quantum-mechanical systems in more than one dimension using matrix product states and the more general tensor product states. Contrary to classical systems, quantum many-body systems possess an exponentially larger number of degrees of freedom, thereby significantly complicating their numerical treatment on a classical computer. For this thesis two different representations of quantum many-body states were employed. The first, the so-called matrix product states (MPS) form the basis for the extremely successful density matrix renormalization group (DMRG) algorithm. While originally conceived for one-dimensional systems, MPS are in principle capable of describing arbitrary quantum many-body states. Using concepts from quantum information theory it is possible to show that MPS provide a representation of one-dimensional quantum systems that scales polynomially in the number of particles, therefore allowing an efficient simulation of one-dimensional systems on a classical computer. One of the key results of this thesis is that MPS representations are indeed efficient enough to describe even large systems in two dimensions, thereby enabling the simulation of such systems using DMRG. As a demonstration of the power of the DMRG algorithm, it is applied to the Heisenberg antiferromagnet with spin S=1/2 on the kagome lattice. This model's ground state has long been under debate, with proposals ranging from static spin configurations to so-called quantum spin liquids, states where quantum fluctuations destroy conventional order and give rise to exotic quantum orders. Using a fully SU(2)-symmetric implementation allowed us to handle the exponential growth of entanglement and to perform a large-scale study of this system, finding the ground state for cylinders of up to 700 sites. Despite employing a one-dimensional algorithm for a two-dimensional system, we were able to compute the spin gap (i.e. the energy gap to the first spinful excitation) and study the ground state properties, such as the decay of correlation functions, the static spin structure factors, and the structure and distribution of the nearest-neighbor spin-spin correlations. Additionally, by applying a new tool from quantum information theory, the topological entanglement entropy, we could also with high confidence demonstrate the ground state of this model to be the elusive gapped Z2 quantum spin liquid with topological order. To complement this study, we also considered the extension of MPS to higher dimensions, known as tensor product states (TPS). We implemented an optimization algorithm exploiting symmetries for this class of states and applied it to the bilinear-biquadratic-bicubic Heisenberg model with spin S=3/2 on the z=3 Bethe lattice. By carefully analyzing the simulation data we were able to determine the presence of both conventional and symmetry-protected topological order in this model, thereby demonstrating the analytically predicted existence of the Haldane phase in higher dimensions within an extended region of the phase diagram. Key properties of this symmetry-protected topological order include a doubling of the levels in the entanglement spectrum and the presence of edge spins, both of which were confirmed in our simulations. This finding simultaneously validated the applicability of the novel TPS algorithms to the search for exotic order.
This thesis is focused on the application and development of numerical methods for studying quantum many- body spin systems. Part I sets the stage for research projects. In chapter 1, I go over major milestones in the field of frustrated magnetism, how it fits into a larger category of condensed matter physics, and recent developments that provide context for the research projects presented in part II. Chapters 2 3 provide the reader with necessary concepts in physics of phase transitions and numerical methods that are often employed when studying ground states of quantum many body systems. Part II focuses on the research projects carried out during my PhD research. In chapter 4 I present a machine-learning inspired method for solving quantum many body problems. This method addresses issues of representation of variational wave functions and their optimization. In chapter 5, I present a numerical study on the embedding of an exactly solvable point in a phase diagram of an extended Heisenberg model on a kagome lattice. Chapter 6 studies a quantum Heisenberg model on the stuffed honeycomb lattice. Phase transitions in quantum many body systems manifest themselves as a reorganization of low energy spectra of the system. They come in a variety of forms ranging from symmetry breaking transitions, where a local order parameter is formed and to topological phase transitions where the change of the ground state can be attributed to non-local changes in the entanglement patterns. In chapter 2, I review classic models studied in the context of frustrated magnetism and briefly go over the formulation of Landau symmetry breaking theory and common order parameters. The breakdown of Landau symmetry-breaking theory for topological phase transition is discussed using the example of Resonance Valence Bond (RVB) state, which serves as a qualitative exemplification of spin liquid phases. I discuss means of detecting phase transitions with a focus on methods that are applicable in a numerical setting. Chapter 3 presents elements of the numerical machinery that the rest of the thesis employs and builds upon. I discuss traditional trade-offs and provide rough metrics used for analysis of strengths and drawbacks of different approaches. Then I focus on Exact Diagonalization (ED) and Variational Monte Carlo (VMC) methods that are used in later chapters. In chapter 4, I describe the application of machine learning techniques to the variational approach to the quantum many body problem. I introduce Computational Graph States (CGS), a new class of variational wave functions based on the dataflow programming paradigm. The use of computational graphs makes composing new wave functions easy due to its modular structure. Automatic differentiation allows us to obtain efficient derivatives with respect to variational parameters automatically. We also develop a novel optimization scheme, Supervised Wave-function Optimization. It combines ideas of supervised learning and imaginary time evolution. Our method resolves a number of technical challenges and allows seamless optimization of complicated variational models. During the exploration process we constructed a new class of system size invariant model wave functions that can be trained on one system size and generalized to other system sizes with high accuracy. In chapter 5, we study the phase diagram of the kagome antiferromagnet in the J2, Jz plane, which corresponds to second nearest neighbor interactions and spin anisotropy. We report an exactly solvable point with macroscopic degeneracy at J2 = 0, Jz = −0.5 and show that many phases emanate from around this point. In addition, we find a previously unknown transition within the spin-liquid region suggesting a possibility of two distinct spin liquids. In chapter 6, we consider a Heisenberg model on the stuffed honeycomb lattice which interpolated between the triangular and honeycomb lattices. Using exact diagonalization we compute the phase diagram of this model. We find that for this model, the triangular lattice J1-J2 spin liquid region is extended to a larger domain that has lower symmetries than that of the triangular lattice. We perform a variational study of projected BCS wave functions which we directly compare to exact ground states on 36 site cluster as well as consider variational energies from larger system sizes. We argue that construction of symmetric states is essential to obtain good variational energies and give evidence that Dirac spin-liquid state is likely to characterize the spin liquid phase.
One of the most important predictions of quantum mechanics is the existence of quantum entanglement. Since the inception of quantum theory, the study of entanglement has become ubiquitous in several branches of physics, from high energy to condensed matter physics and quantum information. To this avail, different measures of entanglement in quantum systems have been introduced. The most successful ones are the von Neumann and the Rényi entanglement entropies, which have been heavily studied both in lattice models and in quantum field theories (QFTs). A different, more complete characterization of entanglement is provided by the entanglement (or modular) Hamiltonian, defined as the logarithm of the reduced density matrix. One of the key property of this operator is that, under generic conditions, in 1+1-dimensional conformal field theories (CFTs) presents a local and few-body structure. Beyond just being of theoretical interest, this property has been recently used to efficiently construct the ground state of spin chains in quantum simulators. A recently introduced generalisation of the entanglement Hamiltonian is the negativity Hamiltonian, a non-Hermitian operator which is more useful for mixed states and tripartite entanglement. In this thesis we study these quantifications in many-body quantum systems and in QFTs. In the first part of the thesis, we focus on fermionic systems at equilibrium, and we obtain analytical expressions for the entanglement and the negativity Hamiltonians in the massless Dirac fermion QFT, both in presence of boundaries and at finite temperature. We compare these predictions in the continuum QFT with numerical calculations on the lattice, showing that a proper comparison requires a non-trivial continuum limit of the lattice expression, which we detail. We also investigate a lattice realisation of a non-unitary CFT, the bc-ghost theory, showing that the entanglement Hamiltonian presents an additional term with respect to the unitary case. In the second part, we study the entanglement entropies in an integrable QFT, the massless renormalisation group flow from the tricritical to the critical Ising CFT, demonstrating that the entropy is monotonically decreasing under the renormalisation group flow, as expected from the entropic c-theorem. In the third part we consider out-of-equilibrium problems. We study the melting of a domain wall configuration, computing analytically the entanglement Hamiltonian and showing that the insertion of a defect dramatically modifies the entanglement properties, changing the growth of the entropies from logarithmic to linear in time. Finally, we study generic quantum quenches in free fermionic theories, using the celebrated quasiparticle picture to obtain a prediction for the evolution of the entanglement Hamiltonian during the quench.