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/9911071 · Physical review. B, Condensed matter

Green’s-function Monte Carlo for lattice fermions: Application to the t−J model

C. Stephen Hellberg, Efstratios Manousakis

Published May 1, 2000• 46 citations
arXivHTMLPDF

Abstract

We develop a general numerical method to study the zero-temperature properties of strongly correlated electron models on large lattices. The technique, which resembles Green's-function Monte Carlo, projects the ground-state component from a trial wave function with no approximations. We use this method to determine the phase diagram of the two-dimensional $t\ensuremath{-}J$ model, using the Maxwell construction to investigate electronic phase separation. The shell effects of fermions on finite-sized periodic lattices are minimized by keeping the number of electrons fixed at a closed-shell configuration and varying the size of the lattice. Results obtained for various electron numbers corresponding to different closed shells indicate that the finite-size effects in our calculation are small. For any value of interaction strength, we find that there is always a value of the electron density above which the system can lower its energy by forming a two-component phase separated state. Our results are compared with other calculations on the $t\ensuremath{-}J$ model. We find that the most accurate results are consistent with phase separation at all interaction strengths.

DOI
10.1103/physrevb.61.11787

Related papers

Phase separation in the two-dimensional Hubbard model: A fixed-node quantum Monte Carlo study

Andrea C. Cosentini, Massimo Capone, Leonardo Guidoni +1 more

Fixed-node Green's-function Monte Carlo calculations have been performed for very large 16\ifmmode\times\else\texttimes\fi{}16 two-dimensional Hubbard lattices, large interaction strengths $U=10,20,$ and $40$, and many $(15\ensuremath{\sim}20)$ densities between empty and half-filling. The nodes were fixed by a simple Slater-Gutzwiller trial wave function. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of $U$. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape, and boundary-condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions.

arXiv:cond-mat/9801299•Physical review. B, Condensed matter•Dec 1, 1998•35 citations

Phase Separation at all Interaction Strengths in thet-JModel

C. Stephen Hellberg, Efstratios Manousakis

We investigate the phase diagram of the two-dimensional t-J model using a recently developed Green's function Monte Carlo method for lattice fermions. We use the technique to calculate exact ground-state energies of the model on large lattices. In contrast to many previous studies, we find the model phase separates for all values of $J/t$. In particular, it is unstable at the hole dopings and interaction strengths at which the model was thought to describe the cuprate superconductors.

arXiv:cond-mat/9611195•Physical Review Letters•Jun 16, 1997•137 citations

Phase separation in the two-dimensional Hubbard model: A fixed-node quantum Monte Carlo study

Andrea C. Cosentini, Guidoni L, Giovanni B. Bachelet +1 more

Fixed-node Green's-function Monte Carlo calculations have been performed for very large 16x16 two-dimensional Hubbard lattices, large interaction strengths U=10,20, and 40, and many (15 similar to 20) densities between empty and half-filling. The nodes were fixed by a simple Slater-Gutzwiller trial wave function. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of U. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape, and boundary-condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions. RI Capone, Massimo/A-7762-2008; Guidoni, Leonardo/G-3759-2011

•Jan 1, 1998•0 citations

Fixed-node Monte Carlo study of the two-dimensional Hubbard model

Guozhong An, J. M. J. van Leeuwen

Abstract This study extended the fixed-node QMC method to a lattice model for hopping fermions using Green’s function techniques. In general, the modifications of the procedures for continuous systems to adapt them to a lattice system were minor, and there appeared to be no surprises. Applications were made to the two-dimensional Hubbard model with different on-site interactions, lattice dimensions, boundary conditions, node specifications, and trial wavefunctions for importance sampling. Results were compared with those of Hartree-Fock and variational QMC calculations. For the HF calculations with zero interaction, the exact solutions are known, and the Green’s function QMC results agreed with those within statistical error. Trial wavefunctions specifying the nodes were obtained directly as those of HF determinantal solutions for zero interaction, or as variationally optimized determinantal functions with correlation terms. The procedure for fixed nodes was to eliminate walkers crossing the boundaries of a nodal region, a procedure not quite correct, but probably a reasonable approximation. The results in all cases were consistent with the HF and variational QMC calculations, giving lower energies than either for nonzero interaction energies. The existence of a domain-wall phase suggested by the HF and variational QMC results was confirmed by the fixed-node Green’s function results. It was noted that the fixed-node method would be applicable to systems other than the Hubbard model.

