Density Functional Theory Correlation Energy Research Articles

Embracing local suppression and enhancement of dynamic correlation effects in a CASΠDFT method for efficient description of excited states.

The recently proposed CASΠDFT method combines the reliable description of nondynamic electron correlation with the complete active space (CAS) wavefunction and the efficient treatment of dynamic correlation by density functional theory (DFT). This marriage is accomplished by adopting the DFT correlation energy functional modified with the local correction function of the on-top pair density (Π). The role of the correction function is to sensitize the correlation functional to local effects of suppression and enhancement of dynamic correlation and to account for an adequate amount of dynamic correlation energy. In this work we show that the presence of covalent and ionic configurations in a wavefunction gives rise to spatial regions where the effects of suppression and enhancement of correlation energy, respectively, dominate. The results obtained for the potential energy curves of the excited states of the hydrogen molecule prove that CASΠDFT is reliable for states that change their character along the dissociation curve. The method is also applied to the lowest excited states of six-membered heterocyclic nitrogen compounds such as pyridine, pyrazine, pyrimidine, and pyridazine. The obtained excitation energies for the n → π* and π → π* excitations confirm the good performance of CASΠDFT for excited states. The absolute average error of the method is 0.1 eV lower than that of the CCSD method and higher by the same amount than that of the more expansive CC3 variant. Compared with the coupled cluster methods, this encouraging performance of CASΠDFT is achieved at the negligible computational cost of obtaining the correlation energy.

Orbital-dependent correlation energy in density-functional theory based on a second-order perturbation approach: Success and failure

We have developed a second-order perturbation theory (PT) energy functional within density-functional theory (DFT). Based on PT with the Kohn-Sham (KS) determinant as a reference, this new ab initio exchange-correlation functional includes an exact exchange (EXX) energy in the first order and a correlation energy including all single and double excitations from the KS reference in the second order. The explicit dependence of the exchange and correlation energy on the KS orbitals in the functional fits well into our direct minimization approach for the optimized effective potential, which is a very efficient method to perform fully self-consistent calculations for any orbital-dependent functionals. To investigate the quality of the correlation functional, we have applied the method to selected atoms and molecules. For two-electron atoms and small molecules described with small basis sets, this new method provides excellent results, improving both second-order Moller-Plesset expression and any conventional DFT results significantly. For larger systems, however, it performs poorly, converging to very low unphysical total energies. The failure of PT based energy functionals is analyzed, and its origin is traced back to near degeneracy problems due to the orbital- and eigenvalue-dependent algebraic structure of the correlation functional. The failure emerges in the self-consistent approach but not in perturbative post-EXX calculations, emphasizing the crucial importance of self-consistency in testing new orbital-dependent energy functionals.

Adiabatic integration formula for the correlation energy functional of the Hartree-Fock density

An adiabatic integration formula for the quantum chemistry correlation energy functional of the Hartree–Fock density, EcQC[n], is presented. The functional EcQC[n] is meant to be added to the completed Hartree–Fock energy to produce the exact ground-state energy of the system under consideration. The initial slope of the integrand in this connection formula is identified as a second-order energy and an explicit expression for the initial slope of the integrand is presented. Our expression should be useful for arriving at new improved approximations to EcQC[n]. Previous numerical results by Huang and Umrigar (1997) Phys Rev A 56:290, for two-electron densities are proved, and a generalization to more than two electrons is presented. Results obtained by means of the present density functional theory correlation energy functionals, when used to approximate the initial slope in our adiabatic integration formula for EcQC[n], are compared against exact numbers.

Closed-form expression relating the second-order component of the density functional theory correlation energy to its functional derivative

For helping to improve approximations to the density-functional exchange-correlation energy, Exc[n], and its functional derivative, the difference between the second-order component of the correlation energy, Ec(2)[n], and the integral ∫dr vc(2)([n];r)n(r), involving its functional derivative, vc(2)([n];r), is given in terms of only the occupied Kohn–Sham orbitals and the exchange potential. The quantity 2Ec(2)[n] is especially significant because it is the initial slope in the adiabatic connection formula for Exc[n]. The analytic expression for 2Ec(2)[n]−∫dr vc(2)([n];r)n(r) is obtained for any spherically symmetric two-electron test density. Numerical examples are presented.