Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:0806.1761

Exact Ground State Energy of Hubbard Rings in the Atomic Limit

William Hodge, Natalie A. W. Holzwarth, Kerr, W. C.

Published Jun 10, 2008• 1 citation
arXivHTMLPDF

Abstract

Using a straightforward extension of the analysis of Lieb and Wu, we derive a simple analytic form for the ground state energy of a one-dimensional Hubbard ring in the atomic limit. This result is valid for an \textit{arbitrary} number of lattice sites $L$ and electrons $N \leq L$. Furthermore, our analysis, including an application of the theory of stochastic matrices, provides insight into the degeneracy and spin properties of the ground states in the atomic limit. We give numerical results which illustrate how the atomic limit is approached.

DOI
10.48550/arxiv.0806.1761

Related papers

Exact and variational investigations of Hubbard rings

William Hodge

We study the one-dimensional Hubbard model for interacting electrons in a crystal lattice, using two methods: solution of the Lieb-Wu equations in the atomic limit of no intersite hopping, and variational techniques to obtain a lower bound for the ground state energy. From the Lieb-Wu equations we derive an exact expression for the ground state energy in the atomic limit. We obtain our result by exploiting the spin symmetry of the Hubbard Hamiltonian. Our energy is lower than that obtained by previous published results, and we are able to explain why our results differ. Our method also gives information about the magnetic ordering of the ground state as a function of the number of electrons. Our results are valid for arbitrary electron density values. The variational techniques we use are applied to the two-particle reduced-density matrix (2-RDM), because this quantity suffices to obtain expectation values of all one- and two-particle operators, and because it is less computationally intensive than obtaining the complete N -particle wave function. Because the 2-RDM must be positive semidefinite, we use an efficient programming algorithm known as the Semidefinite Programming Algorithm. A naive search through the space of two-particle matrices encounters a difficult problem, known as the " N -representability problem", which is to assure that the trial density matrices actually come from a wave function that has the required exchange symmetry. This problem is ameliorated by imposing constraints on the variational search that eliminate trial matrices with the incorrect exchange symmetry. These constraints express known physical properties of the ground state (for example its total spin and its translational symmetry) and the so-called D, Q , and G two-positivity conditions that express the required exchange symmetry. To calculate matrix elements, we use the two-particle eigenstates of the reduced Hubbard Hamiltonian; this choice of basis significantly decreases the number of matrix elements needed for energy calculations. After an optimal two-particle density matrix has been obtained, we calculate additional properties of the ground state, including spin correlation functions, chemical potentials, and the momentum distribution. We compare all of these properties to values obtained by exact numerical solution of the Schrödinger equation.

•Jan 1, 2008•0 citations

GROUND-STATE PROPERTIES AND ELEMENTARY EXCITATIONS OF THE ONE-DIMENSIONAL HUBBARD MODEL

J. M. P. Carmelo, D. Baeriswyl

The ground state of the one-dimensional Hubbard Hamiltonian is discussed on the basis of the Bethe ansatz solution of Lieb and Wu. A simple analytical representation of the distribution functions for charge and spin degrees of freedom is derived. The elementary excitations are associated with electron-hole pairs of two pseudo-fermion systems, one for the charge, the other for the magnetic excitations. Spinons and holons correspond to ground state configurations of odd-numbered rings. It is argued that holons will form bound pairs in weakly coupled chains.

International Journal of Modern Physics B•Oct 1, 1988•7 citations

Phase transitions and exact ground-state properties of the one-dimensional Hubbard model in a magnetic field

Chia‐Wei Yang, ARMEN N. KOCHARIAN, Y. L. Chiang

The exact phase diagrams of the one-dimensional Hubbard model, both attractive and repulsive, in the presence of an arbitrary magnetic field h for various electron concentrations n and on-site interaction strengths U <0 or U >0 are calculated and investigated. The exact ground-state properties, namely, the ground-state energy, the average spin (magnetization), the concentration of the doubly occupied sites, the kinetic energy, the chemical potential, the spin (magnetic) susceptibility and the charge compressibility, are calculated and examined over a wide range of interaction strengths U for various h and n . It is found that the spin susceptibility at half-filling is non-analytic and changes discontinuously as U →0. The exact theory shows the absence of a charge energy gap in the U ⩽0 region for all n and provides, for the chemical potential, the rigorous upper and lower bounds for half-filled and empty bands respectively. The analytical results derived for the weak-coupling limit and asymptotic expansions for strong coupling are in full agreement with the numerical calculations.

