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

Low-energy effective Hamiltonians for correlated electron systems beyond density functional theory

Motoaki Hirayama, Takashi Miyake, Masatoshi Imada, Silke Biermann

Published Aug 1, 2017• 30 citations
arXivHTMLPDF

Abstract

We propose a refined scheme of deriving an effective low-energy Hamiltonian for materials with strong electronic Coulomb correlations beyond density functional theory (DFT). By tracing out the electronic states away from the target degrees of freedom in a controlled way by a perturbative scheme, we construct an effective Hamiltonian for a restricted low-energy target space incorporating the effects of high-energy degrees of freedom in an effective manner. The resulting effective Hamiltonian can afterwards be solved by accurate many-body solvers. We improve this ``multiscale ab initio scheme for correlated electrons'' (MACE) primarily in two directions by elaborating and combining two frameworks developed by Hirayama et al. [M. Hirayama, T. Miyake, and M. Imada, Phys. Rev. B 87, 195144 (2013)] and Casula et al. [M. Casula, P. Werner, L. Vaugier, F. Aryasetiawan, T. Miyake, A. J. Millis, and S. Biermann, Phys. Rev. Lett. 109, 126408 (2012)]: (1) Double counting of electronic correlations between the DFT and the low-energy solver is avoided by using the constrained $GW$ scheme; and (2) the frequency dependent interactions emerging from the partial trace summation are successfully separated into a nonlocal part that is treated following ideas by Hirayama et al. and a local part treated nonperturbatively in the spirit of Casula et al. and are incorporated into the renormalization of the low-energy dispersion. The scheme is favorably tested on the example of ${\mathrm{SrVO}}_{3}$.

DOI
10.1103/physrevb.96.075102

Related papers

Diagrammatics in the Dual Space, or There and Back Again

E. A. Stepanov

Accurately describing many-body effects in multi-orbital systems remains a major challenge in theoretical condensed matter physics. At present, there is a significant methodological gap between the numerical tools used in ab initio computational materials science and those developed to study strong electronic correlations. The former can treat realistic, large-scale systems but typically neglect many-body effects, while the latter focus on simplified models with only a few degrees of freedom, as only such models can be solved accurately in the presence of strong interactions. The purpose of this thesis is to bridge these two approaches and establish a systematic theoretical framework for realistic correlated electronic materials. This involves a full-cycle methodology that begins with constructing ab initio interacting models from density-functional theory, solving them using dynamical mean-field theory to capture local correlations, and extending beyond to incorporate non-local collective electronic fluctuations. To this end, we introduce the "dual" approach to strong correlations, which includes the dual fermion, dual boson, and dual triply irreducible local expansion methods. The central idea of the dual theories is to shift the reference point of the conventional Feynman diagrammatic expansion from a non-interacting electronic system to an interacting but exactly solvable one. Integrating out this reference system recasts the expansion in an effective dual space, where all diagrammatic building blocks are renormalized by the corresponding impurity quantities. This procedure transforms a non-perturbative expansion in the original variables into a perturbative expansion in terms of dual fermionic and bosonic fields, exact in both the weak- and strong-coupling limits. In this thesis, we collect and systematize the major developments of dual techniques achieved to date.

arXiv:2512.12389•Dec 13, 2025•0 citations

Self-consistent many-body perturbation theory in range-separated density-functional theory: A one-electron reduced-density-matrix-based formulation

Emmanuel Fromager, Hans Jørgen Aa. Jensen

In many cases, density-functional theory (DFT) with current standard approximate functionals offers a relatively accurate and computationally cheap description of the short-range dynamic electron correlation effects. However, in general, standard DFT does not treat the dispersion interaction effects adequately which, on the other hand, can be described by many-body perturbation theory (MBPT). It is therefore of interest to develop a hybrid model which combines the best of both the MBPT and DFT approaches. This can be achieved by splitting the two-electron interaction into long-range and short-range parts; the long-range part is then treated by MBPT and the short-range part by DFT. This work deals with the formulation of a general MBPT-DFT model (i.e., valid for any type of zeroth-order Hamiltonian) based on such a range separation. Applying the Rayleigh-Schr\"odinger formalism in this context, one finds that the generalized Bloch equation becomes self-consistent at each order of perturbation. This complication arises because the short-range part of the energy is a functional of the exact electron density, which is expanded in a perturbation series. In order to address this ``self-consistency problem'' and provide computable orbital-based expressions for any order of perturbation, a general one-electron reduced-density-matrix-based formalism is proposed. Two applications of our general formalism are presented: The derivation of a hybrid second-order M\o{}ller-Plesset-DFT model and the formulation of a range-separated optimized effective potential method based on a hybrid G\"orling-Levy-type MBPT-DFT model.

Physical Review A•Aug 7, 2008•63 citations

Strong Correlation in Kohn-Sham Density Functional Theory

Francesc Malet, Paola Gori‐Giorgi

We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.

arXiv:1207.2775•Physical Review Letters•Dec 10, 2012•118 citations

Electronic Structure Calculation by First Principles for Strongly Correlated Electron Systems

