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:2003.05215 · Physical Review Research

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

Marcel Klett, Nils Wentzell, Thomas Schäfer, Fedor Šimkovic, Olivier Parcollet, Sabine Andergassen, P. Hansmann

Published Sep 23, 2020• 29 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevresearch.2.033476

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

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

Mott transition and pseudogap of the square-lattice Hubbard model: Results from center-focused cellular dynamical mean-field theory

Michael Meixner, Henri Menke, Marcel Klett +4 more

The recently proposed center-focused post-processing procedure [Phys. Rev. Res. 2, 033476 (2020)] of cellular dynamical mean-field theory suggests that central sites of large impurity clusters are closer to the exact solution of the Hubbard model than the edge sites. In this paper, we systematically investigate results in the spirit of this center-focused scheme for several cluster sizes up to 8×8 in and out of particle-hole symmetry. First we analyze the metal-insulator crossovers and transitions of the half-filled Hubbard model on a simple square lattice. We find that the critical interaction of the crossover is reduced with increasing cluster sizes and the critical temperature abruptly drops for the 4×4 cluster. Second, for this cluster size, we apply the center-focused scheme to a system with more realistic tight-binding parameters, investigating its pseudogap regime as a function of temperature and doping, where we find doping dependent metal-insulator crossovers, Lifshitz transitions and a strongly renormalized Fermi-liquid regime. Additionally to diagnosing the real space origin of the suppressed antinodal spectral weight in the pseudogap regime, we can infer hints towards underlying charge ordering tendencies.

arXiv:2310.17302•SciPost Physics•Feb 27, 2024•14 citations

Quantum Monte Carlo study of strongly correlated electrons: Cellular dynamical mean-field theory

Bumsoo Kyung, Gabriel Kotliar, A.–M. S. Tremblay

We study the Hubbard model using the cellular dynamical mean-field theory (CDMFT) with quantum Monte Carlo (QMC) simulations. We present the algorithmic details of CDMFT with the Hirsch-Fye QMC method for the solution of the self-consistently embedded quantum cluster problem. We use the one- and two-dimensional half filled Hubbard model to gauge the performance of $\mathrm{CDMFT}+\mathrm{QMC}$ particularly for small clusters by comparing with the exact results and also with other quantum cluster methods. We calculate single-particle Green's functions and self-energies on small clusters to study their size dependence in one and two dimensions. It is shown that in one dimension, CDMFT with two sites in the cluster is already able to describe with high accuracy the evolution of the density as a function of the chemical potential and the compressibility divergence at the Mott transition, in good agreement with the exact Bethe ansatz result. With increasing $U$ the result on small clusters rapidly approaches that of the infinite size cluster. Large scattering rate and a positive slope in the real part of the self-energy in one dimension suggest that the system is a non-Fermi liquid for all the parameters studied here. In two dimensions, at intermediate to strong coupling, even the smallest cluster $({N}_{c}=2\ifmmode\times\else\texttimes\fi{}2)$ accounts for more than 95% of the correlation effect of the infinite-size cluster in the single particle spectrum, suggesting that some of the important problems in strongly correlated electron systems may be studied highly accurately with a reasonable computational effort. Finally, as an application that is sensitive to details of correlations, we show that $\mathrm{CDMFT}+\mathrm{QMC}$ can describe spin-charge separated Luttinger liquid physics in one dimension. The spinon and holon branches appear only for sufficiently large system sizes.

arXiv:cond-mat/0601271•Physical Review B•May 10, 2006•42 citations

Real-space renormalized dynamical mean field theory

Dai Kubota, 志朗 酒井, Masatoshi Imada

We propose real-space renormalized dynamical mean field theory (rr-DMFT) to deal with large clusters in the framework of a cluster extension of the DMFT. In the rr-DMFT, large clusters are decomposed into multiple smaller clusters through a real-space renormalization. In this work, the renormalization effect is taken into account only at the lowest order with respect to the intercluster coupling, which nonetheless reproduces exactly both the noninteracting and atomic limits. Our method allows us large cluster-size calculations which are intractable with the conventional cluster extensions of the DMFT with impurity solvers, such as the continuous-time quantum Monte Carlo and exact diagonalization methods. We benchmark the rr-DMFT for the two-dimensional Hubbard model on a square lattice at and away from half filling, where the spatial correlations play important roles. Our results on the spin structure factor indicate that the growth of the antiferromagnetic spin correlation is taken into account beyond the decomposed cluster size. We also show that the self-energy obtained from the large-cluster solver is reproduced by our method better than the solution obtained directly for the smaller cluster. When applied to the Mott metal-insulator transition, the rr-DMFT is able to reproduce the reduced critical value for the Coulomb interaction comparable to the large cluster result.