•Jun 18, 2007•0 citations

Finite-size scaling with modified boundary conditions

Sandro Sorella

An efficient scheme is introduced for a fast and smooth convergence to the thermodynamic limit of electronic properties obtained with finite-size calculations on correlated Hamiltonians. This is obtained by modifying the energy levels of the free electron part of the Hamiltonian in a way consistent with the corresponding one-particle density of states in the thermodynamic limit. After this modification, free electron ground state energies, exact in the thermodynamic limit, are obtained with finite-size calculations and for all the particular fillings that satisfy the so called ``closed-shell condition.'' For those fillings the auxiliary field quantum Monte Carlo technique is particularly efficient and, by combining it with the present method, we provide strong numerical evidence that phase separation occurs in the low doping region and moderate $U\ensuremath{\lesssim}4t$ regime of this model.

arXiv:1505.02560•Physical Review B•Jun 29, 2015•26 citations

Constrained path Monte Carlo method for fermion ground states

Shiwei Zhang, Joseph Carlson, James E. Gubernatis

We describe and discuss a recently proposed quantum Monte Carlo algorithm to compute the ground-state properties of various systems of interacting fermions. In this method, the ground state is projected from an initial wave function by a branching random walk in an overcomplete basis of Slater determinants. By constraining the determinants according to a trial wave function |${\mathrm{\ensuremath{\psi}}}_{\mathrm{T}}$〉, we remove the exponential decay of signal-to-noise ratio characteristic of the sign problem. The method is variational and is exact if |${\mathrm{\ensuremath{\psi}}}_{\mathrm{T}}$〉 is exact. We illustrate the method by describing in detail its implementation for the two-dimensional one-band Hubbard model. We show results for lattice sizes up to 16\ifmmode\times\else\texttimes\fi{}16 and for various electron fillings and interaction strengths. With simple single-determinant wave functions as |${\mathrm{\ensuremath{\psi}}}_{\mathrm{T}}$〉, the method yields accurate (often to within a few percent) estimates of the ground-state energy as well as correlation functions, such as those for electron pairing. We conclude by discussing possible extensions of the algorithm.

