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

Worm-improved estimators in continuous-time quantum Monte Carlo

Patrik Gunacker, Markus Wallerberger, T. Ribic, Andreas Hausoel, Giorgio Sangiovanni, Karsten Held

Published Sep 30, 2016• 49 citations
arXivHTMLPDF

Abstract

We derive the improved estimators for general interactions and employ these for the continuous-time quantum Monte Carlo method. Using a worm algorithm we show how measuring higher-ordered correlators leads to an improved high-frequency behavior in irreducible quantities such as the one-particle self-energy or the irreducible two-particle vertex for non-density-density interactions. A good knowledge of the asymptotics of the two-particle vertex is essential for calculating nonlocal electronic correlations using diagrammatic extensions to the dynamical mean field theory as well as for calculating susceptibilities. We test our algorithm against analytic results for the multiorbital atomic limit and the Falicov-Kimball model.

DOI
10.1103/physrevb.94.125153

Related papers

Continuous-time quantum Monte Carlo calculation of multiorbital vertex asymptotics

Josef Kaufmann, Patrik Gunacker, Karsten Held

We derive the equations for calculating the high-frequency asymptotics of the local two-particle vertex function for a multiorbital impurity model. These relate the asymptotics for a general local interaction to equal-time two-particle Green's functions, which we sample using continuous-time quantum Monte Carlo simulations with a worm algorithm. As specific examples we study the single-orbital Hubbard model and the three ${t}_{2g}$ orbitals of ${\mathrm{SrVO}}_{3}$ within dynamical mean-field theory (DMFT). We demonstrate how the knowledge of the high-frequency asymptotics reduces the statistical uncertainties of the vertex and further eliminates finite-box-size effects. The proposed method benefits the calculation of nonlocal susceptibilities in DMFT and diagrammatic extensions of DMFT.

arXiv:1703.09407•Physical review. B./Physical review. B•Jul 7, 2017•40 citations

Continuous-time quantum Monte Carlo using worm sampling

Patrik Gunacker, Markus Wallerberger, Emanuel Gull +3 more

We present a worm sampling method for calculating one- and two-particle Green's functions using continuous-time quantum Monte Carlo simulations in the hybridization expansion (CT-HYB). Instead of measuring Green's functions by removing hybridization lines from partition function configurations, as in conventional CT-HYB, the worm algorithm directly samples the Green's function. We show that worm sampling is necessary to obtain general two-particle Green's functions which are not of density-density type and that it improves the sampling efficiency when approaching the atomic limit. Such two-particle Green's functions are needed to compute off-diagonal elements of susceptibilities and occur in diagrammatic extensions of the dynamical mean-field theory and in efficient estimators for the single-particle self-energy.

arXiv:1506.01905•Physical Review B•Oct 1, 2015•51 citations

Improved estimators for the self-energy and vertex function in hybridization-expansion continuous-time quantum Monte Carlo simulations

Hartmut Hafermann, Kelly R. Patton, Philipp Werner

We propose efficient measurement procedures for the self-energy and vertex function of the Anderson impurity model within the hybridization expansion continuous-time quantum Monte Carlo algorithm. The method is based on the measurement of higher-order correlation functions related to the quantities being sought through the equation of motion, a technique previously introduced in the numerical renormalization-group context. For the case of interactions of density-density type, the additional correlators can be obtained at essentially no additional computational cost. In combination with a recently introduced method for filtering the Monte Carlo noise using a representation in terms of orthogonal polynomials, we obtain data with unprecedented accuracy. This leads to an enhanced stability in analytical continuations of the self-energy or in two-particle-based theories such as the dual fermion approach. As an illustration of the method we reexamine the previously reported spin-freezing and high-spin to low-spin transitions in a two-orbital model with density-density interactions. In both cases, the vertex function undergoes significant changes, which suggests significant corrections to the dynamical mean-field solutions in dual fermion calculations.

arXiv:1108.1936•Physical Review B•May 3, 2012•118 citations

Symmetric improved estimators for continuous-time quantum Monte Carlo

Josef Kaufmann, Patrik Gunacker, Alexander Kowalski +2 more

We derive equations of motion for Green's functions of the multiorbital Anderson impurity model by differentiating symmetrically with respect to all time arguments. The resulting equations relate the one- and two-particle Green's function to correlators of up to six particles at four times. As an application we consider continuous-time quantum Monte Carlo simulations in the hybridization expansion, which hitherto suffered from notoriously high noise levels at large Matsubara frequencies. Employing the derived symmetric improved estimators overcomes this problem.

