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:cond-mat/0701060 · Annalen der Physik

Rigorous solution of a Hubbard model extended by nearest-neighbour Coulomb and exchange interaction on a triangle and tetrahedron

Rolf Schumann

Published Feb 7, 2008• 19 citations
arXivHTMLPDF

Abstract

The Hubbard model extended by either nearest-neighbour Coulomb correlation and/or nearest neighbour Heisenberg exchange is solved analytically for a triangle and tetrahedron.All eigenvalues and eigenvectors are given as functions of the model parameters in a closed form.The groundstate crossings and degeneracies are discussed both for the canonical and grand-canonical energy levels.The grand canonical potential Ω(µ, T, h) and the electron occupation N (µ, T, h) of the related cluster gases were calculated for arbitrary values (attractive and repulsive) of the three interaction constants.In the pure Hubbard model we found various steps in N (µ, T = 0, h) higher than one.It is shown that the various degeneracies of the grand-canonical energy levels are partially lifted by an antiferromagnetic exchange interaction, whereas a moderate ferromagnetic exchange modifies only slightly the results of the pure Hubbard model.A repulsive nn Coulomb correlation lifts these degeneracies completely.The relation of the cluster gas results to extended systems is discussed.

DOI
10.1002/andp.200710281

Related papers

Rigorous solution of a Hubbard model extended by nearest‐neighbour Coulomb and exchange interaction on a triangle and tetrahedron

Rolf Schumann

Abstract The Hubbard model extended by both nearest‐neighbour (nn) Coulomb correlation and nearest‐neighbour Heisenberg exchange is solved rigorously for a triangle and tetrahedron. All eigenvalues and eigenvectors are given as functions of the model parameters in a closed analytical form. For fixed electron numbers we found a multitude of level crossings, both in the ground state and in the excited states in dependence on the various model parameters. By coupling an ensemble of clusters to an electron bath we get the cluster gas model or the cluster gas approximation, if an extended array of weak‐interacting clusters is considered. The grand‐canonical potential Ω (μ, T , h ) and the electron occupation N (μ, T , h ) of the related cluster gases were calculated for arbitrary values (attractive and repulsive) of the three interaction constants. For the cluster gases without the additional interactions we found various steps in N (μ, T = 0, h = 0) higher than one. The reason is the degeneration of ground states differing in their electron occupation by more than one electron. For the triangular cluster gas we have one such degeneration point. For the tetrahedral cluster gas two. As a consequence, we do not find areas with one electron in the μ‐ U ground‐state phase diagram of the triangular cluster gas or with one, two and five electrons in the case of the tetrahedral cluster gas. The degeneration point of the triangular cluster gas can not be destroyed by an applied magnetic field. This holds also for the lower degeneration point of the tetrahedral cluster gas. Otherwise, the upper degeneration point breaks down at a critical magnetic field h c . The dependence of h c on U shows a maximum for strong on‐site correlation. The influence of nn‐exchange and nn‐Coulomb correlation on the ground‐state phase diagrams is calculated. Whereas antiferromagnetic nn‐exchange breaks the degeneration points of the tetrahedral cluster gas partially only, a repulsive nn‐Coulomb correlation lifts the underlying degeneracies completely. Otherwise both ferromagnetic nn‐exchange and attractive nn‐Coulomb interaction stabilise the degeneration points. The consequences of the cluster gas results for extended cluster arrays are discussed.

Annalen der Physik•Mar 17, 2008•5 citations

The Hubbard model extended by nearest‐neighbor Coulomb and exchange interaction on a cubic cluster – rigorous and exact results

Rolf Schumann, David Zwicker

