Computational Investigation of RO2 + HO2 and RO2 + RO2 Reactions of Monoterpene Derived First-Generation Peroxy Radicals Leading to Radical Recycling.

The oxidation of biogenically emitted volatile organic compounds (BVOC) plays an important role in the formation of secondary organic aerosols (SOA) in the atmosphere. Peroxy radicals (RO2) are central intermediates in the BVOC oxidation process. Under clean (low-NO x) conditions, the main bimolecular sink reactions for RO2 are with the hydroperoxy radical (HO2) and with other RO2 radicals. Especially for small RO2, the RO2 + HO2 reaction mainly leads to closed-shell hydroperoxide products. However, there exist other known RO2 + HO2 and RO2 + RO2 reaction channels that can recycle radicals and oxidants in the atmosphere, potentially leading to lower-volatility products and enhancing SOA formation. In this work, we present a thermodynamic overview of two such reactions: (a) RO2 + HO2 → RO + OH + O2 and (b) R'O2 + RO2 → R'O + RO + O2 for selected monoterpene + oxidant derived peroxy radicals. The monoterpenes considered are α-pinene, ß-pinene, limonene, trans-ß-ocimene, and Δ3-carene. The oxidants considered are the hydroxyl radical (OH), the nitrate radical (NO3), and ozone (O3). The reaction Gibbs energies were calculated at the DLPNO-CCSD(T)/def2-QZVPP//ωB97X-D/aug-cc-pVTZ level of theory. All reactions studied here were found to be exergonic in terms of Gibbs energy. On the basis of a comparison with previous mechanistic studies, we predict that reaction a and reaction b are likely to be most important for first-generation peroxy radicals from O3 oxidation (especially for ß-pinene), while being less so for most first-generation peroxy radicals from OH and NO3 oxidation. This is because both reactions are comparatively more exergonic for the O3 oxidized systems than their OH and NO3 oxidized counterparts. Our results indicate that bimolecular reactions of certain complex RO2 may contribute to an increase in radical and oxidant recycling under high HO2 conditions in the atmosphere, which can potentially enhance SOA formation.

Benchmarking ab initio binding energies of hydrogen-bonded molecular clusters based on FTIR spectroscopy.

Models of formation and growth of atmospheric aerosols are highly dependent on accurate cluster binding energies. These are most often calculated by ab initio electronic structure methods but remain associated with significant uncertainties. We present a computational benchmarking study of the Gibbs free binding energies in molecular complexes and clusters based on gas phase FTIR spectroscopy. The acetonitrile-HCl molecular complex is identified via its redshifted H-Cl stretching vibrational mode. We determine the Gibbs free binding energy, ΔG°295 K, to between 4.8 and 7.9 kJ mol(-1) and compare this range to predictions from several widely used electronic structure methods, including five density functionals, Møller-Plesset perturbation theory, and five coupled cluster methods up to CCSDT quality, considering also the D3 dispersion correctional scheme. With some exceptions, we find that most electronic structure methods overestimate ΔG°295 K. The effects of vibrational anharmonicity is approximated using scaling factors, reducing ΔG°295 K by ca. 1.8 kJ mol(-1), whereby ΔG°295 K predictions well within the experimental range can be obtained.