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

Journal of the Physical Society of Japan

Two-Dimensional Hubbard Model –Metal Insulator Transition Studied by Monte Carlo Calculation–

Nobuo Furukawa, Masatoshi Imada

Published Sep 15, 1992• 218 citations
Publisher

Abstract

The Hubbard model on a square lattice in the ground state is investigated. Various aspects of Mott transition at half-filling are clarified from the quantum Monte Carlo study. Critical exponents of the transition for charge and spin correlations are estimated. For the doping concentration δ, the charge susceptibility is proportional to δ -1 indicating the divergence of charge mass for δ→0, while the spin susceptibility is always finite and not singular near δ=0. Several theoretical consequences of the above results are discussed. Incommensurate peak amplitude of spin correlation scales as (δ-δ c ) -1 with δ c <0.01 in disagreement with RPA results ln (δ-δ c ). The antiferromagnetic order at the half-filling also shows strong correlation character and is not sensitive to the shape of the fermi surface, while the incommensurate peak observed away from the half-filling is sensitively suppressed by the loss of partial nesting consistently with the weak coupling picture. Comparisons with experimental indications of high- T c oxides are made.

DOI
10.1143/jpsj.61.3331

Related papers

Charge Gap, Charge Susceptibility and Spin Correlationin the Hubbard Model on a Square Lattice

Nobuo Furukawa, Masatoshi Imada

Quantum Monte Carlo study of the Hubbard model on a square lattice in the ground state is reported. For the on-site interaction U =4 scaled by the transfer, the charge excitation gap at the half-filling is estimated to be Δ c =0.58±0.08. Near half-filling, the charge compressibility κ follows the form κ∝δ -1 for the doping concentration δ, indicating a divergent behavior as the system approaches half-filling, while κ=0 at the half-filling. The incommensurate spin correlation length ξ scales as ξ∝1/\(\sqrt{\delta}\).

Journal of the Physical Society of Japan•Nov 15, 1991•63 citations

Mott metal-insulator transition in the half-filled Hubbard model on the triangular lattice

Massimo Capone, Luca Capriotti, Federico Becca +1 more

We investigate the metal-insulator transition in the half-filled Hubbard model on the two-dimensional triangular lattice using both the Kotliar-Ruckenstein slave-boson technique and an exact numerical diagonalization of finite clusters. Contrary to the case of a square lattice, where a perfect nesting of the Fermi surface leads to a metal-insulator transition at arbitrarily small values of U, always accompanied by antiferromagnetic ordering, on a triangular lattice, due to the lack of perfect nesting, the transition takes place at a finite value of U, and frustration induces a nontrivial competition among different magnetic phases. Indeed, within the mean-field approximation in the slave-boson approach, as the interaction grows the paramagnetic metal turns into a metallic phase with incommensurate spiral ordering. Increasing the interaction further, a linear spin-density wave is stabilized, and finally for strong coupling the latter phase undergoes a first-order transition toward an antiferromagnetic insulator. No trace of the intermediate phases is seen in the exact diagonalization results, indicating a transition between a paramagnetic metal and an antiferromagnetic insulator.

arXiv:cond-mat/0006437•Physical review. B, Condensed matter•Feb 2, 2001•51 citations

Ground-state phase diagram of the half-filled bilayer Hubbard model

Michael Golor, Timo Reckling, Laura Classen +2 more

Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account of the ground-state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite on-site repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonic finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure.

arXiv:1409.1103•Physical Review B•Nov 18, 2014•23 citations

Charge Susceptibility of the Two-Dimensional Hubbard Model near the Half-Filled Band -A Monte Carlo Study-

Hiromi Otsuka

Using the quantum Monte Carlo simulation technique for the lattice fermion systems, we study the charge susceptibility of the two-dimensional Hubbard near the half-filled band. Lattice sizes up to 6×6 and inverse temperatures up to 6 t ( t : hopping parameter) are studied. We find that the behavior of the charge susceptibility is strongly dependent on the electron concentration at low temperature: the charge susceptibility is suppressed at the half-filled band, while it is considerably enhanced in non-half-filled band cases by the electron correlation. The above behavior of the charge susceptibility is quite similar to the recent Bethe-Ansatz result for the one-dimensional Hubbard model.

Journal of the Physical Society of Japan•Aug 15, 1990•10 citations

Filling-Control Metal-Insulator Transition in the Hubbard Model Studied by the Operator Projection Method

Shigeki Onoda, Masatoshi Imada

