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:2602.16351

Liouvillian interpolation of the self-energy of cluster dynamical mean-field theories

Mathias Pelz, Jan von Delft, Andreas Gleis

Published Feb 18, 2026• 0 citations
arXivHTMLPDF

Abstract

Two widely-used non-local extensions of dynamical mean-field theory (DMFT), namely cellular DMFT (CDMFT) and the dynamical cluster approximation (DCA), both yield self-energies marred by having some unphysical properties: CDMFT yields real-space self-energies that are not translationally invariant, and DCA yields momentum-space self-energies with discontinuities in their momentum dependence. It is often desirable to remove these flaws by post-processing cluster DMFT results, using strategies called periodization for CDMFT and interpolation for DCA---for brevity, we refer to both cases as \textit{interpolation}. However, traditional interpolation approaches struggle to reliably capture intricate structures such as hole pockets in the hole-doped square-lattice Hubbard model. Further, these approaches interpolate frequency-dependent functions, which may lead to causality violations. Here, we propose \textit{Liouvillian interpolation}, a novel, intuitive, and robust scheme for interpolating cluster DMFT results. Our key idea is to interpolate frequency-\textit{independent} matrix elements of the single-particle irreducible part of the Liouvillian, obtained from a continued-fraction expansion of the cluster DMFT self-energy. We demonstrate that the ingredients of such an expansion possess a more local Fourier expansion than the functions involved in traditional interpolation schemes, and that Liouvillian interpolation conserves causality by construction. We illustrate our method for the one-dimensional Hubbard model using CDMFT, and for the two-dimensional Hubbard model using four-site CDMFT and four-patch DCA. For the latter, $\mathcal{L}$-interpolation can (depending on doping) yield Fermi and Luttinger arcs which together seem to form a closed surface.

DOI
10.48550/arxiv.2602.16351

Related papers

Cellular Dynamical Mean-field Theory in large Impurity Clusters

Marcel Klett

Even after decades of intense research, the single band Hubbard model representing the fundamental model for interacting quantum systems and electronic correlations remains an unsolved cornerstone paradigm in theoretical solid state physics. Within the non-perturbative dynamical mean-field theory (DMFT), the lattice problem is mapped to a self-consistent auxiliary quantum impurity model leading to a local approximation of the self-energy. Nonlocal correlations can be included by the real-space cluster extension of the DMFT (CDMFT). Enlarging the unit cell makes the calculation of a numerically exact solution in CDMFT very challenging. We here use a quantum Monte-Carlo approach in the imaginary-time space in order to solve the self-consistency equations for different cluster sizes. In order to restore the broken translational symmetry re-periodization schemes for the Green function, the self-energy, or its cumulants have been introduced. However, these suffer from ambiguity and may even lead to convergence problems when attempted inside the self-consistency loop. The comparison to numerically exact diagrammatic quantum Monte Carlo calculations shows that introducing a so-called centerfocused extrapolation (CFE) to approximate the lattice self-energy yields very accurate results. Moreover, the CFE converges faster with the cluster size than previous periodization schemes. We here perform a detailed CDMFT analysis for the single-band Hubbard model on the 2D square lattice, reaching real-space cluster sizes of up to 9x9 sites. In addition to spectral properties, we also compute two-particle correlation functions which are not accessible in DMFT. Using benchmarks against diagrammatic Monte Carlo at high temperature, we show that the cluster spin susceptibility can be extrapolated to the exact result at large cluster size. In particular, the exponential decay of spin-spin correlations is very well captured by CDMFT calculations, even when correlations extend beyond the size of the cluster. We further present results at lower temperature T and larger $U$ than the range currently accessible with diagrammatic Monte Carlo methods techniques. The CDMFT+CFE represents therefore a powerful computational tool to access the physics of non-local correlations beyond dynamical mean-field theory.

•Dec 29, 2020•0 citations

Dynamical cluster approximation with continuous lattice self-energy

Peter Staar, Thomas Maier, T. C. Schulthess