Journal of Physics Condensed Matter•Aug 2, 2000•24 citations

Perturbation expansion of the ground‐state energy for the one‐dimensional cyclic Hubbard system in the Hückel limit

M. Takahashi, Paul Bracken, Jakub Čı́žek +1 more

Abstract The perturbation expansion coefficients for the ground‐state energy of the half‐filled one‐dimensional Hubbard model with N = 4ν + 2, (ν = 1, 2,…) sites and satisfying cyclic boundary conditions were calculated in the Hückel limit (U/β → 0), where β designates the one‐electron hopping or resonance integral, and U, the two‐electron on‐site Coulomb integral. This was achieved by solving—order by order—the Lieb–Wu equations, a system of transcendental equations that determines the set of quasi‐momenta {ki} and spin variables {τα} for this model. The exact values for these quantities were found for the N = 6 member ring up to the 20th order in terms of the coupling constant B = U/2β, as well as numerically for 10 ⩽ N ⩽ 50, and the N = 6 Lieb–Wu system was reduced to a system of sextic algebraic equations. These results are compared with those of the infinite system derived analytically by Misurkin and Ovchinnikov [Teor. Mat. Fiz. 11, 127 (1972)]. It is further shown how the results of this article can be used for the interpolation by the root of a polynomial. It is demonstrated that a polynomial of relatively small degree provides a very good approximation for the energy in the intermediate coupling region. © 1995 John Wiley & Sons, Inc.

International Journal of Quantum Chemistry•Mar 5, 1995•10 citations

Solution of the one-dimensional Hubbard model for arbitrary electron density and largeU

J. M. P. Carmelo, D. Baeriswyl

The one-dimensional Hubbard model is studied for general electron densities and large U. On the basis of the coupled integral equations of Lieb and Wu exact asymptotic expansions are derived to fifth order in ${U}^{\mathrm{\ensuremath{-}}1}$ for the ground-state energy, the kinetic energy, and the number of doubly occupied sites. A closed-form expression for the distribution functions of charge and spin degrees of freedom is proposed which yields the same large-U expansion as that derived directly from the Lieb-Wu equations. Our analytical results are in excellent agreement with previous numerical studies.

Physical review. B, Condensed matter•May 1, 1988•77 citations

Falicov-Kimball model and its relation to the Hubbard model: Studies on clusters

Romuald L yżwa, Zbigniew Domański

Ground-state properties of the asymmetric Hubbard model are studied on one-dimensional clusters (rings) with N=6, 8, and 10 sites. The ground-state energy, correlation functions, and phase diagrams are determined for various sets of model parameters. Computations are performed both by exact diagonalization (for N=6) and by an approximate method by which correlation effects can be examined on larger clusters than exact diagonalization allows. In the limiting cases (one limit corresponds to the spinless Falicov-Kimball model and the other to the Hubbard model) our results agree quite well with those obtained analytically for infinite systems.