Thermodynamic and dynamical properties of filling-control metal-insulator transition (MIT) in the Hubbard model are studied by the operator projection method, especially in two dimensions. This is a non-perturbative analytic approach to many-body systems. The present theory incorporates the Mott-Hubbard, Brinkman-Rice and Slater pictures of the MIT into a unified framework, together with reproducing low-energy narrow band arising from spin-charge fluctuations. At half filling, single-particle spectra $A(ω, k)$ show formation of two Hubbard bands and their antiferromagnetic shadows separated by a Mott gap in the plane of energy $ω$ and momentum $k$ with lowering temperatures. These four bands produce splitting to two low-energy narrow bands and two SDW-like bands in the dispersion. Near half filling, the low-energy narrow band persists at low temperatures. This narrow band has a particularly weak dispersion and large weights around $(π, 0)$ and $(0, π)$ momenta. The velocity of these low-energy excitations is shown to vanish towards the MIT, indicating the mass divergence as in the Brinkman-Rice picture, but most prominently around $(π,0)$ and $(0,π)$ with strong momentum dependence. This reflects the suppression of the coherence near the MIT. Main structures in $A(ω, k)$ show remarkable similarities to quantum Monte-Carlo results in two dimensions as well as in the one-dimensional Hubbard model. The charge compressibility appears to diverge with decreasing doping concentration in both one and two dimensions in agreement with the exact and quantum Monte-Carlo results. We also discuss implications of the flat dispersion formed near the Fermi level to the observations in high-$\Tc$ cuprate superconductors.

arXiv:cond-mat/0107291•Journal of the Physical Society of Japan•Nov 1, 2001•33 citations

Magnetic susceptibility of the two-dimensional Hubbard model

Adriana Moreo

Using quantum Monte Carlo techniques we study the magnetic susceptibility of the Hubbard model on a two-dimensional square lattice for values of the coupling U/t ranging from 4 to 20, and as a function of temperature and doping. For U/t=4 the magnetic susceptibility has a maximum at half filling and decreases with doping, in qualitative agreement with recent random-phase-approximation calculations. However, for U/t=10 the susceptibility increases with hole doping near half filling resembling the experimental results obtained for ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$${\mathrm{Sr}}_{\mathit{x}}$${\mathrm{CuO}}_{4}$. At low electronic density the temperature dependence of the susceptibility is compatible with that of a Fermi liquid.

