This PhD thesis is concerned with the investigation of Extended Hubbard Models on one-dimensional lattices (chains); the approach used to such investigation is that of determining exact results. The extended Hubbard models are models of correlated electrons on lattices; they are the extensions of a model proposed by Hubbard in \\cite{HUB} to explain the metal-insulator transition in
d-transition compounds as the effect of electronic correlations, usually neglected in band theories. Later on, a need for extensions of such model was first felt within the context of high-temperature superconductivity; indeed, the different behaviour of doped cuprates with respect to conventional superconductors suggested that a microscopic theory of this phenomenon could be formulated in terms of electron-electron interaction, without explicitly invoking the phonon coupling, like in BCS theory. Although an exhaustive microscopic formulation for cuprate superconductivity has not been achieved yet, the Extended Hubbard Models have been in the meanwhile applied to other problems in condensed matter; they are actually considered the reference paradigm for treating all those systems in which electronic correlations cannot be neglected. \\\\As observed above, this thesis focuses in particular on 1-D chains: indeed in last years the progress in technological synthesization and small-scale observations has shown that there are many classes of materials exhibiting a strong one-dimensional character, such as organic (super)conductors, linear-chain polymers, quantum wires, quantum-dot chains, one-dimensional cuprates. Many experimental results confirm that such materials exhibit a quite rich behaviour, with different phases, charge-spin separation, as well as a strong electronic correlations. However, a sufficient understanding of these phenomena on a theoretical point of view has not been reached yet. This is due to the fact the usual theoretical methods can be shown to fail in D=1, and theoreticians are thus forced to develop new schemes for approaching the problem of many interacting fermions in such dimension. The latter therefore represents a fascinating challenge. Some progress has been made and many aspects are nowadays more clear; nevertheless, due to the lack of traditional techniques, the need for a comparison with exact results is perhaps more urgent in D=1 than in higher dimensions. The aim of this thesis is to provide some results to the research in this field. The present work is organized as follows: the two first chapters are introductory, while chapters 3 and 4 mostly contain the author's original contribution. More explicitly, {\\bf chapter 1} is an introduction to the Hubbard models (both the ordinary and the Extended) and is meant for the reader who is non-familiar with these models. We tried to emphasize the motivation that historically led to their introduction and the physical contexts in which they are applied, not necessarily in one-dimension. \\\\{\\bf Chapter 2} is instead devoted to general problematics of one dimensional physics. We have tried to stress that there is a number of physically interesting materials that can be regarded as actual 1-D systems; we give a brief account of the difficulties leading the usual theoretical approaches, such as mean-field and Fermi Liquid, to fail in D=1. While briefly describing the important results of the proposed alternative formulations (Luttinger Liquid, Renormalization Group), we also point out the points that are still unsatisfactory, and emphasize the importance that exactly integrable models have in this field. \\\\Matching the subjects of chapter 1 and 2, the subsequent ones contain the results of our investigation of exactly solvable Extended Hubbard model on 1-D chains. \\\\We wish to emphasize to the reader non familiar with integrable systems that the literature on this subject is often strongly concerned with the mathematical aspects of the involved methods, leaving the determination of the physical features of the derived integrable models as a secondary issue. On the contrary, we believe that both aspects of the research deserve a relevant role in the study. For these reasons, while chapter 3 contains new results about the general formulation of integrability for fermionic systems, chapter 4 is devoted to the investigation of the physical features for some models whose integrability has been proved in the previous chapter. \\\\More explicitly, in {\\bf chapter 3} (apart from a preliminary introduction of sec.3.1) the theory of integrable {\\it fermionic} systems has been developed within the frame of the Quantum Inverse Scattering Method (QISM). The latter is a technique for obtaining integrable models, traditionally formulated with success for spin systems; its more recent application to electron systems was leaving many unsatisfactory aspects. To overcome this weaknesses, in sec.3.2 we both clarify the theory of matrix representation for fermionic operators, and describe a method for obtaining integrable models of fermions (straightforwardly readable in terms of fermion operators), starting from the the solutions of a
mathbbC-number equation (known as Yang-Baxter Equation). In sec.3.3 we first work out a general method of search for solutions of such equation; secondly we explicitly find 96 solutions, corresponding to 96 integrable Extended Hubbard models (most of which are new) for which a set of mutually commuting operators is proved to exist. In sec.3.4 we formulate a general theorem concerning the symmetries of the integrable models in QISM, and widely discuss the symmetries and supersymmetries for the class of extended Hubbard Models.\\\\ {\\bf Chapter 4} is devoted to the investigation of the physical properties of the integrable Extended Hubbard models found in chapter 3. We introduce a technique, called Sutherland's species technique, that, when applicable, yields an effective reduction of the Hilbert space and considerably simplifies the problem of diagonalizing the Hamiltonian. In fact, such technique is proved to be applicable to all the 96 models found above, and is used to determine the exact ground state phase diagram and the thermodynamics of the models. In particular we both recover in a quite straightforward way the phase diagram of a known model, and that of two new models (Model 2 and Model 3) that were not known before. \\\\The phase diagram exhibit structured forms with respect to the filling (density of electrons) and the on-site Coulomb repulsion (
U); the pressure and the compressibility in the ground state are also calculated. Model 2 exhibits a superconducting phase particularly enhanced towards positive values of
U, as well as a band-controlled insulator-superconductor transition at half filling; model 3 is the first example of exactly solved model with an asymmetric shape with respect to half-filling, due to the absence of particle-hole invariance of the Hamiltonian; a filling controlled metal-insulator transition is found to occur for this model. Generalizations to the other cases are also discussed. \\\\As to the thermodynamics, we compute the exact partition function of Model 3, which we believe to be particular worth of interest due to the above asymmetry. We then derive the exact specific heat of the model as a function of the temperature, which exhibits an interesting two-peak behaviour as well as a nearly universal crossing point. These results are discussed and compared with similar ones present in the literature. The chapters are of course linked each to another according to the above scheme; nevertheless, it has been a purpose of the author to make them as independently readable as possible. In particular, the reader who is more interested in the physical aspects of the models rather than in the general formalism of fermionic QISM can skip chapter 3 (at least in a first-sight reading) and pass directly to chapter 4; apart from some references to previous formulas the understanding should not be affected.\\\\{\\it \\noindent \\underline{Notational Remark:}} The customary symbol to denote the electron filling is
n; this has in fact been used throughout the two first introductory chapters. Nevertheless, in order to avoid possible misunderstanding with other extensive quantities, in chapter 3~and~4~the~symbol~has been changed into
rhoQuantum many body physics is an exciting research area, involving novel phases of matter with fundamentally new properties, but is also notoriously hard due to complexity of interacting quantum systems. Some popular approaches involve approximation techniques and numerical simulations, which are known to fail in several important cases. In contrast, rigorous mathematical tools, such as exact solutions and operator inequalities, have a narrower range of applicability, but provide guaranteed results and insights into the underlying physical mechanisms. This research aims to develop new techniques in this direction and use them to explore novel phases of matter. My first direction is to construct toy models where exact solutions are possible. Such models are important as they prove that certain physical phenomena are theoretically possible in nature, and often lead to the discovery of new phases. Such insights are provided by three exactly solvable models I have discovered: (1) a family of 1D quantum spin models hosting free parastatistical quasiparticles (an exotic type of identical particles beyond fermions, bosons, and anyons), proving for the first time that parastatistics is theoretically possible as an emergent phenomenon; (2) a 3D classical Ising model whose phases are characterized by topological features of certain loop observables, suggesting existence of previously unknown classical phases and phase transitions with topological order parameters; and (3) a family of models with exact p-wave superconducting ground states demonstrating the existence of Majorana quasiparticles and non-Abelian statistics in particle number-conserving systems. My second direction is to derive rigorous bounds and exact constraints on physical observables, which are applicable to large families of quantum many-body systems. I present three directions of progress: (1) a method that dramatically improves the upper bounds on the speed of information propagation in locally-interacting systems, which significantly extends the scope of these bounds and enables new applications; (2) bounds on finite-size errors in numerical simulations of lattice systems, including quench dynamics and gapped ground states; and (3) a locality bound on gapped ground states of power-law interacting systems, which leads to a generalization of the aforementioned error bounds to such systems. These error bounds have important theoretical implications such as proving the existence of the thermodynamic limit and stability of phases, and are practically useful in determining the validity of finite-size numerical simulations.
In this thesis I present contributions to the development, analysis and application of tensor network state methods for numerically exact quantum dynamics in one and two-dimensional lattice systems. The setting of numerically exact quantum dynamics is introduced in Chapter 2. This includes a discussion of exact diagonalization approaches and massively parallel implementations thereof as well as a brief introduction of tensor network states. In Chapter 3, I perform a detailed analysis of the performance of n-ary tree tensor network states for simulating the dynamics of two-dimensional lattices. This constitutes the first application of this class of tensor network to dynamics in two spatial dimensions, a long-standing challenge, and the method is found to perform on par with existing state-of-the-art approaches. Chapter 4 showcases the efficacy of a novel tensor network format I developed, tailored to electron-phonon coupled problems in their single-electron sector, through an application to the Holstein model. The applicability of the approach to a broad range of parameters of the model allows to reveal the strong influence of the spread of the electron distribution on the initial state of the phonons at the site where the electron is introduced, for which a simple physical picture is offered. I depart from method development in Chapter 5 and analyse the prospects of using tensor network states evolved using the time-dependent variational principle as an approximate approach to determine asymptotic transport properties with a finite, moderate computational effort. The method is shown to not yield the correct asymptotics in a clean, non-integrable system and can thus not be expected to work in generic systems, outside of finely tuned parameter regimes of certain models. Chapters 6 and 7 are concerned with studies of spin transport in long-range interacting systems using tensor network state methods. For the clean case, discussed in Chapter 6, we find that for sufficiently short-ranged interactions, the spreading of the bulk of the excitation is diffusive and thus dominated by the local part of the interaction, while the tail of the excitation decays with a powerlaw that is twice as large as the powerlaw of the interaction. Similarly, in the disordered case, analysed in Chapter 7, we find subdiffusive transport of spin and sub-linear growth of entanglement entropy. This behaviour is in agreement with the behaviour of systems with local interactions at intermediate disorder strength, but provides evidence against the phenomelogical Griffith picture of rare, strongly disordered insulating regions. We generalize the latter to long-ranged interactions and show that it predicts to diffusion, in contrast to the local case where it results in subdiffusive behaviour.
Introduction to integrable many-body systems I This is the first volume of a three-volume introductory course about integrable (exactly solvable) systems of interacting bodies. The aim of the course is to derive and analyze, on an elementary mathematical and physical level, the Bethe ansatz solutions, ground-state properties and the thermodynamics of integrable many-body systems in many domains of physics: Nonrelativistic one-dimensional continuum Fermi and Bose gases; One-dimensional quantum models of condensed matter physics like the Heisenberg, Hubbard and Kondo models; Relativistic models of the (1+1)-dimensional Quantum Field Theory like the Luttinger model, the sine-Gordon model and its fermionic analog the Thirring model; Two-dimensional classical models, especially the symmetric Coulomb gas. In the first part of this volume, we deal with nonrelativistic one-dimensional continuum Fermi and Bose quantum gases of spinless (identical) particles with specific types of pairwise interactions like the short-range δ-function and hard-core interactions, and the long-range 1/
This dissertation concerns the quantum many-body problem, which is the problem of predicting the properties of systems of several quantum particles from the first principles of quantum mechanics. Included under this umbrella are various problems of fundamental importance in quantum chemistry, condensed matter physics, and materials science. Of particular interest is the electronic structure problem, the problem of determining the state of the electrons in a system with fixed atomic nuclei. Since direct numerical solution of the many-body Schrödinger equation is intractable even for systems of moderate size, a diverse array of approximate methods has been developed. The broad goals of this dissertation are to improve the mathematical understanding of certain widely-used approximations, as well as to propose new methods. Roughly speaking, we consider three (overlapping) categories of methods: Green's function methods, embedding methods, and variational methods.One can understand Green's function methods in terms of many-body perturbation theory, which computes series expansions of physical quantities about a non-interacting reference system. These expansions can be expressed graphically in terms of Feynman diagrams, which can in turn be reorganized, in some cases, into an expansion in terms of so-called bold diagrams. Green's function methods can be specified by choosing a subset of bold diagrams to approximate the sum. At the same time, such methods can be understood in terms of an object known as the Luttinger-Ward (LW) functional, which admits a representation in terms of the bold diagrams. Many aspects of these constructions are purely formal, and indeed the existence of the fermionic LW functional as a single-valued functional has recently been called into question. To contribute to the understanding of these issues, we provide rigorous proofs of the combinatorial construction and analytic interpretation of the bold diagrams in the simplified setting of a classical field theory. In this setting we also provide a rigorous non-perturbative construction of the LW functional via convex duality and prove several key properties, including continuity up to the boundary of its domain and asymptotics in the limit of large interaction.Quantum embedding methods, meanwhile, view a large system as being composed of smaller fragments that are treated with high accuracy and embedded in the larger system in a mutually consistent way. Inspired by a connection between the boundary analysis of the LW functional and embedding, we perform similar analysis for the 1-RDM theory for fermionic systems, which is also developed via convex duality, illustrating a relation to fermionic embedding methods such as the density matrix embedding theory (DMET).Another embedding method of note is the dynamical mean-field theory (DMFT), which is at the same time a Green's function method that can be understood in terms of the LW functional. DMFT relies on the solution of impurity problems, which specify the embedding of an interacting system into a non-interacting bath. Underlying DMFT is a result about the sparsity pattern of the self-energy matrix for impurity problems, which to our knowledge has not been proved in the literature. We provide a rigorous proof of this result in various classical and quantum settings. We go on to investigate the fermionic DMFT in depth, identifying the key mathematical structures that appear in the algorithmic loop for solving it and using these to prove the well-posedness of this loop, in a certain sense.Finally, we introduce a suite of new approaches to the quantum many-body problem that provide variational lower bounds to the ground-state energy. These methods, which combine the themes of convexity and embedding, are based on novel convex relaxations of the variational principles for the ground-state energies of many-body systems. To begin, we recover a second-quantized version of the formalism of strictly correlated electrons (SCE), which yields an exact expression for the exchange-correlation functional in Kohn-Sham density functional theory in the limit of infinite Coulomb repulsion in terms of the solution of a multi-marginal optimal transport problem. We introduce a semidefinite relaxation method for approximately solving this problem and obtaining a lower bound for the ground-state energy. The ideas underlying this relaxation are generalized considerably, outside the context of SCE, to yield much tighter lower bounds, which we validate numerically for both quantum spin systems and fermionic systems. We also describe how these relaxation methods can be interpreted as embedding methods via convex duality.
THESIS SUMMARYTEXT:This thesis is made of two parts. In the first one, the issue of entanglement in many body systems is addressed. The concept of entanglement and some of the recent progress on the study of entropy of entanglement in many body quantum systems are reviewed. Emphasis is placed on the scaling properties of entropy for one-dimensional models at quantum phase transitions. Then, we focus on the area-law scaling of the entanglement entropy. An explicit computation in arbitrary dimensions of the entanglement entropy of the ground state of a discretized scalar free field theory that shows the expected area law result is also presented. For this system, it is shown that area law scaling is a manifestation of a deeper reordering of the vacuum produced by majorization relations.To finish this first part, the issue of how simple can a quantum system be such as to give a highly entangled ground state is addressed. In particular, we propose a Hamiltonian of a XX model with a ground state whose entropy scales linearly with the size of the block. It provides a simple example of a one dimensional system of spin-1/2 particles with nearest neighbour interactions that violates area-law for the entanglement entropy.The second part of this thesis deals with the problem of simulating quantum mechanics for highly entangled systems. Two different approaches to this issue are considered. One consists of using ultra-cold atoms systems as quantum simulators. With this aim, some experimental techniques related to cold atoms that allow to simulate strongly correlated many body quantum systems are reviewed an explicit example of simulation is presented. In particular, we analyze how to achieve a Mott state of Laughlin wave functions in an optical lattice and study the consequences of considering anharmonic corrections to each single site potential expansion that were not taken into account until now.Finally, a different approach to simulate strongly correlated systems is considered: to use small quantum computers to simulate them. An explicit quantum algorithm that creates the Laughlin state for an arbitrary number of particles n in the case of falling fraction equal to one is presented. We further prove the optimality of the circuit using permutation theory arguments and we compute exactly how entanglement develops along the action of each gate. We also discuss its experimental feasibility decomposing the qudits and the gates in terms of qubits and two qubit-gates as well as the generalization to arbitrary falling fraction.KEYWORDS: Entanglement, Many body quantum systems, Quantum Information Condensed Matter, Cold atoms, Spin chains, Quantum simulator, Quantum computation.