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

Adaptively truncated Hilbert space based impurity solver for dynamical mean-field theory

Ara Go, Andrew J. Millis

Published Aug 28, 2017• 36 citations
arXivHTMLPDF

Abstract

We present an impurity solver based on adaptively truncated Hilbert spaces. The solver is particularly suitable for dynamical mean-field theory in circumstances where quantum Monte Carlo approaches are ineffective. It exploits the sparsity structure of quantum impurity models, in which the interactions couple only a small subset of the degrees of freedom. We further introduce an adaptive truncation of the particle or hole excited spaces, which enables computations of Green functions with an accuracy needed to avoid unphysical (sign change of imaginary part) self-energies. The method is benchmarked on the one-dimensional Hubbard model.

DOI
10.1103/physrevb.96.085139

Related papers

Quantum Monte Carlo impurity solvers for multi-orbital problems and frequency-dependent interactions

Hiroshi Shinaoka, Fakher F. Assaad, Nils Blümer +1 more

The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.

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

Next Generation Dynamical Mean-Field Theory Simulations with the Adaptive Sampling Configuration Interaction Method

Carlos Mejuto-Zaera, Norm M. Tubman, K. Birgitta Whaley

In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. In this Letter, we present the ASCI-DMFT method and use the one- and two-dimensional Hubbard models to exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground state DMFT simulations. Our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.

•Nov 13, 2017•1 citation

Dynamical mean field theory simulations with the adaptive sampling configuration interaction method

Carlos Mejuto-Zaera, Norm M. Tubman, K. Birgitta Whaley

In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero-temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. Here, we present the ASCI-DMFT method and example calculations on the one- and two-dimensional Hubbard models that exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground-state DMFT simulations, which allows us to study the bath discretization error in simulations with small clusters, as well as to address cluster sizes beyond the current state of the art. Our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.

arXiv:1711.04771•Physical review. B./Physical review. B•Sep 30, 2019•31 citations

Neural-Network Quantum Embedding Solvers for Correlated Materials

Agnes Valenti, Ina Park, Antoine Georges +2 more

Quantum impurity solvers are the computational bottleneck of quantum embedding approaches to correlated materials, such as dynamical mean-field theory (DMFT). We show that neural networks trained on synthetic, material-agnostic data learn the impurity mapping from hybridization functions and local interactions to Green's functions with quantitative accuracy for both model systems and real materials, providing fast solvers for single- and multi-orbital models. Benchmarks against numerically controlled quantum Monte Carlo show that the method reproduces the Mott transition, multi-orbital phase diagrams of Hubbard-Kanamori models, and the electronic properties of SrVO$_3$ and SrMnO$_3$. The learned solvers achieve orders-of-magnitude speedup and can initialize controlled calculations, dramatically accelerating DMFT while preserving accuracy.

arXiv:2603.15741•Mar 16, 2026•0 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

A Scalable Configuration-Interaction Impurity Solver via Active Learning

Jeongmoo Lee, Ara Go

Finite-Hamiltonian impurity solvers provide direct real-frequency spectra and a natural route to enlarged impurity Hamiltonians, but their applicability is limited by the rapid Hilbert-space growth with the number of bath or other added one-particle orbitals. We introduce an active-learning extension of adaptive-truncation configuration interaction (AL-ATCI) that identifies the determinant manifold relevant to the correlated state. The approximation is systematically controlled by the query size N_query, which also provides an internal convergence parameter when no external benchmark is available. Over the benchmark range studied here, the computational cost grows only weakly with bath size, because enlarging the bath mainly expands the combinatorial determinant space rather than the physically relevant manifold. In dynamical mean-field-theory benchmarks for the one-dimensional Hubbard model, AL-ATCI reproduces exact-diagonalization accuracy and extends cellular calculations to clusters as large as N_c = 10. For a three-orbital rotationally invariant Sr2RuO4 impurity problem, we demonstrate systematic convergence of dynamical quantities and a highly compressed configuration space as N_b is increased from 9 to 18. These results substantially alleviate the bath-discretization bottleneck of exact-diagonalization- and configuration-interaction-based impurity solvers and make large-bath and enlarged-orbital calculations more practical.