Physical review. B, Condensed matter•Oct 15, 1994•14 citations
Using a straightforward extension of the analysis of Lieb and Wu, we derive a simple analytic form for the ground state energy of a one-dimensional Hubbard ring in the atomic limit. This result is valid for an \textit{arbitrary} number of lattice sites LLL and electrons N≤LN \leq LN≤L. Furthermore, our analysis, including an application of the theory of stochastic matrices, provides insight into the degeneracy and spin properties of the ground states in the atomic limit. We give numerical results which illustrate how the atomic limit is approached.
We study the one-dimensional Hubbard model for interacting electrons in a crystal lattice, using two methods: solution of the Lieb-Wu equations in the atomic limit of no intersite hopping, and variational techniques to obtain a lower bound for the ground state energy. From the Lieb-Wu equations we derive an exact expression for the ground state energy in the atomic limit. We obtain our result by exploiting the spin symmetry of the Hubbard Hamiltonian. Our energy is lower than that obtained by previous published results, and we are able to explain why our results differ. Our method also gives information about the magnetic ordering of the ground state as a function of the number of electrons. Our results are valid for arbitrary electron density values. The variational techniques we use are applied to the two-particle reduced-density matrix (2-RDM), because this quantity suffices to obtain expectation values of all one- and two-particle operators, and because it is less computationally intensive than obtaining the complete N -particle wave function. Because the 2-RDM must be positive semidefinite, we use an efficient programming algorithm known as the Semidefinite Programming Algorithm. A naive search through the space of two-particle matrices encounters a difficult problem, known as the " N -representability problem", which is to assure that the trial density matrices actually come from a wave function that has the required exchange symmetry. This problem is ameliorated by imposing constraints on the variational search that eliminate trial matrices with the incorrect exchange symmetry. These constraints express known physical properties of the ground state (for example its total spin and its translational symmetry) and the so-called D, Q , and G two-positivity conditions that express the required exchange symmetry. To calculate matrix elements, we use the two-particle eigenstates of the reduced Hubbard Hamiltonian; this choice of basis significantly decreases the number of matrix elements needed for energy calculations. After an optimal two-particle density matrix has been obtained, we calculate additional properties of the ground state, including spin correlation functions, chemical potentials, and the momentum distribution. We compare all of these properties to values obtained by exact numerical solution of the Schrödinger equation.
The ground state of the one-dimensional Hubbard Hamiltonian is discussed on the basis of the Bethe ansatz solution of Lieb and Wu. A simple analytical representation of the distribution functions for charge and spin degrees of freedom is derived. The elementary excitations are associated with electron-hole pairs of two pseudo-fermion systems, one for the charge, the other for the magnetic excitations. Spinons and holons correspond to ground state configurations of odd-numbered rings. It is argued that holons will form bound pairs in weakly coupled chains.
The exact phase diagrams of the one-dimensional Hubbard model, both attractive and repulsive, in the presence of an arbitrary magnetic field h for various electron concentrations n and on-site interaction strengths U <0 or U >0 are calculated and investigated. The exact ground-state properties, namely, the ground-state energy, the average spin (magnetization), the concentration of the doubly occupied sites, the kinetic energy, the chemical potential, the spin (magnetic) susceptibility and the charge compressibility, are calculated and examined over a wide range of interaction strengths U for various h and n . It is found that the spin susceptibility at half-filling is non-analytic and changes discontinuously as U →0. The exact theory shows the absence of a charge energy gap in the U ⩽0 region for all n and provides, for the chemical potential, the rigorous upper and lower bounds for half-filled and empty bands respectively. The analytical results derived for the weak-coupling limit and asymptotic expansions for strong coupling are in full agreement with the numerical calculations.
Abstract The perturbation expansion coefficients for the ground‐state energy of the half‐filled one‐dimensional Hubbard model with N = 4ν + 2, (ν = 1, 2,…) sites and satisfying cyclic boundary conditions were calculated in the Hückel limit (U/β → 0), where β designates the one‐electron hopping or resonance integral, and U, the two‐electron on‐site Coulomb integral. This was achieved by solving—order by order—the Lieb–Wu equations, a system of transcendental equations that determines the set of quasi‐momenta {ki} and spin variables {τα} for this model. The exact values for these quantities were found for the N = 6 member ring up to the 20th order in terms of the coupling constant B = U/2β, as well as numerically for 10 ⩽ N ⩽ 50, and the N = 6 Lieb–Wu system was reduced to a system of sextic algebraic equations. These results are compared with those of the infinite system derived analytically by Misurkin and Ovchinnikov [Teor. Mat. Fiz. 11, 127 (1972)]. It is further shown how the results of this article can be used for the interpolation by the root of a polynomial. It is demonstrated that a polynomial of relatively small degree provides a very good approximation for the energy in the intermediate coupling region. © 1995 John Wiley & Sons, Inc.
The one-dimensional Hubbard model is studied for general electron densities and large U. On the basis of the coupled integral equations of Lieb and Wu exact asymptotic expansions are derived to fifth order in ${U}^{\mathrm{\ensuremath{-}}1}$ for the ground-state energy, the kinetic energy, and the number of doubly occupied sites. A closed-form expression for the distribution functions of charge and spin degrees of freedom is proposed which yields the same large-U expansion as that derived directly from the Lieb-Wu equations. Our analytical results are in excellent agreement with previous numerical studies.
Ground-state properties of the asymmetric Hubbard model are studied on one-dimensional clusters (rings) with N=6, 8, and 10 sites. The ground-state energy, correlation functions, and phase diagrams are determined for various sets of model parameters. Computations are performed both by exact diagonalization (for N=6) and by an approximate method by which correlation effects can be examined on larger clusters than exact diagonalization allows. In the limiting cases (one limit corresponds to the spinless Falicov-Kimball model and the other to the Hubbard model) our results agree quite well with those obtained analytically for infinite systems.