Insight into organic reactions from the direct random phase approximation and its corrections

Self-Interaction-Corrected Random Phase Approximation

The short-range correlation energy of the random phase approximation (RPA) is too negative and is often corrected by local or nonlocal methods. These beyond-RPA corrections usually lead to a mixed performance for thermodynamics and dissociation properties. RPA+ is an additive correction based on density functional approximations that often gives realistic total energies for atoms or solids. RPA+ adds a moderate correction to the ionization energies/electron affinities of RPA but does not yield an improvement beyond RPA for atomization energies of molecules. This incompleteness results in severely underestimated atomization energies just like in RPA. Exchange-correlation kernels within the Dyson equation could simultaneously improve atomization, ionization energies, and electron affinities, but their implementation is computationally less feasible in localized basis set codes. In preceding work ( Phys. Rev. A 100, 2019022515), two of the authors proposed a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities. gRPA+ was found to yield a large improvement of ionization energies and electron affinities of light atoms over RPA, and a smaller improvement over RPA+. Within this work, we investigate to what extent this improvement transfers to atomization energies, ionization energies, and electron affinities of molecules, using a modified gRPA+ (mgRPA+) method that can be applied in codes with localized basis functions. We thereby aim to understand the applicability of beyond-RPA corrections based on density functional approximations.

Range-Separated Density-Functional Theory in Combination with the Random Phase Approximation: An Accuracy Benchmark

A formulation of range-separated random phase approximation (RPA) based on our efficient ω-CDGD-RI-RPA [J. Chem. Theory Comput. 2018, 14, 2505] method and a large scale benchmark study are presented. By application to the GMTKN55 data set, we obtain a comprehensive picture of the performance of range-separated RPA in general main group thermochemistry, kinetics, and noncovalent interactions. The results show that range-separated RPA performs stably over the broad range of molecular chemistry included in the GMTKN55 set. It improves significantly over semilocal DFT but it is still less accurate than modern dispersion corrected double-hybrid functionals. Furthermore, range-separated RPA shows a faster basis set convergence compared to standard full-range RPA making it a promising applicable approach with only one empirical parameter.

Shrinking Self-Interaction Errors with the Fermi-Löwdin Orbital Self-Interaction-Corrected Density Functional Approximation

The self-interaction error (SIE) is one of the major drawbacks of practical exchange-correlation functionals for Kohn-Sham density functional theory. Despite this, the use of methods that explicitly remove SIE from approximate density functionals is scarce in the literature due to their relatively high computational cost and lack of consistent improvement over standard modern functionals. In this article we assess the performance of a novel approach recently proposed by Pederson, Ruzsinszky, and Perdew [ J. Chem. Phys. 2014, 140, 121103] for performing self-interaction free calculations in density functional theory based on Fermi orbitals. To this end, we employ test sets consisting of reaction energies that are considered particularly sensitive to SIE. We found that the parameter-free Fermi-Löwdin orbital self-interaction correction method combined with the standard local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof (PBE) functionals gives a much better estimate of reaction energies compared to their parent LSDA and PBE functionals for most of the reactions in these two sets. They also perform on par with the global PBE0 and range-separated LC-ωPBE hybrids, which partially eliminate the SIE by including Hartree-Fock exchange. This shows the potential of the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for practical density functional calculations without SIE.

Performance and Scope of Perturbative Corrections to Random-Phase Approximation Energies

It has been suspected since the early days of the random-phase approximation (RPA) that corrections to RPA correlation energies result mostly from short-range correlation effects and are thus amenable to perturbation theory. Here we test this hypothesis by analyzing formal and numerical results for the most common beyond-RPA perturbative corrections, including the bare second-order exchange (SOX), second-order screened exchange (SOSEX), and approximate exchange kernel (AXK) methods. Our analysis is facilitated by efficient and robust algorithms based on the resolution-of-the-identity (RI) approximation and numerical frequency integration, which enable benchmark beyond-RPA calculations on medium- and large-size molecules with size-independent accuracy. The AXK method systematically improves upon RPA, SOX, and SOSEX for reaction barrier heights, reaction energies, and noncovalent interaction energies of main-group compounds. The improved accuracy of AXK compared with SOX and SOSEX is attributed to stronger screening of bare SOX in AXK. For reactions involving transition-metal compounds, particularly 3d transition-metal dimers, the AXK correction is too small and can even have the wrong sign. These observations are rationalized by a measure α̅ of the effective coupling strength for beyond-RPA correlation. When the effective coupling strength increases beyond a critical α̅ value of approximately 0.5, the RPA errors increase rapidly and perturbative corrections become unreliable. Thus, perturbation theory can systematically correct RPA but only for systems and properties qualitatively well captured by RPA, as indicated by small α̅ values.

Beyond energies: geometry predictions with the XYG3 type of doubly hybrid density functionals

The scaled mean absolute deviations (s-MADs) of the optimized geometric parameters for covalent bondings (the CCse set), nonbonded interactions (the S22G30 set) and the transition state structures (the TSG36 set), with Tot referring to the averaged s-MAD for general performances.