Abstract The Hubbard model on a cube was revisited and extended by both nearest‐neighbor Coulomb correlation W and nearest‐neighbor Heisenberg exchange J. The complete eigensystem was computed exactly for all electron occupancies and all model parameters ranging from minus infinity to plus infinity. For two electrons on the cluster the eigensystem is given in analytical form. For six electrons and infinite on‐site correlation U we determinded the groundstate and the groundstate energy of the pure Hubbard model analytically. For fixed electron numbers we found a multitude of ground state level crossings depending on the various model parameters. Furthermore the groundstates of the pure Hubbard model in dependence on a magnetic field h coupled to the spins are shown for the complete U‐h plane. The critical magnetic field, where the zero spin groundstate breaks down is given for four and six electrons. Suprisingly we found parameter regions, where the ground state spin does not depend monotonously on J in the extended model. For the cubic cluster gas, i.e. an ensemble of clusters coupled to an electron bath, we calculated the density n (μ, T, h) and the thermodynamical density of states from the grand potential. The ground states and the various spin‐spin correlation functions are studied for both attractive and repulsive values of the three interaction constants. We determined the various anomalous degeneration lines, where n (μ, T = 0, h = 0) shows steps higher than one, since in this parameter regions exotic phenomena as phase separation are to expect in extended models. For the cases where these lines end in triple points, i.e. groundstates of three different occupation numbers are degenerated, we give the related parameter values. Regarding the influence of the nn‐exchange and the nn‐Coulomb correlation onto the anomalous degeneration we find both lifting and inducing of degeneracies depending on the parameter values.

arXiv:0906.3589•Annalen der Physik•Apr 30, 2010•12 citations

Rigorous solution of a Hubbard model extended by nearest‐neighbour Coulomb and isotropic exchange interaction on a triangle and tetrahedron

Rolf Schumann

In connection with the study of the extended Hubbard model on an isosceles triangular cluster, we discovered a factor two error in the coding of the Ising part of the Heisenberg term in the Hamiltonian. In consequence the rigorous results remain valid, but not for the isotropic Heisenberg model, but for an anisotropic Heisenberg model.

Annalen der Physik•Apr 15, 2009•1 citation

Cluster approach to the extended Hubbard model

W. N. Mei, Y. C. Lee

