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

Physical review. B, Condensed matter

Monte Carlo simulations of one-dimensional fermion systems

J. E. Hirsch, Robert L. Sugar, Douglas J. Scalapino, R. Blankenbecler

Published Nov 1, 1982• 487 citations
Publisher

Abstract

We discuss a new method to perform numerical simulations of one-dimensional systems with fermion and boson degrees of freedom. The method is based on a direct-space, imaginary-time representation of the fermion field. It is fast so that systems having up to 100 sites can easily be simulated. In addition, the method provides an intuitive physical "picture" of the ground state of a one-dimensional many-body system. We discuss in detail how to implement the method and how to compute various physical quantities. In particular, we show how to extend the method to study averages of off-diagonal quantities in an occupation-number representation. To assess the accuracy of our procedure, we apply it to free fermions in one dimension and compare with exact results. We then study a model of spinless interacting fermions and obtain the expected phase structure and behavior of correlation functions. We also consider the extended Hubbard model at various points in its phase diagram and study the behavior of spin-density, charge-density, and pairing correlation functions. We then study the Gross-Neveu model and show how the behavior depends on the number of fermion flavors. Finally, we consider an electron-phonon model and study its behavior both in the one-particle polaron sector and in the half-filled-band case. Along the way we show pictures of the ground-state configurations that give physical insight into the properties of the systems, like charge-density-wave, spin-density-wave, and superconducting states, "fractional charges," and solitons. We conclude by comparing our method with other methods and discuss the possibility of extending it to higher dimensions.

DOI
10.1103/physrevb.26.5033

Related papers

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 study of the one-dimensional Holstein model of spinless fermions

Ross H. McKenzie, Chris J. Hamer, D. W. Murray

The Holstein model of spinless fermions interacting with dispersionless phonons in one dimension is studied by a Green's function Monte Carlo technique. The ground-state energy, first fermionic excited state, density wave correlations, and mean lattice displacements are calculated for lattices of up to 16 sites, for one fermion per two sites, i.e., a half-filled band. Results are obtained for values of the fermion hopping parameter of t=0.1\ensuremath{\omega}, \ensuremath{\omega}, and 10\ensuremath{\omega}, where \ensuremath{\omega} is the phonon frequency. At a finite fermion-phonon coupling g there is a transition from a metallic phase to an insulating phase in which there is charge-density-wave order. Finite-size scaling is found to hold in the metallic phase and is used to extract the coupling dependence of the Luttinger liquid parameters, ${\mathit{u}}_{\mathrm{\ensuremath{\rho}}}$ and ${\mathit{K}}_{\mathrm{\ensuremath{\rho}}}$, the velocity of charge excitations, and the correlation exponent, respectively. For free fermions (g=0) and for strong coupling (${\mathit{g}}^{2}$\ensuremath{\gg}t\ensuremath{\omega}) our results agree well with known analytic results. \textcopyright{} 1996 The American Physical Society.

arXiv:cond-mat/9512046•Physical review. B, Condensed matter•Apr 15, 1996•62 citations

Computational Studies of Correlated Electrons in One-Dimension

Pinaki Sengupta

"We have studied the extended Hubbard and the Peierls-extended Hubbard models in one dimension using primarily computational methods. We discuss the derivation of the models and the computational methods used to study them, focusing primarily on our extensions over existing numerical methods. In particular, we present the construction of ""operator-loop"" updates for spin-1/2 fermions and the inclusion of electron-phonon interactions in the Stochastic Series Expansion Quantum Monte Carlo method and the use of Quantum Parallel Tempering to improve the quality of the numerical data. For the extended Hubbard model at half-filling, we have confirmed the existence of a recently reported bond-order wave phase for the ground state in a small region of parameter space between the well-known charge- and spin-density wave phases. We have studied in detail the nature of the associated quantum phase transitions between all the phases. Upon doping away from half-filling, the doubly-degenerate, spin and charge gapped long-range-ordered ground states of the model are shown to evolve to a Luther-Emery liquid with no charge gap but a finite spin gap. For the Peierls-extended Hubbard model, which includes both electron-phonon and electron-electron interactions, we have studied the dimerization transition of the half-filled system in the presence of fully quantum mechanical phonons. We show that the transition to a dimerized state occurs only when the electron-phonon coupling exceeds a finite critical value that depends on the bare phonon frequency. This is in contrast to previous studies in similar models with adiabatic phonons that had predicted a dimerized state for arbitrarily small values of the electron-phonon coupling."

