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

Simulating the Euclidean time Schroedinger equations using an Intel iPSC/860 hypercube: Application to the t-J model of high-{Tc} superconductivity

Mitchel D. Kovarik, Ted Barnes

Published Oct 1, 1993• 0 citations
OpenAlex

Abstract

We describe a Monte Carlo simulation of a dynamical fermion problem in two spatial dimensions on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the ``minus sign problem`` of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. Generalization to physical hopping parameter values wig only require use of an improved trial wavefunction for importance sampling.

Related papers

A Monte Carlo study of the minus sign problem in the t-J model using an Intel iPSC/860 hypercube

Mitchel D. Kovarik, Ted Barnes

We describe a Monte Carlo simulation of the 2-dimensional t-J model on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the ``minus sign problem`` of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. We demonstrate that a qualitative picture at higher hopping parameters may be found by extrapolating weak hopping results where the minus sign problem is less severe. Generalization to physical hopping parameter values will only require use of an improved trial wavefunction for importance sampling.

•Jan 1, 1993•0 citations

Simulating the Euclidean Time — Schrödinger Equation Using an Intel iPSC/860 Hypercube: Application to the t — J Model of High-T c Superconductivity

Mitchel D. Kovarik, Ted Barnes

Springer proceedings in physics•Jan 1, 1993•0 citations

MONTE CARLO DIAGONALIZATION OF VERY LARGE MATRICES: APPLICATION TO FERMION SYSTEMS

Hans A. De Raedt, Wolfgang von der Linden

All known Quantum-Monte-Carlo algorithms for fermions suffer from the so-called “minus-sign-problem” which is detrimental to the application of these simulation methods to fermion systems at very low temperatures and/or of very many lattice sites. We identify the origin of this fundamental problem, demonstrate that it is a very general feature, not necessarily related to the presence of fermionic degrees of freedom. We describe an novel algorithm which does not suffer from the minus-sign problem. Illustrative results for the two-dimensional Hubbard model are presented

International Journal of Modern Physics C•Feb 1, 1992•6 citations

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

Perturbative solution of fermionic sign problem in quantum Monte Carlo computations

Sergei Iskakov, Mikhail I. Katsnelson, Alexander I. Lichtenstein

Abstract We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1 t , which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter $${t}^{{\prime} }$$ t ′ for high-temperature superconducting cuprates.

arXiv:2303.01607•npj Computational Materials•Feb 15, 2024•6 citations

Flat-histogram methods in quantum Monte Carlo simulations: Application to the t-J model

Nikolaos G. Diamantis, Efstratios Manousakis

We discuss that flat-histogram techniques can be appropriately applied in the sampling of quantum Monte Carlo simulation in order to improve the statistical quality of the results at long imaginary time or low excitation energy. Typical imaginary-time correlation functions calculated in quantum Monte Carlo are subject to exponentially growing errors as the range of imaginary time grows and this smears the information on the low energy excitations. We show that we can extract the low energy physics by modifying the Monte Carlo sampling technique to one in which configurations which contribute to making the histogram of certain quantities flat are promoted. We apply the diagrammatic Monte Carlo (diag-MC) method to the motion of a single hole in the t-J model and we show that the implementation of flat-histogram techniques allows us to calculate the Green's function in a wide range of imaginary-time. In addition, we show that applying the flat-histogram technique alleviates the "sign"-problem associated with the simulation of the single-hole Green's function at long imaginary time.

Journal of Physics Conference Series•Sep 1, 2016•0 citations
We describe a Monte Carlo simulation of a dynamical fermion problem in two spatial dimensions on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the minus sign problem of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. Generalization to physical hopping parameter values wig only require use of an improved trial wavefunction for importance sampling.
We describe a Monte Carlo simulation of the 2-dimensional t-J model on an Intel iPSC/860 hypercube. The problem studied is the determination of the dispersion relation of a dynamical hole in the t-J model of the high temperature superconductors. Since this problem involves the motion of many fermions in more than one spatial dimensions, it is representative of the class of systems that suffer from the ``minus sign problem`` of dynamical fermions which has made Monte Carlo simulation very difficult. We demonstrate that for small values of the hole hopping parameter one can extract the entire hole dispersion relation using the GRW Monte Carlo algorithm, which is a simulation of the Euclidean time Schroedinger equation, and present results on 4 {times} 4 and 6 {times} 6 lattices. We demonstrate that a qualitative picture at higher hopping parameters may be found by extrapolating weak hopping results where the minus sign problem is less severe. Generalization to physical hopping parameter values will only require use of an improved trial wavefunction for importance sampling.
All known Quantum-Monte-Carlo algorithms for fermions suffer from the so-called “minus-sign-problem” which is detrimental to the application of these simulation methods to fermion systems at very low temperatures and/or of very many lattice sites. We identify the origin of this fundamental problem, demonstrate that it is a very general feature, not necessarily related to the presence of fermionic degrees of freedom. We describe an novel algorithm which does not suffer from the minus-sign problem. Illustrative results for the two-dimensional Hubbard model are presented
Abstract We have developed a strong-coupling perturbation scheme for a general doped Hubbard model around a particle-hole-symmetric reference system, which is free from the fermionic sign problem. Our approach is based on the lattice determinantal Quantum Monte Carlo (QMC) method in both continuous and discrete time versions for large periodic clusters in a fermionic bath. By considering the first-order perturbation in the shift of the chemical potential and the second-neighbor hopping, we are able to obtain an accurate electronic spectral function for a range of parameters that correspond to optimally doped cuprate systems at temperatures of up to T = 0.1 t , which are challenging to access using straightforward lattice QMC calculations. We also discuss the formation of the pseudogap and the nodal-antinodal dichotomy for a doped Hubbard system with the interaction parameter U equal to the bandwidth and the optimal value of the next-nearest-neighbor hopping parameter t′{t}^{{\prime} }t′ t ′ for high-temperature superconducting cuprates.
We discuss that flat-histogram techniques can be appropriately applied in the sampling of quantum Monte Carlo simulation in order to improve the statistical quality of the results at long imaginary time or low excitation energy. Typical imaginary-time correlation functions calculated in quantum Monte Carlo are subject to exponentially growing errors as the range of imaginary time grows and this smears the information on the low energy excitations. We show that we can extract the low energy physics by modifying the Monte Carlo sampling technique to one in which configurations which contribute to making the histogram of certain quantities flat are promoted. We apply the diagrammatic Monte Carlo (diag-MC) method to the motion of a single hole in the t-J model and we show that the implementation of flat-histogram techniques allows us to calculate the Green's function in a wide range of imaginary-time. In addition, we show that applying the flat-histogram technique alleviates the "sign"-problem associated with the simulation of the single-hole Green's function at long imaginary time.