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:0708.1749 · Physical Review B

Efficiency of quantum Monte Carlo impurity solvers for the dynamical mean-field theory

Nils Blümer

Published Nov 28, 2007• 31 citations
arXivHTMLPDF

Abstract

Since the inception of the dynamical mean-field theory, numerous numerical studies have relied on the Hirsch-Fye quantum Monte Carlo (HF-QMC) method for solving the associated impurity problem. Recently developed continuous-time algorithms (CT-QMC) avoid the Trotter discretization error and allow for faster configuration updates, which makes them candidates for replacing HF-QMC. We demonstrate, however, that a state-of-the-art implementation of HF-QMC (with extrapolation of discretization $\ensuremath{\Delta}\ensuremath{\tau}\ensuremath{\rightarrow}0$) is competitive with CT-QMC. A quantitative analysis of Trotter errors in HF-QMC estimates and of appropriate $\ensuremath{\Delta}\ensuremath{\tau}$ values is included.

DOI
10.1103/physrevb.76.205120

Related papers

Multigrid Hirsch-Fye quantum Monte Carlo method for dynamical mean-field theory

Nils Blümer

We present a new algorithm which allows for direct numerically exact solutions within dynamical mean-field theory (DMFT). It is based on the established Hirsch-Fye quantum Monte Carlo (HF-QMC) method. However, the DMFT impurity model is solved not at fixed imaginary-time discretization Delta_tau, but for a range of discretization grids; by extrapolation, unbiased Green functions are obtained in each DMFT iteration. In contrast to conventional HF-QMC, the multigrid algorithm converges to the exact DMFT fixed points. It extends the useful range of Delta_tau, is precise and reliable even in the immediate vicinity of phase transitions and is more efficient, also in comparison to continuous-time methods. Using this algorithm, we show that the spectral weight transfer at the Mott transition has been overestimated in a recent density matrix renormalization group study.

arXiv:0801.1222•Jan 8, 2008•7 citations

Continuous-time auxiliary field Monte Carlo for quantum impurity models

Emanuel Gull, Olivier Parcollet, Matthias Troyer

We present a continuous-time Monte Carlo method for quantum impurity models, which combines a weak-coupling expansion with an auxiliary-field decomposition. The method is considerably more efficient than Hirsch-Fye and free of time discretization errors, and is particularly useful as impurity solver in large cluster dynamical mean-field theory (DMFT) calculations.

•Nov 20, 2012•122 citations

Continuous-time auxiliary-field Monte Carlo for quantum impurity models

Emanuel Gull, Philipp Werner, Olivier Parcollet +1 more

We present a continuous-time Monte Carlo method for quantum impurity models, which combines a weak-coupling expansion with an auxiliary-field decomposition. The method is considerably more efficient than Hirsch-Fye and free of time discretization errors, and is particularly useful as impurity solver in large cluster dynamical mean-field theory (DMFT) calculations.

arXiv:0802.3222•Europhysics Letters (EPL)•May 27, 2008•265 citations

Continuous-time quantum Monte Carlo solver for dynamical mean field theory in the compact Legendre representation

Evan Sheridan, Cédric Weber, Evgeny A. Plekhanov +1 more

Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to $d$- and $f$-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for ${\mathrm{SrVO}}_{3}$ properties are presented. The computational efficiency of this method is also demonstrated.

arXiv:1811.12869•Physical review. B./Physical review. B•May 30, 2019•14 citations

Numerically exact Green functions from Hirsch-Fye quantum Monte Carlo simulations

Nils Blümer

We present a new method for extracting numerically exact imaginary-time Green functions from standard Hirsch-Fye quantum Monte Carlo (HF-QMC) simulations within dynamical mean-field theory. By analytic continuation, angular resolved spectra are obtained without the discretization bias previously associated with HF-QMC results. The method is shown to be accurate even at very low temperatures (T=W/800 for bandwidth W) in the strongly correlated regime.

arXiv:0712.1290•Dec 8, 2007•1 citation

Quasi-continuous-time impurity solver for the dynamical mean-field theory with linear scaling in the inverse temperature

D. Rost, Fakher F. Assaad, Nils Blümer

We present an algorithm for solving the self-consistency equations of the dynamical mean-field theory (DMFT) with high precision and efficiency at low temperatures. In each DMFT iteration, the impurity problem is mapped to an auxiliary Hamiltonian, for which the Green function is computed by combining determinantal quantum Monte Carlo (BSS-QMC) calculations with a multigrid extrapolation procedure. The method is numerically exact, i.e., yields results which are free of significant Trotter errors, but retains the BSS advantage, compared to direct QMC impurity solvers, of linear (instead of cubic) scaling with the inverse temperature. The new algorithm is applied to the half-filled Hubbard model close to the Mott transition; detailed comparisons with exact diagonalization, Hirsch-Fye QMC, and continuous-time QMC are provided.

