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

Diagrammatic determinantal quantum Monte Carlo methods: Projective schemes and applications to the Hubbard-Holstein model

Fakher F. Assaad, Thomas C. Lang

Published Jul 23, 2007• 126 citations
arXivHTMLPDF

Abstract

We extend the weak-coupling diagrammatic determinantal algorithm to projective schemes as well as to the inclusion of phonon degrees of freedom. The projective approach provides a very efficient algorithm to access zero-temperature properties. To implement phonons, we integrate them out in favor of a retarded density-density interaction and simulate the resulting purely electronic action with the weak-coupling diagrammatic determinantal algorithm. Both extensions are tested within the dynamical mean-field approximation for the Hubbard and Hubbard-Holstein models.

DOI
10.1103/physrevb.76.035116

Related papers

Wirkungsbasierte Quanten-Monte-Carlo-Methoden für Fermion-Boson-Modelle

Weber, Manuel

This work deals with the development and application of novel quantum Monte Carlo methods to simulate fermion-boson models. Our developments are based on the path-integral formalism, where the bosonic degrees of freedom are integrated out exactly to obtain a retarded fermionic interaction. We give an overview of three methods that can be used to simulate retarded interactions. In particular, we develop a novel quantum Monte Carlo method with global directed-loop updates that solves the autocorrelation problem of previous approaches and scales linearly with system size. We demonstrate its efficiency for the Peierls transition in the Holstein model and discuss extensions to other fermion-boson models as well as spin-boson models. Furthermore, we show how with the help of generating functionals bosonic observables can be recovered directly from the Monte Carlo configurations. This includes estimators for the boson propagator, the fidelity susceptibility, and the specific heat of the Holstein model. The algorithmic developments of this work allow us to study the specific heat of the spinless Holstein model covering its entire parameter range. Its key features are explained from the single-particle spectral functions of electrons and phonons. In the adiabatic limit, the spectral properties are calculated exactly as a function of temperature using a classical Monte Carlo method and compared to results for the Su-Schrieffer-Heeger model.

•Jan 1, 2019•0 citations

Determinant quantum Monte Carlo study of the two-dimensional single-band Hubbard-Holstein model

Steven Johnston, Elizabeth Nowadnick, Y. F. Kung +3 more

We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.

arXiv:1306.2968•Physical Review B•Jun 24, 2013•82 citations

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

Dynamical vertex approximation for retarded interactions: Application to the Hubbard-Holstein model

Emin Moghadas, Samuele Giuli, Massimo Capone +1 more

We extend the dynamical vertex approximation (D$Γ$A), in its most commonly used ladder-version, to systems with retarded electronic interactions. Such an extended scheme, which we name "dynamical-$U$ D$Γ$A", allows us to study strongly interacting electrons coupled to bosons including non-local spatial correlations on all length scales on top of the purely local ones captured by dynamical mean-field theory. We demonstrate the applicability of the method in the specific case of the coupling to a single dispersionless phononic mode, by studying the Hubbard-Holstein model on the square and cubic lattices, focusing on regimes of sizable Hubbard interactions where strong electronic correlations compete with phonon-mediated interactions. As illustrative examples, we consider the antiferromagnetic phase in the three-dimensional model and the d-wave superconductor in two dimensions, analyzing the effect of the electron-phonon coupling, including its frequency dependence, on both instabilities.

arXiv:2610.08512•Oct 6, 2026•0 citations

Efficient Dynamical Mean Field Simulation of the Holstein-Hubbard Model

Philipp Werner, Andrew J. Millis

We present a method for solving impurity models with electron-phonon coupling, which treats the phonons efficiently and without approximations. The algorithm is applied to the Holstein-Hubbard model in the dynamical mean field approximation, where it allows access to strong interactions, very low temperatures, and arbitrary fillings. We show that a renormalized Migdal-Eliashberg theory provides a reasonable description of the phonon contribution to the electronic self-energy in strongly doped systems, but fails if the quasiparticle energy becomes of order of the phonon frequency.