arXiv:cond-mat/9607062•Physical review. B, Condensed matter•Mar 15, 1997•388 citations
We develop a general numerical method to study the zero-temperature properties of strongly correlated electron models on large lattices. The technique, which resembles Green's-function Monte Carlo, projects the ground-state component from a trial wave function with no approximations. We use this method to determine the phase diagram of the two-dimensional t\ensuremath−Jt\ensuremath{-}Jt\ensuremath−J model, using the Maxwell construction to investigate electronic phase separation. The shell effects of fermions on finite-sized periodic lattices are minimized by keeping the number of electrons fixed at a closed-shell configuration and varying the size of the lattice. Results obtained for various electron numbers corresponding to different closed shells indicate that the finite-size effects in our calculation are small. For any value of interaction strength, we find that there is always a value of the electron density above which the system can lower its energy by forming a two-component phase separated state. Our results are compared with other calculations on the t\ensuremath−Jt\ensuremath{-}Jt\ensuremath−J model. We find that the most accurate results are consistent with phase separation at all interaction strengths.
Fixed-node Green's-function Monte Carlo calculations have been performed for very large 16\ifmmode\times\else\texttimes\fi{}16 two-dimensional Hubbard lattices, large interaction strengths U=10,20,U=10,20,U=10,20, and 404040, and many $(15\ensuremath{\sim}20)$ densities between empty and half-filling. The nodes were fixed by a simple Slater-Gutzwiller trial wave function. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of UUU. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape, and boundary-condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions.
We investigate the phase diagram of the two-dimensional t-J model using a recently developed Green's function Monte Carlo method for lattice fermions. We use the technique to calculate exact ground-state energies of the model on large lattices. In contrast to many previous studies, we find the model phase separates for all values of J/tJ/tJ/t. In particular, it is unstable at the hole dopings and interaction strengths at which the model was thought to describe the cuprate superconductors.
Fixed-node Green's-function Monte Carlo calculations have been performed for very large 16x16 two-dimensional Hubbard lattices, large interaction strengths U=10,20, and 40, and many (15 similar to 20) densities between empty and half-filling. The nodes were fixed by a simple Slater-Gutzwiller trial wave function. For each value of U we obtained a sequence of ground-state energies which is consistent with the possibility of a phase separation close to half-filling, with a hole density in the hole-rich phase which is a decreasing function of U. The energies suffer, however, from a fixed-node bias: more accurate nodes are needed to confirm this picture. Our extensive numerical results and their test against size, shell, shape, and boundary-condition effects also suggest that phase separation is quite a delicate issue, on which simulations based on smaller lattices than considered here are unlikely to give reliable predictions. RI Capone, Massimo/A-7762-2008; Guidoni, Leonardo/G-3759-2011
Abstract This study extended the fixed-node QMC method to a lattice model for hopping fermions using Green’s function techniques. In general, the modifications of the procedures for continuous systems to adapt them to a lattice system were minor, and there appeared to be no surprises. Applications were made to the two-dimensional Hubbard model with different on-site interactions, lattice dimensions, boundary conditions, node specifications, and trial wavefunctions for importance sampling. Results were compared with those of Hartree-Fock and variational QMC calculations. For the HF calculations with zero interaction, the exact solutions are known, and the Green’s function QMC results agreed with those within statistical error. Trial wavefunctions specifying the nodes were obtained directly as those of HF determinantal solutions for zero interaction, or as variationally optimized determinantal functions with correlation terms. The procedure for fixed nodes was to eliminate walkers crossing the boundaries of a nodal region, a procedure not quite correct, but probably a reasonable approximation. The results in all cases were consistent with the HF and variational QMC calculations, giving lower energies than either for nonzero interaction energies. The existence of a domain-wall phase suggested by the HF and variational QMC results was confirmed by the fixed-node Green’s function results. It was noted that the fixed-node method would be applicable to systems other than the Hubbard model.
An efficient scheme is introduced for a fast and smooth convergence to the thermodynamic limit of electronic properties obtained with finite-size calculations on correlated Hamiltonians. This is obtained by modifying the energy levels of the free electron part of the Hamiltonian in a way consistent with the corresponding one-particle density of states in the thermodynamic limit. After this modification, free electron ground state energies, exact in the thermodynamic limit, are obtained with finite-size calculations and for all the particular fillings that satisfy the so called ``closed-shell condition.'' For those fillings the auxiliary field quantum Monte Carlo technique is particularly efficient and, by combining it with the present method, we provide strong numerical evidence that phase separation occurs in the low doping region and moderate $U\ensuremath{\lesssim}4t$ regime of this model.
We describe and discuss a recently proposed quantum Monte Carlo algorithm to compute the ground-state properties of various systems of interacting fermions. In this method, the ground state is projected from an initial wave function by a branching random walk in an overcomplete basis of Slater determinants. By constraining the determinants according to a trial wave function |${\mathrm{\ensuremath{\psi}}}_{\mathrm{T}}$〉, we remove the exponential decay of signal-to-noise ratio characteristic of the sign problem. The method is variational and is exact if |${\mathrm{\ensuremath{\psi}}}_{\mathrm{T}}$〉 is exact. We illustrate the method by describing in detail its implementation for the two-dimensional one-band Hubbard model. We show results for lattice sizes up to 16\ifmmode\times\else\texttimes\fi{}16 and for various electron fillings and interaction strengths. With simple single-determinant wave functions as |${\mathrm{\ensuremath{\psi}}}_{\mathrm{T}}$〉, the method yields accurate (often to within a few percent) estimates of the ground-state energy as well as correlation functions, such as those for electron pairing. We conclude by discussing possible extensions of the algorithm.