Masatoshi Imada, Takashi Miyake

Recent trends of ab initio studies and progress in methodologies for electronic structure calculations of strongly correlated electron systems are discussed. The interest for developing efficient methods is motivated by recent discoveries and characterizations of strongly correlated electron materials and by requirements for understanding mechanisms of intriguing phenomena beyond a single-particle picture. A three-stage scheme is developed as renormalized multi-scale solvers (RMS) utilizing the hierarchical electronic structure in the energy space. It provides us with an ab initio downfolding of the global band structure into low-energy effective models followed by low-energy solvers for the models. The RMS method is illustrated with examples of several materials. In particular, we overview cases such as dynamics of semiconductors, transition metals and its compounds including iron-based superconductors and perovskite oxides, as well as organic conductors of kappa-ET type.

arXiv:1009.3851•Journal of the Physical Society of Japan•Nov 10, 2010•140 citations

Systematically improvable multiscale solver for correlated electron systems

Alexei A. Kananenka, Emanuel Gull, Dominika Zgid

The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. The general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.

arXiv:1410.5118•Physical Review B•Mar 23, 2015•91 citations

Ab initio procedure for constructing effective models of correlated materials with entangled band structure

Takashi Miyake, F. Aryasetiawan, Masatoshi Imada

In a previous work [Phys. Rev. B 77, 085122 (2008)], a procedure for constructing low-energy models of electrons in solids was proposed. The procedure starts with dividing the Hilbert space into two subspaces: the low-energy part (``$d$ space'') and the rest of the space (``$r$ space''). The low-energy model is constructed for the $d$ space by eliminating the degrees of freedom of the $r$ space. The thus derived model contains the strength of electron correlation expressed by a partially screened Coulomb interaction, calculated in the constrained random-phase approximation (cRPA), where screening channels within the $d$ space, ${P}_{d}$, are subtracted. One conceptual problem of this established downfolding method is that for entangled bands it is not clear how to cut out the $d$ space and how to distinguish ${P}_{d}$ from the total polarization. Here, we propose a simple procedure to overcome this difficulty. The $d$ space is defined to be an isolated set of bands generated from a set of maximally localized Wannier basis, which consequently defines ${P}_{d}$. The $r$ subspace is constructed as the complementary space orthogonal to the $d$ subspace, resulting in two sets of completely disentangled bands. Using these disentangled bands, the effective parameters of the $d$ space are uniquely determined by the cRPA method. The method is successfully applied to $3d$ transition metals.