The dynamical cluster approximation (DCA) is a systematic extension beyond the single-site approximation in dynamical mean field theory, to include spatially nonlocal correlations in quantum many-body simulations of strongly correlated systems. We extend the DCA with a continuous lattice self-energy in order to achieve better convergence with cluster size. This method, which we call DCA${}^{+}$, cures the cluster-shape dependence problems of the DCA, without suffering from causality violations of previous attempts to interpolate the cluster self-energy. A practical approach based on standard inference techniques is given to deduce the continuous lattice self-energy from an interpolated cluster self-energy. We study the pseudogap region of a hole-doped two-dimensional Hubbard model and find that, in the DCA${}^{+}\phantom{\rule{0.222222em}{0ex}}$algorithm, the self-energy and pseudogap temperature ${T}^{*}$ converge monotonously with cluster size. Introduction of a continuous lattice self-energy eliminates artificial long-range correlations and thus significantly reduces the sign problem of the quantum Monte Carlo cluster solver in the DCA${}^{+}\phantom{\rule{0.222222em}{0ex}}$algorithm compared to the normal DCA. Simulations with much larger cluster sizes thus become feasible, which, along with the improved convergence in cluster size, raises hope that precise extrapolations to the exact infinite cluster size limit can be reached for other physical quantities as well.

arXiv:1304.3624•Physical Review B•Sep 3, 2013•36 citations

Cluster-size dependence in cellular dynamical mean-field theory

志朗 酒井, Giorgio Sangiovanni, Marcello Civelli +3 more

We examine the cluster-size dependence of the cellular dynamical mean-field theory (CDMFT) applied to the two-dimensional Hubbard model. Employing the continuous-time quantum Monte Carlo method as the solver for the effective cluster model, we obtain CDMFT solutions for 4-, 8-, 12-, and 16-site clusters at a low temperature. Comparing various periodization schemes, which are used to construct the infinite-lattice quantities from the cluster results, we find that the cumulant periodization yields the fastest convergence for the hole-doped Mott insulator where the most severe size dependence is expected. We also find that the convergence is much faster around $(0,0)$ and $(\frac{\ensuremath{\pi}}{2},\frac{\ensuremath{\pi}}{2})$ than around $(\ensuremath{\pi},0)$ and $(\ensuremath{\pi},\ensuremath{\pi})$. The cumulant-periodized self-energy seems to be close to its thermodynamic limit already for a 16-site cluster in the range of parameters studied. The 4-site results remarkably agree well with the 16-site results, indicating that the previous studies based on the 4-site cluster capture the essence of the physics of doped Mott insulators.

arXiv:1112.3227•Physical Review B•Jan 4, 2012•69 citations

Continuous momentum dependence in the dynamical cluster approximation

Urs R. Hähner, Thomas Maier, T. C. Schulthess

The dynamical cluster approximation (DCA) is a quantum cluster extension to the single-site dynamical mean-field theory that incorporates spatially nonlocal dynamic correlations systematically and nonperturbatively. The ${\mathrm{DCA}}^{+}$ algorithm addresses the cluster shape dependence of the DCA and improves the convergence with cluster size by introducing a lattice self-energy with continuous momentum dependence. However, we show that the ${\mathrm{DCA}}^{+}$ algorithm is plagued by a fundamental problem when its self-consistency equations are formulated using the bare Green's function of the cluster. This problem is most severe in the strongly correlated regime at low doping, where the ${\mathrm{DCA}}^{+}$ self-energy becomes overly metallic and local, and persists to cluster sizes where the standard DCA has long converged. In view of the failure of the ${\mathrm{DCA}}^{+}$ algorithm, we propose to complement DCA simulations with a post-interpolation procedure for single-particle and two-particle correlation functions to preserve continuous momentum dependence and the associated benefits in the DCA. We demonstrate the effectiveness of this practical approach with results for the half-filled and hole-doped two-dimensional Hubbard model.

arXiv:2002.06866•Physical review. B./Physical review. B•May 11, 2020•4 citations

Dynamical Mean Field Theory: Application to the Strongly Correlated Materials

Prosenjit Haldar

In this chapter, we present a bird’s-eye view of dynamical mean field theory (DMFT) and its extensions to highly correlated systems. We first demonstrate the elementary introduction to DMFT on a single site with detailed derivations and then a solution of the simple interacting models: the Hubbard model and the Falicov–Kimball model on infinite dimensional Bethe lattice using DMFT. Thereafter, we focus on two kinds of non-local extensions of single-site DMFT - (a) cluster dynamical mean field theory (CDMFT) on real space and (b) dynamical cluster approximation (DCA) on momentum space. Finally, we conclude the chapter with the application of DMFT combining with density functional theory (DFT) in the real-material systems (LDA+DMFT).

River Publishers eBooks•Sep 30, 2025•0 citations

Real-space cluster dynamical mean-field theory: Center-focused extrapolation on the one- and two particle-levels

