This thesis focuses on the development and application of tensor networks methods for simulating one-dimensional quantum many-body systems. The general aim is to unveil numerical evidence for a new universality class of 1D interacting fermion systems referred to as quasi-Fermi liquids (qFLs). Realizing such systems beyond a purely analytical context and elucidating their peculiar characteristics is the principle value of this thesis.Chapter 1 begins by introducing the essential physics of one-dimensional systems, known as Luttinger liquids relative to their higher dimensional counterparts. Then it introduces non-linear Luttinger liquids juxtaposed against the defining quasi-Fermi liquid (qFL) field theory that motivates the following lattice studies. Chapter 2 provides an in-depth review of matrix product states (MPS), the main methodology used to explore quasi-Fermi liquids. The chapter serves as a pedagogical introduction to the field and details those novel algorithmic developments necessary for the study of fermionic systems. The central result of this chapter is the section on finite entanglement scaling techniques applied to dynamic correlation functions. Chapter 3 begins the study of lattice models that exhibit quasi-Fermi liquid characteristics. Here we explore the ground states exclusively by computing their static correlation functions. The central results include Luttinger parameter phase space diagrams which demonstrate the existence of the quasi-Fermi liquid across multiple lattice models. Chapter 4 studies the excitations of the qFL lattice models. We begin with results from non-MPS based methods which were essential for benchmarking our novel application of MPS to fermionic spectral functions. The central results of the chapter are the power-law exponents that demonstrate the low-energy universality of the qFL and our analysis of the dynamic structure factor. Chapter 5 concludes with a recapitulation of the main results and an outlook on future work. The Appendix contains those tensor diagrams necessary for a numerical implementation of those algorithms introduced in Chapter 2.--Author's abstract
This thesis focuses on the development and application of tensor networks methods for simulating one-dimensional quantum many-body systems. The general aim is to unveil numerical evidence for a new universality class of 1D interacting fermion systems referred to as quasi-Fermi liquids (qFLs). Realizing such systems beyond a purely analytical context and elucidating their peculiar characteristics is the principle value of this thesis.Chapter 1 begins by introducing the essential physics of one-dimensional systems, known as Luttinger liquids relative to their higher dimensional counterparts. Then it introduces non-linear Luttinger liquids juxtaposed against the defining quasi-Fermi liquid (qFL) field theory that motivates the following lattice studies. Chapter 2 provides an in-depth review of matrix product states (MPS), the main methodology used to explore quasi-Fermi liquids. The chapter serves as a pedagogical introduction to the field and details those novel algorithmic developments necessary for the study of fermionic systems. The central result of this chapter is the section on finite entanglement scaling techniques applied to dynamic correlation functions. Chapter 3 begins the study of lattice models that exhibit quasi-Fermi liquid characteristics. Here we explore the ground states exclusively by computing their static correlation functions. The central results include Luttinger parameter phase space diagrams which demonstrate the existence of the quasi-Fermi liquid across multiple lattice models. Chapter 4 studies the excitations of the qFL lattice models. We begin with results from non-MPS based methods which were essential for benchmarking our novel application of MPS to fermionic spectral functions. The central results of the chapter are the power-law exponents that demonstrate the low-energy universality of the qFL and our analysis of the dynamic structure factor. Chapter 5 concludes with a recapitulation of the main results and an outlook on future work. The Appendix contains those tensor diagrams necessary for a numerical implementation of those algorithms introduced in Chapter 2.--Author's abstract
We present a class of one-dimensional, generic, spinless fermion lattice Hamiltonians that express quasi-Fermi liquid physics, manifesting both Luttinger and Fermi liquid features solely due to irrelevant interactions. Using infinite matrix product state techniques, we unveil its universal structure by calculating static and dynamic responses. Key features include a finite discontinuity in the momentum distribution at the Fermi level, despite power-law singularities in the spectral function protected by particle-hole symmetry. Away from half filling, Landau quasiparticles emerge. Charge dynamics show either high-energy bound states or a concentration of spectral weight within the continuum for attractive or repulsive interactions, respectively. These universal features are realized across multiple models and energy scales, thus reifying the quasi-Fermi liquid as a unique paradigm for one-dimensional fermions.
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.
The study of low-dimensional, strongly-interacting quantum systems has been of interest in the past few decades due to new and exotic physics that show great potential in the advancement of fundamental physics as well as technological applications. In order to investigate such systems, a vast number of techniques have been developed throughout the years, all of which have their respective strengths and weaknesses. In the world of one-dimensional (1D) quantum systems, the density-matrix renormalization group (DMRG) algorithm and matrix product states (MPS) have proven to be powerful and reliable tools that provide invaluable information while providing a wide range of data that support both theoretical and experimental studies. Coupled with the development and refinement of tensor networks - an ansatz for quantum many-body wavefunctions, in particular MPS, the DMRG algorithm has seen a spike of improvements that has pushed the boundaries of its application into more sophisticated and complex systems.This thesis broadly covers the study of several 1D many-body quantum systems represented by infinite MPS (iMPS) - a class of ansatz that represent translationally-invariant, 1D many-body quantum states. More specifically, it is divided into two main parts.The first part investigates the 1D topological Kondo insulator (TKI). This is an effective model designed to investigate how strong interactions between conduction electrons and localized moments result in a symmetry-protected topological (SPT) phase, specifically the Haldane phase. By studying the effect of the electron repulsion as well as the competition between the local and non-local couplings between the electrons and the local moments, it is found that a topological phase transition from the SPT phase to a topologically trivial phase occurs with the assistance of electron interaction. The properties of both topologically trivial and nontrivial groundstates were revealed via means of quantum information theoretic quantities such as the entanglement entropy, entanglement spectrum and non-local order parameters. These quantities unveil the nature, structure and the effect that the different interactions have on the groundstate properties. Besides the phases, the properties of the critical point was studied and classified. By obtaining the critical exponents and central charge, the universality class of this transition was found to be that of 1D free fermions/bosons. Furthermore, the presence of electron interactions was found not to alter this universality class.The second part of this thesis is a project that stems from the necessity to overcome the inherent limitation of an iMPS in representing a system at a cusp of a quantum phase transition (QPT), specifically, in determining the critical point and exponents. In the thermodynamic limit, the laws of statistical mechanics dictates that the system’s correlation length diverges at the critical point. However, when representing a critical state with an iMPS, the finite bond dimension implies that the correlation length of the iMPS remains finite, albeit large. This presents a fundamental challenge in extracting critical data since the critical exponents are defined for systems in the thermodynamic limit. To overcome this shortcoming, finite-entanglement scaling (FES) - the analog of finite-size scaling, has previously been used to locate the critical point and extract the critical exponents for systems represented by finite-bond dimensional iMPS. While FES serves its purpose sufficiently well, the computational cost is high for just about any moderately complex, real-world system. To overcome such difficulties, the second part of this thesis combines new ideas of using higher-order cumulants of the order parameter with previously-known iMPS techniques to determine the critical point and extract the exponents, thus relating statistical mechanics of QPTs and tensor networks of critical states. Some of the main methods used and developed here are FES, the application of the Binder cumulant in iMPS and the new cumulant exponent relation. Putting these tools together provides a scheme that can be used to efficiently locate the critical point and extract the critical exponents at a lower computational cost than previously-known methods. This is shown through comparisons of several exemplary models such as the 1D and quasi-1D transverse field Ising model, the 1D TKI, the S = 1 Heisenberg chain with single-ion anisotropy and the 1D Bose-Hubbard model.
The properties of fermions in one dimension are markedly different from those in higher dimensions. In two or three dimensions, the many-body behaviour of fermions is described by the Fermi liquid model, whose most prominent result is that, instead of studying the behaviour of individual particles, it is better to describe the emergent phenomena of the whole system, known as quasiparticles. These are collective modes of excitations born out of particle-particle interactions, which resemble the constituent single particles, albeit with renormalised mass. In one dimension, however, the formation of such a correlated state is unfavourable and the quasiparticles are completely absent. Instead, 1D fermions are described by a model known as the Tomonaga-Luttinger liquid. Central to this model is the process of bosonisation, a mathematical technique which maps the excitation of fermions to bosonic modes of spin-density and charge-density waves. As a prerequisite of the bosonisation technique, the fermionic single-particle dispersion is approximated as a linear function. Consequently, the Luttinger liquid model is limited in its scope of applicability. Specifically, it is valid only in the vicinity of the Fermi points, where the excitation energy is low. Recent advancement in the theoretical field has expanded our understanding of 1D interacting fermions beyond the limit of linearity. The mobile-impurity model, proposed by L. Glazman and his co-authors, is a generalised model of 1D fermions which extends the scope of applicability of the Tomanaga-Luttinger model, accounting for the curvature of single-particle dispersion relations. It extends the region of validity of the tradition TLL theory from the neighbourhood of the Fermi points to a band which follows the single-particle dispersion. This thesis presents the results of an experimental work which provides direct evidence of non-linear behaviour predicted by the new model. It revolves around the measurement of the spectral function of electrons in surface-gate-defined quantum wires on a double-quantum-well, GaAs/AlGaAs device, and the comparison and fitting of the result to the non-linear theory. The spectral function is a theoretical concept frequently used in discussions on many-body problems. Its physical interpretation makes it closely related to experimental observables, such as the tunnelling current between the double quantum wells in this work. This thesis consists of three parts. Part one concerns the theories of the various many-body models and the experimental technique: Chapter 1 briefly introduces the concept of Fermi liquid, before focusing on a detailed account of the Luttinger-liquid model. The mobile-impurity model will be introduced in the framework of the Luttinger liquid model afterwards, followed by a review of experimental works reported in the literature on the study of TLL physics. In Chapter 2, the principles of the experimental tunnelling spectroscopy technique are described. The second part of the thesis details the structure of the test device used in the experiment and its fabrication. In particular, we will discuss the intricacies of the three-dimensional EBL fabrication process developed specifically for this project. The setup and running of the measurement circuit will also be discussed in this part of the thesis, in the context of a prerequisite sample characterisation measurement. In the final part of the work, data from the main experiment will be presented and analysed, in search for the electron-electron interaction effects which are predicted by the theoretical discussions in the first chapter.
For many years, the Luttinger liquid theory has served as a useful paradigm for the description of one-dimensional (1D) quantum fluids in the limit of low energies. This theory is based on a linearization of the dispersion relation of the particles constituting the fluid. Recent progress in understanding 1D quantum fluids beyond the low-energy limit is reviewed, where the nonlinearity of the dispersion relation becomes essential. The novel methods which have been developed to tackle such systems combine phenomenology built on the ideas of the Fermi-edge singularity and the Fermi-liquid theory, perturbation theory in the interaction strength, and new ways of treating finite-size properties of integrable models. These methods can be applied to a wide variety of 1D fluids, from 1D spin liquids to electrons in quantum wires to cold atoms confined by 1D traps. Existing results for various dynamic correlation functions are reviewed, in particular, the dynamic structure factor and the spectral function. Moreover, it is shown how a dispersion nonlinearity leads to finite particle lifetimes and its impact on the transport properties of 1D systems at finite temperatures is discussed. The conventional Luttinger liquid theory is a special limit of the new theory, and the relation between the two is explained.
The research work in this thesis is based on strongly interacting quantum lattice systems. The biggest part of research was conducted using state-of-the-art tensor network simulations. Tensor networks provide efficient and highly accurate representations of quantum states when any simply connected patch of the quantum state is slightly entangled. Matrix Product States (MPS) is a tensor network representation for quantum states which is quasi-exact for one-dimensional systems when entanglement entropy of any bipartition of the state follows ``area-law". Projected Entangled Pair States (PEPS) is the extension of MPS for higher dimensional systems, where entanglement entropy area-law is non-trivial. These formalisms and their relevant ground state optimization techniques, Density Matrix Renormalization Group (DMRG) for MPS and Simple Update (SU) with Tensor Renormalization Group (TRG) for PEPS are thoroughly analysed. Entanglement entropy area-law is usually obeyed by the ground state of local Hamiltonians, while generic highly excited states follow a volume-law (they span a finite part of the Hilbert space). Recently, a class of interacting system was shown to undergo a dynamical phase transition (Many Body Localization) where the entropy of every eigenstate follows an area-law. This transition is achieved when the system is highly disordered and the quantum many body state becomes localized similarly to the free particle Anderson localization. The area-law property makes highly excited eigenstates of Many Body Localized (MBL) systems efficiently represented by the MPS ansatz. We develop a highly optimized algorithm (eDMRG) which goes beyond ground state optimization and successfully target eigenstates in any part of the spectrum. This algorithm is used together with analytical calculations, based on local integrability of MBL systems, to identify the universal behaviour of the entanglement spectrum of highly excited eigenstates in MBL systems. In the second part we study interacting bosons in 2D optical lattices. Strongly interacting bosons are simulated using the Bose-Hubbard (BH) model when the interactions are strictly local and the Extended Bose-Hubbard (EBH) model, when additional dipolar interactions are present. The ground states of BH/EBH Hamiltonians in a hexagonal lattice are studied. The phase spectrum includes various insulating and critical phases which are studied in detail using infinite-PEPS ansatz. Additional results on entanglement entropy scaling with respect to the filling of the lattice are presented. Finally, a strongly interacting Harper-Hofstadter Hamiltonian is realized by combining synthetic fields with the BH model. The homogeneity of the system is then broken by a parabolic trapping potential, similar to the ones used in cold atom experiments. Using time-dependent Gutzwiller ansatz (GA) the expansion dynamics of the cloud in large square lattices are studied. In contrast to the expansion dynamics of the BH model, it is found that the synthetic fields generate a self-trapping effect. Using phenomenology and simulations the dynamics are studied in the hard-core and soft-core limit.