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:1609.02026 · Physical review. B./Physical review. B

Continuous-time quantum Monte Carlo for fermion-boson lattice models: Improved bosonic estimators and application to the Holstein model

Manuel Weber, Fakher F. Assaad, Martin Hohenadler

Published Dec 27, 2016• 15 citations
arXivHTMLPDF

Abstract

We extend the continuous-time interaction-expansion quantum Monte Carlo method with respect to measuring observables for fermion-boson lattice models. Using generating functionals, we express expectation values involving boson operators, which are not directly accessible because simulations are done in terms of a purely fermionic action, as integrals over fermionic correlation functions. We also demonstrate that certain observables can be inferred directly from the vertex distribution, and present efficient estimators for the total energy and the phonon propagator of the Holstein model. Furthermore, we generalize the covariance estimator of the fidelity susceptibility, an unbiased diagnostic for phase transitions, to the case of retarded interactions. The new estimators are applied to half-filled spinless and spinful Holstein models in one dimension. The observed renormalization of the phonon mode across the Peierls transition in the spinless model suggests a soft-mode transition in the adiabatic regime. The critical point is associated with a minimum in the phonon kinetic energy and a maximum in the fidelity susceptibility.

DOI
10.1103/physrevb.94.245138

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

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

Directed-Loop Quantum Monte Carlo Method for Retarded Interactions

Manuel Weber, Fakher F. Assaad, Martin Hohenadler

The directed-loop quantum Monte Carlo method is generalized to the case of retarded interactions. Using the path integral, fermion-boson or spin-boson models are mapped to actions with retarded interactions by analytically integrating out the bosons. This yields an exact algorithm that combines the highly efficient loop updates available in the stochastic series expansion representation with the advantages of avoiding a direct sampling of the bosons. The application to electron-phonon models reveals that the method overcomes the previously detrimental issues of long autocorrelation times and exponentially decreasing acceptance rates. For example, the resulting dramatic speedup allows us to investigate the Peierls quantum phase transition on chains of up to 1282 sites.

arXiv:1704.07913•Physical Review Letters•Aug 31, 2017•19 citations

Phonon spectral function of the one-dimensional Holstein-Hubbard model

Manuel Weber, Fakher F. Assaad, Martin Hohenadler

We use the continuous-time interaction expansion (CT-INT) quantum Monte Carlo method to calculate the phonon spectral function of the one-dimensional Holstein-Hubbard model at half-filling. Our results are consistent with a soft-mode Peierls transition in the adiabatic regime, and the existence of a central peak related to long-range order in the Peierls phase. We explain a previously observed feature at small momenta in terms of a hybridization of charge and phonon excitations. Tuning the system from a Peierls to a metallic phase with a nonzero Hubbard interaction suppresses the central peak, but a significant renormalization of the phonon dispersion remains. In contrast, the dispersion is only weakly modified in the Mott phase. We discuss finite-size effects, the relation to the dynamic charge structure factor, as well as additional sum rules and their implications. Finally, we reveal the existence of a discrete symmetry in a continuum field theory of the Holstein model, which is spontaneously broken in the Peierls phase.

arXiv:1504.07446•Physical Review B•Jun 29, 2015•22 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

Effect of Electron-Phonon Interaction Range for a Half-Filled Band in One Dimension

Martin Hohenadler, Fakher F. Assaad, Holger Fehske

We demonstrate that fermion-boson models with nonlocal interactions can be simulated at finite band filling with the continuous-time quantum Monte Carlo method. We apply this method to explore the influence of the electron-phonon interaction range for a half-filled band in one dimension, covering the full range from the Holstein to the Fröhlich regime. The phase diagram contains metallic, Peierls, and phase-separated regions. Nonlocal interactions suppress the Peierls instability, and thereby lead to almost degenerate power-law exponents for charge and pairing correlations.

arXiv:1205.0612•Physical Review Letters•Sep 13, 2012•32 citations
We extend the continuous-time interaction-expansion quantum Monte Carlo method with respect to measuring observables for fermion-boson lattice models. Using generating functionals, we express expectation values involving boson operators, which are not directly accessible because simulations are done in terms of a purely fermionic action, as integrals over fermionic correlation functions. We also demonstrate that certain observables can be inferred directly from the vertex distribution, and present efficient estimators for the total energy and the phonon propagator of the Holstein model. Furthermore, we generalize the covariance estimator of the fidelity susceptibility, an unbiased diagnostic for phase transitions, to the case of retarded interactions. The new estimators are applied to half-filled spinless and spinful Holstein models in one dimension. The observed renormalization of the phonon mode across the Peierls transition in the spinless model suggests a soft-mode transition in the adiabatic regime. The critical point is associated with a minimum in the phonon kinetic energy and a maximum in the fidelity susceptibility.
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.
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.
The directed-loop quantum Monte Carlo method is generalized to the case of retarded interactions. Using the path integral, fermion-boson or spin-boson models are mapped to actions with retarded interactions by analytically integrating out the bosons. This yields an exact algorithm that combines the highly efficient loop updates available in the stochastic series expansion representation with the advantages of avoiding a direct sampling of the bosons. The application to electron-phonon models reveals that the method overcomes the previously detrimental issues of long autocorrelation times and exponentially decreasing acceptance rates. For example, the resulting dramatic speedup allows us to investigate the Peierls quantum phase transition on chains of up to 1282 sites.
We use the continuous-time interaction expansion (CT-INT) quantum Monte Carlo method to calculate the phonon spectral function of the one-dimensional Holstein-Hubbard model at half-filling. Our results are consistent with a soft-mode Peierls transition in the adiabatic regime, and the existence of a central peak related to long-range order in the Peierls phase. We explain a previously observed feature at small momenta in terms of a hybridization of charge and phonon excitations. Tuning the system from a Peierls to a metallic phase with a nonzero Hubbard interaction suppresses the central peak, but a significant renormalization of the phonon dispersion remains. In contrast, the dispersion is only weakly modified in the Mott phase. We discuss finite-size effects, the relation to the dynamic charge structure factor, as well as additional sum rules and their implications. Finally, we reveal the existence of a discrete symmetry in a continuum field theory of the Holstein model, which is spontaneously broken in the Peierls phase.
"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 demonstrate that fermion-boson models with nonlocal interactions can be simulated at finite band filling with the continuous-time quantum Monte Carlo method. We apply this method to explore the influence of the electron-phonon interaction range for a half-filled band in one dimension, covering the full range from the Holstein to the Fröhlich regime. The phase diagram contains metallic, Peierls, and phase-separated regions. Nonlocal interactions suppress the Peierls instability, and thereby lead to almost degenerate power-law exponents for charge and pairing correlations.