It is common wisdom that a better understanding of how to quantify electron correlation should lead to more efficient treatments of the correlation problem. Here, we not only substantiate this view but turn it into a concrete, systematic framework for ground states. Since prominent correlation measures such as the nonfreeness depend only on the one-particle reduced density matrix, modeling the correlation energy as a function of such a measure naturally places the framework within reduced density matrix functional theory, while keeping the cost at mean-field level. Strikingly, the variable conjugate to the measure emerges as a ``correlation temperature'', thereby connecting correlation measures, functional theory, and effective temperatures. The framework rationalizes methods such as i-DMFT, and we illustrate its key concepts analytically for the Hubbard dimer and numerically for small molecular systems. This brings a fundamental question within reach: how much of the correlation problem can be solved at mean-field cost?
It is common wisdom that a better understanding of how to quantify electron correlation should lead to more efficient treatments of the correlation problem. Here, we not only substantiate this view but turn it into a concrete, systematic framework for ground states. Since prominent correlation measures such as the nonfreeness depend only on the one-particle reduced density matrix, modeling the correlation energy as a function of such a measure naturally places the framework within reduced density matrix functional theory, while keeping the cost at mean-field level. Strikingly, the variable conjugate to the measure emerges as a ``correlation temperature'', thereby connecting correlation measures, functional theory, and effective temperatures. The framework rationalizes methods such as i-DMFT, and we illustrate its key concepts analytically for the Hubbard dimer and numerically for small molecular systems. This brings a fundamental question within reach: how much of the correlation problem can be solved at mean-field cost?
Capturing static and dynamic electron correlation simultaneously at manageable computational cost remains a key challenge in electronic structure theory. Motivated by the recurring role of effective temperature in self-consistent field methods we extend the adiabatic connection mechanism from density functional theory to include a temperature dimension. This gives rise to thermally enhanced hybrid functionals, in which static correlation is captured via fractional occupations generated by an effective "correlation temperature" within thermal Hartree-Fock, while dynamic correlation is described by established density functional approximations. The framework formally unifies information-density matrix functional theory (i-DMFT) and thermally-assisted-occupation (TAO-)DFT as special cases, clarifying their theoretical basis and the respective prescriptions for the temperature parameter. We assess the performance of several thermally enhanced hybrids for different cases where i-DMFT was previously shown to fail. Furthermore, we discuss and test different schemes for determining the correlation temperature and investigate the accuracy of different such choices for a benchmark set derived from biphenyl. The results help to clarify the achievable balance of accuracy, generality, and efficiency, establishing thermally enhanced hybrids as a systematic route toward low-cost treatment of both correlation types.
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.
Electron correlation plays an important role in determining the properties of physical systems, from atoms and molecules to condensed phases. Recent theoretical progress in the field has proven possible using both analytical methods and numerical many-body treatments, for realistic systems as well as for simplified models. Within the models, one may mention the jellium and the Hubbard. The jellium model, while providing a simple, rough approximation to conduction electrons in metals, also constitutes a key ingredient in the treatment of electrons in condensed phases within density-functional formalism. The Hubbard- and the related Heisenberg-model Hamiltonians, on the other hand, are designed to treat situations in which very strong correlations tend to bring about site localization of electrons. The character of the interactions in these lattice models allows for a local treatment of correlations. This is achieved by the use of projection techniques that were first proposed by Gutzwiller for the multicenter problem, being the natural extension of the Coulson-Fischer treatment of the H2 molecule. Much work in this area is analytic or semianalytic and requires approximations. However, a full many-body treatment of both realistic and simplified models is possible by resorting to numerical simulations, i.e., to the so-called quantum Monte Carlo method. This method, which can be implemented in a number of ways, has been applied to atoms, molecules, and solids. In spite of continuing progress, technical problems still remain. Thus one may mention the fermion sign problem and the increase in computational time with the nuclear charge in atomic and related situations. Still, this method provides, to date, one of the most accurate ways to calculate correlation energies, both in atomic and in multicenter problems.
Electron correlation plays an important role in determining the properties of physical systems, from atoms and molecules to condensed phases. Recent theoretical progress in the field has proven possible using both analytical methods and numerical many-body treatments, for realistic systems as well as for simplified models. Within the models, one may mention the jellium and the Hubbard. The jellium model, while providing a simple, rough approximation to conduction electrons in metals, also constitutes a key ingredient in the treatment of electrons in condensed phases within density-functional formalism. The Hubbard- and the related Heisenberg-model Hamiltonians, on the other hand, are designed to treat situations in which very strong correlations tend to bring about site localization of electrons. The character of the interactions in these lattice models allows for a local treatment of correlations. This is achieved by the use of projection techniques that were first proposed by Gutzwiller for the multicenter problem, being the natural extension of the Coulson-Fischer treatment of the H2 molecule. Much work in this area is analytic or semianalytic and requires approximations. However, a full many-body treatment of both realistic and simplified models is possible by resorting to numerical simulations, i.e., to the so-called quantum Monte Carlo method. This method, which can be implemented in a number of ways, has been applied to atoms, molecules, and solids. In spite of continuing progress, technical problems still remain. Thus one may mention the fermion sign problem and the increase in computational time with the nuclear charge in atomic and related situations. Still, this method provides, to date, one of the most accurate ways to calculate correlation energies, both in atomic and in multicenter problems.
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.
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.