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/9705026 · International Journal of Modern Physics C

Comparison of Calculations for the Hubbard Model Obtained with Quantum-Monte-Carlo, Exact, and Stochastic Diagonalization

Thomas Husslein, Werner Fettes, I. Morgenstern

Published Apr 1, 1997• 31 citations
arXivHTMLPDF

Abstract

In this paper we compare numerical results for the ground state of the Hubbard model obtained by Quantum-Monte-Carlo simulations with results from exact and stochastic diagonalizations. We find good agreement for the ground state energy and superconducting correlations for both, the repulsive and attractive Hubbard model. Special emphasis lies on the superconducting correlations in the repulsive Hubbard model, where the small magnitude of the values obtained by Monte-Carlo simulations gives rise to the question, whether these results might be caused by fluctuations or systematic errors of the method. Although we notice that the Quantum-Monte-Carlo method has convergence problems for large interactions, coinciding with a minus sign problem, we confirm the results of the diagonalization techniques for small and moderate interaction strengths. Additionally we investigate the numerical stability and the convergence of the Quantum-Monte-Carlo method in the attractive case, to study the influence of the minus sign problem on convergence. Also here in the absence of a minus sign problem we encounter convergence problems for strong interactions.

DOI
10.1142/s0129183197000333

Related papers

Monte Carlo diagonalization of many-body problems: Application to fermion systems

Hans A. De Raedt, Wolfgang von der Linden

Existing quantum Monte Carlo algorithms suffer the so-called minus-sign problem. We propose a stochastic algorithm to compute ground-state properties of quantum lattice models that is free of minus-sign problems. Illustrative results for the two-dimensional Hubbard model are presented.

Physical review. B, Condensed matter•Apr 15, 1992•31 citations

Some Findings on the Minus-Sign Problem for the Hubbard Model in the Projector Quantum Monte Carlo Method

Hans‐Georg Matuttis, Nobuyasu Ito

The quantum Monte Carlo method for the Hubbard model using the Hubbard-Stratonovich transformation is reexamined. It is observed that the value of the ground state energy for a 3×2 and 4×4 system is estimated correctly by ignoring the minus-sign of the negative sampling weight appearing in the simulation using the projector method. It is also observed that the average of sign is positive and the correlation between the energy and the sign is weak for these systems. The reasons why the ground state energy is estimated correctly is given on a phenomenological base.

Journal of the Physical Society of Japan•Jun 15, 2001•2 citations

Quantum Monte Carlo and exact diagonalization study of a dynamic Hubbard model

J. E. Hirsch

A one-dimensional model of electrons locally coupled to spin-1/2 degrees of freedom is studied by numerical techniques. The model is one in the class of dynamic Hubbard models that describe the relaxation of an atomic orbital upon double-electron occupancy due to electron-electron interactions. We study the parameter regime where pairing occurs in this model by exact diagonalization of small clusters. World-line quantum Monte Carlo simulations support the results of exact diagonalization for larger systems and show that the kinetic energy is lowered when pairing occurs. The qualitative physics of this model and others in its class, obtained through approximate analytic calculations, is that superconductivity occurs through hole undressing even in parameter regimes where the effective on-site interaction is strongly repulsive. Our numerical results confirm the expected qualitative behavior and show that pairing will occur in a substantially larger parameter regime than predicted by the approximate low-energy effective Hamiltonian.

arXiv:cond-mat/0201005•Physical review. B, Condensed matter•May 22, 2002•27 citations

Finite Size Scaling in the tt'-Hubbard Model and in BCS-Reduced Hubbard Models

Werner Fettes, I. Morgenstern

With the projector quantum Monte Carlo algorithm and the stochastic diagonalization it is possible to calculate the ground state of the Hubbard model for small finite clusters. Nevertheless the usual finite size scaling of the Hubbard model has problems of deducing the behavior of the infinite system correctly from the numerical data of small system sizes. Therefore we study the finite size scaling of the superconducting correlation functions in superconducting BCS-reduced Hubbard models to analyze the finite size behavior in small finite clusters. The ground state of the BCS-reduced Hubbard models is calculated with the stochastic diagonalization without any approximations. As result of these analyses we propose a new finite size scaling ansatz for the Hubbard model, which is able describe the finite size effects in a consistent way taking the corrections to scaling into account, which are dominant for weak interaction strength and small clusters. With this new finite size scaling ansatz it is possible to give evidence for superconductivity for all interaction strengths for both the attractive tt'-Hubbard model (with s-wave symmetry) and the repulsive tt'-Hubbard model (with dx2-y2-wave symmetry).

International Journal of Modern Physics C•Oct 1, 1998•4 citations

