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

Projection-based DMRG-in-DFT embedding corrected by non-additive exchange-correlation

Enzo Monino, Daria Drwal, Pavel Beran, Michał Hapka, Libor Veis, Katarzyna Pernal

Published Nov 17, 2025• 0 citations
arXivHTMLPDF

Abstract

The projection-based wave function (WF)-in-DFT embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [Beran, P. et al. J. Phys. Chem. Lett. 2023, 14, 3, 716-722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H20 chain and the cleavage of a triple CN bond in propionitrile.

DOI
10.48550/arxiv.2511.13346

Related papers

Projection-Based DMRG-in-DFT Embedding Corrected by Nonadditive Exchange-Correlation

Enzo Monino, Daria Drwal, Pavel Beran +3 more

The projection-based wave function in density functional theory (WF-in-DFT) embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [ Beran, P. J. Phys. Chem. Lett. 2023, 14, 716–722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H 20 chain and the cleavage of a triple CN bond in propionitrile.

Journal of Chemical Theory and Computation•Feb 9, 2026•3 citations

Projection-Based Density Matrix Renormalization Group in Density Functional Theory Embedding

Pavel Beran, Katarzyna Pernal, Fabijan Pavošević +1 more

The density matrix renormalization group (DMRG) method has already proved itself as a very efficient and accurate computational method, which can treat large active spaces and capture the major part of strong correlation. Its application on larger molecules is, however, limited by its own computational scaling as well as demands of methods for treatment of the missing dynamical electron correlation. In this work, we present the first step in the direction of combining DMRG with density functional theory (DFT), one of the most employed quantum chemical methods with favorable scaling, by means of the projection-based wave function (WF)-in-DFT embedding. On two proof-of-concept but important molecular examples, we demonstrate that the developed DMRG-in-DFT approach provides a very accurate description of molecules with a strongly correlated fragment.

arXiv:2210.16289•The Journal of Physical Chemistry Letters•Jan 17, 2023•8 citations

Why projection-based WF-in-DFT cannot be exact, even with the exact exchange–correlation functional. Formal and practical sources of errors

Enzo Monino, Daria Drwal, Michał Hapka +2 more

We establish theoretical foundations for embedding a correlated wavefunction in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation that equates two, in principle, distinct kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-density functional theory (DFT) formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based density matrix renormalization group-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (pair-density functional theory) instead of a semilocal generalized gradient approximation (GGA) does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.

The Journal of Chemical Physics•Jun 9, 2026•2 citations

Why Projection-Based DMRG-in-DFT Cannot Be Exact, Even with the Exact Exchange-Correlation Functional

Enzo Monino, Daria Drwal, Michał Hapka +2 more

We establish the theoretical foundations for embedding a correlated wave function in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation which equates two, in principle distinct, kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-DFT formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based DMRG-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (PDFT) instead of a semilocal GGA does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.

arXiv:2603.05303•Mar 5, 2026•0 citations

Density Matrix Renormalization Group with Dynamical Correlation via Adiabatic Connection

Pavel Beran, Mikuláš Matoušek, Michał Hapka +2 more

The quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as one of the methods of choice for calculations of strongly correlated molecular systems. Despite its great ability to capture strong electronic correlation in large active spaces, it is not suitable for computations of dynamical electron correlation. In this work, we present a new approach to the electronic structure problem of strongly correlated molecules, in which DMRG is responsible for a proper description of the strong correlation, whereas dynamical correlation is computed via the recently developed adiabatic connection (AC) technique which requires only up to two-body active space reduced density matrices. We report the encouraging results of this approach on typical candidates for DMRG computations, namely, n -acenes ( n = 2 → 7), Fe(II)–porphyrin, and the Fe 3 S 4 cluster.

arXiv:2108.12803•Journal of Chemical Theory and Computation•Nov 11, 2021•51 citations

A Stochastic Cluster Expansion for Electronic Correlation in Large Systems

Annabelle Canestraight, Anthony J. Dominic, Andrés Montoya−Castillo +2 more

Accurate many-body treatments of condensed-phase systems are challenging because correlated solvers such as full configuration interaction (FCI) and the density matrix renormalization group (DMRG) scale exponentially with system size. Downfolding and embedding approaches mitigate this cost but typically require prior selection of a correlated subspace, which can be difficult to determine in heterogeneous or extended systems. Here, we introduce a stochastic cluster expansion framework for efficiently recovering the total correlation energy of large systems with near-DMRG accuracy, without the need to select an active space a priori. By combining correlation contributions from randomly sampled environment orbitals with an exactly treated subspace of interest, the method reproduces total energies for non-reacting and reactive systems while drastically reducing computational cost. The approach also provides a quantitative diagnostic for molecule-solvent correlation, guiding principled embedding decisions. This framework enables systematically improvable many-body calculations in extended systems, opening the door to high-accuracy studies of chemical processes in condensed phase environments.

arXiv:2602.12254•Feb 12, 2026•0 citations
The projection-based wave function (WF)-in-DFT embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [Beran, P. et al. J. Phys. Chem. Lett. 2023, 14, 3, 716-722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H20 chain and the cleavage of a triple CN bond in propionitrile.
The projection-based wave function in density functional theory (WF-in-DFT) embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [ Beran, P. J. Phys. Chem. Lett. 2023, 14, 716–722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H 20 chain and the cleavage of a triple CN bond in propionitrile.
The density matrix renormalization group (DMRG) method has already proved itself as a very efficient and accurate computational method, which can treat large active spaces and capture the major part of strong correlation. Its application on larger molecules is, however, limited by its own computational scaling as well as demands of methods for treatment of the missing dynamical electron correlation. In this work, we present the first step in the direction of combining DMRG with density functional theory (DFT), one of the most employed quantum chemical methods with favorable scaling, by means of the projection-based wave function (WF)-in-DFT embedding. On two proof-of-concept but important molecular examples, we demonstrate that the developed DMRG-in-DFT approach provides a very accurate description of molecules with a strongly correlated fragment.
We establish theoretical foundations for embedding a correlated wavefunction in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation that equates two, in principle, distinct kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-density functional theory (DFT) formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based density matrix renormalization group-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (pair-density functional theory) instead of a semilocal generalized gradient approximation (GGA) does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.
We establish the theoretical foundations for embedding a correlated wave function in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation which equates two, in principle distinct, kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-DFT formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based DMRG-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (PDFT) instead of a semilocal GGA does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.
The quantum chemical version of the density matrix renormalization group (DMRG) method has established itself as one of the methods of choice for calculations of strongly correlated molecular systems. Despite its great ability to capture strong electronic correlation in large active spaces, it is not suitable for computations of dynamical electron correlation. In this work, we present a new approach to the electronic structure problem of strongly correlated molecules, in which DMRG is responsible for a proper description of the strong correlation, whereas dynamical correlation is computed via the recently developed adiabatic connection (AC) technique which requires only up to two-body active space reduced density matrices. We report the encouraging results of this approach on typical candidates for DMRG computations, namely, n -acenes ( n = 2 → 7), Fe(II)–porphyrin, and the Fe 3 S 4 cluster.
Accurate many-body treatments of condensed-phase systems are challenging because correlated solvers such as full configuration interaction (FCI) and the density matrix renormalization group (DMRG) scale exponentially with system size. Downfolding and embedding approaches mitigate this cost but typically require prior selection of a correlated subspace, which can be difficult to determine in heterogeneous or extended systems. Here, we introduce a stochastic cluster expansion framework for efficiently recovering the total correlation energy of large systems with near-DMRG accuracy, without the need to select an active space a priori. By combining correlation contributions from randomly sampled environment orbitals with an exactly treated subspace of interest, the method reproduces total energies for non-reacting and reactive systems while drastically reducing computational cost. The approach also provides a quantitative diagnostic for molecule-solvent correlation, guiding principled embedding decisions. This framework enables systematically improvable many-body calculations in extended systems, opening the door to high-accuracy studies of chemical processes in condensed phase environments.