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:1305.6798 · Physical Review B

Equation of state of the fermionic two-dimensional Hubbard model

J. P. F. LeBlanc, Emanuel Gull

Published Oct 7, 2013• 55 citations
arXivHTMLPDF

Abstract

We present results for the equation of state of the two-dimensional Hubbard model on an isotropic square lattice as obtained from a controlled and numerically exact large-cluster dynamical mean field simulation. Our results are obtained for large but finite systems and are extrapolated to infinite system size using a known finite-size scaling relation, and are supplemented by reliable error bars accounting for all sources of errors. We establish the importance of examining the decay of spatial spin correlations to determine whether a sufficiently large cluster has been used and with this in mind we present the energy, entropy, double occupancy, and nearest-neighbor spin correlations extrapolated to the thermodynamic limit. We discuss the implications of these calculations on pseudogap physics of the 2D Hubbard model away from half filling, where we find a strong behavioral shift in energy below a temperature ${T}^{*}$ which becomes more pronounced for larger clusters. Finally, we provide reference calculations and tables for the equation of state for values of doping away from half filling which are of interest to cold-atom experiments.

DOI
10.1103/physrevb.88.155108

Related papers

Two-dimensional Hubbard model: Numerical simulation study

J. E. Hirsch

We have studied the two-dimensional Hubbard model on a square lattice with nearest-neighbor hopping. We first discuss the properties of the model within the mean-field approximation: Because of the form of the band structure, some peculiar features are found. We then discuss the simulation algorithm used and compare simulation results with exact results for 6-site chains to test the reliability of the approach. We present results for thermodynamic properties and correlation functions for lattices up to 8\ifmmode\times\else\texttimes\fi{}8 in spatial size. The system is found to be an antiferromagnetic insulator for all values of the coupling constant at zero temperature in the half-filled-band case, but the long-range order is much smaller than predicted by mean-field theory. We perform a finite-size-scaling analysis to determine the character of the transition at zero coupling. For non-half-filled-band cases, our results suggest that the system is always paramagnetic, in contradiction with Hartree-Fock predictions. The system does not show tendency to ferromagnetism nor triplet superconductivity in the parameter range studied. We also discuss some properties of the attractive Hubbard model in the half-filled-band case.

Physical review. B, Condensed matter•Apr 1, 1985•1,082 citations

Thermodynamics of the quantum critical point at finite doping in the two-dimensional Hubbard model studied via the dynamical cluster approximation

Karlis Mikelsons, Ehsan Khatami, Dimitrios Galanakis +3 more

We study the thermodynamics of the two-dimensional Hubbard model within the dynamical cluster approximation. We use continuous time quantum Monte Carlo as a cluster solver to avoid the systematic error which complicates the calculation of the entropy and potential energy (double occupancy). We find that at a critical filling, there is a pronounced peak in the entropy divided by temperature, S/T, and in the normalized double occupancy as a function of doping. At this filling, we find that specific heat divided by temperature, C/T, increases strongly with decreasing temperature and kinetic and potential energies vary like T2 lnT. These are all characteristics of quantum critical behavior. © 2009 The American Physical Society.

Civil War Book Review•Mar 17, 2010•0 citations

Thermodynamics of the quantum critical point at finite doping in the two-dimensional Hubbard model studied via the dynamical cluster approximation

Karlis Mikelsons, Ehsan Khatami, D. Galanakis +3 more

We study the thermodynamics of the two-dimensional Hubbard model within the dynamical cluster approximation. We use continuous time quantum Monte Carlo as a cluster solver to avoid the systematic error which complicates the calculation of the entropy and potential energy (double occupancy). We find that at a critical filling, there is a pronounced peak in the entropy divided by temperature, $S/T$, and in the normalized double occupancy as a function of doping. At this filling, we find that specific heat divided by temperature, $C/T$, increases strongly with decreasing temperature and kinetic and potential energies vary like ${T}^{2}\text{ }\text{ln}\text{ }T$. These are all characteristics of quantum critical behavior.

arXiv:0909.0498•Physical Review B•Oct 8, 2009•39 citations

Benchmark study of the two-dimensional Hubbard model with auxiliary-field quantum Monte Carlo method

Mingpu Qin, Hao Shi, Shiwei Zhang

Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of $U/t$ from 2 to 8, for a range of densities including half-filling and $n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75$, and $0.875$. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.

arXiv:1605.09421•Physical review. B./Physical review. B•Aug 4, 2016•100 citations

Pseudogaps in the 2D Hubbard Model

Carey Huscroft, Mark Jarrell, Thomas Maier +2 more

We study the pseudogaps in the spectra of the 2D Hubbard model using both finite-size and dynamical cluster approximation (DCA) quantum Monte Carlo calculations. At half-filling, a charge pseudogap, accompanied by non-Fermi-liquid behavior in the self-energy, is shown to persist in the thermodynamic limit. The DCA (finite-size) method systematically underestimates (overestimates) the width of the pseudogap. A spin pseudogap is not seen at half-filling. At finite doping, a divergent d-wave pair susceptibility is observed.

Physical Review Letters•Jan 1, 2001•132 citations

The Mott metal-insulator transition in half-filled two-dimensional Hubbard models

P Sahebsara

We study the Mott transition in the two dimensional Hubbard model by using the variational cluster approximation. The transition potential obtained is roughly Uc ≈ 2 and 6 for square and triangular lattices, respectively. A comparison between results of this approximation and other quantum cluster methods is presented. Our zero-temperature calculation at strong coupling show that the transition on the triangular and square lattices occur at lower values of compared with other numerical techniques such as DMFT, CDMFT, and DCA. We also study the thermodynamic limit by an extrapolation to infinite size.

