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

Performance analysis of a physically constructed orthogonal representation of imaginary-time Green's function

Naoya Chikano, Junya Otsuki, Hiroshi Shinaoka

Published Jul 3, 2018• 41 citations
arXivHTMLPDF

Abstract

The imaginary-time Green's function is a building block of various numerical methods for correlated electron systems. Recently, it was shown that a model-independent compact orthogonal representation of the Green's function can be constructed by decomposing its spectral representation. We investigate the performance of this so-called intermediate representation (IR) from several points of view. First, we develop an efficient algorithm for computing the IR basis functions of arbitrary high degree. Second, for two simple models, we study the number of coefficients required to represent the Green's function within a given tolerance. We show that the number of coefficients grows only as $O(ln\ensuremath{\beta})$ for fermions and converges to a constant for bosons as temperature $T=1/\ensuremath{\beta}$ decreases. Third, we show that this remarkable feature is ascribed to the properties of the physically constructed basis functions. The fermionic basis functions on the real-frequency axis have features whose width is scaled as $O(T)$, which is consistent with the low-$T$ properties of quasiparticles in a Fermi liquid state. On the other hand, the properties of the bosonic basis functions are consistent with those of spin/orbital susceptibilities at low $T$. These results demonstrate the potential wide applications of the IR to calculations of correlated systems.

DOI
10.1103/physrevb.98.035104

Related papers

Discrete Lehmann representation of imaginary time Green's functions

Jason Kaye, Kun Chen, Olivier Parcollet

We present an efficient basis for imaginary time Green's functions based on a low-rank decomposition of the spectral Lehmann representation. The basis functions are simply a set of well-chosen exponentials, so the corresponding expansion may be thought of as a discrete form of the Lehmann representation using an effective spectral density which is a sum of $\ensuremath{\delta}$ functions. The basis is determined only by an upper bound on the product $\ensuremath{\beta}{\ensuremath{\omega}}_{max}$, with $\ensuremath{\beta}$ the inverse temperature and ${\ensuremath{\omega}}_{max}$ an energy cutoff, and a user-defined error tolerance $\ensuremath{\epsilon}$. The number $r$ of basis functions scales as $O(log(\ensuremath{\beta}{\ensuremath{\omega}}_{max})log(1/\ensuremath{\epsilon}))$. The discrete Lehmann representation of a particular imaginary time Green's function can be recovered by interpolation at a set of $r$ imaginary time nodes. Both the basis functions and the interpolation nodes can be obtained rapidly using standard numerical linear algebra routines. Due to the simple form of the basis, the discrete Lehmann representation of a Green's function can be explicitly transformed to the Matsubara frequency domain, or obtained directly by interpolation on a Matsubara frequency grid. We benchmark the efficiency of the representation on simple cases, and with a high-precision solution of the Sachdev-Ye-Kitaev equation at low temperature. We compare our approach with the related intermediate representation method, and introduce an improved algorithm to build the intermediate representation basis and a corresponding sampling grid.

arXiv:2107.13094•Physical review. B./Physical review. B•Jun 13, 2022•61 citations

Compressing Green's function using intermediate representation between imaginary-time and real-frequency domains

Hiroshi Shinaoka, Junya Otsuki, Masayuki Ohzeki +1 more

Model-independent compact representations of imaginary-time data are presented in terms of the intermediate representation (IR) of analytical continuation. We demonstrate the efficiency of the IR through continuous-time quantum Monte Carlo calculations of an Anderson impurity model. We find that the IR yields a significantly compact form of various types of correlation functions. This allows the direct quantum Monte Carlo measurement of Green's functions in a compressed form, which considerably reduces the computational cost and memory usage. Furthermore, the present framework will provide general ways to boost the power of cutting-edge diagrammatic/quantum Monte Carlo treatments of many-body systems.

arXiv:1702.03054•Physical review. B./Physical review. B•Jul 25, 2017•171 citations

Analytic Origin of Green-Function Compression in the Intermediate Representation

Takahiro Misawa

