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:2606.17933

Three-dimensional repulsive Hubbard model

Alexei Sherman

Published Jun 16, 2026• 0 citations
arXivHTMLPDF

Abstract

The three-dimensional repulsive Hubbard model is investigated using the strong-coupling diagram technique. For half-filling, the boundary between paramagnetic and antiferromagnetic states is determined for the range of the Hubbard repulsion $4t\leq U\leq12t$, where $t$ is the hopping integral between neighboring sites. Along this boundary, the density of states is calculated, and it demonstrates the Mott transition at $U\approx9t$. For $U\geq6t$ and half-filling, the density of states has the shape inherent in the strong electron repulsion, while for $U=4t$, its shape points to weak coupling. The dependence of the Néel temperature $T_{\rm N}$ on the electron concentration $\bar{n}$ is investigated for the cases $U=4t$ and $12t$. In the former case, $T_{\rm N}$ decreases monotonously with $\bar{n}$, while in the latter case, there is a plateau in the dependence near $\bar{n}=0.87$. The plateau is connected to a reconstruction of the density of states caused by an effective weakening of electron coupling due to electron depopulation. For $U=12t$ and half-filling, the magnetic critical exponent $γ\approx1.4$, which is close to the value in the Heisenberg model. Some features resembling the first-order phase transition, revealing themselves in a finite crystal, are discussed.

DOI
10.48550/arxiv.2606.17933

Related papers

Properties of the half-filled Hubbard model investigated by the strong coupling diagram technique

A. Sherman

The equation for the electron Green's function of the fermionic Hubbard model, derived using the strong coupling diagram technique, is solved self-consistently for the near-neighbor form of the kinetic energy and for half-filling. In this case the Mott transition occurs at the Hubbard repulsion Uc ≈ 6.96t, where t is the hopping constant. The calculated spectral functions, density of states (DOS) and momentum distribution are compared with results of Monte Carlo simulations. A satisfactory agreement was found for U > Uc and for temperatures, at which magnetic ordering and spin correlations are suppressed. For U < Uc and lower temperatures the theory describes qualitatively correctly the positions and widths of spectral continua, variations of spectral shapes and occupation numbers with changing wave vector and repulsion. The locations of spectral maxima turn out to be close to the positions of δ-function peaks in the Hubbard-I approximation.

arXiv:1410.8289•International Journal of Modern Physics B•Apr 15, 2015•8 citations

Two-dimensional extended Hubbard model at half-filling

A. Sherman

Abstract We consider the extended Hubbard model on a two-dimensional square lattice at half-filling. The model is investigated using the strong coupling diagram technique. We sum infinite series of ladder diagrams allowing for full-scale charge and spin fluctuations and the actual short-range antiferromagnetic order for nonzero temperatures. In agreement with earlier results, we find the first-order phase transition in the charge subsystem occurring at v = v c ≳ U/4 with v and U the intersite and on-site Coulomb repulsion constants. The transition reveals itself in an abrupt sign change of a sharp maximum in the zero-frequency charge susceptibility at the corner of the Brillouin. States arising at the transition have alternating deviations of electron occupations from the mean value on neighboring sites. Due to fluctuations, these alternating occupation deviations have short-range order. For the considered parameters, such behavior is found for U ≲ 5t with t the hopping constant. For the insulating case U ≳ 6t, in which the transition is not observed, we find a continuous growth of the Mott gap with v. The evolution of the electron density of states with increasing v is also considered.

arXiv:2302.01062•Physica Scripta•May 4, 2023•3 citations

The t − t ′− t ″− U Hubbard model and Fermi-level peak

A Sherman

