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
rhoIn this Thesis, we develop and apply an accurate method to determine the ground-state properties of strongly-correlated electrons in the framework of lattice model Hamiltonians. In lattice density-functional theory (LDFT) the basic variable is the one-particle density matrix γ. From a generalized HK theorem, the ground state energy Egs[γgs] = minγ E[γ] is obtained by minimizing the energy-functional E[γ] over all the physical or representable γ. The energy functional can be divided into two contributions: the kinetic-energy functional T [γ], which linear dependence on γ is exactly known, and the correlation-energy functional W [γ], which dependence on γ is not explicitly known. Finding accurate approximations for W [γ] constitutes the actual challenge of this thesis. Part of this work is based on previous studies which derive an approximation of W [γ] for the Hubbard Hamiltonian based on scaling hypothesis and an exact analytical results for the Hubbard dimer. However, this approach is limited to spin independent and homogeneous systems. In order to extend the scope of LDFT, we have developed three different approaches to derive W [γ] in order to study broken symmetry problems. First we extend the previous scaling functional to system having charge transfer. By a systematic study of the dependence of W [γ] on the charge distribution we find similar scaling properties as for the homogeneous case. An extension to the Hubbard model on bipartite lattices is then derived and applied to finite and infinite systems with repulsive and attractive interactions. The high accuracy of this functional is demonstrated. Nevertheless, this approach is difficult to transfer to more complex systems, since the calculation of W [γ] applies to the system as a whole. To overcome this problem, we derive another approximation based on local scaling properties. This functional is formulated on-site so that it can be applied to any ordered or disordered Hamiltonian with on-site interaction. As applications we study the metal-insulator transition in the ionic Hubbard model in oneand two-dimensions, as well as the one-dimensional Hubbard chain with first and second nearest neighbor hoppings. Finally, we develop a numerical approach of W [γ], on the basis on exact diagonalizations of an effective many-body Hamiltonian corresponding to a cluster surrounded by an effective medium. This effective Hamiltonian depends on the density matrix γ and allows to derive approximation to W [γ] which improve systematically with increasing cluster size. The formulation is spin-dependent and allows a straightforward generalization to multi-orbital correlated systems such spd-Hamiltonians. In addition it take into account the effects of short-range charge and spin fluctuations in the functional. The accuracy of the method is demonstrated for the Hubbard model by comparison with the Bethe-Ansatz solution (1D) and quantum Monte Carlo simulations (2D). Finally, an outlook or relevant future developments of the theory is provided.
Free-fermion solvability in quantum spin systems is increasingly understood to be governed by a graph Clifford algebra defined from the frustration graph of the Hamiltonian. When the frustration graph belongs to certain classes, such as the even-hole-free and claw-free (ECF) class, the Hamiltonian is solvable by hidden free fermions: it admits a free-fermion solution although it does not reduce to a Majorana bilinear under the Jordan-Wigner transformation. However, unlike in the Jordan-Wigner case, where each mode is a linear combination of single Majorana fermions, the explicit operator structure of the hidden free-fermion modes---and that of the local conserved charges---has remained obscure. In this work, we derive a path-product expansion that expresses each free-fermion mode as a linear combination of products along induced paths in the extended frustration graph. The expansion follows from the Krylov generating function and yields the modes directly, without using the transfer matrix or nonlocal conserved charges; the resulting mode decomposition also computes infinite-temperature dynamical correlation functions. We further obtain explicit expressions for local conserved charges as linear combinations of path products along induced paths; these charges apply beyond the free-fermion (ECF) class to more general claw-free frustration graphs. We also identify a family of generalized conserved charges containing both the known nonlocal conserved charges and these local charges. For the homogeneous periodic Fendley model, the local conserved charges exhibit the same Catalan-tree pattern as those of the spin-
1/2 XXX chain.
Integrable models are characterized by the existence of an extensive number of local conserved quantities.Their existence is closely related to various exact methods, such as the Bethe ansatz and quantum inverse scattering.They also play an important role in nonequilibrium dynamics in the integrable models.Their steady states and transport phenomena exhibit exotic properties, which are also related to this point.Local conserved quantities in a quantum integrable model in one dimension are generated from the transfer matrix of the corresponding classical two-dimensional model.However, it is still non-trivial to obtain their explicit expression.The meaning of "local" here is that each conserved quantity consists of the sum of observables acting on a finite range of sites called its support.We note that the local conserved quantities in this context are not O (1)-body observables in general.Their explicit expressions enable us to analyze theories of nonequilibrium dynamics using the local conserved quantities in detail, for example, the generalized Gibbs ensemble, the Mazur bound, and the generalized hydrodynamics.In this thesis, we have studied the local conserved quantities in the XYZ spin-1/2 chain, which is an integrable model corresponding to the classical eight-vertex model.We have obtained their explicit expression for arbitrary coupling constants of the XYZ chain.The model at the isotropic point called the Heisenberg or XXX chain is an exceptional case in which the expression had been known.We have checked that our results reproduce the known expression for the XXX chain as a special case.To obtain the explicit expression, we have employed the notation of doubling-product operators, which was introduced by Shiraishi to study the non-integrability of the XYZ spin-1/2 chain in a magnetic field.Using this notation, we have found recursion relations of the coefficients among different terms in the conserved quantities in the XYZ chain.Then, we have solved the recursion relations and obtained the explicit expression.We have shown that the local conserved quantity Q k for a fixed support k is essentially unique i ii ABSTRACT if 2 ≤ k ≤ L/2 with L being the number of sites.We have obtained the explicit expression for all {Q k } for 1 ≤ k ≤ L/2 and one Q k for each k with L/2 < k ≤ L.We have also studied the currents of these conserved quantities based on the respective continuity equations.We have obtained their explicit expression and numerically investigated their time evolution.In particular, we have studied the effects of an integrabilitybreaking perturbation on the time evolution of the currents and found that this term drastically changes their long-time behaviors as follows.As for the initial state, we have chosen several energy eigenstates of the XXX chain with a nonzero energy current, which is also known to be conserved.In particular, we have studied the time evolution with the Hamiltonian of the XXZ chain with and without a next-nearest-neighbor interaction and seen the effects of this integrability-breaking perturbation.We have found that, for these initial states, the long-time average of the total currents is nonzero for the integrable cases while zero for the non-integrable cases.We have also studied this behavior analytically and shown that the degeneracy of the energy eigenstates is essential.Our results for the XYZ spin-1/2 chain have revealed the structure of the conserved quantities and currents in detail.They provide materials for the development of algebraic studies in integrable models and will be useful for detailed analyses of nonequilibrium dynamics in these models.