Information compression plays a central role in diverse fields of modern science and technology, from communication theory to machine learning. In condensed-matter physics, the intermediate representation (IR) basis has recently been developed as an efficient method for compressing imaginary-time Green functions, which are fundamental quantities for describing quantum many-body systems. This compression relies on the rapid decay of the singular values with the basis index and the unusually weak growth of the effective rank with inverse temperature. Because of these useful features, the IR basis is now widely used as a standard method in quantum many-body calculations. However, the analytic origin of its compression capability has remained unclear. Here we uncover a finite-Laplace-transform structure underlying the IR kernel, which reveals that the eigenfunctions of the IR kernel admit a natural expansion in terms of classical special functions, the oblate spheroidal wave functions. This finite-Laplace-transform structure also enables us to analytically clarify the compression mechanism of the IR basis. Our results provide a mathematical foundation for the compression of imaginary-time Green functions, connecting quantum many-body physics with theories of information compression and finite integral transforms.

arXiv:2605.24814•May 24, 2026•0 citations

Decomposing Imaginary-Time Feynman Diagrams Using Separable Basis Functions: Anderson Impurity Model Strong-Coupling Expansion

Jason Kaye, Zhen Huang, Hugo U. R. Strand +1 more

We present a deterministic algorithm for the efficient evaluation of imaginary-time diagrams based on the recently introduced discrete Lehmann representation (DLR) of imaginary-time Green’s functions. In addition to the efficient discretization of diagrammatic integrals afforded by its approximation properties, the DLR basis is separable in imaginary-time, allowing us to decompose diagrams into linear combinations of nested sequences of one-dimensional products and convolutions. Focusing on the strong-coupling bold-line expansion of generalized Anderson impurity models, we show that our strategy reduces the computational complexity of evaluating an M th-order diagram at inverse temperature β and spectral width ω max from O ( ( β ω max ) 2 M − 1 ) for a direct quadrature to O ( M ( log ( β ω max ) ) M + 1 ) , with controllable high-order accuracy. We benchmark our algorithm using third-order expansions for multiband impurity problems with off-diagonal hybridization and spin-orbit coupling, presenting comparisons with exact diagonalization and quantum Monte Carlo approaches. In particular, we perform a self-consistent dynamical mean-field theory calculation for a three-band Hubbard model with strong spin-orbit coupling representing a minimal model of Ca 2 RuO 4 , demonstrating the promise of the method for modeling realistic strongly correlated multiband materials. For both strong and weak coupling expansions of low and intermediate order, in which diagrams can be enumerated, our method provides an efficient, straightforward, and robust blackbox evaluation procedure. In this sense, it fills a gap between diagrammatic approximations of the lowest order, which are simple and inexpensive but inaccurate, and those based on Monte Carlo sampling of high-order diagrams. Published by the American Physical Society 2024

arXiv:2307.08566•Physical Review X•Aug 26, 2024•12 citations

Efficient ab initio many-body calculations based on sparse modeling of Matsubara Green's function

Hiroshi Shinaoka, Naoya Chikano, Emanuel Gull +6 more

This lecture note reviews recently proposed sparse-modeling approaches for efficient ab initio many-body calculations based on the data compression of Green's functions. The sparse-modeling techniques are based on a compact orthogonal basis, an intermediate representation (IR) basis, for imaginary-time and Matsubara Green's functions. A sparse sampling method based on the IR basis enables solving diagrammatic equations efficiently. We describe the basic properties of the IR basis, the sparse sampling method and its applications to ab initio calculations based on the GW approximation and the Migdal--Eliashberg theory. We also describe a numerical library for the IR basis and the sparse sampling method, sparse-ir, and provide its sample codes. This lecture note follows the Japanese review article [H. Shinaoka et al., Solid State Physics 56(6), 301 (2021)].

arXiv:2106.12685•SciPost Physics Lecture Notes•Sep 22, 2022•30 citations

Overcomplete compact representation of two-particle Green's functions