Abstract Using the strong-coupling diagram technique, low-temperature spectral properties of the two-dimensional fermionic Hubbard model are considered for strong and moderate Hubbard repulsions U . The electron hopping to the nearest, second and third neighbors is taken into account with hopping constants t , t ′ = − 0.3 t and t ″ = 0.2 t , respectively. The nonzero values of t ′ and t ″ lead to strong asymmetry in magnetic properties with respect to the hole and electron doping—for U = 8 t strong antiferromagnetic correlations are retained up to the electron concentration n ¯ ≈ 1.25 , while they are destroyed completely at n ¯ ≈ 0.87 . When the temperature is decreased to T ≲ 0.1 t , in a wide range of electron concentrations there appear narrow and intensive peaks at the Fermi level (FL) in densities of states. For U ≲ 6 t the peaks are seen even at half-filling, while for larger U they arise as the FL leaves the Mott gap. The peaks are connected with a narrow band emerging at low temperatures. We identify states forming the band with spin-polaron excitations—bound states of correlated electrons and mobile spin excitations. Obtained low-temperature spectral functions are used for interpreting the peak-dip-hump structure observed in the photoemission of Nd 2− x Ce x CuO 4 . In the case of hole doping, the calculated Fermi contour contains arcs near nodal points with pseudogaps near antinodal points, while for electron doping the spectral intensity is suppressed at hot spots, in agreement with experimental observations in cuprate perovskites.

arXiv:1810.04967•Physica Scripta•Feb 18, 2019•10 citations

Ground-state phase diagram of the triangular lattice Hubbard model by the density-matrix renormalization group method

Tomonori Shirakawa, Takami Tohyama, J. Kokalj +2 more

Two-dimensional density-matrix renormalization group method is employed to examine the ground-state phase diagram of the Hubbard model on the triangular lattice at half-filling. The calculation reveals two discontinuities in the double occupancy with increasing the repulsive Hubbard interaction $U$ at ${U}_{\mathrm{c}1}\ensuremath{\sim}7.8t$ and ${U}_{\mathrm{c}2}\ensuremath{\sim}9.9t$ ($t$ being the hopping integral), indicating that there are three phases separated by first-order transitions. The absence of any singularity in physical quantities for $0\ensuremath{\le}U<{U}_{\mathrm{c}1}$ implies a metallic phase in this regime. For $U>{U}_{\mathrm{c}2}$, the local spin density induced by an applied pinning magnetic field exhibits a three sublattice feature, which is compatible with the ${120}^{\ensuremath{\circ}}$ N\'eel-ordered state realized in the limit of $U\ensuremath{\rightarrow}\ensuremath{\infty}$. For ${U}_{\mathrm{c}1}<U<{U}_{\mathrm{c}2}$, a response to the applied pinning magnetic field is comparable to that in the metallic phase with a relatively large spin correlation length, but showing neither valence bond nor chiral magnetic order, which therefore resembles gapless spin liquid. However, the spin structure factor for the intermediate phase exhibits the maximum at the $\mathrm{K}$ and ${\mathrm{K}}^{\ensuremath{'}}$ points in the momentum space, which is not compatible to spin liquid with a large spinon Fermi surface. The calculation also finds that the pairing correlation function monotonically decreases with increasing $U$ and thus the superconductivity is unlikely in the intermediate phase.

arXiv:1606.06814•Physical review. B./Physical review. B•Nov 15, 2017•88 citations

Two-dimensional extended Hubbard model: doping, next-nearest neighbor hopping and phase diagrams

A. Sherman

Abstract Using the strong coupling diagram technique, we investigate the extended Hubbard model on a two-dimensional square lattice. This approach allows for charge and spin fluctuations and a short-range antiferromagnetic order at nonzero temperatures. The model features the first-order phase transition to states with alternating site occupations (SAO) at the intersite repulsion v = v c . In this work, we show that doping decreases v c . For a nonzero next-nearest neighbor hopping t ′ , less mobile carriers produce a stronger fall in v c . For half-filling and t ′ = 0 , we consider phase diagrams for a fixed temperature T, on-site U, and intersite repulsions. The diagrams contain regions of SAO, Mott insulator, and several metallic states distinguished by their densities of states. Two of them are characterized by a dip and a peak at the Fermi level. The dip originates from the Slater mechanism for itinerant electrons, while the peak to bound states of electrons with localized magnetic moments. The existence of these two metallic regions in the phase diagram is a manifestation of the Pomeranchuk effect. The boundary between these regions reveals itself as a kink in the curve v c (T) and as a maximum in the T dependence of the double occupancy D. Our calculated D for different values of v, U, and T are in semiquantitative agreement with Monte Carlo results.

arXiv:2304.14683•Physica Scripta•Oct 4, 2023•8 citations

Ground-State Phase Diagram of a Half-Filled One-Dimensional Extended Hubbard Model

Eric Jeckelmann