•Jan 1, 2001•0 citations

Quantum Monte Carlo Simulations of Hubbard and Anderson Models.

Zhang, Yi

We present, in this dissertation, a numerical simulation method to study interacting fermion systems. The general simulation procedures are discussed in connection with a description of a Quantum Monte Carlo simulation algorithm for interacting electrons on lattices. The algorithm presented here has been used to simulate interacting electrons on lattices, and it makes possible the study of substantially larger systems than can be studied by other numerical methods. As long as certain limits of applicability are respected, model Hamiltonians of interacting electrons can be studied without resort to uncontrolled approximations. The method then provides nearly exact solutions to model Hamiltonians of many-body systems, in the sense that the degree of error can be controlled. We also discuss some results obtained from simulations of the extended Hubbard model and and the periodic Anderson model. For the extended Hubbard model in two dimensions, different regions are identified where electron correlation produces antiferromagnetism, charge-density-waves, and singlet pairing superconducting behaviors. We also find regions where transitions from one type of correlation to another occur. For the symmetric periodic Anderson model, we observe the formation of local spin moments at high temperatures and their quenching at low temperatures, as in the single impurity Anderson model. In addition, we find antiferromagnetic interaction between f local moments at low temperatures which are not present in a single impurity problem.

Civil War Book Review•Jan 1, 1990•0 citations

Quantum Monte Carlo simulations of the one-dimensional extended Hubbard model

W. R. Somsky, James E. Gubernatis

We report preliminary results of an investigation of the thermodynamic properties of the extended Hubbard model in one- dimension, calculated with the world-line Monte Carlo method described by Hirsch et al. With strictly continuous world-lines, we are able to measure the expectation of operators that conserve fermion number locally, such as the energy and (spatial) occupation number. By permitting the world-lines to be broken'' stochastically, we may also measure the expectation of operators that conserve fermion number only globally, such as the single-particle Green's function. For a 32 site lattice we present preliminary calculations of the average electron occupancy as a function of wavenumber when U = 4, V = 0 and {beta} = 16. For a half-filled band we find no indications of a Fermi surface. Slightly away from half-filling, we find Fermi-surface-like behavior similar to that found in other numerical investigations. 8 refs., 3 figs.

