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

Physical review. B, Condensed matter

Stability of the split-band solution and energy gap in the narrow-band region of the Hubbard model

Tadashi Arai, Morrel H. Cohen

Published Apr 15, 1980• 5 citations
Publisher

Abstract

By inserting quasielectron energies $\ensuremath{\omega}$ calculated from the fully renormalized Green's function of the Hubbard model obtained in the preceding paper into the exact expression of Galitskii and Migdal, the ground-state energy, the chemical potential, and the dynamic- and thermodynamic-stability conditions are calculated in the narrow-band region. The results show that as long as the interaction energy $I$ is finite, electrons in the narrow-band region do not obey the Landau theory of Fermi liquids, and a gap appears between the lowest quasielectron energy $\ensuremath{\omega}$ and the chemical potential $\ensuremath{\mu}$ for any occupation $n$, regardless of whether the lower band is exactly filled or not. This unusual behavior is possible because, when an electron is added to the system of $N$ electrons, the whole system relaxes due to the strong interaction, introducing a relaxation energy difference between the two quantities. We also show that all previous solutions which exhibit the split-band structure, including Hubbard's work, yield the same conclusion that electrons do not behave like Landau quasiparticles. However, the energy gap is calculated to be negative at least for some occupations $n$, demonstrating the dynamic instability of those solutions. They also exhibit thermodynamic instability for certain occupations, while the fully renormalized solution, having sufficient electron correlations built in, satisfies the dynamic and thermodynamic stability conditions for all occupations. When the lower band is nearly filled, the nature of the solution is shown to change, making the coherent motion of electrons with fixed $k$ values more difficult. In the pathological limit where $I=\ensuremath{\infty}$, however, the gap vanishes, yielding a metallic state.

DOI
10.1103/physrevb.21.3309

Related papers

Quasiparticle excitations, the relaxation, and transport properties in the narrow-band region of the Hubbard model

Tadashi Arai