arXiv:1303.2004•Physical Review E•May 29, 2013•16 citations
Since the inception of the dynamical mean-field theory, numerous numerical studies have relied on the Hirsch-Fye quantum Monte Carlo (HF-QMC) method for solving the associated impurity problem. Recently developed continuous-time algorithms (CT-QMC) avoid the Trotter discretization error and allow for faster configuration updates, which makes them candidates for replacing HF-QMC. We demonstrate, however, that a state-of-the-art implementation of HF-QMC (with extrapolation of discretization \ensuremathΔ\ensuremathτ\ensuremath→0\ensuremath{\Delta}\ensuremath{\tau}\ensuremath{\rightarrow}0\ensuremathΔ\ensuremathτ\ensuremath→0) is competitive with CT-QMC. A quantitative analysis of Trotter errors in HF-QMC estimates and of appropriate \ensuremathΔ\ensuremathτ\ensuremath{\Delta}\ensuremath{\tau}\ensuremathΔ\ensuremathτ values is included.
We present a new algorithm which allows for direct numerically exact solutions within dynamical mean-field theory (DMFT). It is based on the established Hirsch-Fye quantum Monte Carlo (HF-QMC) method. However, the DMFT impurity model is solved not at fixed imaginary-time discretization Delta_tau, but for a range of discretization grids; by extrapolation, unbiased Green functions are obtained in each DMFT iteration. In contrast to conventional HF-QMC, the multigrid algorithm converges to the exact DMFT fixed points. It extends the useful range of Delta_tau, is precise and reliable even in the immediate vicinity of phase transitions and is more efficient, also in comparison to continuous-time methods. Using this algorithm, we show that the spectral weight transfer at the Mott transition has been overestimated in a recent density matrix renormalization group study.
Dynamical mean-field theory (DMFT) is one of the most widely-used methods to treat accurately electron correlation effects in ab-initio real material calculations. Many modern large-scale implementations of DMFT in electronic structure codes involve solving a quantum impurity model with a continuous-time quantum Monte Carlo (CT-QMC) solver [Rubtsov et al., Phys. Rev. B 72, 035122 (2005); Werner et al., Phys. Rev. Lett. 97, 076405 (2006); Werner and Millis, Phys. Rev. B 74, 155107 (2006); Gull et al., Rev. Mod. Phys. 83, 349 (2011)]. The main advantage of CT-QMC is that, unlike standard quantum Monte Carlo approaches, it is able to generate the local Green's functions $G(\ensuremath{\tau})$ of the correlated system on an arbitrarily fine imaginary time $\ensuremath{\tau}$ grid, and is free of any systematic errors. In this work, we extend a hybrid QMC solver proposed by Khatami et al. [Phys. Rev. E 81, 056703 (2010)] and Rost et al. [Phys. Rev. E 87, 053305 (2013)] to a multiorbital context. This has the advantage of enabling impurity solver QMC calculations to scale linearly with inverse temperature $\ensuremath{\beta}$, and permit its application to ddd- and fff-band materials. In addition, we present a Green's-function processing scheme which generates accurate quasicontinuous imaginary time solutions of the impurity problem which overcome errors inherent to standard QMC approaches. This solver and processing scheme are incorporated into a full DFT+DMFT calculation using the castep DFT code [Clark et al., Z. Kristallogr. 220, 567 (2005)]. Benchmark calculations for SrVO3{\mathrm{SrVO}}_{3}SrVO3​ properties are presented. The computational efficiency of this method is also demonstrated.
We present an algorithm for solving the self-consistency equations of the dynamical mean-field theory (DMFT) with high precision and efficiency at low temperatures. In each DMFT iteration, the impurity problem is mapped to an auxiliary Hamiltonian, for which the Green function is computed by combining determinantal quantum Monte Carlo (BSS-QMC) calculations with a multigrid extrapolation procedure. The method is numerically exact, i.e., yields results which are free of significant Trotter errors, but retains the BSS advantage, compared to direct QMC impurity solvers, of linear (instead of cubic) scaling with the inverse temperature. The new algorithm is applied to the half-filled Hubbard model close to the Mott transition; detailed comparisons with exact diagonalization, Hirsch-Fye QMC, and continuous-time QMC are provided.