arXiv:0906.1344•Physical Review B•Oct 27, 2009•214 citations
We propose a refined scheme of deriving an effective low-energy Hamiltonian for materials with strong electronic Coulomb correlations beyond density functional theory (DFT). By tracing out the electronic states away from the target degrees of freedom in a controlled way by a perturbative scheme, we construct an effective Hamiltonian for a restricted low-energy target space incorporating the effects of high-energy degrees of freedom in an effective manner. The resulting effective Hamiltonian can afterwards be solved by accurate many-body solvers. We improve this ``multiscale ab initio scheme for correlated electrons'' (MACE) primarily in two directions by elaborating and combining two frameworks developed by Hirayama et al. [M. Hirayama, T. Miyake, and M. Imada, Phys. Rev. B 87, 195144 (2013)] and Casula et al. [M. Casula, P. Werner, L. Vaugier, F. Aryasetiawan, T. Miyake, A. J. Millis, and S. Biermann, Phys. Rev. Lett. 109, 126408 (2012)]: (1) Double counting of electronic correlations between the DFT and the low-energy solver is avoided by using the constrained GWGWGW scheme; and (2) the frequency dependent interactions emerging from the partial trace summation are successfully separated into a nonlocal part that is treated following ideas by Hirayama et al. and a local part treated nonperturbatively in the spirit of Casula et al. and are incorporated into the renormalization of the low-energy dispersion. The scheme is favorably tested on the example of SrVO3{\mathrm{SrVO}}_{3}SrVO3​.
Accurately describing many-body effects in multi-orbital systems remains a major challenge in theoretical condensed matter physics. At present, there is a significant methodological gap between the numerical tools used in ab initio computational materials science and those developed to study strong electronic correlations. The former can treat realistic, large-scale systems but typically neglect many-body effects, while the latter focus on simplified models with only a few degrees of freedom, as only such models can be solved accurately in the presence of strong interactions. The purpose of this thesis is to bridge these two approaches and establish a systematic theoretical framework for realistic correlated electronic materials. This involves a full-cycle methodology that begins with constructing ab initio interacting models from density-functional theory, solving them using dynamical mean-field theory to capture local correlations, and extending beyond to incorporate non-local collective electronic fluctuations. To this end, we introduce the "dual" approach to strong correlations, which includes the dual fermion, dual boson, and dual triply irreducible local expansion methods. The central idea of the dual theories is to shift the reference point of the conventional Feynman diagrammatic expansion from a non-interacting electronic system to an interacting but exactly solvable one. Integrating out this reference system recasts the expansion in an effective dual space, where all diagrammatic building blocks are renormalized by the corresponding impurity quantities. This procedure transforms a non-perturbative expansion in the original variables into a perturbative expansion in terms of dual fermionic and bosonic fields, exact in both the weak- and strong-coupling limits. In this thesis, we collect and systematize the major developments of dual techniques achieved to date.
In many cases, density-functional theory (DFT) with current standard approximate functionals offers a relatively accurate and computationally cheap description of the short-range dynamic electron correlation effects. However, in general, standard DFT does not treat the dispersion interaction effects adequately which, on the other hand, can be described by many-body perturbation theory (MBPT). It is therefore of interest to develop a hybrid model which combines the best of both the MBPT and DFT approaches. This can be achieved by splitting the two-electron interaction into long-range and short-range parts; the long-range part is then treated by MBPT and the short-range part by DFT. This work deals with the formulation of a general MBPT-DFT model (i.e., valid for any type of zeroth-order Hamiltonian) based on such a range separation. Applying the Rayleigh-Schr\"odinger formalism in this context, one finds that the generalized Bloch equation becomes self-consistent at each order of perturbation. This complication arises because the short-range part of the energy is a functional of the exact electron density, which is expanded in a perturbation series. In order to address this ``self-consistency problem'' and provide computable orbital-based expressions for any order of perturbation, a general one-electron reduced-density-matrix-based formalism is proposed. Two applications of our general formalism are presented: The derivation of a hybrid second-order M\o{}ller-Plesset-DFT model and the formulation of a range-separated optimized effective potential method based on a hybrid G\"orling-Levy-type MBPT-DFT model.
We use the exact strong-interaction limit of the Hohenberg-Kohn energy density functional to approximate the exchange-correlation energy of the restricted Kohn-Sham scheme. Our approximation corresponds to a highly nonlocal density functional whose functional derivative can be easily constructed, thus transforming exactly, in a physically transparent way, an important part of the electron-electron interaction into an effective local one-body potential. We test our approach on quasi-one-dimensional systems, showing that it captures essential features of strong correlation that restricted Kohn-Sham calculations using the currently available approximations cannot describe.
Recent trends of ab initio studies and progress in methodologies for electronic structure calculations of strongly correlated electron systems are discussed. The interest for developing efficient methods is motivated by recent discoveries and characterizations of strongly correlated electron materials and by requirements for understanding mechanisms of intriguing phenomena beyond a single-particle picture. A three-stage scheme is developed as renormalized multi-scale solvers (RMS) utilizing the hierarchical electronic structure in the energy space. It provides us with an ab initio downfolding of the global band structure into low-energy effective models followed by low-energy solvers for the models. The RMS method is illustrated with examples of several materials. In particular, we overview cases such as dynamics of semiconductors, transition metals and its compounds including iron-based superconductors and perovskite oxides, as well as organic conductors of kappa-ET type.
The development of numerical methods capable of simulating realistic materials with strongly correlated electrons, with controllable errors, is a central challenge in quantum many-body physics. Here we describe a framework for a general multiscale method based on embedding a self-energy of a strongly correlated subsystem into a self-energy generated by a method able to treat large weakly correlated systems approximately. As an example, we present the embedding of an exact diagonalization self-energy into a self-energy generated from self-consistent second-order perturbation theory. Using a quantum impurity model, generated from a cluster dynamical mean field approximation to the two-dimensional Hubbard model, as a benchmark, we illustrate that our method allows us to obtain accurate results at a fraction of the cost of typical Monte Carlo calculations. We test the method in multiple regimes of interaction strengths and dopings of the model. The general embedding framework we present avoids difficulties such as double counting corrections, frequency-dependent interactions, or vertex functions. As it is solely formulated at the level of the single-particle Green's function, it provides a promising route for the simulation of realistic materials that are currently difficult to study with other methods.
In a previous work [Phys. Rev. B 77, 085122 (2008)], a procedure for constructing low-energy models of electrons in solids was proposed. The procedure starts with dividing the Hilbert space into two subspaces: the low-energy part (``ddd space'') and the rest of the space (``rrr space''). The low-energy model is constructed for the ddd space by eliminating the degrees of freedom of the rrr space. The thus derived model contains the strength of electron correlation expressed by a partially screened Coulomb interaction, calculated in the constrained random-phase approximation (cRPA), where screening channels within the ddd space, Pd{P}_{d}Pd​, are subtracted. One conceptual problem of this established downfolding method is that for entangled bands it is not clear how to cut out the ddd space and how to distinguish Pd{P}_{d}Pd​ from the total polarization. Here, we propose a simple procedure to overcome this difficulty. The ddd space is defined to be an isolated set of bands generated from a set of maximally localized Wannier basis, which consequently defines Pd{P}_{d}Pd​. The rrr subspace is constructed as the complementary space orthogonal to the ddd subspace, resulting in two sets of completely disentangled bands. Using these disentangled bands, the effective parameters of the ddd space are uniquely determined by the cRPA method. The method is successfully applied to 3d3d3d transition metals.