Numerical Studies on the Hubbard Model and thet-JModel in One- and Two-Dimensions

Masatoshi Imada, Yasuhiro Hatsugai

Numerical results obtained from Monte Carlo simulation in the ground state and at finite temperatures as well as the exact diagonalization and the transfer matrix method are reported. Efficiency of the Monte Carlo method in the ground state is examined. Spin, charge and superconducting correlations are investigated for the Hubbard and the t - J model in one and two dimensions. The momentum distribution in the one-dimensional Hubbard model shows fermi-liquid-like behavior at least in two orders of magnitude smaller energy scale than the band width. The short-range incommensurate spin correlation is observed and analyzed in the doped Hubbard and the t - J model both in one and two dimensions. In the ground state, the superconducting correlation shows the absence of system size dependence in two dimensions at the filling smaller than 0.8 and in one dimension at any filling. The binding energy of two fermions and the spin and charge distortion around an itinerant fermion are also discussed in the t - J and the coupled spin-fermion model.

Journal of the Physical Society of Japan•Oct 15, 1989•227 citations

Quantum Monte Carlo studies of strongly correlated systems for quantum simulators

Stephan Humeniuk

Many strongly correlated quantum systems are difficult to study numerically because of the exponential growth of the Hilbert space combined with the failure of numerical methods in important parameter regimes. This includes the long-standing problems of high-temperature superconductivity and frustrated magnetism. A promising solution to this dilemma are quantum simulators, well-controlled (synthetic) quantum systems, emulating model Hamiltonians not only of condensed matter systems, where parameters can be tuned at will and model complexity can be added step by step. This dissertation presents numerical simulation studies that are motivated by quantum simulators. Using large-scale quantum Monte Carlo computer simulations and mean-field theory the semiquantitative ground state phase diagram of a Penning ion trap quantum simulator, realizing a two-dimensional model of quantum magnetism with long-range interactions, is obtained. We study also thermal and quantum phase transitions in one-dimensional transverse-field Ising models with long-range interactions, which exhibit a Kosterlitz-Thouless transition and universality classes with continuously varying critical exponents. The second part of this thesis is concerned with the Fermi-Hubbard model, which has been realized recently with ultracold atoms in optical lattices, both for repulsive and attractive interactions. Equipped with a quantum gas microscope, which allows single-site single-atom detection, these experimental setups can collect histograms of the particle number from repeated projective measurements. This dissertation presents a method to compute such full counting statistics for quadratic operators such as particle number or magnetization in the framework of determinantal quantum Monte Carlo. From an analysis of the full counting statistics of the particle number for different subsystem sizes the size of preformed pairs or Cooper pairs in the attractive Hubbard model is inferred. Furthermore, excellent agreement of the numerically computed probability distribution of the staggered magnetization with a recent quantum gas experiment demonstrates that current experiments are capable of resolving the differences between the Hubbard model and its limiting Heisenberg model. An extensive method section reviews the stochastic series expansion quantum Monte Carlo method for Ising models in a transverse field, which is free of the sign problem for arbitrary interactions, with special emphasis on the efficient simulation of long-range interactions. Also a self-contained introduction to the determinantal quantum Monte Carlo method is provided.

