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:0712.1290

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

Nils Blümer

Published Dec 8, 2007• 1 citation
arXivHTMLPDF

Abstract

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.

DOI
10.48550/arxiv.0712.1290

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

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

Nils Blümer

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.

arXiv:0708.1749•Physical Review B•Nov 28, 2007•31 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

Cluster solver for dynamical mean-field theory with linear scaling in inverse temperature

Ehsan Khatami, C. R. Lee, Zhaojun Bai +2 more

Dynamical mean-field theory and its cluster extensions provide a very useful approach for examining phase transitions in model Hamiltonians and, in combination with electronic structure theory, constitute powerful methods to treat strongly correlated materials. The key advantage to the technique is that, unlike competing real-space methods, the sign problem is well controlled in the Hirsch-Fye (HF) quantum Monte Carlo used as an exact cluster solver. However, an important computational bottleneck remains; the HF method scales as the cube of the inverse temperature, β . This often makes simulations at low temperatures extremely challenging. We present here a method based on determinant quantum Monte Carlo which scales linearly in β , with a quadratic term that comes in to play for the number of time slices larger than hundred, and demonstrate that the sign problem is identical to HF.

arXiv:0904.1239•Physical Review E•May 12, 2010•18 citations

Flat histogram quantum Monte Carlo for analytic continuation to real time

Diamantis, Nikolaos G., Efstratios Manousakis

The Quantum Monte Carlo (QMC) method can yield the imaginary-time dependence of a correlation function $C(τ)$ of an operator $\hat O$. The analytic continuation to real-time proceeds by means of a "numerical inversion" of these data to find the response function or spectral density $A(ω)$ corresponding to $\hat O$. Such a technique is very sensitive to the statistical errors in $C(τ)$ especially for large values of $τ$, when we are interested in the low-energy excitations. In this paper, we find that if we use the flat histogram technique in the QMC method, in such a way to make the {\it histogram of} $C(τ)$ flat, the results of the analytic continuation for low-energy excitations improve using the same amount of computational time. To demonstrate the idea we select an exactly soluble version of the single-hole motion in the $t-J$ model and the diagrammatic Monte Carlo technique.

arXiv:1411.5984•Nov 21, 2014•0 citations

Green's functions from real-time bold-line Monte Carlo

Guy Cohen, David R. Reichman, Andrew J. Millis +1 more

We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation.

arXiv:1401.0526•Physical Review B•Mar 31, 2014•77 citations
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.
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.
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.
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.
Dynamical mean-field theory and its cluster extensions provide a very useful approach for examining phase transitions in model Hamiltonians and, in combination with electronic structure theory, constitute powerful methods to treat strongly correlated materials. The key advantage to the technique is that, unlike competing real-space methods, the sign problem is well controlled in the Hirsch-Fye (HF) quantum Monte Carlo used as an exact cluster solver. However, an important computational bottleneck remains; the HF method scales as the cube of the inverse temperature, β . This often makes simulations at low temperatures extremely challenging. We present here a method based on determinant quantum Monte Carlo which scales linearly in β , with a quadratic term that comes in to play for the number of time slices larger than hundred, and demonstrate that the sign problem is identical to HF.
The Quantum Monte Carlo (QMC) method can yield the imaginary-time dependence of a correlation function C(τ)C(τ)C(τ) of an operator O^\hat OO^. The analytic continuation to real-time proceeds by means of a "numerical inversion" of these data to find the response function or spectral density A(ω)A(ω)A(ω) corresponding to O^\hat OO^. Such a technique is very sensitive to the statistical errors in C(τ)C(τ)C(τ) especially for large values of τττ, when we are interested in the low-energy excitations. In this paper, we find that if we use the flat histogram technique in the QMC method, in such a way to make the {\it histogram of} C(τ)C(τ)C(τ) flat, the results of the analytic continuation for low-energy excitations improve using the same amount of computational time. To demonstrate the idea we select an exactly soluble version of the single-hole motion in the t−Jt-Jt−J model and the diagrammatic Monte Carlo technique.
We present two methods for computing two-time correlation functions or Green's functions from real-time bold-line continuous-time quantum Monte Carlo. One method is a formally exact generalized auxiliary lead formalism by which spectral properties may be obtained from single-time observables. The other involves the evaluation of diagrams contributing to two-time observables directly on the Keldysh contour. Additionally, we provide a detailed description of the bold-line Monte Carlo method. Our methods are general and numerically exact, and able to reliably resolve high-energy features such as band edges. We compare the spectral functions obtained from real-time methods to analytically continued spectral functions obtained from imaginary-time Monte Carlo, thus probing the limits of analytic continuation.