The density-matrix renormalization group is used to study the phase diagram of the one-dimensional half-filled Hubbard model with on-site (U) and nearest-neighbor (V) repulsion and hopping t. A critical line V(c)(U) approximately U/2 separates a Mott insulating phase from a charge-density-wave phase. The formation of bound charge excitations for V>2t changes the phase transition from continuous to first-order at a tricritical point U(t) approximately 3.7t, V(t)=2t. A frustrating effective antiferromagnetic spin coupling induces a bond-order-wave phase on the critical line V(c)(U) for U(t)<U less, similar 7t.

arXiv:cond-mat/0204244•Physical Review Letters•Nov 15, 2002•133 citations
The three-dimensional repulsive Hubbard model is investigated using the strong-coupling diagram technique. For half-filling, the boundary between paramagnetic and antiferromagnetic states is determined for the range of the Hubbard repulsion 4t≤U≤12t4t\leq U\leq12t4t≤U≤12t, where ttt is the hopping integral between neighboring sites. Along this boundary, the density of states is calculated, and it demonstrates the Mott transition at U≈9tU\approx9tU≈9t. For U≥6tU\geq6tU≥6t and half-filling, the density of states has the shape inherent in the strong electron repulsion, while for U=4tU=4tU=4t, its shape points to weak coupling. The dependence of the Néel temperature TNT_{\rm N}TN​ on the electron concentration nˉ\bar{n}nˉ is investigated for the cases U=4tU=4tU=4t and 12t12t12t. In the former case, TNT_{\rm N}TN​ decreases monotonously with nˉ\bar{n}nˉ, while in the latter case, there is a plateau in the dependence near nˉ=0.87\bar{n}=0.87nˉ=0.87. The plateau is connected to a reconstruction of the density of states caused by an effective weakening of electron coupling due to electron depopulation. For U=12tU=12tU=12t and half-filling, the magnetic critical exponent γ≈1.4γ\approx1.4γ≈1.4, which is close to the value in the Heisenberg model. Some features resembling the first-order phase transition, revealing themselves in a finite crystal, are discussed.
The equation for the electron Green's function of the fermionic Hubbard model, derived using the strong coupling diagram technique, is solved self-consistently for the near-neighbor form of the kinetic energy and for half-filling. In this case the Mott transition occurs at the Hubbard repulsion Uc ≈ 6.96t, where t is the hopping constant. The calculated spectral functions, density of states (DOS) and momentum distribution are compared with results of Monte Carlo simulations. A satisfactory agreement was found for U > Uc and for temperatures, at which magnetic ordering and spin correlations are suppressed. For U < Uc and lower temperatures the theory describes qualitatively correctly the positions and widths of spectral continua, variations of spectral shapes and occupation numbers with changing wave vector and repulsion. The locations of spectral maxima turn out to be close to the positions of δ-function peaks in the Hubbard-I approximation.
Abstract We consider the extended Hubbard model on a two-dimensional square lattice at half-filling. The model is investigated using the strong coupling diagram technique. We sum infinite series of ladder diagrams allowing for full-scale charge and spin fluctuations and the actual short-range antiferromagnetic order for nonzero temperatures. In agreement with earlier results, we find the first-order phase transition in the charge subsystem occurring at v = v c ≳ U/4 with v and U the intersite and on-site Coulomb repulsion constants. The transition reveals itself in an abrupt sign change of a sharp maximum in the zero-frequency charge susceptibility at the corner of the Brillouin. States arising at the transition have alternating deviations of electron occupations from the mean value on neighboring sites. Due to fluctuations, these alternating occupation deviations have short-range order. For the considered parameters, such behavior is found for U ≲ 5t with t the hopping constant. For the insulating case U ≳ 6t, in which the transition is not observed, we find a continuous growth of the Mott gap with v. The evolution of the electron density of states with increasing v is also considered.
Abstract Using the strong-coupling diagram technique, low-temperature spectral properties of the two-dimensional fermionic Hubbard model are considered for strong and moderate Hubbard repulsions U . The electron hopping to the nearest, second and third neighbors is taken into account with hopping constants t , t ′ = − 0.3 t and t ″ = 0.2 t , respectively. The nonzero values of t ′ and t ″ lead to strong asymmetry in magnetic properties with respect to the hole and electron doping—for U = 8 t strong antiferromagnetic correlations are retained up to the electron concentration n ¯ ≈ 1.25 , while they are destroyed completely at n ¯ ≈ 0.87 . When the temperature is decreased to T ≲ 0.1 t , in a wide range of electron concentrations there appear narrow and intensive peaks at the Fermi level (FL) in densities of states. For U ≲ 6 t the peaks are seen even at half-filling, while for larger U they arise as the FL leaves the Mott gap. The peaks are connected with a narrow band emerging at low temperatures. We identify states forming the band with spin-polaron excitations—bound states of correlated electrons and mobile spin excitations. Obtained low-temperature spectral functions are used for interpreting the peak-dip-hump structure observed in the photoemission of Nd 2− x Ce x CuO 4 . In the case of hole doping, the calculated Fermi contour contains arcs near nodal points with pseudogaps near antinodal points, while for electron doping the spectral intensity is suppressed at hot spots, in agreement with experimental observations in cuprate perovskites.
Two-dimensional density-matrix renormalization group method is employed to examine the ground-state phase diagram of the Hubbard model on the triangular lattice at half-filling. The calculation reveals two discontinuities in the double occupancy with increasing the repulsive Hubbard interaction UUU at ${U}_{\mathrm{c}1}\ensuremath{\sim}7.8t$ and ${U}_{\mathrm{c}2}\ensuremath{\sim}9.9t$ (ttt being the hopping integral), indicating that there are three phases separated by first-order transitions. The absence of any singularity in physical quantities for $0\ensuremath{\le}U<{U}_{\mathrm{c}1}$ implies a metallic phase in this regime. For U>Uc2U>{U}_{\mathrm{c}2}U>Uc2​, the local spin density induced by an applied pinning magnetic field exhibits a three sublattice feature, which is compatible with the ${120}^{\ensuremath{\circ}}$ N\'eel-ordered state realized in the limit of $U\ensuremath{\rightarrow}\ensuremath{\infty}$. For Uc1<U<Uc2{U}_{\mathrm{c}1}<U<{U}_{\mathrm{c}2}Uc1​<U<Uc2​, a response to the applied pinning magnetic field is comparable to that in the metallic phase with a relatively large spin correlation length, but showing neither valence bond nor chiral magnetic order, which therefore resembles gapless spin liquid. However, the spin structure factor for the intermediate phase exhibits the maximum at the K\mathrm{K}K and ${\mathrm{K}}^{\ensuremath{'}}$ points in the momentum space, which is not compatible to spin liquid with a large spinon Fermi surface. The calculation also finds that the pairing correlation function monotonically decreases with increasing UUU and thus the superconductivity is unlikely in the intermediate phase.
Abstract Using the strong coupling diagram technique, we investigate the extended Hubbard model on a two-dimensional square lattice. This approach allows for charge and spin fluctuations and a short-range antiferromagnetic order at nonzero temperatures. The model features the first-order phase transition to states with alternating site occupations (SAO) at the intersite repulsion v = v c . In this work, we show that doping decreases v c . For a nonzero next-nearest neighbor hopping t ′ , less mobile carriers produce a stronger fall in v c . For half-filling and t ′ = 0 , we consider phase diagrams for a fixed temperature T, on-site U, and intersite repulsions. The diagrams contain regions of SAO, Mott insulator, and several metallic states distinguished by their densities of states. Two of them are characterized by a dip and a peak at the Fermi level. The dip originates from the Slater mechanism for itinerant electrons, while the peak to bound states of electrons with localized magnetic moments. The existence of these two metallic regions in the phase diagram is a manifestation of the Pomeranchuk effect. The boundary between these regions reveals itself as a kink in the curve v c (T) and as a maximum in the T dependence of the double occupancy D. Our calculated D for different values of v, U, and T are in semiquantitative agreement with Monte Carlo results.
The density-matrix renormalization group is used to study the phase diagram of the one-dimensional half-filled Hubbard model with on-site (U) and nearest-neighbor (V) repulsion and hopping t. A critical line V(c)(U) approximately U/2 separates a Mott insulating phase from a charge-density-wave phase. The formation of bound charge excitations for V>2t changes the phase transition from continuous to first-order at a tricritical point U(t) approximately 3.7t, V(t)=2t. A frustrating effective antiferromagnetic spin coupling induces a bond-order-wave phase on the critical line V(c)(U) for U(t)<U less, similar 7t.