On the accuracy of explicitly correlated coupled-cluster interaction energies — have orbital results been beaten yet?

Cross Second Virial Coefficients of the H2O-H2 and H2S-H2 Systems from First-Principles

The cross second virial coefficients B12 for the interactions of water (H2O) with molecular hydrogen (H2) and of hydrogen sulfide (H2S) with H2 were obtained at temperatures in the range from 150 to 2000 K from new intermolecular potential energy surfaces (PESs) for the respective molecule pairs. The PESs are based on interaction energies determined for about 12 000 configurations of each molecule pair employing different high-level quantum-chemical ab initio methods up to coupled cluster with single, double, triple, and perturbative quadruple excitations [CCSDT(Q)]. Furthermore, the interaction energies were corrected for scalar relativistic effects. Both classical and semiclassical values for B12 were extracted from the PESs using the Mayer-sampling Monte Carlo approach. While our results for the H2O-H2 system validate the older first-principles results of Hodges et al. [J. Chem. Phys. 2004, 120, 710-720], B12 for the H2S-H2 system was, to the best of our knowledge, hitherto neither measured experimentally nor predicted from first principles.

Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12

We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol^-1, and this maximum error is reduced to below 0.2 kcal mol^-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol^-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol^-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.

Explicitly Correlated Dispersion and Exchange Dispersion Energies in Symmetry-Adapted Perturbation Theory

The individual interaction energy terms in symmetry-adapted perturbation theory (SAPT) not only have different physical interpretations but also converge to their complete basis set (CBS) limit values at quite different rates. Dispersion energy is notoriously the slowest converging interaction energy contribution, and exchange dispersion energy, while smaller in absolute value, converges just as poorly in relative terms. To speed up the basis set convergence of the lowest-order SAPT dispersion and exchange dispersion energies, we borrow the techniques from explicitly correlated (F12) electronic structure theory and develop practical expressions for the closed-shell E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 contributions. While the latter term has been derived and implemented for the first time, the former correction was recently proposed by Przybytek [ J. Chem. Theory Comput. 2018 , 14 , 5105 - 5117 ] using an Ansatz with a full optimization of the explicitly correlated amplitudes. In addition to reimplementing the fully optimized variant of E_disp⁽²⁰⁾-F12, we propose three approximate Ansätze that substantially improve the scaling of the method and at the same time avoid the numerical instabilities of the unrestricted optimization. The performance of all four resulting flavors of E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 is first tested on helium, neon, argon, water, and methane dimers, with orbital and auxiliary basis sets up to aug-cc-pV5Z and aug-cc-pV5Z-RI, respectively. The double- and triple-ζ basis set calculations are then extended to the entire A24 database of noncovalent interaction energies and compared with CBS estimates for E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾ computed using conventional SAPT with basis sets up to aug-cc-pV6Z with midbond functions. It is shown that the F12 treatment is highly successful in improving the basis set convergence of the SAPT terms, with the F12 calculations in an X-tuple ζ basis about as accurate as conventional calculations in bases with cardinal numbers (X + 2) for E_disp⁽²⁰⁾ and either (X + 1) or (X + 2) for E_exch-disp⁽²⁰⁾. While the full amplitude optimization affords the highest accuracy for both corrections, the much simpler and numerically stable optimized diagonal Ansatz is a very close second. We have also tested the performance of the simple F12 correction based on the second-order Møller-Plesset perturbation theory, SAPT-F12(MP2) [ Frey , J. A. ; Chem. Rev. 2016 , 116 , 5614 - 5641 ] and observed that it is also quite successful in speeding up the basis set convergence of conventional E_disp⁽²⁰⁾ + E_exch-disp⁽²⁰⁾, albeit with some outliers.

Ab initio studies of the ground and first excited states of the Sr-H2 and Yb-H2 complexes

Accurate intermolecular potential-energy surfaces (IPESs) for the ground and first excited states of the Sr-H₂ and Yb-H₂ complexes were calculated. After an extensive methodological study, the coupled cluster with single, double, and non-iterative triple excitation method with the Douglas-Kroll-Hess Hamiltonian and correlation-consistent basis sets of triple-ζ quality extended with 2 sets of diffuse functions and a set of midbond functions were chosen. The obtained ground-state IPESs are similar in both complexes, being relatively isotropic with two minima and two transition states (equivalent by symmetry). The global minima correspond to the collinear geometries with R = 5.45 and 5.10 Å and energies of -27.7 and -31.7 cm^-1 for the Sr-H₂ and Yb-H₂ systems, respectively. The calculated surfaces for the Sr(³P)-H₂ and Yb(³P)-H₂ states are deeper and more anisotropic, and they exhibit similar patterns within both complexes. The deepest surfaces, where the singly occupied p-orbital of the metal atom is perpendicular to the intermolecular axis, are characterised by the global minima of ca. -2053 and -2260 cm^-1 in the T-shape geometries at R = 2.41 and 2.29 Å for Sr-H₂ and Yb-H₂, respectively. Additional calculations for the complexes of Sr and Yb with the He atom revealed a similar, strong dependence of the interaction energy on the orientation of the p-orbital in the Sr(³P)-He and Yb(³P)-He states.

Development of a potential energy surface for the O3–Ar system: rovibrational states of the complex

Collisional stabilization is an important step in the process of atmospheric formation of ozone.

Ab initio study of the CO–N2 complex: a new highly accurate intermolecular potential energy surface and rovibrational spectrum

We present a highly accurate ab initio intermolecular potential-energy surface and rovibrational spectrum for the CO–N₂ complex.

CBS extrapolation in electronic structure pushed to the end: a revival of minimal and sub-minimal basis sets

The complete basis set (CBS) limit is secluded in calculations of electronic structure, and hence CBS extrapolation draws immediate attention.

A theoretical benchmark study of the spectroscopic constants of the very heavy rare gas dimers

The binding energy of the superheavy dimer Uuo₂ is considerably larger than that of its lighter homologues, despite a 40% reduction due to spin-other orbit interaction.

An accurate benchmark description of the interactions between carbon dioxide and polyheterocyclic aromatic compounds containing nitrogen

We assessed the performance of a large variety of modern density functional theory approaches for the adsorption of carbon dioxide on molecular models of pyridinic N-doped graphene.