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:1205.2791 · Chinese Physics B

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

Li Huang

Published Nov 1, 2016• 4 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1088/1674-1056/25/11/117101

Related papers

Orthogonal polynomial representation of imaginary-time Green’s functions

Lewin Boehnke, Hartmut Hafermann, Michel Ferrero +2 more

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.

arXiv:1104.3215•Physical Review B•Aug 12, 2011•215 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

Hybrid quantum-classical algorithm for computing imaginary-time correlation functions

Rihito Sakurai, Wataru Mizukami, Hiroshi Shinaoka

Quantitative descriptions of strongly correlated materials pose a considerable challenge in condensed matter physics and chemistry. A promising approach to address this problem is quantum embedding methods. In particular, the dynamical mean-field theory (DMFT) maps the original system to an effective quantum impurity model comprising correlated orbitals embedded in an electron bath. The biggest bottleneck in DMFT calculations is numerically solving the quantum impurity model, i.e., computing the Green's function. Past studies have proposed theoretical methods to compute the Green's function of a quantum impurity model in polynomial time using a quantum computer. So far, however, efficient methods for computing the imaginary-time Green's functions have not been established despite the advantages of the imaginary-time formulation. We propose a quantum-classical hybrid algorithm for computing imaginary-time Green's functions on quantum devices with limited hardware resources by applying the variational quantum simulation. Using a quantum circuit simulator, we verified this algorithm by computing Green's functions for a dimer model as well as a four-site impurity model obtained by DMFT calculations of the single-band Hubbard model, although our method can be applied to general imaginary-time correlation functions.

arXiv:2112.02764•Physical Review Research•Jun 16, 2022•30 citations

Improved Hubbard-I approximation impurity solver for quantum impurity models

Huanhuan Qiu, Jianing Zhuang, Li Huang +2 more

We develop a fast impurity solver which is based on the combination of Hubbard-I approximation and hybridization expansion continuous-time quantum Monte Carlo algorithm. This solver inherits the advantages of both algorithms. In order to demonstrate the power and usefulness of this solver, we use it to study the magnetic phase transitions of single-band and two-band Hubbard models in the framework of single-site dynamical mean-field theory. The calculated results agree well with those obtained by hybridization expansion quantum impurity solver. It is suggested that this solver is very suitable to solve the multi-orbital quantum impurity models efficiently and can be used to study more realistic systems with magnetic long-range order in the future.

Journal of Physics Condensed Matter•Nov 6, 2018•0 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

Fast and Accurate Scheme for Green Functions and Eigenvectors: Regulated Polynomial Expansion without Gibbs Oscillation

Shigetoshi Sota, Masaki Itoh

An oscillation-free Fourier expansion is proposed and applied to the polynomial expansion of resolvent operator, as an alternative treatment for Gibbs phenomenon to the kernel polynomial method (KPM) by Voter et al. [Phys. Rev. B 53 (1996) 12733]. Adopting regulated Legendre polynomial as basis functions, the recurrence algorithm underlying the KPM is practically preserved, while it eliminates unphysical oscillations entirely. The resolution is uniform in the frequency domain, and improved without limitation by extending the truncation range of the expansion. The method can also be used to calculate the eigenvector with remarkably high precision. It is further extended to the time domain, being able to deal with the quantum evolution in the same algorithm. A numerical test was performed against a lattice dynamics of 19 3 atoms, as well as of a computer generated Lennard–Jones glass of 1000 atoms. Not only it serves as a most efficient simulation scheme for a bulk, confirmed was also that the calculation converges to exact results within realistic computation time.

Journal of the Physical Society of Japan•May 11, 2007•15 citations
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.
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.
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.
Quantitative descriptions of strongly correlated materials pose a considerable challenge in condensed matter physics and chemistry. A promising approach to address this problem is quantum embedding methods. In particular, the dynamical mean-field theory (DMFT) maps the original system to an effective quantum impurity model comprising correlated orbitals embedded in an electron bath. The biggest bottleneck in DMFT calculations is numerically solving the quantum impurity model, i.e., computing the Green's function. Past studies have proposed theoretical methods to compute the Green's function of a quantum impurity model in polynomial time using a quantum computer. So far, however, efficient methods for computing the imaginary-time Green's functions have not been established despite the advantages of the imaginary-time formulation. We propose a quantum-classical hybrid algorithm for computing imaginary-time Green's functions on quantum devices with limited hardware resources by applying the variational quantum simulation. Using a quantum circuit simulator, we verified this algorithm by computing Green's functions for a dimer model as well as a four-site impurity model obtained by DMFT calculations of the single-band Hubbard model, although our method can be applied to general imaginary-time correlation functions.
We develop a fast impurity solver which is based on the combination of Hubbard-I approximation and hybridization expansion continuous-time quantum Monte Carlo algorithm. This solver inherits the advantages of both algorithms. In order to demonstrate the power and usefulness of this solver, we use it to study the magnetic phase transitions of single-band and two-band Hubbard models in the framework of single-site dynamical mean-field theory. The calculated results agree well with those obtained by hybridization expansion quantum impurity solver. It is suggested that this solver is very suitable to solve the multi-orbital quantum impurity models efficiently and can be used to study more realistic systems with magnetic long-range order in the future.
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.
An oscillation-free Fourier expansion is proposed and applied to the polynomial expansion of resolvent operator, as an alternative treatment for Gibbs phenomenon to the kernel polynomial method (KPM) by Voter et al. [Phys. Rev. B 53 (1996) 12733]. Adopting regulated Legendre polynomial as basis functions, the recurrence algorithm underlying the KPM is practically preserved, while it eliminates unphysical oscillations entirely. The resolution is uniform in the frequency domain, and improved without limitation by extending the truncation range of the expansion. The method can also be used to calculate the eigenvector with remarkably high precision. It is further extended to the time domain, being able to deal with the quantum evolution in the same algorithm. A numerical test was performed against a lattice dynamics of 19 3 atoms, as well as of a computer generated Lennard–Jones glass of 1000 atoms. Not only it serves as a most efficient simulation scheme for a bulk, confirmed was also that the calculation converges to exact results within realistic computation time.