arXiv:1906.00880•Physical review. B./Physical review. B•Aug 8, 2019•37 citations

Local and Nonlocal Correlations in Interacting Electron Systems

Josef Kaufmann

Strong electronic correlations are the root of many fascinating phenomena of condensed matter physics. However, they also make calculations notoriously difficult. This thesis is concerned with the development and improvement of computational methods for strongly correlated electron systems. After a brief recapitulation of the Green’s function formalism,which is the theoretical basis of most such methods, we discuss the analytical continuation from imaginary to real frequencies. Especially the maximum entropy method is explained in detail and an implementation there of is presented.The remainder of this thesis deals with methods to compute Green’s functions on thei maginary frequency axis. First, the Anderson impurity model (AIM) is considered. It can be solved numerically – a task for which several methods exist. We focus on continuous-time quantum Monte Carlo with worm sampling, through which various kinds of correlation functions can be calculated. On this basis we derive and implement symmetric improved estimators that allow for the computation of the self-energy and vertex functions where the level of stochastic noise is largely reduced with respect to conventional estimators.As a next step the dynamical mean-field theory (DMFT) is introduced as a method to capture local correlations in the Hubbard model (HM) by self-consistently mapping it onto the AIM. We study the convergence behavior of this iterative procedure and the improvement of convergence due to mixing. Then, we compare electron scattering rates from DMFT in the weak-coupling regime of the square-lattice HM to results from Boltzmann scattering theory.While local correlations account for many interesting phenomena, such as the Mott-Hubbard metal-insulator transition, they are not sufficient to correctly describe, e.g., the physics of phase transitions in two- or three-dimensional systems. It is thus the final step in this thesis to also include nonlocal correlations by means of the dynamical vertex approximation (DΓA). For this we develop a new self-consistency scheme, which is then applied to several systems: We study the square-lattice HM, where we compare to benchmark methods in the weak-coupling regime, but also present results for intermediate and strong coupling. Furthermore, a two-band model with doping and a three-band Wannier model for strontium vanadate are studied. As an application of our new self-consistent DΓA, we study the effects of geometrical frustration in the kagome lattice, where we also compare our results to determinant quantum Monte Carlo.

•May 1, 2021•1 citation

Continuous-time quantum Monte Carlo method for fermions: Beyond auxiliary field framework

Alexey N. Rubtsov, Alexander I. Lichtenstein

A numerically exact continuous-time quantum Monte Carlo algorithm for finite fermionic systems with nonlocal interactions is proposed. The scheme is particularly applicable for general multiband time-dependent correlations, since it does not invoke Hubbard-Stratonovich transformation. The present determinantal grand-canonical method is based on a stochastic series expansion for the partition function in the interaction representation. The results for the Green function and for the time-dependent susceptibility of multiorbital supersymmetric impurity model with a spin-flip interaction are presented.