•Jan 1, 1989•0 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 discuss a new method to perform numerical simulations of one-dimensional systems with fermion and boson degrees of freedom. The method is based on a direct-space, imaginary-time representation of the fermion field. It is fast so that systems having up to 100 sites can easily be simulated. In addition, the method provides an intuitive physical "picture" of the ground state of a one-dimensional many-body system. We discuss in detail how to implement the method and how to compute various physical quantities. In particular, we show how to extend the method to study averages of off-diagonal quantities in an occupation-number representation. To assess the accuracy of our procedure, we apply it to free fermions in one dimension and compare with exact results. We then study a model of spinless interacting fermions and obtain the expected phase structure and behavior of correlation functions. We also consider the extended Hubbard model at various points in its phase diagram and study the behavior of spin-density, charge-density, and pairing correlation functions. We then study the Gross-Neveu model and show how the behavior depends on the number of fermion flavors. Finally, we consider an electron-phonon model and study its behavior both in the one-particle polaron sector and in the half-filled-band case. Along the way we show pictures of the ground-state configurations that give physical insight into the properties of the systems, like charge-density-wave, spin-density-wave, and superconducting states, "fractional charges," and solitons. We conclude by comparing our method with other methods and discuss the possibility of extending it to higher dimensions.
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.
The Holstein model of spinless fermions interacting with dispersionless phonons in one dimension is studied by a Green's function Monte Carlo technique. The ground-state energy, first fermionic excited state, density wave correlations, and mean lattice displacements are calculated for lattices of up to 16 sites, for one fermion per two sites, i.e., a half-filled band. Results are obtained for values of the fermion hopping parameter of t=0.1\ensuremath{\omega}, \ensuremath{\omega}, and 10\ensuremath{\omega}, where \ensuremath{\omega} is the phonon frequency. At a finite fermion-phonon coupling g there is a transition from a metallic phase to an insulating phase in which there is charge-density-wave order. Finite-size scaling is found to hold in the metallic phase and is used to extract the coupling dependence of the Luttinger liquid parameters, ${\mathit{u}}_{\mathrm{\ensuremath{\rho}}}$ and ${\mathit{K}}_{\mathrm{\ensuremath{\rho}}}$, the velocity of charge excitations, and the correlation exponent, respectively. For free fermions (g=0) and for strong coupling (g2{\mathit{g}}^{2}g2\ensuremath{\gg}t\ensuremath{\omega}) our results agree well with known analytic results. \textcopyright{} 1996 The American Physical Society.
"We have studied the extended Hubbard and the Peierls-extended Hubbard models in one dimension using primarily computational methods. We discuss the derivation of the models and the computational methods used to study them, focusing primarily on our extensions over existing numerical methods. In particular, we present the construction of ""operator-loop"" updates for spin-1/2 fermions and the inclusion of electron-phonon interactions in the Stochastic Series Expansion Quantum Monte Carlo method and the use of Quantum Parallel Tempering to improve the quality of the numerical data. For the extended Hubbard model at half-filling, we have confirmed the existence of a recently reported bond-order wave phase for the ground state in a small region of parameter space between the well-known charge- and spin-density wave phases. We have studied in detail the nature of the associated quantum phase transitions between all the phases. Upon doping away from half-filling, the doubly-degenerate, spin and charge gapped long-range-ordered ground states of the model are shown to evolve to a Luther-Emery liquid with no charge gap but a finite spin gap. For the Peierls-extended Hubbard model, which includes both electron-phonon and electron-electron interactions, we have studied the dimerization transition of the half-filled system in the presence of fully quantum mechanical phonons. We show that the transition to a dimerized state occurs only when the electron-phonon coupling exceeds a finite critical value that depends on the bare phonon frequency. This is in contrast to previous studies in similar models with adiabatic phonons that had predicted a dimerized state for arbitrarily small values of the electron-phonon coupling."
We present, in this dissertation, a numerical simulation method to study interacting fermion systems. The general simulation procedures are discussed in connection with a description of a Quantum Monte Carlo simulation algorithm for interacting electrons on lattices. The algorithm presented here has been used to simulate interacting electrons on lattices, and it makes possible the study of substantially larger systems than can be studied by other numerical methods. As long as certain limits of applicability are respected, model Hamiltonians of interacting electrons can be studied without resort to uncontrolled approximations. The method then provides nearly exact solutions to model Hamiltonians of many-body systems, in the sense that the degree of error can be controlled. We also discuss some results obtained from simulations of the extended Hubbard model and and the periodic Anderson model. For the extended Hubbard model in two dimensions, different regions are identified where electron correlation produces antiferromagnetism, charge-density-waves, and singlet pairing superconducting behaviors. We also find regions where transitions from one type of correlation to another occur. For the symmetric periodic Anderson model, we observe the formation of local spin moments at high temperatures and their quenching at low temperatures, as in the single impurity Anderson model. In addition, we find antiferromagnetic interaction between f local moments at low temperatures which are not present in a single impurity problem.
We report preliminary results of an investigation of the thermodynamic properties of the extended Hubbard model in one- dimension, calculated with the world-line Monte Carlo method described by Hirsch et al. With strictly continuous world-lines, we are able to measure the expectation of operators that conserve fermion number locally, such as the energy and (spatial) occupation number. By permitting the world-lines to be broken'' stochastically, we may also measure the expectation of operators that conserve fermion number only globally, such as the single-particle Green's function. For a 32 site lattice we present preliminary calculations of the average electron occupancy as a function of wavenumber when U = 4, V = 0 and {beta} = 16. For a half-filled band we find no indications of a Fermi surface. Slightly away from half-filling, we find Fermi-surface-like behavior similar to that found in other numerical investigations. 8 refs., 3 figs.
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.