arXiv:cond-mat/0701730•Physical Review Letters•Oct 4, 2007•180 citations

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
We extend the weak-coupling diagrammatic determinantal algorithm to projective schemes as well as to the inclusion of phonon degrees of freedom. The projective approach provides a very efficient algorithm to access zero-temperature properties. To implement phonons, we integrate them out in favor of a retarded density-density interaction and simulate the resulting purely electronic action with the weak-coupling diagrammatic determinantal algorithm. Both extensions are tested within the dynamical mean-field approximation for the Hubbard and Hubbard-Holstein models.
This work deals with the development and application of novel quantum Monte Carlo methods to simulate fermion-boson models. Our developments are based on the path-integral formalism, where the bosonic degrees of freedom are integrated out exactly to obtain a retarded fermionic interaction. We give an overview of three methods that can be used to simulate retarded interactions. In particular, we develop a novel quantum Monte Carlo method with global directed-loop updates that solves the autocorrelation problem of previous approaches and scales linearly with system size. We demonstrate its efficiency for the Peierls transition in the Holstein model and discuss extensions to other fermion-boson models as well as spin-boson models. Furthermore, we show how with the help of generating functionals bosonic observables can be recovered directly from the Monte Carlo configurations. This includes estimators for the boson propagator, the fidelity susceptibility, and the specific heat of the Holstein model. The algorithmic developments of this work allow us to study the specific heat of the spinless Holstein model covering its entire parameter range. Its key features are explained from the single-particle spectral functions of electrons and phonons. In the adiabatic limit, the spectral properties are calculated exactly as a function of temperature using a classical Monte Carlo method and compared to results for the Su-Schrieffer-Heeger model.
We have performed numerical studies of the Hubbard-Holstein model in two dimensions using determinant quantum Monte Carlo (DQMC). Here, we present details of the method, emphasizing the treatment of the lattice degrees of freedom, and then study the filling and behavior of the fermion sign as a function of model parameters. We find a region of parameter space with large Holstein coupling where the fermion sign recovers despite large values of the Hubbard interaction. This indicates that studies of correlated polarons at finite carrier concentrations are likely accessible to DQMC simulations. We then restrict ourselves to the half-filled model and examine the evolution of the antiferromagnetic structure factor, other metrics for antiferromagnetic and charge-density-wave order, and energetics of the electronic and lattice degrees of freedom as a function of electron-phonon coupling. From this we find further evidence for a competition between charge-density-wave and antiferromagnetic orders at half-filling.
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.
We extend the dynamical vertex approximation (DΓΓΓA), in its most commonly used ladder-version, to systems with retarded electronic interactions. Such an extended scheme, which we name "dynamical-UUU DΓΓΓA", allows us to study strongly interacting electrons coupled to bosons including non-local spatial correlations on all length scales on top of the purely local ones captured by dynamical mean-field theory. We demonstrate the applicability of the method in the specific case of the coupling to a single dispersionless phononic mode, by studying the Hubbard-Holstein model on the square and cubic lattices, focusing on regimes of sizable Hubbard interactions where strong electronic correlations compete with phonon-mediated interactions. As illustrative examples, we consider the antiferromagnetic phase in the three-dimensional model and the d-wave superconductor in two dimensions, analyzing the effect of the electron-phonon coupling, including its frequency dependence, on both instabilities.
We present a method for solving impurity models with electron-phonon coupling, which treats the phonons efficiently and without approximations. The algorithm is applied to the Holstein-Hubbard model in the dynamical mean field approximation, where it allows access to strong interactions, very low temperatures, and arbitrary fillings. We show that a renormalized Migdal-Eliashberg theory provides a reasonable description of the phonon contribution to the electronic self-energy in strongly doped systems, but fails if the quasiparticle energy becomes of order of the phonon frequency.
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.