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

Comparative study of numerical techniques for the simulation of a one-dimensional spinless fermion system

G.M. Buend�a

Published Nov 1, 1990• 0 citations
Publisher

Abstract

A comparative study of the efficiency of different algorithms (Metropolis, Langevin, and Hybrid) for the simulation of a one-dimensional spinless fermion model with a half-filled band is presented. The Metropolis updating of the fields is performed with an algorithm introduced by Blankenbecler, Scalapino, and Sugar. The density-density correlation function, the order parameter, and the autocorrelation are calculated. A study of the dependence of the systematic errors on the temperature and the step size for the Langevin and hybrid algorithms is presented. The relation between the accuracy and the speed of the algorithms is studied.

DOI
10.1103/physrevb.42.8163

Related papers

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

Comparative study of the discrete and the continuous Hubbard-Stratonovich transformation for a one-dimensional spinless fermion model

G. M. Buendía

A comparative study of the discrete and the continuous Hubbard-Stratonovich transformation is presented. With the use of the Metropolis algorithm, the Monte Carlo time correlation functions for a one-dimensional spinless fermion model are calculated for both methods. The results show the discrete method to be more efficient.

Physical review. B, Condensed matter•Mar 1, 1986•26 citations

Monte Carlo simulations of one-dimensional fermion systems

J. E. Hirsch, Robert L. Sugar, Douglas J. Scalapino +1 more

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.

Physical review. B, Condensed matter•Nov 1, 1982•487 citations

Global Monte Carlo algorithms for many-fermion systems

Michael J. Creutz

I discuss algorithms for simulating many-fermion systems via global updatings of auxiliary fields followed by an accept-reject stage which eliminates finite-step-size errors. When the system size is larger than the correlation length, these procedures should require computer time growing only slightly faster than linearly with the system volume V. A corrected Langevin scheme should asymptotically display a ${V}^{4/3}$ behavior, while the hybrid Monte Carlo scheme can behave as ${V}^{5/4}$. I present some tests of the latter algorithm on a simple model of interacting electrons on a two-dimensional lattice.

Physical review. D. Particles, fields, gravitation, and cosmology/Physical review. D. Particles and fields•Aug 15, 1988•151 citations

A New Hybrid Algorithm for Fermion Simulations Applied to the Hubbard Model

Masatoshi Imada

A new hybrid-type algorithm is proposed for simulations of fermion lattices. The direct calculation of the matrix inversion operation is replaced with a stochastic iteration procedure. The computation time grows linearly with the increase of spatial system size. Applications to the Hubbard model on the 6 2 , 12 2 and 16 2 square lattices are illustrated with no indication of enhancement for superconducting correlation in the non-half-filled band.

Journal of the Physical Society of Japan•Jan 15, 1988•16 citations

Stochastic series expansion simulation of the t−V model

Lei Wang, Ye-Hua Liu, Matthias Troyer

We present an algorithm for the efficient simulation of the half-filled spinless $t\text{\ensuremath{-}}V$ model on bipartite lattices, which combines the stochastic series expansion method with determinantal quantum Monte Carlo techniques widely used in fermionic simulations. The algorithm scales linearly in the inverse temperature, cubically with the system size, and is free from the time-discretization error. We use it to map out the finite-temperature phase diagram of the spinless $t\text{\ensuremath{-}}V$ model on the honeycomb lattice and observe a suppression of the critical temperature of the charge-density-wave phase in the vicinity of a fermionic quantum critical point.

arXiv:1602.02095•Physical review. B./Physical review. B•Apr 11, 2016•28 citations
A comparative study of the efficiency of different algorithms (Metropolis, Langevin, and Hybrid) for the simulation of a one-dimensional spinless fermion model with a half-filled band is presented. The Metropolis updating of the fields is performed with an algorithm introduced by Blankenbecler, Scalapino, and Sugar. The density-density correlation function, the order parameter, and the autocorrelation are calculated. A study of the dependence of the systematic errors on the temperature and the step size for the Langevin and hybrid algorithms is presented. The relation between the accuracy and the speed of the algorithms is studied.
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 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.
I discuss algorithms for simulating many-fermion systems via global updatings of auxiliary fields followed by an accept-reject stage which eliminates finite-step-size errors. When the system size is larger than the correlation length, these procedures should require computer time growing only slightly faster than linearly with the system volume V. A corrected Langevin scheme should asymptotically display a V4/3{V}^{4/3}V4/3 behavior, while the hybrid Monte Carlo scheme can behave as V5/4{V}^{5/4}V5/4. I present some tests of the latter algorithm on a simple model of interacting electrons on a two-dimensional lattice.
We present an algorithm for the efficient simulation of the half-filled spinless $t\text{\ensuremath{-}}V$ model on bipartite lattices, which combines the stochastic series expansion method with determinantal quantum Monte Carlo techniques widely used in fermionic simulations. The algorithm scales linearly in the inverse temperature, cubically with the system size, and is free from the time-discretization error. We use it to map out the finite-temperature phase diagram of the spinless $t\text{\ensuremath{-}}V$ model on the honeycomb lattice and observe a suppression of the critical temperature of the charge-density-wave phase in the vicinity of a fermionic quantum critical point.