Application of a convergent, composite coupled cluster approach to bound state, adiabatic electron affinities in atoms and small molecules

Sub 20 cm-1 computational prediction of the CH bond energy - a case of systematic error in computational thermochemistry

The bond dissociation energy of methylidyne, D0(CH), is studied using an improved version of the High-Accuracy Extrapolated ab initio Thermochemistry (HEAT) approach as well as the Feller-Peterson-Dixon (FPD) model chemistry. These calculations, which include basis sets up to nonuple (aug-cc-pCV9Z) quality, are expected to be capable of providing results substantially more accurate than the ca. 1 kJ mol-1 level that is characteristic of standard high-accuracy protocols for computational thermochemistry. The calculated 0 K CH bond energy (27 954 ± 15 cm-1 for HEAT and 27 956 ± 15 cm-1 for FPD), along with equivalent treatments of the CH ionization energy and the CH+ dissociation energy (85 829 ± 15 cm-1 and 32 946 ± 15 cm-1, respectively), were compared to the existing benchmarks from Active Thermochemical Tables (ATcT), uncovering an unexpected difference for D0(CH). This has prompted a detailed reexamination of the provenance of the corresponding ATcT benchmark, allowing the discovery and subsequent correction of a systematic error present in several published high-level calculations, ultimately yielding an amended ATcT benchmark for D0(CH). Finally, the current theoretical results were added to the ATcT Thermochemical Network, producing refined ATcT estimates of 27 957.3 ± 6.0 cm-1 for D0(CH), 32 946.7 ± 0.6 cm-1 for D0(CH+), and 85 831.0 ± 6.0 cm-1 for IE(CH).

Elaborated thermochemical treatment of HF, CO, N2, and H2O: Insight into HEAT and its extensions

Empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N₂, and H₂O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm^-1 (0.12 kJ mol^-1), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol-where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure-this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

New Insights into the Remarkable Difference between CH5- and SiH5. J Phys Chem A 2021;125:7414-7424. [PMID: 34424705 DOI: 10.1021/acs.jpca.1c05357]

It has long been known that there is a fundamental difference in the electronic structures of CH₅^- and SiH₅^-, two isoelectronic molecules. The former is a saddle point for the S_N2 exchange reaction H^- + CH₄ → [CH₅^-]^‡ → CH₄ + H^-, while the latter is a stable molecule that is bound relative to SiH₄ + H^-. SCGVB calculations indicate that this difference is the result of a dramatic difference in the nature of the axial electron pairs in these anions. In SiH₅^-, the axial pairs represent two stable bonds-a weak recoupled pair bond dyad. In CH₅^-, the axial electron pairs represent an intermediate transition between the electron pairs in the reactants and those in the products. Furthermore, the axial orbitals at the saddle point in CH₅^- are highly overlapping, giving rise to strong Pauli repulsion and a high barrier for the S_N2 exchange reaction.

Ab initio dynamics of hydrogen abstraction from N2H4 by OH radicals: an RRKM-based master equation study

The detailed kinetic mechanism of the N₂H₄ + OH reaction is comprehensively reported for a wide condition range of conditions (i.e., 200–3000 K & 1–7600 Torr) using the CCSD(T)/CBS//M06-2X/6-311++G(3df,2p) level and the RRKM-based master equation rate model.

EOM-CC guide to Fock-space travel: the C2 edition

Electronic structure calculations for C₂, C₂⁻, and C₂²⁻ using the CC/EOM-CC family of methods. Results illustrate that EOM-CCSD provides an attractive alternative to MR approaches.

A theoretical study of the adiabatic and vertical ionization potentials of water

Theoretical predictions of the three lowest adiabatic and vertical ionization potentials of water were obtained from the Feller-Peterson-Dixon approach. This approach combines multiple levels of coupled cluster theory with basis sets as large as aug-cc-pV8Z in some cases and various corrections up to and including full configuration interaction theory. While agreement with experiment for the adiabatic ionization potential of the lowest energy ²B₁ state was excellent, differences for other states were much larger, sometimes exceeding 10 kcal/mol (0.43 eV). Errors of this magnitude are inconsistent with previous benchmark work on 52 adiabatic ionization potentials, where a root mean square of 0.20 kcal/mol (0.009 eV) was found. Difficulties in direct comparisons between theory and experiment for vertical ionization potentials are discussed. With regard to the differences found for the ²A₁/²Π_u and ²B₂ adiabatic ionization potentials, a reinterpretation of the experimental spectrum appears justified.