arXiv:1601.05916•Physical review. B./Physical review. B•May 12, 2016•3 citations

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

Mathias Pelz, Jan von Delft, Andreas Gleis

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.

arXiv:2602.16351•Feb 18, 2026•0 citations
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.
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.
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 recently proposed center-focused post-processing procedure [Phys. Rev. Res. 2, 033476 (2020)] of cellular dynamical mean-field theory suggests that central sites of large impurity clusters are closer to the exact solution of the Hubbard model than the edge sites. In this paper, we systematically investigate results in the spirit of this center-focused scheme for several cluster sizes up to 8×8 in and out of particle-hole symmetry. First we analyze the metal-insulator crossovers and transitions of the half-filled Hubbard model on a simple square lattice. We find that the critical interaction of the crossover is reduced with increasing cluster sizes and the critical temperature abruptly drops for the 4×4 cluster. Second, for this cluster size, we apply the center-focused scheme to a system with more realistic tight-binding parameters, investigating its pseudogap regime as a function of temperature and doping, where we find doping dependent metal-insulator crossovers, Lifshitz transitions and a strongly renormalized Fermi-liquid regime. Additionally to diagnosing the real space origin of the suppressed antinodal spectral weight in the pseudogap regime, we can infer hints towards underlying charge ordering tendencies.
We study the Hubbard model using the cellular dynamical mean-field theory (CDMFT) with quantum Monte Carlo (QMC) simulations. We present the algorithmic details of CDMFT with the Hirsch-Fye QMC method for the solution of the self-consistently embedded quantum cluster problem. We use the one- and two-dimensional half filled Hubbard model to gauge the performance of CDMFT+QMC\mathrm{CDMFT}+\mathrm{QMC}CDMFT+QMC particularly for small clusters by comparing with the exact results and also with other quantum cluster methods. We calculate single-particle Green's functions and self-energies on small clusters to study their size dependence in one and two dimensions. It is shown that in one dimension, CDMFT with two sites in the cluster is already able to describe with high accuracy the evolution of the density as a function of the chemical potential and the compressibility divergence at the Mott transition, in good agreement with the exact Bethe ansatz result. With increasing UUU the result on small clusters rapidly approaches that of the infinite size cluster. Large scattering rate and a positive slope in the real part of the self-energy in one dimension suggest that the system is a non-Fermi liquid for all the parameters studied here. In two dimensions, at intermediate to strong coupling, even the smallest cluster $({N}_{c}=2\ifmmode\times\else\texttimes\fi{}2)$ accounts for more than 95% of the correlation effect of the infinite-size cluster in the single particle spectrum, suggesting that some of the important problems in strongly correlated electron systems may be studied highly accurately with a reasonable computational effort. Finally, as an application that is sensitive to details of correlations, we show that CDMFT+QMC\mathrm{CDMFT}+\mathrm{QMC}CDMFT+QMC can describe spin-charge separated Luttinger liquid physics in one dimension. The spinon and holon branches appear only for sufficiently large system sizes.
We propose real-space renormalized dynamical mean field theory (rr-DMFT) to deal with large clusters in the framework of a cluster extension of the DMFT. In the rr-DMFT, large clusters are decomposed into multiple smaller clusters through a real-space renormalization. In this work, the renormalization effect is taken into account only at the lowest order with respect to the intercluster coupling, which nonetheless reproduces exactly both the noninteracting and atomic limits. Our method allows us large cluster-size calculations which are intractable with the conventional cluster extensions of the DMFT with impurity solvers, such as the continuous-time quantum Monte Carlo and exact diagonalization methods. We benchmark the rr-DMFT for the two-dimensional Hubbard model on a square lattice at and away from half filling, where the spatial correlations play important roles. Our results on the spin structure factor indicate that the growth of the antiferromagnetic spin correlation is taken into account beyond the decomposed cluster size. We also show that the self-energy obtained from the large-cluster solver is reproduced by our method better than the solution obtained directly for the smaller cluster. When applied to the Mott metal-insulator transition, the rr-DMFT is able to reproduce the reduced critical value for the Coulomb interaction comparable to the large cluster result.
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.