Marcel Klett, Nils Wentzell, Thomas Schäfer +4 more

We revisit the cellular dynamical mean-field theory (CDMFT) for the single-band Hubbard model on the square lattice at half filling, reaching real-space cluster sizes of up to 9 × 9 sites.Using benchmarks against direct lattice diagrammatic Monte Carlo at high temperature, we show that the self-energy obtained from a cluster center-focused extrapolation converges faster with the cluster size than the periodization schemes previously introduced in the literature.The same benchmark also shows that the cluster spin susceptibility can be extrapolated to the exact result at large cluster size, even though its spatial extension is larger than the cluster size.

arXiv:2003.05215•Physical Review Research•Sep 23, 2020•29 citations
Two widely-used non-local extensions of dynamical mean-field theory (DMFT), namely cellular DMFT (CDMFT) and the dynamical cluster approximation (DCA), both yield self-energies marred by having some unphysical properties: CDMFT yields real-space self-energies that are not translationally invariant, and DCA yields momentum-space self-energies with discontinuities in their momentum dependence. It is often desirable to remove these flaws by post-processing cluster DMFT results, using strategies called periodization for CDMFT and interpolation for DCA---for brevity, we refer to both cases as \textit{interpolation}. However, traditional interpolation approaches struggle to reliably capture intricate structures such as hole pockets in the hole-doped square-lattice Hubbard model. Further, these approaches interpolate frequency-dependent functions, which may lead to causality violations. Here, we propose \textit{Liouvillian interpolation}, a novel, intuitive, and robust scheme for interpolating cluster DMFT results. Our key idea is to interpolate frequency-\textit{independent} matrix elements of the single-particle irreducible part of the Liouvillian, obtained from a continued-fraction expansion of the cluster DMFT self-energy. We demonstrate that the ingredients of such an expansion possess a more local Fourier expansion than the functions involved in traditional interpolation schemes, and that Liouvillian interpolation conserves causality by construction. We illustrate our method for the one-dimensional Hubbard model using CDMFT, and for the two-dimensional Hubbard model using four-site CDMFT and four-patch DCA. For the latter, L\mathcal{L}L-interpolation can (depending on doping) yield Fermi and Luttinger arcs which together seem to form a closed surface.
Even after decades of intense research, the single band Hubbard model representing the fundamental model for interacting quantum systems and electronic correlations remains an unsolved cornerstone paradigm in theoretical solid state physics. Within the non-perturbative dynamical mean-field theory (DMFT), the lattice problem is mapped to a self-consistent auxiliary quantum impurity model leading to a local approximation of the self-energy. Nonlocal correlations can be included by the real-space cluster extension of the DMFT (CDMFT). Enlarging the unit cell makes the calculation of a numerically exact solution in CDMFT very challenging. We here use a quantum Monte-Carlo approach in the imaginary-time space in order to solve the self-consistency equations for different cluster sizes. In order to restore the broken translational symmetry re-periodization schemes for the Green function, the self-energy, or its cumulants have been introduced. However, these suffer from ambiguity and may even lead to convergence problems when attempted inside the self-consistency loop. The comparison to numerically exact diagrammatic quantum Monte Carlo calculations shows that introducing a so-called centerfocused extrapolation (CFE) to approximate the lattice self-energy yields very accurate results. Moreover, the CFE converges faster with the cluster size than previous periodization schemes. We here perform a detailed CDMFT analysis for the single-band Hubbard model on the 2D square lattice, reaching real-space cluster sizes of up to 9x9 sites. In addition to spectral properties, we also compute two-particle correlation functions which are not accessible in DMFT. Using benchmarks against diagrammatic Monte Carlo at high temperature, we show that the cluster spin susceptibility can be extrapolated to the exact result at large cluster size. In particular, the exponential decay of spin-spin correlations is very well captured by CDMFT calculations, even when correlations extend beyond the size of the cluster. We further present results at lower temperature T and larger UUU than the range currently accessible with diagrammatic Monte Carlo methods techniques. The CDMFT+CFE represents therefore a powerful computational tool to access the physics of non-local correlations beyond dynamical mean-field theory.
The dynamical cluster approximation (DCA) is a systematic extension beyond the single-site approximation in dynamical mean field theory, to include spatially nonlocal correlations in quantum many-body simulations of strongly correlated systems. We extend the DCA with a continuous lattice self-energy in order to achieve better convergence with cluster size. This method, which we call DCA+{}^{+}+, cures the cluster-shape dependence problems of the DCA, without suffering from causality violations of previous attempts to interpolate the cluster self-energy. A practical approach based on standard inference techniques is given to deduce the continuous lattice self-energy from an interpolated cluster self-energy. We study the pseudogap region of a hole-doped two-dimensional Hubbard model and find that, in the DCA+{}^{+}\phantom{\rule{0.222222em}{0ex}}+algorithm, the self-energy and pseudogap temperature T∗{T}^{*}T∗ converge monotonously with cluster size. Introduction of a continuous lattice self-energy eliminates artificial long-range correlations and thus significantly reduces the sign problem of the quantum Monte Carlo cluster solver in the DCA+{}^{+}\phantom{\rule{0.222222em}{0ex}}+algorithm compared to the normal DCA. Simulations with much larger cluster sizes thus become feasible, which, along with the improved convergence in cluster size, raises hope that precise extrapolations to the exact infinite cluster size limit can be reached for other physical quantities as well.
We examine the cluster-size dependence of the cellular dynamical mean-field theory (CDMFT) applied to the two-dimensional Hubbard model. Employing the continuous-time quantum Monte Carlo method as the solver for the effective cluster model, we obtain CDMFT solutions for 4-, 8-, 12-, and 16-site clusters at a low temperature. Comparing various periodization schemes, which are used to construct the infinite-lattice quantities from the cluster results, we find that the cumulant periodization yields the fastest convergence for the hole-doped Mott insulator where the most severe size dependence is expected. We also find that the convergence is much faster around (0,0)(0,0)(0,0) and $(\frac{\ensuremath{\pi}}{2},\frac{\ensuremath{\pi}}{2})$ than around $(\ensuremath{\pi},0)$ and $(\ensuremath{\pi},\ensuremath{\pi})$. The cumulant-periodized self-energy seems to be close to its thermodynamic limit already for a 16-site cluster in the range of parameters studied. The 4-site results remarkably agree well with the 16-site results, indicating that the previous studies based on the 4-site cluster capture the essence of the physics of doped Mott insulators.
The dynamical cluster approximation (DCA) is a quantum cluster extension to the single-site dynamical mean-field theory that incorporates spatially nonlocal dynamic correlations systematically and nonperturbatively. The DCA+{\mathrm{DCA}}^{+}DCA+ algorithm addresses the cluster shape dependence of the DCA and improves the convergence with cluster size by introducing a lattice self-energy with continuous momentum dependence. However, we show that the DCA+{\mathrm{DCA}}^{+}DCA+ algorithm is plagued by a fundamental problem when its self-consistency equations are formulated using the bare Green's function of the cluster. This problem is most severe in the strongly correlated regime at low doping, where the DCA+{\mathrm{DCA}}^{+}DCA+ self-energy becomes overly metallic and local, and persists to cluster sizes where the standard DCA has long converged. In view of the failure of the DCA+{\mathrm{DCA}}^{+}DCA+ algorithm, we propose to complement DCA simulations with a post-interpolation procedure for single-particle and two-particle correlation functions to preserve continuous momentum dependence and the associated benefits in the DCA. We demonstrate the effectiveness of this practical approach with results for the half-filled and hole-doped two-dimensional Hubbard model.
In this chapter, we present a bird’s-eye view of dynamical mean field theory (DMFT) and its extensions to highly correlated systems. We first demonstrate the elementary introduction to DMFT on a single site with detailed derivations and then a solution of the simple interacting models: the Hubbard model and the Falicov–Kimball model on infinite dimensional Bethe lattice using DMFT. Thereafter, we focus on two kinds of non-local extensions of single-site DMFT - (a) cluster dynamical mean field theory (CDMFT) on real space and (b) dynamical cluster approximation (DCA) on momentum space. Finally, we conclude the chapter with the application of DMFT combining with density functional theory (DFT) in the real-material systems (LDA+DMFT).
We revisit the cellular dynamical mean-field theory (CDMFT) for the single-band Hubbard model on the square lattice at half filling, reaching real-space cluster sizes of up to 9 × 9 sites.Using benchmarks against direct lattice diagrammatic Monte Carlo at high temperature, we show that the self-energy obtained from a cluster center-focused extrapolation converges faster with the cluster size than the periodization schemes previously introduced in the literature.The same benchmark also shows that the cluster spin susceptibility can be extrapolated to the exact result at large cluster size, even though its spatial extension is larger than the cluster size.