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/0305387 · Physical Review B

Quantum Monte Carlo and variational approaches to the Holstein model

Martin Hohenadler, Hans Gerd Evertz, Wolfgang von der Linden

Published Jan 15, 2004• 86 citations
arXivHTMLPDF

Abstract

Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principal component representation of the phonon degrees of freedom allows to sample completely uncorrelated phonon configurations. The combination of these elements enables us to perform efficient simulations for a wide range of temperature, phonon frequency, and electron-phonon coupling on clusters large enough to avoid finite-size effects. The algorithm is tested in one dimension and the data are compared with exact-diagonalization results and with existing work. Moreover, the ideas presented here can also be applied to the many-electron case. In the one-electron case considered here, the physics of the Holstein model can be described by a simple variational approach.

DOI
10.1103/physrevb.69.024301

Related papers

New QMC and Polaron Approaches to the Holstein model

Martin Hohenadler, Hans Gerd Evertz, Wolfgang von der Linden

Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient Quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principle component representation of the phonon degrees of freedom allows to sample completely uncorrelated phonon configurations. The combination of these elements enables us to perform efficient simulations for a wide range of temperature, phonon frequency and electron-phonon coupling on clusters large enough to avoid finite-size effects. The algorithm is tested in one dimension and the data are compared with exact-diagonalization results and with existing work. Moreover, the new ideas presented here can also be applied to the many-electron case. In the one-electron limit the physics of the Holstein model can be described by a simple variational approach.

Physical Review B•Jan 1, 2003•0 citations

Variational Monte Carlo method for electron-phonon coupled systems

Takahiro Ohgoe, Masatoshi Imada

We develop a variational Monte Carlo (VMC) method for electron-phonon coupled systems. The VMC method has been extensively used for investigating strongly correlated electrons over the last decades. However, its applications to electron-phonon coupled systems have been severely restricted because of its large Hilbert space. Here, we propose a variational wave function with a large number of variational parameters, which is suitable and tractable for systems with electron-phonon coupling. In the proposed wave function, we implement an unexplored electron-phonon correlation factor, which takes into account the effect of the entanglement between electrons and phonons. The method is applied to systems with diagonal electron-phonon interactions, i.e., interactions between charge densities and lattice displacements (phonons). As benchmarks, we compare VMC results with previous results obtained by the exact diagonalization, the Green function Monte Carlo method and the density matrix renormalization group for the Holstein and Holstein-Hubbard model. From these benchmarks, we show that the present method offers an efficient way to treat strongly coupled electron-phonon systems.

arXiv:1403.4150•Physical Review B•May 29, 2014•28 citations

Fast and scalable quantum Monte Carlo simulations of electron-phonon models

Benjamin Cohen-Stead, Owen Bradley, Cole Miles +3 more

We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and we report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but it incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.

arXiv:2203.01291•Physical review. E•Jun 15, 2022•35 citations

Temperature and quantum phonon effects on Holstein-Hubbard bipolarons

Martin Hohenadler, Wolfgang von der Linden

The one-dimensional Holstein-Hubbard model with two electrons of opposite spin is studied using an extension of a recently developed quantum Monte Carlo method and a very simple yet rewarding variational approach, both based on a canonically transformed Hamiltonian. The quantum Monte Carlo method yields very accurate results in the regime of small but finite phonon frequencies, characteristic of many strongly correlated materials, e.g., the cuprates and the manganites. The influence of electron-electron repulsion, phonon frequency, and temperature on the bipolaron state is investigated. Thermal dissociation of the intersite bipolaron is observed at high temperatures, and its relation to an existing theory of the manganites is discussed.

arXiv:cond-mat/0409573•Physical Review B•May 31, 2005•28 citations

Variational and auxiliary field Monte Carlo for the Hubbard and Hubbard-Holstein models: an accurate finite-size scaling and a "sign problem" solution

Seher Karakuzu

