Random-Phase Approximation in Many-Body Noncovalent Systems: Methane in a Dodecahedral Water Cage

Contributions beyond direct random-phase approximation in the binding energy of solid ethane, ethylene, and acetylene

The random-phase approximation (RPA) includes a subset of higher than second-order correlation-energy contributions, but stays in the same complexity class as the second-order Møller-Plesset perturbation theory (MP2) in both Gaussian-orbital and plane-wave codes. This makes RPA a promising ab initio electronic structure approach for the binding energies of molecular crystals. Still, some issues stand out in practical applications of RPA. Notably, compact clusters of nonpolar molecules are poorly described, and the interaction energies strongly depend on the reference single-determinant state. Using the many-body expansion of the binding energy of a crystal, we investigate those issues and the effect of beyond-RPA corrections. We find the beneficial effect of quartic-scaling exchange and non-ring coupled-cluster doubles corrections. The nonadditive interactions in compact trimers of molecules are improved by using the self-consistent Hartree-Fock orbitals instead of the usual Kohn-Sham states, but this kind of orbital input also leads to underestimated dimer energies. Overall, a substantial improvement over the RPA with a renormalized singles approach is possible at a modest quartic-scaling cost, which encourages further research into additional RPA corrections.

Post-Kohn-Sham Random-Phase Approximation and Correction Terms in the Expectation-Value Coupled-Cluster Formulation

Using expectation-value coupled-cluster theory and many-body perturbation theory (MBPT), we formulate a series of corrections to the post-Kohn-Sham (post-KS) random-phase approximation (RPA) energy. The beyond-RPA terms are of two types: those accounting for the non-Hartree-Fock reference and those introducing the coupled-cluster doubles non-ring contractions. The contributions of the former type, introduced via the semicanonical orbital basis, drastically reduce the binding strength in noncovalent systems. The good accuracy is recovered by the attractive third-order doubles correction referred to as Ec2g. The existing RPA approaches based on KS orbitals neglect most of the proposed corrections but can perform well thanks to error cancellation. The proposed method accounts for every contribution in the state-of-the-art renormalized second-order perturbation theory (rPT2) approach but adds additional terms which initially contribute in the third order of MBPT. The cost of energy evaluation scales as noniterative O ( N 4 ) in the implementation with low-rank tensor decomposition. The numerical tests of the proposed approach demonstrate accurate results for noncovalent dimers of polar molecules and for the challenging many-body noncovalent cluster of CH4···(H2O)20.

Efficient Calculation of the Dispersion Energy for Multireference Systems with Cholesky Decomposition: Application to Excited-State Interactions

Accurate and efficient prediction of dispersion interactions in excited-state complexes poses a challenge due to the complex nature of electron correlation effects that need to be simultaneously considered. We propose an algorithm for computing the dispersion energy in nondegenerate ground- or excited-state complexes with arbitrary spin. The algorithm scales with the fifth power of the system size due to employing Cholesky decomposition of Coulomb integrals and a recently developed recursive formula for density response functions of the monomers. As a numerical illustration, we apply the new algorithm in the framework of multiconfigurational symmetry adapted perturbation theory, SAPT(MC), to study interactions in dimers with localized excitons. The SAPT(MC) analysis reveals that the dispersion energy may be the main force stabilizing excited-state dimers.

Efficient Adiabatic Connection Approach for Strongly Correlated Systems: Application to Singlet-Triplet Gaps of Biradicals

Strong electron correlation can be captured with multireference wave function methods, but an accurate description of the electronic structure requires accounting for the dynamic correlation, which they miss. In this work, a new approach for the correlation energy based on the adiabatic connection (AC) is proposed. The AC_n method accounts for terms up to order n in the coupling constant, and it is size-consistent and free from instabilities. It employs the multireference random phase approximation and the Cholesky decomposition technique, leading to a computational cost growing with the fifth power of the system size. Because of the dependence on only one- and two-electron reduced density matrices, AC_n is more efficient than existing ab initio multireference dynamic correlation methods. AC_n affords excellent results for singlet-triplet gaps of challenging organic biradicals. The development presented in this work opens new perspectives for accurate calculations of systems with dozens of strongly correlated electrons.

Energies, structures, and harmonic frequencies of small water clusters from the direct random phase approximation

The binding energies, structures, and vibrational frequencies of water clusters up to 20 molecules are computed at the direct random phase approximation (RPA) level of theory and compared to theoretical benchmarks. Binding energies of the WATER27 set, which includes neutral and positively and negatively charged clusters, are predicted to be too low in the complete basis set limit by an average of 7 kcal/mol (9%) and are worse than the results from the best density functional theory methods or from the Møller-Plesset theory. The RPA shows significant basis set size dependence for binding energies. The order of the relative energies of the water hexamer and dodecamer isomers is predicted correctly by the RPA. The mean absolute deviation for angles and distances for neutral clusters up to the water hexamer are 0.2° and 0.6 pm, respectively, using quintuple-ζ basis sets. The relative energetic order of the hexamer isomers is preserved upon optimization. Vibrational frequencies for these systems are underestimated by several tens of wavenumbers for large basis sets, and deviations increase with the basis set size. Overall, the direct RPA method yields accurate structural parameters but systematically underestimates binding energies and shows strong basis set size dependence.