arXiv:cond-mat/0408009•Journal of Experimental and Theoretical Physics Letters•Jul 1, 2004•101 citations
We derive the improved estimators for general interactions and employ these for the continuous-time quantum Monte Carlo method. Using a worm algorithm we show how measuring higher-ordered correlators leads to an improved high-frequency behavior in irreducible quantities such as the one-particle self-energy or the irreducible two-particle vertex for non-density-density interactions. A good knowledge of the asymptotics of the two-particle vertex is essential for calculating nonlocal electronic correlations using diagrammatic extensions to the dynamical mean field theory as well as for calculating susceptibilities. We test our algorithm against analytic results for the multiorbital atomic limit and the Falicov-Kimball model.
We derive the equations for calculating the high-frequency asymptotics of the local two-particle vertex function for a multiorbital impurity model. These relate the asymptotics for a general local interaction to equal-time two-particle Green's functions, which we sample using continuous-time quantum Monte Carlo simulations with a worm algorithm. As specific examples we study the single-orbital Hubbard model and the three t2g{t}_{2g}t2g​ orbitals of SrVO3{\mathrm{SrVO}}_{3}SrVO3​ within dynamical mean-field theory (DMFT). We demonstrate how the knowledge of the high-frequency asymptotics reduces the statistical uncertainties of the vertex and further eliminates finite-box-size effects. The proposed method benefits the calculation of nonlocal susceptibilities in DMFT and diagrammatic extensions of DMFT.
We present a worm sampling method for calculating one- and two-particle Green's functions using continuous-time quantum Monte Carlo simulations in the hybridization expansion (CT-HYB). Instead of measuring Green's functions by removing hybridization lines from partition function configurations, as in conventional CT-HYB, the worm algorithm directly samples the Green's function. We show that worm sampling is necessary to obtain general two-particle Green's functions which are not of density-density type and that it improves the sampling efficiency when approaching the atomic limit. Such two-particle Green's functions are needed to compute off-diagonal elements of susceptibilities and occur in diagrammatic extensions of the dynamical mean-field theory and in efficient estimators for the single-particle self-energy.
We propose efficient measurement procedures for the self-energy and vertex function of the Anderson impurity model within the hybridization expansion continuous-time quantum Monte Carlo algorithm. The method is based on the measurement of higher-order correlation functions related to the quantities being sought through the equation of motion, a technique previously introduced in the numerical renormalization-group context. For the case of interactions of density-density type, the additional correlators can be obtained at essentially no additional computational cost. In combination with a recently introduced method for filtering the Monte Carlo noise using a representation in terms of orthogonal polynomials, we obtain data with unprecedented accuracy. This leads to an enhanced stability in analytical continuations of the self-energy or in two-particle-based theories such as the dual fermion approach. As an illustration of the method we reexamine the previously reported spin-freezing and high-spin to low-spin transitions in a two-orbital model with density-density interactions. In both cases, the vertex function undergoes significant changes, which suggests significant corrections to the dynamical mean-field solutions in dual fermion calculations.
We derive equations of motion for Green's functions of the multiorbital Anderson impurity model by differentiating symmetrically with respect to all time arguments. The resulting equations relate the one- and two-particle Green's function to correlators of up to six particles at four times. As an application we consider continuous-time quantum Monte Carlo simulations in the hybridization expansion, which hitherto suffered from notoriously high noise levels at large Matsubara frequencies. Employing the derived symmetric improved estimators overcomes this problem.
Strong electronic correlations are the root of many fascinating phenomena of condensed matter physics. However, they also make calculations notoriously difficult. This thesis is concerned with the development and improvement of computational methods for strongly correlated electron systems. After a brief recapitulation of the Green’s function formalism,which is the theoretical basis of most such methods, we discuss the analytical continuation from imaginary to real frequencies. Especially the maximum entropy method is explained in detail and an implementation there of is presented.The remainder of this thesis deals with methods to compute Green’s functions on thei maginary frequency axis. First, the Anderson impurity model (AIM) is considered. It can be solved numerically – a task for which several methods exist. We focus on continuous-time quantum Monte Carlo with worm sampling, through which various kinds of correlation functions can be calculated. On this basis we derive and implement symmetric improved estimators that allow for the computation of the self-energy and vertex functions where the level of stochastic noise is largely reduced with respect to conventional estimators.As a next step the dynamical mean-field theory (DMFT) is introduced as a method to capture local correlations in the Hubbard model (HM) by self-consistently mapping it onto the AIM. We study the convergence behavior of this iterative procedure and the improvement of convergence due to mixing. Then, we compare electron scattering rates from DMFT in the weak-coupling regime of the square-lattice HM to results from Boltzmann scattering theory.While local correlations account for many interesting phenomena, such as the Mott-Hubbard metal-insulator transition, they are not sufficient to correctly describe, e.g., the physics of phase transitions in two- or three-dimensional systems. It is thus the final step in this thesis to also include nonlocal correlations by means of the dynamical vertex approximation (DΓA). For this we develop a new self-consistency scheme, which is then applied to several systems: We study the square-lattice HM, where we compare to benchmark methods in the weak-coupling regime, but also present results for intermediate and strong coupling. Furthermore, a two-band model with doping and a three-band Wannier model for strontium vanadate are studied. As an application of our new self-consistent DΓA, we study the effects of geometrical frustration in the kagome lattice, where we also compare our results to determinant quantum Monte Carlo.
A numerically exact continuous-time quantum Monte Carlo algorithm for finite fermionic systems with nonlocal interactions is proposed. The scheme is particularly applicable for general multiband time-dependent correlations, since it does not invoke Hubbard-Stratonovich transformation. The present determinantal grand-canonical method is based on a stochastic series expansion for the partition function in the interaction representation. The results for the Green function and for the time-dependent susceptibility of multiorbital supersymmetric impurity model with a spin-flip interaction are presented.