arXiv:2604.10230•Apr 11, 2026•0 citations
We present an impurity solver based on adaptively truncated Hilbert spaces. The solver is particularly suitable for dynamical mean-field theory in circumstances where quantum Monte Carlo approaches are ineffective. It exploits the sparsity structure of quantum impurity models, in which the interactions couple only a small subset of the degrees of freedom. We further introduce an adaptive truncation of the particle or hole excited spaces, which enables computations of Green functions with an accuracy needed to avoid unphysical (sign change of imaginary part) self-energies. The method is benchmarked on the one-dimensional Hubbard model.
The solution of an auxiliary quantum impurity system is the computationally expensive step in dynamical mean field theory simulations of lattice models and materials. In this review, we discuss Monte Carlo based impurity solvers, which are suitable for a wide range of applications. In particular, we present an efficient implementation of the hybridization expansion approach, which enables the simulation of multiorbital impurity problems with off-diagonal and complex hybridizations, and dynamically screened (retarded) density-density interactions. As a complementary approach, we discuss an impurity solver based on the determinant Monte Carlo method, which scales favorably with inverse temperature and hence provides access to the very low temperature regime. The usefulness of these state-of-the-art impurity solvers is demonstrated with applications to the downfolding problem, i.e., the systematic derivation of dynamically screened interactions for low-energy effective models, and to pyrochlore iridates, where the spin-orbit coupling leads to complex hybridization functions in a multi-orbital system. As a benchmark for cluster extensions of dynamical mean field theory, we also present results from lattice Monte Carlo simulations for the momentum dependence of the pseudo-gap in the half-filled two-dimensional Hubbard model.
In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. In this Letter, we present the ASCI-DMFT method and use the one- and two-dimensional Hubbard models to exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground state DMFT simulations. Our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.
In the pursuit of accurate descriptions of strongly correlated quantum many-body systems, dynamical mean field theory (DMFT) has been an invaluable tool for elucidating the spectral properties and quantum phases of both phenomenological models and ab initio descriptions of real materials. Key to the DMFT process is the self-consistent map of the original system into an Anderson impurity model, the ground state of which is computed using an impurity solver. The power of the method is thus limited by the complexity of the impurity model the solver can handle. Simulating realistic systems generally requires many correlated sites. By adapting the recently proposed adaptive sampling configuration interaction (ASCI) method as an impurity solver, we enable much more efficient zero-temperature DMFT simulations. The key feature of the ASCI method is that it selects only the most relevant Hilbert space degrees of freedom to describe the ground state. This reduces the numerical complexity of the calculation, which will allow us to pursue future DMFT simulations with more correlated impurity sites than in previous works. Here, we present the ASCI-DMFT method and example calculations on the one- and two-dimensional Hubbard models that exemplify its efficient convergence and timing properties. We show that the ASCI approach is several orders of magnitude faster than the current best published ground-state DMFT simulations, which allows us to study the bath discretization error in simulations with small clusters, as well as to address cluster sizes beyond the current state of the art. Our approach can also be adapted for other embedding methods such as density matrix embedding theory and self-energy embedding theory.
Quantum impurity solvers are the computational bottleneck of quantum embedding approaches to correlated materials, such as dynamical mean-field theory (DMFT). We show that neural networks trained on synthetic, material-agnostic data learn the impurity mapping from hybridization functions and local interactions to Green's functions with quantitative accuracy for both model systems and real materials, providing fast solvers for single- and multi-orbital models. Benchmarks against numerically controlled quantum Monte Carlo show that the method reproduces the Mott transition, multi-orbital phase diagrams of Hubbard-Kanamori models, and the electronic properties of SrVO3_33​ and SrMnO3_33​. The learned solvers achieve orders-of-magnitude speedup and can initialize controlled calculations, dramatically accelerating DMFT while preserving accuracy.
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.
Finite-Hamiltonian impurity solvers provide direct real-frequency spectra and a natural route to enlarged impurity Hamiltonians, but their applicability is limited by the rapid Hilbert-space growth with the number of bath or other added one-particle orbitals. We introduce an active-learning extension of adaptive-truncation configuration interaction (AL-ATCI) that identifies the determinant manifold relevant to the correlated state. The approximation is systematically controlled by the query size N_query, which also provides an internal convergence parameter when no external benchmark is available. Over the benchmark range studied here, the computational cost grows only weakly with bath size, because enlarging the bath mainly expands the combinatorial determinant space rather than the physically relevant manifold. In dynamical mean-field-theory benchmarks for the one-dimensional Hubbard model, AL-ATCI reproduces exact-diagonalization accuracy and extends cellular calculations to clusters as large as N_c = 10. For a three-orbital rotationally invariant Sr2RuO4 impurity problem, we demonstrate systematic convergence of dynamical quantities and a highly compressed configuration space as N_b is increased from 9 to 18. These results substantially alleviate the bath-discretization bottleneck of exact-diagonalization- and configuration-interaction-based impurity solvers and make large-bath and enlarged-orbital calculations more practical.