Basis set converged weak interaction energies from conventional and explicitly correlated coupled-cluster approach

Efficient Density-Fitted Explicitly Correlated Dispersion and Exchange Dispersion Energies

The leading-order dispersion and exchange-dispersion terms in symmetry-adapted perturbation theory (SAPT), E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾, suffer from slow convergence to the complete basis set limit. To alleviate this problem, explicitly correlated variants of these corrections, E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12, have been proposed recently. However, the original formalism (M., Kodrycka , J. Chem. Theory Comput. 2019, 15, 5965-5986), while highly successful in terms of improving convergence, was not competitive to conventional orbital-based SAPT in terms of computational efficiency due to the need to manipulate several kinds of two-electron integrals. In this work, we eliminate this need by decomposing all types of two-electron integrals using robust density fitting. We demonstrate that the error of the density fitting approximation is negligible when standard auxiliary bases such as aug-cc-pVXZ/MP2FIT are employed. The new implementation allowed us to study all complexes in the A24 database in basis sets up to aug-cc-pV5Z, and the E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 values exhibit vastly improved basis set convergence over their conventional counterparts. The well-converged E_disp⁽²⁰⁾-F12 and E_exch-disp⁽²⁰⁾-F12 numbers can be substituted for conventional E_disp⁽²⁰⁾ and E_exch-disp⁽²⁰⁾ ones in a calculation of the total SAPT interaction energy at any level (SAPT0, SAPT2+3, ...). We show that the addition of F12 terms does not improve the accuracy of low-level SAPT treatments. However, when the theory errors are minimized in high-level SAPT approaches such as SAPT2+3(CCD)δMP2, the reduction of basis set incompleteness errors thanks to the F12 treatment substantially improves the accuracy of small-basis calculations.

Computational Investigations of Dispersion Interactions between Small Molecules and Graphene-like Flakes

We investigate dispersion interactions in a selection of atomic, molecular, and molecule-surface systems, comparing high-level correlated methods with empirically corrected density functional theory (DFT). We assess the efficacy of functionals commonly used for surface-based calculations, with and without the D3 correction of Grimme. We find that the inclusion of the correction is essential to get meaningful results, but there is otherwise little to distinguish between the functionals. We also present coupled-cluster quality interaction curves for H₂, NO₂, H₂O, and Ar interacting with large carbon flakes, acting as models for graphene surfaces, using novel absolutely localized molecular orbital based methods. These calculations demonstrate that the problems with empirically corrected DFT when investigating dispersion appear to compound as the system size increases, with important implications for future computational studies of molecule-surface interactions.

Surprisingly broad applicability of the cc-pVnZ-F12 basis set for ground and excited states

Excellent convergence properties for the (aug-)cc-pVnZ-F12 basis set family, purpose-made for explicitly correlated calculations, are demonstrated with conventional wave function methods and Kohn-Sham density functional theory for various ground and excited-state calculations. Among the ground-state properties studied are dipole moments, covalent bond lengths, and interaction and reaction energies. For excited states, we looked at vertical excitation energies, UV absorption, and excited-state absorption spectra. Convergence is compared against the basis sets cc-pVnZ, def2-nVD, aug-pcseg-n, and nZaPa-NR. It is established that the cc-pVnZ-F12 family consistently yields results of n + 1 quality and better. Especially, the cc-pVDZ-F12 basis set is found to be a basis set of good cost vs performance trade-off.

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.

Thermophysical Properties of Gaseous H2S–N2 Mixtures from First-Principles Calculations

The cross second virial coefficient and three dilute gas transport properties (shear viscosity, thermal conductivity, and binary diffusion coefficient) of mixtures of hydrogen sulfide (H2S) and nitrogen (N2) were determined with high accuracy at temperatures up to 1200 K using statistical thermodynamics and the kinetic theory of molecular gases, respectively. The required intermolecular potential energy surface (PES) for the H2S–N2 interaction is presented in this work, while the H2S–H2S and N2–N2 PESs were reported previously. All three PESs are based on high-level quantum-chemical ab initio (i.e. first-principles) calculations. There is only very limited experimental information available on the second virial coefficients of H2S–N2 mixtures, and there appear to be no experimental data at all for the transport properties. Thus, the present predictions constitute a substantial increase in our knowledge of the thermophysical properties of this system, which are of practical relevance for modeling sour natural gas.

A high-level ab initio study of the N2 + N2 reaction channel

A new six-dimensional (6D) global potential energy surface (PES) is proposed for the full range description of the interaction of the N2(1Σg+)+N2(1Σg+) system governing collisional processes, including N atom exchange. The related potential energy values were determined using high-level ab initio methods. The calculations were performed at a coupled-cluster with single and double and perturbative triple excitations level of theory in order to have a first full range picture of the PES. Subsequently, in order to accurately describe the stretching of the bonds of the two interacting N2 molecules by releasing the constraints of being considered as rigid rotors, for the same molecular geometries higher level of theory multi reference calculations were performed. Out of the calculated values a 6D 4-atoms global PES was produced for use in dynamical calculations. The ab initio calculations were made possible by the combined use of High Throughput Computing and High Performance Computing techniques within the frame of a computing grid empowered molecular simulator.