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

Orthogonal polynomial representation of imaginary-time Green’s functions

Lewin Boehnke, Hartmut Hafermann, Michel Ferrero, Frank Lechermann, Olivier Parcollet

Published Aug 12, 2011• 215 citations
arXivHTMLPDF

Abstract

We study the expansion of single-particle and two-particle imaginary-time Matsubara Green's functions of quantum impurity models in the basis of Legendre orthogonal polynomials. We discuss various applications within the dynamical mean-field theory (DMFT) framework. The method provides a more compact representation of the Green's functions than standard Matsubara frequencies, and therefore significantly reduces the memory-storage size of these quantities. Moreover, it can be used as an efficient noise filter for various physical quantities within the continuous-time quantum Monte Carlo impurity solvers recently developed for DMFT and its extensions. In particular, we show how to use it for the computation of energies in the context of realistic DMFT calculations in combination with the local density approximation to the density functional theory (LDA+DMFT) and for the calculation of lattice susceptibilities from the local irreducible vertex function.

DOI
10.1103/physrevb.84.075145

Related papers

Kernel polynomial representation for imaginary-time Green’s functions in continuous-time quantum Monte Carlo impurity solver

Li Huang

Inspired by the recent proposed Legendre orthogonal polynomial representation of imaginary-time Green's functions, we develop an alternate representation for the Green's functions of quantum impurity models and combine it with the hybridization expansion continuous-time quantum Monte Carlo impurity solver. This representation is based on the kernel polynomial method, which introduces various integral kernels to filter fluctuations caused by the explicit truncations of polynomial expansion series and improve the computational precision significantly. As an illustration of the new representation, we reexamine the imaginary-time Green's functions of single-band Hubbard model in the framework of dynamical mean-field theory. The calculated results suggest that with carefully chosen integral kernels the Gibbs oscillations found in previous orthogonal polynomial representation have been suppressed vastly and remarkable corrections to the measured Green's functions have been obtained.

arXiv:1205.2791•Chinese Physics B•Nov 1, 2016•4 citations

Low-rank Green's function representations applied to dynamical mean-field theory

Nan Sheng, Alexander Hampel, Sophie Beck +4 more

Several recent works have introduced highly compact representations of single-particle Green's functions in the imaginary time and Matsubara frequency domains, as well as efficient interpolation grids used to recover the representations. In particular, the intermediate representation with sparse sampling and the discrete Lehmann representation (DLR) make use of low rank compression techniques to obtain optimal approximations with controllable accuracy. We consider the use of the DLR in dynamical mean-field theory (DMFT) calculations, and in particular show that the standard full Matsubara frequency grid can be replaced by the compact grid of DLR Matsubara frequency nodes. We test the performance of the method for a DMFT calculation of ${\mathrm{Sr}}_{2}{\mathrm{RuO}}_{4}$ at temperature $50\phantom{\rule{0.16em}{0ex}}\mathrm{K}$ using a continuous-time quantum Monte Carlo impurity solver, and demonstrate that Matsubara frequency quantities can be represented on a grid of only 36 nodes with no reduction in accuracy, or increase in the number of self-consistent iterations, despite the presence of significant Monte Carlo noise.

arXiv:2301.07764•Physical review. B./Physical review. B•Jun 13, 2023•8 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

Efficient Temperature-Dependent Green’s Functions Methods for Realistic Systems: Compact Grids for Orthogonal Polynomial Transforms

Alexei A. Kananenka, Jordan J. Phillips, Dominika Zgid

The Matsubara Green's function that is used to describe temperature-dependent behavior is expressed on a numerical grid. While such a grid usually has a couple of hundred points for low-energy model systems, for realistic systems with large basis sets the size of an accurate grid can be tens of thousands of points, constituting a severe computational and memory bottleneck. In this paper, we determine efficient imaginary time grids for the temperature-dependent Matsubara Green's function formalism that can be used for calculations on realistic systems. We show that, because of the use of an orthogonal polynomial transform, we can restrict the imaginary time grid to a few hundred points and reach micro-Hartree accuracy in the electronic energy evaluation. Moreover, we show that only a limited number of orthogonal polynomial expansion coefficients are necessary to preserve accuracy when working with a dual representation of the Green's function or self-energy and transforming between the imaginary time and frequency domain.