•Jun 1, 2008•2 citations
We present results for the equation of state of the two-dimensional Hubbard model on an isotropic square lattice as obtained from a controlled and numerically exact large-cluster dynamical mean field simulation. Our results are obtained for large but finite systems and are extrapolated to infinite system size using a known finite-size scaling relation, and are supplemented by reliable error bars accounting for all sources of errors. We establish the importance of examining the decay of spatial spin correlations to determine whether a sufficiently large cluster has been used and with this in mind we present the energy, entropy, double occupancy, and nearest-neighbor spin correlations extrapolated to the thermodynamic limit. We discuss the implications of these calculations on pseudogap physics of the 2D Hubbard model away from half filling, where we find a strong behavioral shift in energy below a temperature T∗{T}^{*}T∗ which becomes more pronounced for larger clusters. Finally, we provide reference calculations and tables for the equation of state for values of doping away from half filling which are of interest to cold-atom experiments.
We have studied the two-dimensional Hubbard model on a square lattice with nearest-neighbor hopping. We first discuss the properties of the model within the mean-field approximation: Because of the form of the band structure, some peculiar features are found. We then discuss the simulation algorithm used and compare simulation results with exact results for 6-site chains to test the reliability of the approach. We present results for thermodynamic properties and correlation functions for lattices up to 8\ifmmode\times\else\texttimes\fi{}8 in spatial size. The system is found to be an antiferromagnetic insulator for all values of the coupling constant at zero temperature in the half-filled-band case, but the long-range order is much smaller than predicted by mean-field theory. We perform a finite-size-scaling analysis to determine the character of the transition at zero coupling. For non-half-filled-band cases, our results suggest that the system is always paramagnetic, in contradiction with Hartree-Fock predictions. The system does not show tendency to ferromagnetism nor triplet superconductivity in the parameter range studied. We also discuss some properties of the attractive Hubbard model in the half-filled-band case.
We study the thermodynamics of the two-dimensional Hubbard model within the dynamical cluster approximation. We use continuous time quantum Monte Carlo as a cluster solver to avoid the systematic error which complicates the calculation of the entropy and potential energy (double occupancy). We find that at a critical filling, there is a pronounced peak in the entropy divided by temperature, S/T, and in the normalized double occupancy as a function of doping. At this filling, we find that specific heat divided by temperature, C/T, increases strongly with decreasing temperature and kinetic and potential energies vary like T2 lnT. These are all characteristics of quantum critical behavior. © 2009 The American Physical Society.
We study the thermodynamics of the two-dimensional Hubbard model within the dynamical cluster approximation. We use continuous time quantum Monte Carlo as a cluster solver to avoid the systematic error which complicates the calculation of the entropy and potential energy (double occupancy). We find that at a critical filling, there is a pronounced peak in the entropy divided by temperature, S/TS/TS/T, and in the normalized double occupancy as a function of doping. At this filling, we find that specific heat divided by temperature, C/TC/TC/T, increases strongly with decreasing temperature and kinetic and potential energies vary like T2 ln T{T}^{2}\text{ }\text{ln}\text{ }TT2 ln T. These are all characteristics of quantum critical behavior.
Ground-state properties of the Hubbard model on a two-dimensional square lattice are studied by the auxiliary-field quantum Monte Carlo method. Accurate results for energy, double occupancy, effective hopping, magnetization, and momentum distribution are calculated for interaction strengths of U/tU/tU/t from 2 to 8, for a range of densities including half-filling and n=0.3,0.5,0.6,0.75n=0.3,0.5,0.6,\phantom{\rule{0.16em}{0ex}}0.75n=0.3,0.5,0.6,0.75, and 0.8750.8750.875. At half-filling, the results are numerically exact. Away from half-filling, the constrained path Monte Carlo method is employed to control the sign problem. Our results are obtained with several advances in the computational algorithm, which are described in detail. We discuss the advantages of generalized Hartree-Fock trial wave functions and its connection to pairing wave functions, as well as the interplay with different forms of Hubbard-Stratonovich decompositions. We study the use of different twist angle sets when applying the twist averaged boundary conditions. We propose the use of quasirandom sequences, which improves the convergence to the thermodynamic limit over pseudorandom and other sequences. With it and a careful finite size scaling analysis, we are able to obtain accurate values of ground-state properties in the thermodynamic limit. Detailed results for finite-sized systems up to $16\ifmmode\times\else\texttimes\fi{}16$ are also provided for benchmark purposes.
We study the pseudogaps in the spectra of the 2D Hubbard model using both finite-size and dynamical cluster approximation (DCA) quantum Monte Carlo calculations. At half-filling, a charge pseudogap, accompanied by non-Fermi-liquid behavior in the self-energy, is shown to persist in the thermodynamic limit. The DCA (finite-size) method systematically underestimates (overestimates) the width of the pseudogap. A spin pseudogap is not seen at half-filling. At finite doping, a divergent d-wave pair susceptibility is observed.
We study the Mott transition in the two dimensional Hubbard model by using the variational cluster approximation. The transition potential obtained is roughly Uc ≈ 2 and 6 for square and triangular lattices, respectively. A comparison between results of this approximation and other quantum cluster methods is presented. Our zero-temperature calculation at strong coupling show that the transition on the triangular and square lattices occur at lower values of compared with other numerical techniques such as DMFT, CDMFT, and DCA. We also study the thermodynamic limit by an extrapolation to infinite size.