Understanding the properties of strongly correlated materials is a tough work since most of the times analytical tools cannot be applied due to the complicated nature of the problem. Here we apply quantum Monte Carlo techniques in order to tackle with one of the most important and unknown phenomena of many materials: the mechanism of superconductivity. We study the Hubbard-Holstein model which is one of the most simple theoretical models including strong correlation and electron-phonon coupling. By using twist-averaging boundary conditions we eliminate finite-size errors in the Monte Carlo simulations. With this useful tool we first solve the Hubbard-Holstein model by variational Monte Carlo by reporting the phase diagram of the model for different phonon frequencies at half filling. Then we investigate the phase separation away from half filling. Finally we try to attack the problem with an essentially unbiased method based on the auxiliary field quantum Monte Carlo technique and accelerated Langevin dynamics. By curing the sign problem via Cauchy integration in the complex plane where the auxiliary fields have been defined, we report the effect of the electron-phonon coupling on some observables such as magnetization, spin and charge structure factors. We show preliminary results that are already meaningful to understand the nature of the transition between magnetism and charge order in the model at half filling.

•Oct 5, 2018•0 citations

Langevin simulations of a long-range electron-phonon model

G. G. Batrouni, Richard T. Scalettar

We present a quantum Monte Carlo (QMC) study, based on the Langevin equation, of a Hamiltonian describing electrons coupled to phonon degrees of freedom. The bosonic part of the action helps control the variation of the field in imaginary time. As a consequence, the iterative conjugate gradient solution of the fermionic action, which depends on the boson coordinates, converges more rapidly than in the case of electron-electron interactions, such as the Hubbard Hamiltonian. Fourier acceleration is shown to be a crucial ingredient in reducing the equilibration and autocorrelation times. After describing and benchmarking the method, we present results for the phase diagram focusing on the range of the electron-phonon interaction. We delineate the regions of charge density wave formation from those in which the fermion density is inhomogeneous, caused by phase separation. We show that the Langevin approach is more efficient than the determinant QMC method for lattice sizes $N\ensuremath{\gtrsim}8\ifmmode\times\else\texttimes\fi{}8$ and that it therefore opens a potential path to problems including, for example, charge order in the three-dimensional Holstein model.