arXiv:1509.04262•Journal of Chemical Theory and Computation•Jan 6, 2016•49 citations

Exact diagonalization as an impurity solver in dynamical mean field theory

Yi Lu, M. W. Haverkort

The dynamical mean-field theory (DMFT) maps a correlated lattice problem onto an impurity problem of a single correlated site coupled to an uncorrelated bath. Most implementations solve the DMFT equations using quantum Monte-Carlo sampling on the imaginary time and frequency (Matsubara) axis. We will here review alternative methods using exact diagonalization, i.e., representing the many-body ground state of the impurity as a sum over Slater determinants and calculating Green’s functions using iterative Lanczos procedures. The advantage being that these methods have no sign problem, can handle involved multi-orbital Hamiltonians (low crystal symmetry, spin-orbit coupling) and – when working completely on the real axis – do not need a mathematically ill-posed analytical continuation. The disadvantage of traditional implementations of exact diagonalization has been the exponential scaling of the calculation problem as a function of number of bath discretization points. In the last part we will review how recent advances in exact diagonalization can evade the exponential barrier thereby increasing the number of bath discretization points to reach the thermodynamic limit.

The European Physical Journal Special Topics•Jul 1, 2017•17 citations

Sum-over-poles framework for quantum embedding in dynamical mean-field theory

Carbone, Alessandro

This thesis develops a theoretical framework for quantum embedding based on compact sum-over-poles (SOP) representations of Green's functions, with a particular focus on dynamical mean-field theory (DMFT) formulated close to the real-frequency axis. The work starts from the algorithmic inversion method (AIM-SOP) and generalizes it to many-body lattice systems, proving that any interacting fragment embedded in an arbitrary bath can be represented by an auxiliary Hamiltonian with a finite set of non-interacting bath sites. The author embeds this formalism into a DMFT self-consistency loop formulated on complex contours close to the real axis, complemented by a SOP-based fitting strategy and exact-diagonalization (ED) impurity solvers, constructing a DMFT scheme that does not rely on approximate analytic continuation techniques. In particular, the framework enables us to work on compact (SOP) representations of Green's functions and self-energies alone, and leads to a frequency-independent formulation of DMFT, grounded in effective free energy (Baym-Kadanoff and Klein) functionals. Applications to a $1000$-site one-dimensional Hubbard ring demonstrate that compact SOP forms reproduce state-of-the-art ED-DMFT results on the Matsubara frequencies and retain key spectral features as the contour approaches the real axis, while also revealing convergence trends and multiple-solution behavior near challenging correlation regimes. Beyond its numerical validation, this approach paves the way for novel variational implementations of DMFT and realistic applications to materials, harnessing its promising ability to resolve fine spectral features.