Assuming that, in the narrow-band region, the exact Green's function of the Hubbard model essentially retains the two-peak structure proposed by Hubbard, we derive its consequences as rigorously as possible. When a quasielectron is added to the ground state of $N$ electrons in the narrow-band region of the Hubbard model, the whole system is shown to relax due to the strong interaction $I$, thus yielding a relaxation energy $\ensuremath{\Delta}E$ between the quasielectron and the original $N$ electrons. This explains the conclusion in the preceding paper that a gap $\ensuremath{\Delta}E$ analogous to the energy gap in superconductivity appears between the quasielectron spectrum and the chemical potential for any occupation $n$ regardless of whether the lower band is exactly filled or not and that the system cannot be a normal metal. We further found that the ground state is, in fact, a bound state, and electrons involved do not obey the Landau theory of quasiparticles. However, the electron Green's function describes the motion of an electron added to this $N$-electron ground state without the relaxation of the original $N$ electrons, and hence this added electron behaves exactly like a Landau quasiparticle and can be represented by the quasiparticle operator ${\mathrm{A}}_{k\ensuremath{\sigma}}^{\ifmmode\dagger\else\textdagger\fi{}}$ introduced by Luttinger and Nozi\`eres. Only if the original $N$ electrons are allowed to relax is the energy of the system reduced by $\ensuremath{\Delta}E$, yielding the energy gap $\ensuremath{\Delta}E$. When a quasielectron is removed from the system of $N+1$ electrons, the same relaxation energy $\ensuremath{\Delta}E$ is removed, and again the gap $\ensuremath{\Delta}E$ is introduced between the chemical potential $\ensuremath{\mu}$ and the quasihole spectrum, thus creating a gap $2\ensuremath{\Delta}E$ between the quasielectron and quasihole spectra for any occupation $n$. The Fermi level, lying at the center of the gap, never intersects the quasiparticle spectra and therefore, Luttinger's theorem that the volume within the Fermi surface is independent of interactions and equal to the value in the noninteracting limit does not apply to the Hubbard model in the narrow-band region. However, this conclusion does not immediately imply that the system is an insulator. Instead, it is shown that a current-carrying state analogous to the condensate state found by Fr\"ohlich appears under certain conditions. When the number of electrons $N$ is nearly or exactly equal to the number of atoms ${N}_{a}$, quasiparticlelike behavior disappears completely, making the system insulating. The only exception is the pathological limit $I=\ensuremath{\infty}$, where a metallic state with a Fermi volume twice that in the noninteracting limit is found.

Physical review. B, Condensed matter•Apr 15, 1980•8 citations

The ground state problem in the infinite-U Hubbard model

E. V. Kuz’min

The problem of the ground state of the electronic system in the Hubbard model for U=∞ is discussed. The author investigates the normal (singlet or nonmagnetic) N state of the electronic system over the entire range of electron densities n⩽1. It is shown that the energy of the N state ɛ 0 (1) (n) in a one-particle approximation, such as (e.g.) the extended Hartree-Fock approximation, is lower than the energy of the saturated ferromagnetic FM state ɛ FM(n) for all n. The dynamic magnetic susceptibility is calculated in the random phase approximation, and it is shown that the N state is stable over the entire range of electron densities: The static susceptibility (ω=0) does not have a band singularity in the zero-wave vector limit q→0. A formally exact representation is obtained for the mass operator of the one-particle Green’s function, and an approximation of this operator is proposed: M k(E)⋍λF(E), where λ=n(1−n)/(1−n/2)z is the kinematic interaction parameter, z is the number of nearest neighbors, and F(E) is the total single-site Green’s function. For an elliptical density of states the integral equation for F(E) is solved exactly, ad it is shown that the spectral intensity rigorously satisfies the sum rule. The calculated energy of the strongly correlated N state ɛ 0(n)<ɛ FM(n) for all n, and in light of this relationship the author discusses the hypothesis that the ground state of the system is the normal (singlet) state in the thermodynamic limit. The electron distribution function at T=0 differs significantly from the Fermi step; it is “smeared” along the entire energy spectrum, and discontinuities do not occur in the region of the chemical potential m.

Physics of the Solid State•Feb 1, 1997•3 citations

Perturbational treatment of correlation effects in the Hubbard model

G. Bulk, Rainer J. Jelitto

In the general frame of a modified many-body perturbation theory, the single-band Hubbard model is treated in the limit of weak electron-electron interaction. In this modification the decomposition of the Hamiltonian into an effective free system and an effective perturbation is rearranged iteratively in each order of the formalism, resulting in a successive reduction of the interaction part. In comparison with conventional perturbation techniques, this means an extended range of validity. Imposing thermodynamic self-consistency, we calculate the wave-vector-dependent electronic self-energy up to second order and analyze the influence of the Coulomb correlations. Several well-known physical properties of the Hubbard model are recovered in characteristic single-particle quantities. Therefore, the effect of the band occupation n and the electron-electron repulsion U on spectral densities, band structures, densities of states, magnetization, etc. is investigated in detail. Some of the most remarkable results are: first, a splitting of the free Bloch band into two quasiparticle subbands when the coupling strength reaches moderate values. These subbands are linked by a small damping-induced background in the density of states and are separated roughly by an energy of the order of U. In addition, the upper subband is shifted to higher energies in the spectrum proportionally to U while the lower one is nearly fixed. Second, we found a paramagnetic-ferromagnetic phase diagram that has a critical coupling strength ${U}_{c}$ and exhibits strong indications of a lower critical occupation number ${n}_{c}$. There is no spontaneous magnetization in the system for parameters smaller than these values. Compared with the Stoner model, one thus obtains striking improvements. In many cases our conclusions for the weak-coupling limit are qualitatively in good agreement with well-known results in the strong-correlation regime.

Physical review. B, Condensed matter•Jan 1, 1990•34 citations

Functional-derivative study of the Hubbard model. III. Fully renormalized Green's function

Tadashi Arai, Morrel H. Cohen

The functional-derivative method of calculating the Green's function developed earlier for the Hubbard model is generalized and used to obtain a fully renormalized solution. Higher-order functional derivatives operating on the basic Green's functions, $G$ and $\ensuremath{\Gamma}$, are all evaluated explicitly, thus making the solution applicable to the narrow-band region as well as the wide-band region. Correction terms $\ensuremath{\Phi}$ generated from functional derivatives of equal-time Green's functions of the type $\frac{{\ensuremath{\delta}}^{n}〈N〉}{\ensuremath{\delta}{\ensuremath{\epsilon}}^{n}}$, etc., with $n\ensuremath{\ge}2$. It is found that the $\ensuremath{\Phi}'\mathrm{s}$ are, in fact, renormalization factors involved in the self-energy $\ensuremath{\Sigma}$ and that the structure of the $\ensuremath{\Phi}'\mathrm{s}$ resembles that of $\ensuremath{\Sigma}$ and contains the same renormalization factors $\ensuremath{\Phi}$. The renormalization factors $\ensuremath{\Phi}$ are shown to satisfy a set of equations and can be evaluated self-consistently. In the presence of the $\ensuremath{\Phi}'\mathrm{s}$, all difficulties found in the previous results (papers I and II) are removed, and the energy spectrum $\ensuremath{\omega}$ can now be evaluated for all occupations $n$. The Schwinger relation is the only basic relation used in generating this fully self-consistent Green's function, and the Baym-Kadanoff continuity condition is automatically satisfied.

Physical review. B, Condensed matter•Apr 15, 1980•7 citations

Evolution of Fermi-liquid behavior with doping in the Hubbard model: Influence of the band structure

Jungsoo Kim, Dermot Coffey

We calculate the single-particle Green's function for a contact interaction with nearest-neighbor hopping on a square lattice as a function of chemical potential $\ensuremath{\mu}.$ This allows us to investigate the dependence of the leading Fermi-liquid dependencies on the band structure as the Fermi surface evolves from a circle at $\ensuremath{\mu}\ensuremath{\sim}\ensuremath{-}4t$ to a square at $\ensuremath{\mu}=0.$ The form of the single-particle self-energy $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E)$ is determined by the density-density correlation function $\ensuremath{\chi}(\stackrel{\ensuremath{\rightarrow}}{q},\ensuremath{\omega})$ which develops two peaks for $\ensuremath{\mu}\ensuremath{\gtrsim}\ensuremath{-}2.5t$ unlike the parabolic band case. Near half-filling, $\ensuremath{\chi}(\stackrel{\ensuremath{\rightarrow}}{q},\ensuremath{\omega})$ becomes independent of $\ensuremath{\omega},$ one-dimensional behavior, at intermediate values of $\ensuremath{\omega},$ which leads to the one-dimensional behavior in $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E).$ However, with $\ensuremath{\mu}<~\ensuremath{-}0.1t$ there is no influence on the Fermi-liquid dependencies from the spin-density wave instability. We find that throughout the doping region $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E)$ remains qualitatively the same as for the isotropic Fermi surface with quantitative differences. The strong $\stackrel{\ensuremath{\rightarrow}}{p}$ and E dependence of the off-shell self-energy $\ensuremath{\Sigma}(p,E)$ found earlier for the parabolic band is recovered for $\ensuremath{\mu}\ensuremath{\lesssim}\ensuremath{-}t$ but deviates from this develop for $\ensuremath{\mu}\ensuremath{\gtrsim}\ensuremath{-}0.1t.$ The resonance peak width of the spectral function $A(\stackrel{\ensuremath{\rightarrow}}{p}E),$ has a linear dependence in ${\ensuremath{\xi}}_{\stackrel{\ensuremath{\rightarrow}}{p}}$ due to the E dependence of the imaginary part of $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E).$ We point out that an accurate detailed form for $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E)$ would be very difficult to recover from angle-resolved photoemission spectroscopy data for the spectral density. Since the leading corrections are determined by the long wavelength particle-hole excitations, the results found here for the Hubbard model carry over to Hamiltonians with finite range interactions.

Physical review. B, Condensed matter•Aug 15, 2000•2 citations

On correlation effects in a narrow-band model with electron–hole asymmetry

L. Didukh, V. Hankevych

Some correlation effects are studied in a narrow-band model with electron–hole asymmetry. A peculiar feature of the model is that it takes into account the hopping integral associated with electron–electron interaction (correlated hopping). The method of Green’s functions is used to find the quasiparticle energy spectrum, as well as expressions for the upper and lower Hubbard subband widths. It is shown that the specific features of the model lead in a number of cases to consequences differing significantly from those following from the Hubbard model. The obtained results mainly illustrate the peculiarities of the model in the low-temperature region and are used for interpreting the physical properties of systems with narrow energy bands.

Low Temperature Physics•May 1, 1999•3 citations
By inserting quasielectron energies \ensuremathω\ensuremath{\omega}\ensuremathω calculated from the fully renormalized Green's function of the Hubbard model obtained in the preceding paper into the exact expression of Galitskii and Migdal, the ground-state energy, the chemical potential, and the dynamic- and thermodynamic-stability conditions are calculated in the narrow-band region. The results show that as long as the interaction energy III is finite, electrons in the narrow-band region do not obey the Landau theory of Fermi liquids, and a gap appears between the lowest quasielectron energy \ensuremathω\ensuremath{\omega}\ensuremathω and the chemical potential \ensuremathμ\ensuremath{\mu}\ensuremathμ for any occupation nnn, regardless of whether the lower band is exactly filled or not. This unusual behavior is possible because, when an electron is added to the system of NNN electrons, the whole system relaxes due to the strong interaction, introducing a relaxation energy difference between the two quantities. We also show that all previous solutions which exhibit the split-band structure, including Hubbard's work, yield the same conclusion that electrons do not behave like Landau quasiparticles. However, the energy gap is calculated to be negative at least for some occupations nnn, demonstrating the dynamic instability of those solutions. They also exhibit thermodynamic instability for certain occupations, while the fully renormalized solution, having sufficient electron correlations built in, satisfies the dynamic and thermodynamic stability conditions for all occupations. When the lower band is nearly filled, the nature of the solution is shown to change, making the coherent motion of electrons with fixed kkk values more difficult. In the pathological limit where I=\ensuremath∞I=\ensuremath{\infty}I=\ensuremath∞, however, the gap vanishes, yielding a metallic state.
Assuming that, in the narrow-band region, the exact Green's function of the Hubbard model essentially retains the two-peak structure proposed by Hubbard, we derive its consequences as rigorously as possible. When a quasielectron is added to the ground state of NNN electrons in the narrow-band region of the Hubbard model, the whole system is shown to relax due to the strong interaction III, thus yielding a relaxation energy $\ensuremath{\Delta}E$ between the quasielectron and the original NNN electrons. This explains the conclusion in the preceding paper that a gap $\ensuremath{\Delta}E$ analogous to the energy gap in superconductivity appears between the quasielectron spectrum and the chemical potential for any occupation nnn regardless of whether the lower band is exactly filled or not and that the system cannot be a normal metal. We further found that the ground state is, in fact, a bound state, and electrons involved do not obey the Landau theory of quasiparticles. However, the electron Green's function describes the motion of an electron added to this NNN-electron ground state without the relaxation of the original NNN electrons, and hence this added electron behaves exactly like a Landau quasiparticle and can be represented by the quasiparticle operator ${\mathrm{A}}_{k\ensuremath{\sigma}}^{\ifmmode\dagger\else\textdagger\fi{}}$ introduced by Luttinger and Nozi\`eres. Only if the original NNN electrons are allowed to relax is the energy of the system reduced by $\ensuremath{\Delta}E$, yielding the energy gap $\ensuremath{\Delta}E$. When a quasielectron is removed from the system of N+1N+1N+1 electrons, the same relaxation energy $\ensuremath{\Delta}E$ is removed, and again the gap $\ensuremath{\Delta}E$ is introduced between the chemical potential $\ensuremath{\mu}$ and the quasihole spectrum, thus creating a gap $2\ensuremath{\Delta}E$ between the quasielectron and quasihole spectra for any occupation nnn. The Fermi level, lying at the center of the gap, never intersects the quasiparticle spectra and therefore, Luttinger's theorem that the volume within the Fermi surface is independent of interactions and equal to the value in the noninteracting limit does not apply to the Hubbard model in the narrow-band region. However, this conclusion does not immediately imply that the system is an insulator. Instead, it is shown that a current-carrying state analogous to the condensate state found by Fr\"ohlich appears under certain conditions. When the number of electrons NNN is nearly or exactly equal to the number of atoms Na{N}_{a}Na​, quasiparticlelike behavior disappears completely, making the system insulating. The only exception is the pathological limit $I=\ensuremath{\infty}$, where a metallic state with a Fermi volume twice that in the noninteracting limit is found.
The problem of the ground state of the electronic system in the Hubbard model for U=∞ is discussed. The author investigates the normal (singlet or nonmagnetic) N state of the electronic system over the entire range of electron densities n⩽1. It is shown that the energy of the N state ɛ 0 (1) (n) in a one-particle approximation, such as (e.g.) the extended Hartree-Fock approximation, is lower than the energy of the saturated ferromagnetic FM state ɛ FM(n) for all n. The dynamic magnetic susceptibility is calculated in the random phase approximation, and it is shown that the N state is stable over the entire range of electron densities: The static susceptibility (ω=0) does not have a band singularity in the zero-wave vector limit q→0. A formally exact representation is obtained for the mass operator of the one-particle Green’s function, and an approximation of this operator is proposed: M k(E)⋍λF(E), where λ=n(1−n)/(1−n/2)z is the kinematic interaction parameter, z is the number of nearest neighbors, and F(E) is the total single-site Green’s function. For an elliptical density of states the integral equation for F(E) is solved exactly, ad it is shown that the spectral intensity rigorously satisfies the sum rule. The calculated energy of the strongly correlated N state ɛ 0(n)<ɛ FM(n) for all n, and in light of this relationship the author discusses the hypothesis that the ground state of the system is the normal (singlet) state in the thermodynamic limit. The electron distribution function at T=0 differs significantly from the Fermi step; it is “smeared” along the entire energy spectrum, and discontinuities do not occur in the region of the chemical potential m.
In the general frame of a modified many-body perturbation theory, the single-band Hubbard model is treated in the limit of weak electron-electron interaction. In this modification the decomposition of the Hamiltonian into an effective free system and an effective perturbation is rearranged iteratively in each order of the formalism, resulting in a successive reduction of the interaction part. In comparison with conventional perturbation techniques, this means an extended range of validity. Imposing thermodynamic self-consistency, we calculate the wave-vector-dependent electronic self-energy up to second order and analyze the influence of the Coulomb correlations. Several well-known physical properties of the Hubbard model are recovered in characteristic single-particle quantities. Therefore, the effect of the band occupation n and the electron-electron repulsion U on spectral densities, band structures, densities of states, magnetization, etc. is investigated in detail. Some of the most remarkable results are: first, a splitting of the free Bloch band into two quasiparticle subbands when the coupling strength reaches moderate values. These subbands are linked by a small damping-induced background in the density of states and are separated roughly by an energy of the order of U. In addition, the upper subband is shifted to higher energies in the spectrum proportionally to U while the lower one is nearly fixed. Second, we found a paramagnetic-ferromagnetic phase diagram that has a critical coupling strength Uc{U}_{c}Uc​ and exhibits strong indications of a lower critical occupation number nc{n}_{c}nc​. There is no spontaneous magnetization in the system for parameters smaller than these values. Compared with the Stoner model, one thus obtains striking improvements. In many cases our conclusions for the weak-coupling limit are qualitatively in good agreement with well-known results in the strong-correlation regime.
The functional-derivative method of calculating the Green's function developed earlier for the Hubbard model is generalized and used to obtain a fully renormalized solution. Higher-order functional derivatives operating on the basic Green's functions, GGG and $\ensuremath{\Gamma}$, are all evaluated explicitly, thus making the solution applicable to the narrow-band region as well as the wide-band region. Correction terms $\ensuremath{\Phi}$ generated from functional derivatives of equal-time Green's functions of the type $\frac{{\ensuremath{\delta}}^{n}〈N〉}{\ensuremath{\delta}{\ensuremath{\epsilon}}^{n}}$, etc., with $n\ensuremath{\ge}2$. It is found that the $\ensuremath{\Phi}'\mathrm{s}$ are, in fact, renormalization factors involved in the self-energy $\ensuremath{\Sigma}$ and that the structure of the $\ensuremath{\Phi}'\mathrm{s}$ resembles that of $\ensuremath{\Sigma}$ and contains the same renormalization factors $\ensuremath{\Phi}$. The renormalization factors $\ensuremath{\Phi}$ are shown to satisfy a set of equations and can be evaluated self-consistently. In the presence of the $\ensuremath{\Phi}'\mathrm{s}$, all difficulties found in the previous results (papers I and II) are removed, and the energy spectrum $\ensuremath{\omega}$ can now be evaluated for all occupations nnn. The Schwinger relation is the only basic relation used in generating this fully self-consistent Green's function, and the Baym-Kadanoff continuity condition is automatically satisfied.
We calculate the single-particle Green's function for a contact interaction with nearest-neighbor hopping on a square lattice as a function of chemical potential $\ensuremath{\mu}.$ This allows us to investigate the dependence of the leading Fermi-liquid dependencies on the band structure as the Fermi surface evolves from a circle at $\ensuremath{\mu}\ensuremath{\sim}\ensuremath{-}4t$ to a square at $\ensuremath{\mu}=0.$ The form of the single-particle self-energy $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E)$ is determined by the density-density correlation function $\ensuremath{\chi}(\stackrel{\ensuremath{\rightarrow}}{q},\ensuremath{\omega})$ which develops two peaks for $\ensuremath{\mu}\ensuremath{\gtrsim}\ensuremath{-}2.5t$ unlike the parabolic band case. Near half-filling, $\ensuremath{\chi}(\stackrel{\ensuremath{\rightarrow}}{q},\ensuremath{\omega})$ becomes independent of $\ensuremath{\omega},$ one-dimensional behavior, at intermediate values of $\ensuremath{\omega},$ which leads to the one-dimensional behavior in $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E).$ However, with $\ensuremath{\mu}<~\ensuremath{-}0.1t$ there is no influence on the Fermi-liquid dependencies from the spin-density wave instability. We find that throughout the doping region $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E)$ remains qualitatively the same as for the isotropic Fermi surface with quantitative differences. The strong $\stackrel{\ensuremath{\rightarrow}}{p}$ and E dependence of the off-shell self-energy $\ensuremath{\Sigma}(p,E)$ found earlier for the parabolic band is recovered for $\ensuremath{\mu}\ensuremath{\lesssim}\ensuremath{-}t$ but deviates from this develop for $\ensuremath{\mu}\ensuremath{\gtrsim}\ensuremath{-}0.1t.$ The resonance peak width of the spectral function $A(\stackrel{\ensuremath{\rightarrow}}{p}E),$ has a linear dependence in ${\ensuremath{\xi}}_{\stackrel{\ensuremath{\rightarrow}}{p}}$ due to the E dependence of the imaginary part of $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E).$ We point out that an accurate detailed form for $\ensuremath{\Sigma}(\stackrel{\ensuremath{\rightarrow}}{p},E)$ would be very difficult to recover from angle-resolved photoemission spectroscopy data for the spectral density. Since the leading corrections are determined by the long wavelength particle-hole excitations, the results found here for the Hubbard model carry over to Hamiltonians with finite range interactions.
Some correlation effects are studied in a narrow-band model with electron–hole asymmetry. A peculiar feature of the model is that it takes into account the hopping integral associated with electron–electron interaction (correlated hopping). The method of Green’s functions is used to find the quasiparticle energy spectrum, as well as expressions for the upper and lower Hubbard subband widths. It is shown that the specific features of the model lead in a number of cases to consequences differing significantly from those following from the Hubbard model. The obtained results mainly illustrate the peculiarities of the model in the low-temperature region and are used for interpreting the physical properties of systems with narrow energy bands.