The S66 noncovalent interactions benchmark reconsidered using explicitly correlated methods near the basis set limit

The S66 benchmark for noncovalent interactions has been re-evaluated using explicitly correlated methods with basis sets near the one-particle basis set limit. It is found that post-MP2 "high-level corrections" are treated adequately well using a combination of CCSD(F12*) with (aug-)cc-pVTZ-F12 basis sets on the one hand, and (T) extrapolated from conventional CCSD(T)/heavy-aug-cc-pV{D,T}Z on the other hand. Implications for earlier benchmarks on the larger S66x8 problem set in particular, and for accurate calculations on noncovalent interactions in general, are discussed. At a slight cost in accuracy, (T) can be considerably accelerated by using sano-V{D,T}Z+ basis sets, while half-counterpoise CCSD(F12*)(T)/cc-pVDZ-F12 offers the best compromise between accuracy and computational cost.Comment: Australian Journal of Chemistry, in press [Graham S. Chandler special issue

The lowest singlet-triplet excitation energy of BN: a converged coupled cluster perspective

The notoriously small $X ^3\Pi-a ^1\Sigma^+$ excitation energy of the BN diatomic has been calculated using high-order coupled cluster methods. Convergence has been established in both the 1-particle basis set and the coupled cluster expansion. Explicit inclusion of connected quadruple excitations $\hat{T}_4$ is required for even semiquantitative agreement with the limit value, while connected quintuple excitations $\hat{T}_5$ still have an effect of about 60 cm$^{-1}$. Still higher excitations only account for about 10 cm$^{-1}$. Inclusion of inner-shell correlation further reduces $T_e$ by about 60 cm$^{-1}$ at the CCSDT, and 85 cm$^{-1}$ at the CCSDTQ level. Our best estimate, $T_e$=183$\pm$40 cm$^{-1}$, is in excellent agreement with earlier calculations and experiment, albeit with a smaller (and conservative) uncertainty. The dissociation energy of BN($X ^3\Pi$) is $D_e$=105.74$\pm$0.16 kcal/mol and $D_0$=103.57$\pm$0.16 kcal/mol.Comment: J. Chem. Phys., in pres

Prototypical pi-pi dimers re-examined by means of high-level CCSDT(Q) composite ab inito methods

The benzene...ethene and parallel-displaced (PD) benzene...benzene dimers are the most fundamental systems involving p-p stacking interactions. Several high-level ab initio investigations calculated the binding energies of these dimers at the CCSD(T)/CBS level of theory using various approaches such as reduced virtual orbital spaces and/or MP2-based basis set corrections. Here we obtain CCSDT(Q) binding energies using a Weizmann-3-type approach. In particular, we extrapolate the SCF, CCSD, and (T) components using large heavy-atom augmented Gaussian basis sets (namely, SCF/jul-cc-pV{5,6}Z, CCSD/jul-cc-pV{Q,5}Z, and (T)/jul-cc-pV{T,Q}Z). We consider post-CCSD(T) contributions up to CCSDT(Q), inner-shell, scalar-relativistic, and Born-Oppenheimer corrections. Overall, our best relativistic, all-electron CCSDT(Q) binding energies are Delta Ee,all,rel = 1.234 (benzene...ethene) and 2.550 (benzene...benzene PD), Delta H0 = 0.949 (benzene...ethene) and 2.310 (benzene...benzene PD), and Delta H298 = 0.130 (benzene...ethene) and 1.461 (benzene...benzene PD) kcal/mol. Important conclusions are reached regarding the basis set convergence of the SCF, CCSD, (T), and post-CCSD(T) components. Explicitly correlated calculations are used as a sanity check on the conventional binding energies. Overall, post-CCSD(T) contributions are destabilizing by 0.028 (benzene...ethene) and 0.058(benzene...benzene) kcal/mol, thus they cannot be neglected if 0.1 kcal/mol accuracy is sought.Comment: J. Chem. Phys., accepted with minor revisio

Numerical model for dynamic installation of large diameter monopiles

publishedVersio

W4 theory for computational thermochemistry: in pursuit of confident sub-kJ/mol predictions

In an attempt to improve on our earlier W3 theory [J. Chem. Phys. {\bf 120}, 4129 (2004)] we consider such refinements as more accurate estimates for the contribution of connected quadruple excitations ($\hat{T}_4$), inclusion of connected quintuple excitations ($\hat{T}_5$), diagonal Born-Oppenheimer corrections (DBOC), and improved basis set extrapolation procedures. Revised experimental data for validation purposes were obtained from the latest version of the ATcT (Active Thermochemical Tables) Thermochemical Network. We found that the CCSDTQ$-$CCSDT(Q) difference converges quite rapidly with the basis set, and that the formula 1.10[CCSDT(Q)/cc-pVTZ+CCSDTQ/cc-pVDZ$-$CCSDT(Q)/cc-pVDZ] offers a very reliable as well as fairly cost-effective estimate of the basis set limit $\hat{T}_4$ contribution. The largest $\hat{T}_5$ contribution found in the present work is on the order of 0.5 kcal/mol (for ozone). DBOC corrections are significant at the 0.1 kcal/mol level in hydride systems. . Based on the accumulated experience, a new computational thermochemistry protocol for first-and second-row main-group systems, to be known as W4 theory, is proposed. Our W4 atomization energies for a number of key species are in excellent agreement (better than 0.1 kcal/mol on average, 95% confidence intervals narrower than 1 kJ/mol) with the latest experimental data obtained from Active Thermochemical Tables. A simple {\em a priori} estimate for the importance of post-CCSD(T) correlation contributions (and hence a pessimistic estimate for the error in a W2-type calculation) is proposed.Comment: J. Chem. Phys., in press; electronic supporting information available at http://theochem.weizmann.ac.il/web/papers/w4.htm

S66x8 Noncovalent Interactions Revisited: New Benchmark and Performance of Composite Localized Coupled-Cluster Methods

The S66x8 noncovalent interactions benchmark has been re-evaluated at the "sterling silver" level, using explicitly correlated MP2-F12 near the complete basis set limit, CCSD(F12*)/aug-cc-pVTZ-F12, and a (T) correction from conventional CCSD(T)/sano-V{D,T}Z+ calculations. The revised reference value disagrees by 0.1 kcal/mol RMS with the original Hobza benchmark and its revision by Brauer et al, but by only 0.04 kcal/mol variety from the "bronze" level data in Kesharwani et al., Aust. J. Chem. 71, 238-248 (2018). We then used these to assess the performance of localized-orbital coupled cluster approaches with and without counterpoise corrections, such as PNO-LCCSD(T) as implemented in MOLPRO, DLPNO-CCSD (T1) as implemented in ORCA, and LNO-CCSD(T) as implemented in MRCC, for their respective "Normal", "Tight", and "very Tight" settings. We also considered composite approaches combining different basis sets and cutoffs. Furthermore, in order to isolate basis set convergence from domain truncation error, for the aug-cc-pVTZ basis set we compared PNO, DLPNO, and LNO approaches with canonical CCSD(T). We conclude that LNO-CCSD(T) with veryTight criteria performs very well for "raw" (CP-uncorrected), but struggles to reproduce counterpoise-corrected numbers even for veryVeryTight criteria: this means that accurate results can be obtained using either extrapolation from basis sets large enough to quench basis set superposition error (BSSE) such as aug-cc-pV{Q,5}Z, or using a composite scheme such as Tight{T,Q}+1.11[vvTight(T) - Tight(T)]. In contrast, PNO-LCCSD(T) works best with counterpoise, while performance with and without counterpoise is comparable for DLPNO-CCSD(T1). Among more economical methods, the highest accuracies are seen for dRPA75-D3BJ, {\omega}B97M-V, {\omega}B97M(2), revDSD-PBEP86-D4, and DFT(SAPT) with a TDEXX or ATDEXX kernel.Comment: Final published version with CC licens

O(^3P) +CO_2 Collisions at Hyperthermal Energies: Dynamics of Nonreactive Scattering, Oxygen Isotope Exchange, and Oxygen-Atom Abstraction

The dynamics of O(^3P) + CO_2 collisions at hyperthermal energies were investigated experimentally and theoretically. Crossed-molecular-beams experiments at Ecoll = 98.8 kcal mol^(–1) were performed with isotopically labeled ^(12)C^(18)O_2 to distinguish products of nonreactive scattering from those of reactive scattering. The following product channels were observed: elastic and inelastic scattering (^(16)O(^3P) + ^(12)C^(18)O^2), isotope exchange (^(18)O + ^(16)O^(12)C^(18)O), and oxygen-atom abstraction (^(18)O^(16)O + ^(12)C^(18)O). Stationary points on the two lowest triplet potential energy surfaces of the O(^3P) + CO_2 system were characterized at the CCSD(T)/aug-cc-pVTZ level of theory and by means of W4 theory, which represents an approximation to the relativistic basis set limit, full-configuration-interaction (FCI) energy. The calculations predict a planar CO_3(C_(2v),^3A″) intermediate that lies 16.3 kcal mol^(–1) (W4 FCI excluding zero point energy) above reactants and is approached by a C_(2v) transition state with energy 24.08 kcal mol^(–1). Quasi-classical trajectory (QCT) calculations with collision energies in the range 23–150 kcal mol^(–1) were performed at the B3LYP/6-311G(d) and BMK/6-311G(d) levels. Both reactive channels observed in the experiment were predicted by these calculations. In the isotope exchange reaction, the experimental center-of-mass (c.m.) angular distribution, T(θ_(c.m.)), of the ^(16)O^(12)C^(18)O products peaked along the initial CO_2 direction (backward relative to the direction of the reagent O atoms), with a smaller isotropic component. The product translational energy distribution, P(E_T), had a relatively low average of E_T = 35 kcal mol^(–1), indicating that the ^(16)O^(12)C^(18)O products were formed with substantial internal energy. The QCT calculations give c.m. P(E_T) and T(θ_(c.m.)) distributions and a relative product yield that agree qualitatively with the experimental results, and the trajectories indicate that exchange occurs through a short-lived CO_3^* intermediate. A low yield for the abstraction reaction was seen in both the experiment and the theory. Experimentally, a fast and weak ^(16)O^(18)O product signal from an abstraction reaction was observed, which could only be detected in the forward direction. A small number of QCT trajectories leading to abstraction were observed to occur primarily via a transient CO_3 intermediate, albeit only at high collision energies (149 kcal mol^(–1)). The oxygen isotope exchange mechanism for CO_2 in collisions with ground state O atoms is a newly discovered pathway through which oxygen isotopes may be cycled in the upper atmosphere, where O(^3P) atoms with hyperthermal translational energies can be generated by photodissociation of O_3 and O_2

Multichannel interference mitigation methods in radio astronomy

Radio-astronomical observations are increasingly corrupted by RF interference, and online detection and filtering algorithms are becoming essential. To facilitate the introduction of such techniques into radio astronomy, we formulate the astronomical problem in an array signal processing language, and give an introduction to some elementary algorithms from that field. We consider two topics in detail: interference detection by rank estimation of short-term covariance matrices, and spatial filtering by subspace estimation and projection. We discuss experimental data collected at the Westerbork radio telescope, and illustrate the effectiveness of the space-time detection and blanking process on the recovery of a 3C48 absorption line in the presence of GSM mobile telephony interference.Comment: 39 pages, 18 figures.Enhanced figures can be downloaded from http://cas.et.tudelft.nl/~leshem/postscripts/leshem/figs34567.ps.gz To appear in Astrophysical Journal Supplements serie