•Jan 1, 2018•5 citations
In this paper we compare numerical results for the ground state of the Hubbard model obtained by Quantum-Monte-Carlo simulations with results from exact and stochastic diagonalizations. We find good agreement for the ground state energy and superconducting correlations for both, the repulsive and attractive Hubbard model. Special emphasis lies on the superconducting correlations in the repulsive Hubbard model, where the small magnitude of the values obtained by Monte-Carlo simulations gives rise to the question, whether these results might be caused by fluctuations or systematic errors of the method. Although we notice that the Quantum-Monte-Carlo method has convergence problems for large interactions, coinciding with a minus sign problem, we confirm the results of the diagonalization techniques for small and moderate interaction strengths. Additionally we investigate the numerical stability and the convergence of the Quantum-Monte-Carlo method in the attractive case, to study the influence of the minus sign problem on convergence. Also here in the absence of a minus sign problem we encounter convergence problems for strong interactions.
The quantum Monte Carlo method for the Hubbard model using the Hubbard-Stratonovich transformation is reexamined. It is observed that the value of the ground state energy for a 3×2 and 4×4 system is estimated correctly by ignoring the minus-sign of the negative sampling weight appearing in the simulation using the projector method. It is also observed that the average of sign is positive and the correlation between the energy and the sign is weak for these systems. The reasons why the ground state energy is estimated correctly is given on a phenomenological base.
A one-dimensional model of electrons locally coupled to spin-1/2 degrees of freedom is studied by numerical techniques. The model is one in the class of dynamic Hubbard models that describe the relaxation of an atomic orbital upon double-electron occupancy due to electron-electron interactions. We study the parameter regime where pairing occurs in this model by exact diagonalization of small clusters. World-line quantum Monte Carlo simulations support the results of exact diagonalization for larger systems and show that the kinetic energy is lowered when pairing occurs. The qualitative physics of this model and others in its class, obtained through approximate analytic calculations, is that superconductivity occurs through hole undressing even in parameter regimes where the effective on-site interaction is strongly repulsive. Our numerical results confirm the expected qualitative behavior and show that pairing will occur in a substantially larger parameter regime than predicted by the approximate low-energy effective Hamiltonian.
With the projector quantum Monte Carlo algorithm and the stochastic diagonalization it is possible to calculate the ground state of the Hubbard model for small finite clusters. Nevertheless the usual finite size scaling of the Hubbard model has problems of deducing the behavior of the infinite system correctly from the numerical data of small system sizes. Therefore we study the finite size scaling of the superconducting correlation functions in superconducting BCS-reduced Hubbard models to analyze the finite size behavior in small finite clusters. The ground state of the BCS-reduced Hubbard models is calculated with the stochastic diagonalization without any approximations. As result of these analyses we propose a new finite size scaling ansatz for the Hubbard model, which is able describe the finite size effects in a consistent way taking the corrections to scaling into account, which are dominant for weak interaction strength and small clusters. With this new finite size scaling ansatz it is possible to give evidence for superconductivity for all interaction strengths for both the attractive tt'-Hubbard model (with s-wave symmetry) and the repulsive tt'-Hubbard model (with dx2-y2-wave symmetry).
Numerical results obtained from Monte Carlo simulation in the ground state and at finite temperatures as well as the exact diagonalization and the transfer matrix method are reported. Efficiency of the Monte Carlo method in the ground state is examined. Spin, charge and superconducting correlations are investigated for the Hubbard and the t - J model in one and two dimensions. The momentum distribution in the one-dimensional Hubbard model shows fermi-liquid-like behavior at least in two orders of magnitude smaller energy scale than the band width. The short-range incommensurate spin correlation is observed and analyzed in the doped Hubbard and the t - J model both in one and two dimensions. In the ground state, the superconducting correlation shows the absence of system size dependence in two dimensions at the filling smaller than 0.8 and in one dimension at any filling. The binding energy of two fermions and the spin and charge distortion around an itinerant fermion are also discussed in the t - J and the coupled spin-fermion model.
Many strongly correlated quantum systems are difficult to study numerically because of the exponential growth of the Hilbert space combined with the failure of numerical methods in important parameter regimes. This includes the long-standing problems of high-temperature superconductivity and frustrated magnetism. A promising solution to this dilemma are quantum simulators, well-controlled (synthetic) quantum systems, emulating model Hamiltonians not only of condensed matter systems, where parameters can be tuned at will and model complexity can be added step by step. This dissertation presents numerical simulation studies that are motivated by quantum simulators. Using large-scale quantum Monte Carlo computer simulations and mean-field theory the semiquantitative ground state phase diagram of a Penning ion trap quantum simulator, realizing a two-dimensional model of quantum magnetism with long-range interactions, is obtained. We study also thermal and quantum phase transitions in one-dimensional transverse-field Ising models with long-range interactions, which exhibit a Kosterlitz-Thouless transition and universality classes with continuously varying critical exponents. The second part of this thesis is concerned with the Fermi-Hubbard model, which has been realized recently with ultracold atoms in optical lattices, both for repulsive and attractive interactions. Equipped with a quantum gas microscope, which allows single-site single-atom detection, these experimental setups can collect histograms of the particle number from repeated projective measurements. This dissertation presents a method to compute such full counting statistics for quadratic operators such as particle number or magnetization in the framework of determinantal quantum Monte Carlo. From an analysis of the full counting statistics of the particle number for different subsystem sizes the size of preformed pairs or Cooper pairs in the attractive Hubbard model is inferred. Furthermore, excellent agreement of the numerically computed probability distribution of the staggered magnetization with a recent quantum gas experiment demonstrates that current experiments are capable of resolving the differences between the Hubbard model and its limiting Heisenberg model. An extensive method section reviews the stochastic series expansion quantum Monte Carlo method for Ising models in a transverse field, which is free of the sign problem for arbitrary interactions, with special emphasis on the efficient simulation of long-range interactions. Also a self-contained introduction to the determinantal quantum Monte Carlo method is provided.