Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition and the vanishing low-level gap are related to the assumption that the high-level density matrix blocks are non-interacting pure-state $v$ -representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy.Density matrix embedding theory (DMET) formally requires the matching of density matrix blocks obtained from high-level and low-level theories, but this is sometimes not achievable in practical calculations. In such a case, the global band gap of the low-level theory vanishes, and this can require additional numerical considerations. We find that both the violation of the exact matching condition and the vanishing low-level gap are related to the assumption that the high-level density matrix blocks are non-interacting pure-state $v$-representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy.

arXiv:2110.09558Oct 18, 20212 citations

Some mathematical insights on Density Matrix Embedding Theory

Eric Cancès B, Fabian M. Faulstich, Alfred Kirsch +2 more

This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. We provide numerical simulations to support our results and comment on the physical meaning of the assumptions under which they hold true. We show that the violation of these assumptions may yield multiple solutions of the DMET equations. We moreover introduce and discuss a specific N-representability problem inherent to DMET.This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. We provide numerical simulations to support our results and comment on the physical meaning of the assumptions under which they hold true. We show that the violation of these assumptions may yield multiple solutions of the DMET equations. We moreover introduce and discuss a specific N-representability problem inherent to DMET.

arXiv:2305.16472May 25, 20236 citations

Density Matrix Embedding: A Simple Alternative to Dynamical Mean-Field Theory

Gerald Knizia, Garnet Kin‐Lic Chan

We introduce density matrix embedding theory (DMET), a quantum embedding theory for computing frequency-independent quantities, such as ground-state properties, of infinite systems. Like dynamical mean-field theory, DMET maps the bulk interacting system to a simpler impurity model and is exact in the noninteracting and atomic limits. Unlike dynamical mean-field theory, DMET is formulated in terms of the frequency-independent local density matrix, rather than the local Green's function. In addition, it features a finite, algebraically constructible bath of only one bath site per impurity site, with no bath discretization error. Frequency independence and the minimal bath make DMET a computationally simple and efficient method. We test the theory in the one-dimensional and two-dimensional Hubbard models at and away from half filling, and we find that compared to benchmark data, total energies, correlation functions, and metal-insulator transitions are well reproduced, at a tiny computational cost.We introduce density matrix embedding theory (DMET), a quantum embedding theory for computing frequency-independent quantities, such as ground-state properties, of infinite systems. Like dynamical mean-field theory, DMET maps the bulk interacting system to a simpler impurity model and is exact in the noninteracting and atomic limits. Unlike dynamical mean-field theory, DMET is formulated in terms of the frequency-independent local density matrix, rather than the local Green's function. In addition, it features a finite, algebraically constructible bath of only one bath site per impurity site, with no bath discretization error. Frequency independence and the minimal bath make DMET a computationally simple and efficient method. We test the theory in the one-dimensional and two-dimensional Hubbard models at and away from half filling, and we find that compared to benchmark data, total energies, correlation functions, and metal-insulator transitions are well reproduced, at a tiny computational cost.

arXiv:1204.5783Physical Review LettersNov 2, 2012527 citations

Local Potential Functional Embedding Theory: A Self-Consistent Flavor of Density Functional Theory for Lattices without Density Functionals

Sajanthan Sekaran, Matthieu Saubanère, Emmanuel Fromager

Quantum embedding is a divide and conquer strategy that aims at solving the electronic Schrödinger equation of sizeable molecules or extended systems. We establish in the present work a clearer and in-principle-exact connection between density matrix embedding theory (DMET) and density-functional theory (DFT) within the simple but nontrivial one-dimensional Hubbard model. For that purpose, we use our recent reformulation of single-impurity DMET as a Householder transformed density-matrix functional embedding theory (Ht-DMFET). On the basis of well-identified density-functional approximations, a self-consistent local potential functional embedding theory (LPFET) is formulated and implemented. Combining both LPFET and DMET numerical results with our formally exact density-functional embedding theory reveals that a single statically embedded impurity can in principle describe the density-driven Mott–Hubbard transition, provided that a complementary density-functional correlation potential (which is neglected in both DMET and LPFET) exhibits a derivative discontinuity (DD) at half filling. The extension of LPFET to multiple impurities (which would enable to circumvent the modeling of DDs) and its generalization to quantum chemical Hamiltonians are left for future work.Quantum embedding is a divide and conquer strategy that aims at solving the electronic Schrödinger equation of sizeable molecules or extended systems. We establish in the present work a clearer and in-principle-exact connection between density matrix embedding theory (DMET) and density-functional theory (DFT) within the simple but nontrivial one-dimensional Hubbard model. For that purpose, we use our recent reformulation of single-impurity DMET as a Householder transformed density-matrix functional embedding theory (Ht-DMFET). On the basis of well-identified density-functional approximations, a self-consistent local potential functional embedding theory (LPFET) is formulated and implemented. Combining both LPFET and DMET numerical results with our formally exact density-functional embedding theory reveals that a single statically embedded impurity can in principle describe the density-driven Mott–Hubbard transition, provided that a complementary density-functional correlation potential (which is neglected in both DMET and LPFET) exhibits a derivative discontinuity (DD) at half filling. The extension of LPFET to multiple impurities (which would enable to circumvent the modeling of DDs) and its generalization to quantum chemical Hamiltonians are left for future work.