The eigenvalue problem of the two-site and three-site extended Hubbard Hamiltonian, which contains the intraatomic repulsion $U$, the nearest-neighbor repulsion $V$, the isotropic exchange $J$, the anisotropic exchange ${J}^{\ensuremath{'}}$, the hopping term $t$, the semihopping term $S$, and the double hopping term $D$, has been solved exactly and the grand partition function is obtained. The magnetic and thermodynamic quantities such as the magnetic susceptibility and the specific heat are calculated for arbitrary electron density. Several features believed to have wider validity have been deduced from this cluster model.

Physical review. B, Condensed matter•Jul 15, 1981•4 citations

Exact solution of a periodic-cluster Hubbard model with an electron-lattice interaction

César R. Proetto, L. M. Falicov

An exact solution of a modified Hubbard model for a tetrahedral cluster with periodic boundary conditions is presented. It includes a static electron-lattice interaction in which the intersite bond-hopping parameter is a function of the bond electron occupation. It is found that the intrasite Coulomb repulsion and the effective attractive electron-lattice interaction have, generally speaking, opposite (although not identical) effects. They cancel each other exactly along a line of singularities in parameter space. The state vectors are easily obtained analytically. A study is made of magnetic and superconducting static correlation functions as functions of the parameters. Various forms of anisotropic superconducting fluctuations are favored in the various regimes. There is an indication of possible states with coexisting nonuniform, nonisotropic superconductivity and antiferromagnetism caused by a ``compromise'' between the competing forces.

Physical review. B, Condensed matter•Jul 15, 1988•9 citations

Localized fermions on the triangular lattice with Ising-like interactions

Lubomíra Regeciová, Konrad Jerzy Kapcia

The model of localized fermions on the triangular lattice is analyzed in means of the Monte Carlo simulations in the grand canonical ensemble. The Hamiltonian of the system has a form of the extended Hubbard model (at the atomic limit) with nearest-neighbor Ising-like magnetic J interactions and onsite Coulomb U interactions. The model is investigated for both signs of J, arbitrary U interaction, and arbitrary chemical potential μ (or, equivalently, arbitrary particle concentration n). Based on the specific heat capacity and sublattice magnetization analyses, the phase diagrams of the model are determined. For ferromagnetic case (J<0), the transition from the ordered phase (which is a standard ferromagnet and can be stable up to k_{B}T/|J|≈0.61) is found to be second order (for sufficiently large temperatures k_{B}T/|J|≳0.2) or first order (for -1<U/|J|<-0.65 at the half-filling, i.e., n=1). In the case of J>0, the ordered phase occurs in a range of -1/2<U/|J|<0 (for n=1), while for larger U the state with short-range order is also found (also for n≠1). The ordered phase is characterized by an antiferromagnetic arrangement of magnetic moments in two sublattices forming the hexagonal lattice. The transition from this ordered phase, which is found also for μ≠0 (n≠1) and U/|J|>-1/2 is always second order for any model parameters. The ordered phase for J>0 can be stable up to k_{B}T/|J|≈0.06.

arXiv:2412.10808•Physical review. E•Jun 3, 2025•1 citation
The Hubbard model extended by either nearest-neighbour Coulomb correlation and/or nearest neighbour Heisenberg exchange is solved analytically for a triangle and tetrahedron.All eigenvalues and eigenvectors are given as functions of the model parameters in a closed form.The groundstate crossings and degeneracies are discussed both for the canonical and grand-canonical energy levels.The grand canonical potential Ω(µ, T, h) and the electron occupation N (µ, T, h) of the related cluster gases were calculated for arbitrary values (attractive and repulsive) of the three interaction constants.In the pure Hubbard model we found various steps in N (µ, T = 0, h) higher than one.It is shown that the various degeneracies of the grand-canonical energy levels are partially lifted by an antiferromagnetic exchange interaction, whereas a moderate ferromagnetic exchange modifies only slightly the results of the pure Hubbard model.A repulsive nn Coulomb correlation lifts these degeneracies completely.The relation of the cluster gas results to extended systems is discussed.
Abstract The Hubbard model extended by both nearest‐neighbour (nn) Coulomb correlation and nearest‐neighbour Heisenberg exchange is solved rigorously for a triangle and tetrahedron. All eigenvalues and eigenvectors are given as functions of the model parameters in a closed analytical form. For fixed electron numbers we found a multitude of level crossings, both in the ground state and in the excited states in dependence on the various model parameters. By coupling an ensemble of clusters to an electron bath we get the cluster gas model or the cluster gas approximation, if an extended array of weak‐interacting clusters is considered. The grand‐canonical potential Ω (μ, T , h ) and the electron occupation N (μ, T , h ) of the related cluster gases were calculated for arbitrary values (attractive and repulsive) of the three interaction constants. For the cluster gases without the additional interactions we found various steps in N (μ, T = 0, h = 0) higher than one. The reason is the degeneration of ground states differing in their electron occupation by more than one electron. For the triangular cluster gas we have one such degeneration point. For the tetrahedral cluster gas two. As a consequence, we do not find areas with one electron in the μ‐ U ground‐state phase diagram of the triangular cluster gas or with one, two and five electrons in the case of the tetrahedral cluster gas. The degeneration point of the triangular cluster gas can not be destroyed by an applied magnetic field. This holds also for the lower degeneration point of the tetrahedral cluster gas. Otherwise, the upper degeneration point breaks down at a critical magnetic field h c . The dependence of h c on U shows a maximum for strong on‐site correlation. The influence of nn‐exchange and nn‐Coulomb correlation on the ground‐state phase diagrams is calculated. Whereas antiferromagnetic nn‐exchange breaks the degeneration points of the tetrahedral cluster gas partially only, a repulsive nn‐Coulomb correlation lifts the underlying degeneracies completely. Otherwise both ferromagnetic nn‐exchange and attractive nn‐Coulomb interaction stabilise the degeneration points. The consequences of the cluster gas results for extended cluster arrays are discussed.
Abstract The Hubbard model on a cube was revisited and extended by both nearest‐neighbor Coulomb correlation W and nearest‐neighbor Heisenberg exchange J. The complete eigensystem was computed exactly for all electron occupancies and all model parameters ranging from minus infinity to plus infinity. For two electrons on the cluster the eigensystem is given in analytical form. For six electrons and infinite on‐site correlation U we determinded the groundstate and the groundstate energy of the pure Hubbard model analytically. For fixed electron numbers we found a multitude of ground state level crossings depending on the various model parameters. Furthermore the groundstates of the pure Hubbard model in dependence on a magnetic field h coupled to the spins are shown for the complete U‐h plane. The critical magnetic field, where the zero spin groundstate breaks down is given for four and six electrons. Suprisingly we found parameter regions, where the ground state spin does not depend monotonously on J in the extended model. For the cubic cluster gas, i.e. an ensemble of clusters coupled to an electron bath, we calculated the density n (μ, T, h) and the thermodynamical density of states from the grand potential. The ground states and the various spin‐spin correlation functions are studied for both attractive and repulsive values of the three interaction constants. We determined the various anomalous degeneration lines, where n (μ, T = 0, h = 0) shows steps higher than one, since in this parameter regions exotic phenomena as phase separation are to expect in extended models. For the cases where these lines end in triple points, i.e. groundstates of three different occupation numbers are degenerated, we give the related parameter values. Regarding the influence of the nn‐exchange and the nn‐Coulomb correlation onto the anomalous degeneration we find both lifting and inducing of degeneracies depending on the parameter values.
The eigenvalue problem of the two-site and three-site extended Hubbard Hamiltonian, which contains the intraatomic repulsion UUU, the nearest-neighbor repulsion VVV, the isotropic exchange JJJ, the anisotropic exchange ${J}^{\ensuremath{'}}$, the hopping term ttt, the semihopping term SSS, and the double hopping term DDD, has been solved exactly and the grand partition function is obtained. The magnetic and thermodynamic quantities such as the magnetic susceptibility and the specific heat are calculated for arbitrary electron density. Several features believed to have wider validity have been deduced from this cluster model.
An exact solution of a modified Hubbard model for a tetrahedral cluster with periodic boundary conditions is presented. It includes a static electron-lattice interaction in which the intersite bond-hopping parameter is a function of the bond electron occupation. It is found that the intrasite Coulomb repulsion and the effective attractive electron-lattice interaction have, generally speaking, opposite (although not identical) effects. They cancel each other exactly along a line of singularities in parameter space. The state vectors are easily obtained analytically. A study is made of magnetic and superconducting static correlation functions as functions of the parameters. Various forms of anisotropic superconducting fluctuations are favored in the various regimes. There is an indication of possible states with coexisting nonuniform, nonisotropic superconductivity and antiferromagnetism caused by a ``compromise'' between the competing forces.
The model of localized fermions on the triangular lattice is analyzed in means of the Monte Carlo simulations in the grand canonical ensemble. The Hamiltonian of the system has a form of the extended Hubbard model (at the atomic limit) with nearest-neighbor Ising-like magnetic J interactions and onsite Coulomb U interactions. The model is investigated for both signs of J, arbitrary U interaction, and arbitrary chemical potential μ (or, equivalently, arbitrary particle concentration n). Based on the specific heat capacity and sublattice magnetization analyses, the phase diagrams of the model are determined. For ferromagnetic case (J<0), the transition from the ordered phase (which is a standard ferromagnet and can be stable up to k_{B}T/|J|≈0.61) is found to be second order (for sufficiently large temperatures k_{B}T/|J|≳0.2) or first order (for -1<U/|J|<-0.65 at the half-filling, i.e., n=1). In the case of J>0, the ordered phase occurs in a range of -1/2<U/|J|<0 (for n=1), while for larger U the state with short-range order is also found (also for n≠1). The ordered phase is characterized by an antiferromagnetic arrangement of magnetic moments in two sublattices forming the hexagonal lattice. The transition from this ordered phase, which is found also for μ≠0 (n≠1) and U/|J|>-1/2 is always second order for any model parameters. The ordered phase for J>0 can be stable up to k_{B}T/|J|≈0.06.