•Jan 1, 2026•0 citations
We study the expansion of single-particle and two-particle imaginary-time Matsubara Green's functions of quantum impurity models in the basis of Legendre orthogonal polynomials. We discuss various applications within the dynamical mean-field theory (DMFT) framework. The method provides a more compact representation of the Green's functions than standard Matsubara frequencies, and therefore significantly reduces the memory-storage size of these quantities. Moreover, it can be used as an efficient noise filter for various physical quantities within the continuous-time quantum Monte Carlo impurity solvers recently developed for DMFT and its extensions. In particular, we show how to use it for the computation of energies in the context of realistic DMFT calculations in combination with the local density approximation to the density functional theory (LDA+DMFT) and for the calculation of lattice susceptibilities from the local irreducible vertex function.
Inspired by the recent proposed Legendre orthogonal polynomial representation of imaginary-time Green's functions, we develop an alternate representation for the Green's functions of quantum impurity models and combine it with the hybridization expansion continuous-time quantum Monte Carlo impurity solver. This representation is based on the kernel polynomial method, which introduces various integral kernels to filter fluctuations caused by the explicit truncations of polynomial expansion series and improve the computational precision significantly. As an illustration of the new representation, we reexamine the imaginary-time Green's functions of single-band Hubbard model in the framework of dynamical mean-field theory. The calculated results suggest that with carefully chosen integral kernels the Gibbs oscillations found in previous orthogonal polynomial representation have been suppressed vastly and remarkable corrections to the measured Green's functions have been obtained.
Several recent works have introduced highly compact representations of single-particle Green's functions in the imaginary time and Matsubara frequency domains, as well as efficient interpolation grids used to recover the representations. In particular, the intermediate representation with sparse sampling and the discrete Lehmann representation (DLR) make use of low rank compression techniques to obtain optimal approximations with controllable accuracy. We consider the use of the DLR in dynamical mean-field theory (DMFT) calculations, and in particular show that the standard full Matsubara frequency grid can be replaced by the compact grid of DLR Matsubara frequency nodes. We test the performance of the method for a DMFT calculation of Sr2RuO4{\mathrm{Sr}}_{2}{\mathrm{RuO}}_{4}Sr2​RuO4​ at temperature 50K50\phantom{\rule{0.16em}{0ex}}\mathrm{K}50K using a continuous-time quantum Monte Carlo impurity solver, and demonstrate that Matsubara frequency quantities can be represented on a grid of only 36 nodes with no reduction in accuracy, or increase in the number of self-consistent iterations, despite the presence of significant Monte Carlo noise.
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.
The Matsubara Green's function that is used to describe temperature-dependent behavior is expressed on a numerical grid. While such a grid usually has a couple of hundred points for low-energy model systems, for realistic systems with large basis sets the size of an accurate grid can be tens of thousands of points, constituting a severe computational and memory bottleneck. In this paper, we determine efficient imaginary time grids for the temperature-dependent Matsubara Green's function formalism that can be used for calculations on realistic systems. We show that, because of the use of an orthogonal polynomial transform, we can restrict the imaginary time grid to a few hundred points and reach micro-Hartree accuracy in the electronic energy evaluation. Moreover, we show that only a limited number of orthogonal polynomial expansion coefficients are necessary to preserve accuracy when working with a dual representation of the Green's function or self-energy and transforming between the imaginary time and frequency domain.
The dynamical mean-field theory (DMFT) maps a correlated lattice problem onto an impurity problem of a single correlated site coupled to an uncorrelated bath. Most implementations solve the DMFT equations using quantum Monte-Carlo sampling on the imaginary time and frequency (Matsubara) axis. We will here review alternative methods using exact diagonalization, i.e., representing the many-body ground state of the impurity as a sum over Slater determinants and calculating Green’s functions using iterative Lanczos procedures. The advantage being that these methods have no sign problem, can handle involved multi-orbital Hamiltonians (low crystal symmetry, spin-orbit coupling) and – when working completely on the real axis – do not need a mathematically ill-posed analytical continuation. The disadvantage of traditional implementations of exact diagonalization has been the exponential scaling of the calculation problem as a function of number of bath discretization points. In the last part we will review how recent advances in exact diagonalization can evade the exponential barrier thereby increasing the number of bath discretization points to reach the thermodynamic limit.
This thesis develops a theoretical framework for quantum embedding based on compact sum-over-poles (SOP) representations of Green's functions, with a particular focus on dynamical mean-field theory (DMFT) formulated close to the real-frequency axis. The work starts from the algorithmic inversion method (AIM-SOP) and generalizes it to many-body lattice systems, proving that any interacting fragment embedded in an arbitrary bath can be represented by an auxiliary Hamiltonian with a finite set of non-interacting bath sites. The author embeds this formalism into a DMFT self-consistency loop formulated on complex contours close to the real axis, complemented by a SOP-based fitting strategy and exact-diagonalization (ED) impurity solvers, constructing a DMFT scheme that does not rely on approximate analytic continuation techniques. In particular, the framework enables us to work on compact (SOP) representations of Green's functions and self-energies alone, and leads to a frequency-independent formulation of DMFT, grounded in effective free energy (Baym-Kadanoff and Klein) functionals. Applications to a 100010001000-site one-dimensional Hubbard ring demonstrate that compact SOP forms reproduce state-of-the-art ED-DMFT results on the Matsubara frequencies and retain key spectral features as the contour approaches the real axis, while also revealing convergence trends and multiple-solution behavior near challenging correlation regimes. Beyond its numerical validation, this approach paves the way for novel variational implementations of DMFT and realistic applications to materials, harnessing its promising ability to resolve fine spectral features.