arXiv:2202.08071ComputationMar 18, 202213 citations

New perspectives on Density-Matrix Embedding Theory

Alicia Negre, Fabian M. Faulstich, Raehyun Kim +3 more

Quantum embedding methods enable the study of large, strongly correlated quantum systems by (usually self-consistent) decomposition into computationally manageable subproblems, in the spirit of divide-and-conquer methods. Among these, Density Matrix Embedding Theory (DMET) is an efficient approach that enforces self-consistency at the level of one-particle reduced density matrices (1-RDMs), facilitating applications across diverse quantum systems. However, conventional DMET is constrained by the requirement that the global 1-RDM (low-level descriptor) be an orthogonal projector, limiting flexibility in bath construction and potentially impeding accuracy in strongly correlated regimes. In this work, we introduce a generalized DMET framework in which the low-level descriptor can be an arbitrary 1-RDM and the bath construction is based on optimizing a quantitative criterion related to the maximal disentanglement between different fragments. This yields an alternative yet controllable bath space construction for generic 1-RDMs, lifting a key limitation of conventional DMET. We demonstrate its consistency with conventional DMET in appropriate limits and exploring its implications for bath construction, downfolding (impurity Hamiltonian construction), low-level solvers, and adaptive fragmentation. We expect that this more flexible framework, which leads to several new variants of DMET, can improve the robustness and accuracy of DMET.Quantum embedding methods enable the study of large, strongly correlated quantum systems by (usually self-consistent) decomposition into computationally manageable subproblems, in the spirit of divide-and-conquer methods. Among these, Density Matrix Embedding Theory (DMET) is an efficient approach that enforces self-consistency at the level of one-particle reduced density matrices (1-RDMs), facilitating applications across diverse quantum systems. However, conventional DMET is constrained by the requirement that the global 1-RDM (low-level descriptor) be an orthogonal projector, limiting flexibility in bath construction and potentially impeding accuracy in strongly correlated regimes. In this work, we introduce a generalized DMET framework in which the low-level descriptor can be an arbitrary 1-RDM and the bath construction is based on optimizing a quantitative criterion related to the maximal disentanglement between different fragments. This yields an alternative yet controllable bath space construction for generic 1-RDMs, lifting a key limitation of conventional DMET. We demonstrate its consistency with conventional DMET in appropriate limits and exploring its implications for bath construction, downfolding (impurity Hamiltonian construction), low-level solvers, and adaptive fragmentation. We expect that this more flexible framework, which leads to several new variants of DMET, can improve the robustness and accuracy of DMET.

arXiv:2503.09881Mar 12, 20250 citations

A unified density matrix functional construction of quantum baths in density matrix embedding theory beyond the mean-field approximation

Sajanthan Sekaran, Oussama Bindech, Emmanuel Fromager

The equivalence in one-electron quantum baths between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [S. Yalouz et al., J. Chem. Phys. 157, 214112 (2022)]. We prove mathematically that the equivalence remains valid when the density matrix is not idempotent anymore, thus allowing for the construction of correlated (one-electron) quantum baths. A density-matrix functional exactification of DMET is derived within the present unified quantum embedding formalism. Numerical examples reveal that the embedding cluster can be quite sensitive to the level of density-matrix functional approximation used for computing the reference density matrix.The equivalence in one-electron quantum baths between the practical implementation of density matrix embedding theory (DMET) and the more recent Householder-transformed density matrix functional embedding theory has been shown previously in the standard but special case where the reference full-size (one-electron reduced) density matrix, from which the bath is constructed, is idempotent [S. Yalouz et al., J. Chem. Phys. 157, 214112 (2022)]. We prove mathematically that the equivalence remains valid when the density matrix is not idempotent anymore, thus allowing for the construction of correlated (one-electron) quantum baths. A density-matrix functional exactification of DMET is derived within the present unified quantum embedding formalism. Numerical examples reveal that the embedding cluster can be quite sensitive to the level of density-matrix functional approximation used for computing the reference density matrix.

arXiv:2304.14729The Journal of Chemical PhysicsJul 19, 202312 citations

On the pure state

v

-representability of density matrix embedding theory

v

-representable (NI-PS-V), which assumes that the low-level density matrix is constructed following the Aufbau principle. In order to relax the NI-PS-V condition, we develop an augmented Lagrangian method to match the density matrix blocks without referring to the Aufbau principle. Numerical results for 2D Hubbard and hydrogen model systems indicate that in some challenging scenarios, the relaxation of the Aufbau principle directly leads to exact matching of the density matrix blocks, which also yields improved accuracy.