Diagrammatic extensions of dynamical mean field theory (DMFT) such as the dynamical vertex approximation (DΓA) allow us to include non-local correlations beyond DMFT on all length scales and proved their worth for model calculations. Here, we develop and implement an AbinitioDΓA approach for electronic structure calculations of materials. Starting point is the two-particle irreducible vertex in the two particle-hole channels which is approximated by the bare non-local Coulomb interaction and all local vertex corrections. From this we calculate the full non-local vertex and the non-local self-energy through the Bethe-Salpeter equation. The AbinitioDΓA approach naturally generates all local DMFT correlations and all non-local GW contributions, but also further non-local correlations beyond: mixed terms of the former two and non-local spin fluctuations. We apply this new methodology to the prototypical correlated metal SrVO3. arXiv:1610.02998v2 [cond-mat.str-el]

We describe the hybridization-expansion continuous-time quantum Monte Carlo code package "w2dynamics", developed in Wien and Würzburg. We discuss the main features of this multi-orbital quantum impurity solver for the Anderson impurity model, dynamical mean field theory as well as its coupling to density functional theory. The w2dynamics package allows for calculating one-and two-particle quantities; it includes worm and further novel sampling schemes. Details about its download, installation, functioning and the relevant parameters are provided.

In this work, we analyze in detail the occurrence of divergences in the irreducible vertex functions for one of the fundamental models of many-body physics: the Anderson impurity model (AIM). These divergences, a surprising hallmark of the breakdown of many-electron perturbation theory -have been recently observed in several contexts, including the dynamical mean-field solution of the Hubbard model. The numerical calculations for the AIM presented in this work, as well as their comparison with the corresponding results for the Hubbard model, allow us to clarify several open questions about the properties of vertex divergences in a particularly interesting context, the correlated metallic regime at low-temperatures. Specifically, our analysis (i) rules out explicitly the transition to a Mott insulating phase, but not the more general suppression of charge fluctuations (proposed in [O. Gunnarsson et al., Phys. Rev. B 93, 245102 (2016)]), as a necessary condition for the occurrence of vertex divergences, (ii) clarifies their relation with the underlying Kondo physics, and, eventually, (iii) individuates which divergences might also appear on the real frequency axis in the limit of zero temperature, through the discovered scaling properties of the singular eigenvectors.

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 efficient estimators for the single-particle self-energy.

We derive the improved estimators for general interactions and employ these for the continuoustime 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 oneparticle 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 non-local 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 multi-orbital atomic-limit and the Falicov-Kimball model.