Physical review. B, Condensed matter•Aug 1, 1993•39 citations
The Hubbard model on a square lattice in the ground state is investigated. Various aspects of Mott transition at half-filling are clarified from the quantum Monte Carlo study. Critical exponents of the transition for charge and spin correlations are estimated. For the doping concentration δ, the charge susceptibility is proportional to δ -1 indicating the divergence of charge mass for δ→0, while the spin susceptibility is always finite and not singular near δ=0. Several theoretical consequences of the above results are discussed. Incommensurate peak amplitude of spin correlation scales as (δ-δ c ) -1 with δ c <0.01 in disagreement with RPA results ln (δ-δ c ). The antiferromagnetic order at the half-filling also shows strong correlation character and is not sensitive to the shape of the fermi surface, while the incommensurate peak observed away from the half-filling is sensitively suppressed by the loss of partial nesting consistently with the weak coupling picture. Comparisons with experimental indications of high- T c oxides are made.
Quantum Monte Carlo study of the Hubbard model on a square lattice in the ground state is reported. For the on-site interaction U =4 scaled by the transfer, the charge excitation gap at the half-filling is estimated to be Δ c =0.58±0.08. Near half-filling, the charge compressibility κ follows the form κ∝δ -1 for the doping concentration δ, indicating a divergent behavior as the system approaches half-filling, while κ=0 at the half-filling. The incommensurate spin correlation length ξ scales as ξ∝1/δ\sqrt{\delta}δ​.
We investigate the metal-insulator transition in the half-filled Hubbard model on the two-dimensional triangular lattice using both the Kotliar-Ruckenstein slave-boson technique and an exact numerical diagonalization of finite clusters. Contrary to the case of a square lattice, where a perfect nesting of the Fermi surface leads to a metal-insulator transition at arbitrarily small values of U, always accompanied by antiferromagnetic ordering, on a triangular lattice, due to the lack of perfect nesting, the transition takes place at a finite value of U, and frustration induces a nontrivial competition among different magnetic phases. Indeed, within the mean-field approximation in the slave-boson approach, as the interaction grows the paramagnetic metal turns into a metallic phase with incommensurate spiral ordering. Increasing the interaction further, a linear spin-density wave is stabilized, and finally for strong coupling the latter phase undergoes a first-order transition toward an antiferromagnetic insulator. No trace of the intermediate phases is seen in the exact diagonalization results, indicating a transition between a paramagnetic metal and an antiferromagnetic insulator.
Employing a combination of functional renormalization group calculations and projective determinantal quantum Monte Carlo simulations, we examine the Hubbard model on the square lattice bilayer at half filling. From this combined analysis, we obtain a comprehensive account of the ground-state phase diagram with respect to the extent of the system's metallic and (antiferromagnetically ordered) Mott-insulating as well as band-insulating regions. By means of an unbiased functional renormalization group approach, we exhibit the antiferromagnetic Mott-insulating state as the relevant instability of the free metallic state, induced by any weak finite on-site repulsion. Upon performing a careful analysis of the quantum Monte Carlo data, we resolve the difficulty of identifying this antiferromagnetic ground state for finite interlayer hopping in the weak-coupling regime, where nonmonotonic finite-size corrections are shown to relate to the two-sheeted Fermi surface structure of the metallic phase. On the other hand, quantum Monte Carlo simulations are well suited to identify the transition between the Mott-insulating phase and the band insulator in the intermediate-to-strong coupling regime. Here, we compare our numerical findings to indications for the transition region obtained from the functional renormalization group procedure.
Using the quantum Monte Carlo simulation technique for the lattice fermion systems, we study the charge susceptibility of the two-dimensional Hubbard near the half-filled band. Lattice sizes up to 6×6 and inverse temperatures up to 6 t ( t : hopping parameter) are studied. We find that the behavior of the charge susceptibility is strongly dependent on the electron concentration at low temperature: the charge susceptibility is suppressed at the half-filled band, while it is considerably enhanced in non-half-filled band cases by the electron correlation. The above behavior of the charge susceptibility is quite similar to the recent Bethe-Ansatz result for the one-dimensional Hubbard model.
Thermodynamic and dynamical properties of filling-control metal-insulator transition (MIT) in the Hubbard model are studied by the operator projection method, especially in two dimensions. This is a non-perturbative analytic approach to many-body systems. The present theory incorporates the Mott-Hubbard, Brinkman-Rice and Slater pictures of the MIT into a unified framework, together with reproducing low-energy narrow band arising from spin-charge fluctuations. At half filling, single-particle spectra A(ω,k)A(ω, k)A(ω,k) show formation of two Hubbard bands and their antiferromagnetic shadows separated by a Mott gap in the plane of energy ωωω and momentum kkk with lowering temperatures. These four bands produce splitting to two low-energy narrow bands and two SDW-like bands in the dispersion. Near half filling, the low-energy narrow band persists at low temperatures. This narrow band has a particularly weak dispersion and large weights around (π,0)(π, 0)(π,0) and (0,π)(0, π)(0,π) momenta. The velocity of these low-energy excitations is shown to vanish towards the MIT, indicating the mass divergence as in the Brinkman-Rice picture, but most prominently around (π,0)(π,0)(π,0) and (0,π)(0,π)(0,π) with strong momentum dependence. This reflects the suppression of the coherence near the MIT. Main structures in A(ω,k)A(ω, k)A(ω,k) show remarkable similarities to quantum Monte-Carlo results in two dimensions as well as in the one-dimensional Hubbard model. The charge compressibility appears to diverge with decreasing doping concentration in both one and two dimensions in agreement with the exact and quantum Monte-Carlo results. We also discuss implications of the flat dispersion formed near the Fermi level to the observations in high-$\Tc$ cuprate superconductors.
Using quantum Monte Carlo techniques we study the magnetic susceptibility of the Hubbard model on a two-dimensional square lattice for values of the coupling U/t ranging from 4 to 20, and as a function of temperature and doping. For U/t=4 the magnetic susceptibility has a maximum at half filling and decreases with doping, in qualitative agreement with recent random-phase-approximation calculations. However, for U/t=10 the susceptibility increases with hole doping near half filling resembling the experimental results obtained for ${\mathrm{La}}_{2\mathrm{\ensuremath{-}}\mathit{x}}$Srx{\mathrm{Sr}}_{\mathit{x}}Srx​CuO4{\mathrm{CuO}}_{4}CuO4​. At low electronic density the temperature dependence of the susceptibility is compatible with that of a Fermi liquid.