Accurately treating strong electron correlation in quantum chemistry typically requires multireference wave-function methods with steep computational scaling. The recently proposed i-DMFT method promises near configuration-interaction accuracy at mean-field cost by invoking an empirical linear relation between correlation energy and entropy (Collins' conjecture), whose validity remains unclear. We systematically assess this relation across a range of di- and polyatomic molecules, including diverse bond types, third-row elements, different types of geometric distortions, and excited states. We find that the conjectured linearity holds for bond-breaking processes dominated by electron redistribution within orbital pairs, but breaks down for heterolytic dissociation and excited states. In simple molecules, i-DMFT provides a reasonable description of total energies, but does not reliably reproduce reduced density matrices or individual energy components. It further degrades in more complex cases such as ethylene. Based on these results, we formulate criteria for the validity of the conjecture and outline implications for entropy-based reduced density matrix functionals.
Accurately treating strong electron correlation in quantum chemistry typically requires multireference wave-function methods with steep computational scaling. The recently proposed i-DMFT method promises near configuration-interaction accuracy at mean-field cost by invoking an empirical linear relation between correlation energy and entropy (Collins' conjecture), whose validity remains unclear. We systematically assess this relation across a range of di- and polyatomic molecules, including diverse bond types, third-row elements, different types of geometric distortions, and excited states. We find that the conjectured linearity holds for bond-breaking processes dominated by electron redistribution within orbital pairs, but breaks down for heterolytic dissociation and excited states. In simple molecules, i-DMFT provides a reasonable description of total energies, but does not reliably reproduce reduced density matrices or individual energy components. It further degrades in more complex cases such as ethylene. Based on these results, we formulate criteria for the validity of the conjecture and outline implications for entropy-based reduced density matrix functionals.
Accurately describing chemical bond breaking requires capturing strong/static electron correlation, a major challenge for single-reference self-consistent-field methods such as Hartree--Fock and Kohn--Sham density-functional theory (DFT). Fractional occupations substantially improve this situation, as demonstrated by information-density-matrix-functional theory (i-DMFT) and thermally-assisted-occupation DFT (TAO-DFT). Although useful, they are formulated within rather different theoretical frameworks, are somewhat semi-empirical, and their relationship remains unclear. Here we introduce finite-temperature adiabatic-connection reduced density-matrix functional theory (FT-AC-RDMFT), in which i-DMFT and TAO-DFT emerge as different reductions of a unified framework. The effective eigenvalue equation for natural orbitals and the Fermi--Dirac occupation numbers follow directly from the variational principle, while a thermal adiabatic connection relates the empirical entropy coefficient of i-DMFT to a characteristic temperature of the auxiliary non-interacting system. Within this framework, i-DMFT is the exact-exchange-based reduction, whereas TAO-DFT results when exact exchange is replaced by a density-based exchange-correlation functional. Numerical calculations confirm that i-DMFT, apart from a constant, coincides with the exact-exchange realization of the framework, and show that the improved dissociation behavior from finite-temperature occupations extends beyond exact exchange to hybrid and semilocal density functionals. For strongly correlated V2O3, finite-temperature occupations substantially reduce the severely overestimated non-magnetic LDA bulk modulus, bringing it close to experiment. These results establish FT-AC-RDMFT as a common theoretical foundation for finite-temperature occupation-based electronic-structure methods in molecular and solid-state systems.
Ab-initio calculations are dominated by density functional theory (DFT), which provides an efficient and accurate description of the electronic structure for most materials. For materials with strong correlations, however, many of the available density functionals yield poor results. Most well known is the case of transition metal oxides, for which most density functionals produce a qualitatively incorrect description. There is a quest to improve the description by borrowing from methods specifically designed for strongly correlated materials. One of the method which is widely used to study strongly correlated electron system in the thermodynamic limit is the dynamical mean-field theory (DMFT). For calculations pertaining to real materials, hybrid approaches like DFT+DMFT are commonly used. We consider reduced density-matrix functional theory (rDMFT) to be a useful framework for a rigorous formulation of such hybrid theories. The link from rDMFT to many-particle wave functions has been established by Levy's constrained-search algorithm on the one hand. The link to many-body perturbation theory and Green's function, on the other hand, has been provided via the Luttinger-Ward functional. The development in the field of rDMFT proceeded analogously to the development of DFT. In order to avoid the full complexity of an explicit many-body description, most density-matrix functionals are not extracted from the exact expressions i.e. via Levy's constrained-search or Green's function based methods . Rather, one proceeds analogously to the development of density functionals, namely by searching models for the density-matrix functional, that capture the most essential physical effects while having an algebraic dependence on the density matrix. Among such models, the most prominent ones are the Müller functional, Sharma functional and Marques-Lathiotakis functional. One of the major arguments in favor of density-matrix functionals is that one of the most simple functionals, the Müller functional, seems to provide a correct description of the bond-dissociation problem, for which common density functionals fail. Despite the successes, the common density-matrix functionals also reproduce a number of features in a qualitatively incorrect manner. Several such discrepancies of the standard state of the art density-matrix functionals are discussed in details in this thesis. The objective of this thesis is to study the performance of the class of commonly used model density-matrix functionals and to set the stage for the desire to develop the framework for an approximate scheme to evaluate new class of functionals from exact Green's function based method.
In [Phys. Rev. Lett. 128, 013001 (2022)] a novel ground state method was proposed. It has been suggested that this i-DMFT would be a method within one-particle reduced density matrix functional theory (DMFT), capable of describing accurately molecules at various geometries with an information-theoretical nature. We reassess this work and its suggestions from a conceptual and practical point of view, leading to the following conclusions: i) A method which assigns to each molecule M its own functional FM is not a functional theory (striking violation of "universality") ii) even for the simplest systems i-DMFT yields incorrect one-particle reduced density matrices and iii) the use of an information-theoretical concept to describe molecular dissociation limits was not essential. The latter insight may help to fix the deficiency of i-DMFT to not reproduce correctly the smaller occupation numbers and thus to not recover the important dynamic correlations.
The accurate calculation of the (differential) correlation energy is central to the quantum chemical description of bond-formation and bond-dissociation processes. In order to estimate the quality of single- and multireference approaches for this purpose, various diagnostic tools have been developed. In this work, we elaborate on our previous observation [J. Phys. Chem. Lett.2012, 3, 3129] that one- and two-orbital-based entanglement measures provide quantitative means for the assessment and classification of electron correlation effects among molecular orbitals. The dissociation behavior of some prototypical diatomic molecules features all types of correlation effects relevant for chemical bonding. We demonstrate that our entanglement analysis is convenient to dissect these electron correlation effects and to provide a conceptual understanding of bond-forming and bond-breaking processes from the point of view of quantum information theory.
The accuracy of large-scale ab initio calculations is always desired, but they are unfortunately not usually practical in many fields, such as studies of heterogeneous catalysis, photosynthesis, and functional mechanisms of biomolecules. DFT serves as a prevailing method for large calculations due to its compromise between cost and accuracy, but it has a number of failure cases, like bond breaking. On the other hand, the simple theoretical nature of force-field simulations, which often work well, reminds us that large systems are not necessarily exponentially complicated. A recurring theme in the electronic structure of large systems is therefore the use of fragment-based approaches. Existing fragment approaches typically treat the most powerful electron correlations perturbatively, on the same footing with dispersion forces (if at all). In our approach, electrons are not treated individually in the global calculation, but subunits will be chosen. (e.g., amino acid residues or individual water molecules). The most important local states are selected to account for inter-subunit interactions. Most familiar wavefunction methods (SCF, PT, coupled-cluster) can be applied directly to this internally correlated Hamiltonian, and the Hamiltonian itself is systematically improvable to exactitude. Furthermore, different parts of a system can be treated at different accuracy. Our initial investigations focus on the correlation structure of groups of harmonic oscillators, where the coupling inside of groups is stronger than between groups. This mimics discrete molecules while providing us with exact analytical solutions for comparison. Correlations between grouped oscillators will be characterized by both numerical and analytical methods, in order to shed light on building an algorithm for real molecules and systems.
In materials with strong electron correlation the proper treatment of local atomic physics described by orbital occupations is crucial. Reduced density-matrix functional theory is a natural extension of density functional theory for systems that are dominated by orbital physics. We review the current state of reduced density-matrix functional theory (RDMFT). For atomic structure relaxations or ab-initio molecular dynamics the combination of density functional theory (DFT) and dynamical mean-field theory (DMFT) possesses a number of disadvantages, like the cumbersome evaluation of forces. We therefore describe a method, DFT+RDMFT, that combines many-particle effects based on reduced density-matrix functional theory with a density functional-like framework. A recent development is the construction of density-matrix functionals directly from many-particle theory such as methods from quantum chemistry or many-particle Green’s functions. We present the underlying exact theorems and describe current progress towards quantitative functionals.