Hiroshi Shinaoka, Junya Otsuki, Kristjan Haule +4 more

Two-particle Green's functions and the vertex functions play a critical role in theoretical frameworks for describing strongly correlated electron systems. However, numerical calculations at the two-particle level often suffer from large computation time and massive memory consumption. We derive a general expansion formula for the two-particle Green's functions in terms of an overcomplete representation based on the recently proposed ``intermediate representation'' basis. The expansion formula is obtained by decomposing the spectral representation of the two-particle Green's function. We demonstrate that the expansion coefficients decay exponentially, while all high-frequency and long-tail structures in the Matsubara-frequency domain are retained. This representation therefore enables efficient treatment of two-particle quantities and opens a route to the application of modern many-body theories to realistic strongly correlated electron systems.

arXiv:1803.01916•Physical review. B./Physical review. B•May 9, 2018•43 citations
The imaginary-time Green's function is a building block of various numerical methods for correlated electron systems. Recently, it was shown that a model-independent compact orthogonal representation of the Green's function can be constructed by decomposing its spectral representation. We investigate the performance of this so-called intermediate representation (IR) from several points of view. First, we develop an efficient algorithm for computing the IR basis functions of arbitrary high degree. Second, for two simple models, we study the number of coefficients required to represent the Green's function within a given tolerance. We show that the number of coefficients grows only as O(ln\ensuremathβ)O(ln\ensuremath{\beta})O(ln\ensuremathβ) for fermions and converges to a constant for bosons as temperature T=1/\ensuremathβT=1/\ensuremath{\beta}T=1/\ensuremathβ decreases. Third, we show that this remarkable feature is ascribed to the properties of the physically constructed basis functions. The fermionic basis functions on the real-frequency axis have features whose width is scaled as O(T)O(T)O(T), which is consistent with the low-TTT properties of quasiparticles in a Fermi liquid state. On the other hand, the properties of the bosonic basis functions are consistent with those of spin/orbital susceptibilities at low TTT. These results demonstrate the potential wide applications of the IR to calculations of correlated systems.
We present an efficient basis for imaginary time Green's functions based on a low-rank decomposition of the spectral Lehmann representation. The basis functions are simply a set of well-chosen exponentials, so the corresponding expansion may be thought of as a discrete form of the Lehmann representation using an effective spectral density which is a sum of $\ensuremath{\delta}$ functions. The basis is determined only by an upper bound on the product $\ensuremath{\beta}{\ensuremath{\omega}}_{max}$, with $\ensuremath{\beta}$ the inverse temperature and ${\ensuremath{\omega}}_{max}$ an energy cutoff, and a user-defined error tolerance $\ensuremath{\epsilon}$. The number rrr of basis functions scales as $O(log(\ensuremath{\beta}{\ensuremath{\omega}}_{max})log(1/\ensuremath{\epsilon}))$. The discrete Lehmann representation of a particular imaginary time Green's function can be recovered by interpolation at a set of rrr imaginary time nodes. Both the basis functions and the interpolation nodes can be obtained rapidly using standard numerical linear algebra routines. Due to the simple form of the basis, the discrete Lehmann representation of a Green's function can be explicitly transformed to the Matsubara frequency domain, or obtained directly by interpolation on a Matsubara frequency grid. We benchmark the efficiency of the representation on simple cases, and with a high-precision solution of the Sachdev-Ye-Kitaev equation at low temperature. We compare our approach with the related intermediate representation method, and introduce an improved algorithm to build the intermediate representation basis and a corresponding sampling grid.
Model-independent compact representations of imaginary-time data are presented in terms of the intermediate representation (IR) of analytical continuation. We demonstrate the efficiency of the IR through continuous-time quantum Monte Carlo calculations of an Anderson impurity model. We find that the IR yields a significantly compact form of various types of correlation functions. This allows the direct quantum Monte Carlo measurement of Green's functions in a compressed form, which considerably reduces the computational cost and memory usage. Furthermore, the present framework will provide general ways to boost the power of cutting-edge diagrammatic/quantum Monte Carlo treatments of many-body systems.
Information compression plays a central role in diverse fields of modern science and technology, from communication theory to machine learning. In condensed-matter physics, the intermediate representation (IR) basis has recently been developed as an efficient method for compressing imaginary-time Green functions, which are fundamental quantities for describing quantum many-body systems. This compression relies on the rapid decay of the singular values with the basis index and the unusually weak growth of the effective rank with inverse temperature. Because of these useful features, the IR basis is now widely used as a standard method in quantum many-body calculations. However, the analytic origin of its compression capability has remained unclear. Here we uncover a finite-Laplace-transform structure underlying the IR kernel, which reveals that the eigenfunctions of the IR kernel admit a natural expansion in terms of classical special functions, the oblate spheroidal wave functions. This finite-Laplace-transform structure also enables us to analytically clarify the compression mechanism of the IR basis. Our results provide a mathematical foundation for the compression of imaginary-time Green functions, connecting quantum many-body physics with theories of information compression and finite integral transforms.
We present a deterministic algorithm for the efficient evaluation of imaginary-time diagrams based on the recently introduced discrete Lehmann representation (DLR) of imaginary-time Green’s functions. In addition to the efficient discretization of diagrammatic integrals afforded by its approximation properties, the DLR basis is separable in imaginary-time, allowing us to decompose diagrams into linear combinations of nested sequences of one-dimensional products and convolutions. Focusing on the strong-coupling bold-line expansion of generalized Anderson impurity models, we show that our strategy reduces the computational complexity of evaluating an M th-order diagram at inverse temperature β and spectral width ω max from O ( ( β ω max ) 2 M − 1 ) for a direct quadrature to O ( M ( log ( β ω max ) ) M + 1 ) , with controllable high-order accuracy. We benchmark our algorithm using third-order expansions for multiband impurity problems with off-diagonal hybridization and spin-orbit coupling, presenting comparisons with exact diagonalization and quantum Monte Carlo approaches. In particular, we perform a self-consistent dynamical mean-field theory calculation for a three-band Hubbard model with strong spin-orbit coupling representing a minimal model of Ca 2 RuO 4 , demonstrating the promise of the method for modeling realistic strongly correlated multiband materials. For both strong and weak coupling expansions of low and intermediate order, in which diagrams can be enumerated, our method provides an efficient, straightforward, and robust blackbox evaluation procedure. In this sense, it fills a gap between diagrammatic approximations of the lowest order, which are simple and inexpensive but inaccurate, and those based on Monte Carlo sampling of high-order diagrams. Published by the American Physical Society 2024
This lecture note reviews recently proposed sparse-modeling approaches for efficient ab initio many-body calculations based on the data compression of Green's functions. The sparse-modeling techniques are based on a compact orthogonal basis, an intermediate representation (IR) basis, for imaginary-time and Matsubara Green's functions. A sparse sampling method based on the IR basis enables solving diagrammatic equations efficiently. We describe the basic properties of the IR basis, the sparse sampling method and its applications to ab initio calculations based on the GW approximation and the Migdal--Eliashberg theory. We also describe a numerical library for the IR basis and the sparse sampling method, sparse-ir, and provide its sample codes. This lecture note follows the Japanese review article [H. Shinaoka et al., Solid State Physics 56(6), 301 (2021)].
Two-particle Green's functions and the vertex functions play a critical role in theoretical frameworks for describing strongly correlated electron systems. However, numerical calculations at the two-particle level often suffer from large computation time and massive memory consumption. We derive a general expansion formula for the two-particle Green's functions in terms of an overcomplete representation based on the recently proposed ``intermediate representation'' basis. The expansion formula is obtained by decomposing the spectral representation of the two-particle Green's function. We demonstrate that the expansion coefficients decay exponentially, while all high-frequency and long-tail structures in the Matsubara-frequency domain are retained. This representation therefore enables efficient treatment of two-particle quantities and opens a route to the application of modern many-body theories to realistic strongly correlated electron systems.