arXiv:1808.08973•Physical review. B./Physical review. B•Jan 8, 2019•51 citations
Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principal component representation of the phonon degrees of freedom allows to sample completely uncorrelated phonon configurations. The combination of these elements enables us to perform efficient simulations for a wide range of temperature, phonon frequency, and electron-phonon coupling on clusters large enough to avoid finite-size effects. The algorithm is tested in one dimension and the data are compared with exact-diagonalization results and with existing work. Moreover, the ideas presented here can also be applied to the many-electron case. In the one-electron case considered here, the physics of the Holstein model can be described by a simple variational approach.
Based on the canonical Lang-Firsov transformation of the Hamiltonian we develop a very efficient Quantum Monte Carlo algorithm for the Holstein model with one electron. Separation of the fermionic degrees of freedom by a reweighting of the probability distribution leads to a dramatic reduction in computational effort. A principle component representation of the phonon degrees of freedom allows to sample completely uncorrelated phonon configurations. The combination of these elements enables us to perform efficient simulations for a wide range of temperature, phonon frequency and electron-phonon coupling on clusters large enough to avoid finite-size effects. The algorithm is tested in one dimension and the data are compared with exact-diagonalization results and with existing work. Moreover, the new ideas presented here can also be applied to the many-electron case. In the one-electron limit the physics of the Holstein model can be described by a simple variational approach.
We develop a variational Monte Carlo (VMC) method for electron-phonon coupled systems. The VMC method has been extensively used for investigating strongly correlated electrons over the last decades. However, its applications to electron-phonon coupled systems have been severely restricted because of its large Hilbert space. Here, we propose a variational wave function with a large number of variational parameters, which is suitable and tractable for systems with electron-phonon coupling. In the proposed wave function, we implement an unexplored electron-phonon correlation factor, which takes into account the effect of the entanglement between electrons and phonons. The method is applied to systems with diagonal electron-phonon interactions, i.e., interactions between charge densities and lattice displacements (phonons). As benchmarks, we compare VMC results with previous results obtained by the exact diagonalization, the Green function Monte Carlo method and the density matrix renormalization group for the Holstein and Holstein-Hubbard model. From these benchmarks, we show that the present method offers an efficient way to treat strongly coupled electron-phonon systems.
We introduce methodologies for highly scalable quantum Monte Carlo simulations of electron-phonon models, and we report benchmark results for the Holstein model on the square lattice. The determinant quantum Monte Carlo (DQMC) method is a widely used tool for simulating simple electron-phonon models at finite temperatures, but it incurs a computational cost that scales cubically with system size. Alternatively, near-linear scaling with system size can be achieved with the hybrid Monte Carlo (HMC) method and an integral representation of the Fermion determinant. Here, we introduce a collection of methodologies that make such simulations even faster. To combat "stiffness" arising from the bosonic action, we review how Fourier acceleration can be combined with time-step splitting. To overcome phonon sampling barriers associated with strongly bound bipolaron formation, we design global Monte Carlo updates that approximately respect particle-hole symmetry. To accelerate the iterative linear solver, we introduce a preconditioner that becomes exact in the adiabatic limit of infinite atomic mass. Finally, we demonstrate how stochastic measurements can be accelerated using fast Fourier transforms. These methods are all complementary and, combined, may produce multiple orders of magnitude speedup, depending on model details.
The one-dimensional Holstein-Hubbard model with two electrons of opposite spin is studied using an extension of a recently developed quantum Monte Carlo method and a very simple yet rewarding variational approach, both based on a canonically transformed Hamiltonian. The quantum Monte Carlo method yields very accurate results in the regime of small but finite phonon frequencies, characteristic of many strongly correlated materials, e.g., the cuprates and the manganites. The influence of electron-electron repulsion, phonon frequency, and temperature on the bipolaron state is investigated. Thermal dissociation of the intersite bipolaron is observed at high temperatures, and its relation to an existing theory of the manganites is discussed.
Understanding the properties of strongly correlated materials is a tough work since most of the times analytical tools cannot be applied due to the complicated nature of the problem. Here we apply quantum Monte Carlo techniques in order to tackle with one of the most important and unknown phenomena of many materials: the mechanism of superconductivity. We study the Hubbard-Holstein model which is one of the most simple theoretical models including strong correlation and electron-phonon coupling. By using twist-averaging boundary conditions we eliminate finite-size errors in the Monte Carlo simulations. With this useful tool we first solve the Hubbard-Holstein model by variational Monte Carlo by reporting the phase diagram of the model for different phonon frequencies at half filling. Then we investigate the phase separation away from half filling. Finally we try to attack the problem with an essentially unbiased method based on the auxiliary field quantum Monte Carlo technique and accelerated Langevin dynamics. By curing the sign problem via Cauchy integration in the complex plane where the auxiliary fields have been defined, we report the effect of the electron-phonon coupling on some observables such as magnetization, spin and charge structure factors. We show preliminary results that are already meaningful to understand the nature of the transition between magnetism and charge order in the model at half filling.
We present a quantum Monte Carlo (QMC) study, based on the Langevin equation, of a Hamiltonian describing electrons coupled to phonon degrees of freedom. The bosonic part of the action helps control the variation of the field in imaginary time. As a consequence, the iterative conjugate gradient solution of the fermionic action, which depends on the boson coordinates, converges more rapidly than in the case of electron-electron interactions, such as the Hubbard Hamiltonian. Fourier acceleration is shown to be a crucial ingredient in reducing the equilibration and autocorrelation times. After describing and benchmarking the method, we present results for the phase diagram focusing on the range of the electron-phonon interaction. We delineate the regions of charge density wave formation from those in which the fermion density is inhomogeneous, caused by phase separation. We show that the Langevin approach is more efficient than the determinant QMC method for lattice sizes $N\ensuremath{\gtrsim}8\ifmmode\times\else\texttimes\fi{}8$ and that it therefore opens a potential path to problems including, for example, charge order in the three-dimensional Holstein model.