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

New QMC and Polaron Approaches to the Holstein model

Martin Hohenadler, Hans Gerd Evertz, Wolfgang von der Linden

Published Jan 1, 2003• 0 citations
OpenAlex

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 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.

Related papers

Quantum Monte Carlo and variational 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 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.

arXiv:cond-mat/0305387•Physical Review B•Jan 15, 2004•86 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

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

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

Electron-Phonon Systems on a Universal Quantum Computer

Alexandru Macridin, Panagiotis Spentzouris, James Frederick Amundson +1 more

We present an algorithm that extends existing quantum algorithms for simulating fermion systems in quantum chemistry and condensed matter physics to include bosons in general and phonons in particular. We introduce a qubit representation for the low-energy subspace of phonons which allows an efficient simulation of the evolution operator of the electron-phonon systems. As a consequence of the Nyquist-Shannon sampling theorem, the phonons are represented with exponential accuracy on a discretized Hilbert space with a size that increases linearly with the cutoff of the maximum phonon number. The additional number of qubits required by the presence of phonons scales linearly with the size of the system. The additional circuit depth is constant for systems with finite-range electron-phonon and phonon-phonon interactions and linear for long-range electron-phonon interactions. Our algorithm for a Holstein polaron problem was implemented on an Atos quantum learning machine quantum simulator employing the quantum phase estimation method. The energy and the phonon number distribution of the polaron state agree with exact diagonalization results for weak, intermediate, and strong electron-phonon coupling regimes.

arXiv:1802.07347•Physical Review Letters•Sep 12, 2018•124 citations

Polaron band formation in the Holstein model

Gerhard Wellein, Holger Fehske

We present numerical exact results for the polaronic band structure of the Holstein molecular crystal model in one and two dimensions. The use of a direct Lanczos diagonalization technique, preserving the full dynamics and quantum nature of phonons, allows us to analyze in detail the renormalization of both quasiparticle bandwidth and dispersion by the electron-phonon interaction. For the two-dimensional case some of our exact data are compared with the results obtained in the framework of a recently developed finite-cluster strong-coupling perturbation theory.

arXiv:cond-mat/9703041•Physical review. B, Condensed matter•Aug 15, 1997•138 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 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.
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.
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.
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.
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.
We present an algorithm that extends existing quantum algorithms for simulating fermion systems in quantum chemistry and condensed matter physics to include bosons in general and phonons in particular. We introduce a qubit representation for the low-energy subspace of phonons which allows an efficient simulation of the evolution operator of the electron-phonon systems. As a consequence of the Nyquist-Shannon sampling theorem, the phonons are represented with exponential accuracy on a discretized Hilbert space with a size that increases linearly with the cutoff of the maximum phonon number. The additional number of qubits required by the presence of phonons scales linearly with the size of the system. The additional circuit depth is constant for systems with finite-range electron-phonon and phonon-phonon interactions and linear for long-range electron-phonon interactions. Our algorithm for a Holstein polaron problem was implemented on an Atos quantum learning machine quantum simulator employing the quantum phase estimation method. The energy and the phonon number distribution of the polaron state agree with exact diagonalization results for weak, intermediate, and strong electron-phonon coupling regimes.
We present numerical exact results for the polaronic band structure of the Holstein molecular crystal model in one and two dimensions. The use of a direct Lanczos diagonalization technique, preserving the full dynamics and quantum nature of phonons, allows us to analyze in detail the renormalization of both quasiparticle bandwidth and dispersion by the electron-phonon interaction. For the two-dimensional case some of our exact data are compared with the results obtained in the framework of a recently developed finite-cluster strong-coupling perturbation theory.