Methanediol from cloud-processed formaldehyde is only a minor source of atmospheric formic acid.

Atmospheric formic acid is severely underpredicted by models. A recent study proposed that this discrepancy can be resolved by abundant formic acid production from the reaction (1) between hydroxyl radical and methanediol derived from in-cloud formaldehyde processing and provided a chamber-experiment-derived rate constant, k1 = 7.5 × 10-12 cm3 s-1. High-level accuracy coupled cluster calculations in combination with E,J-resolved two-dimensional master equation analyses yield k1 = (2.4 ± 0.5) × 10-12 cm3 s-1 for relevant atmospheric conditions (T = 260-310 K and P = 0-1 atm). We attribute this significant discrepancy to HCOOH formation from other molecules in the chamber experiments. More importantly, we show that reversible aqueous processes result indirectly in the equilibration on a 10 min. time scale of the gas-phase reaction [Formula: see text] (2) with a HOCH2OH to HCHO ratio of only ca. 2%. Although HOCH2OH outgassing upon cloud evaporation typically increases this ratio by a factor of 1.5-5, as determined by numerical simulations, its in-cloud reprocessing is shown using a global model to strongly limit the gas-phase sink and the resulting production of formic acid. Based on the combined findings in this work, we derive a range of 1.2-8.5 Tg/y for the global HCOOH production from cloud-derived HOCH2OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant Keq (2) determined in this work also allows us to estimate the Henry's law constant of methanediol (8.1 × 105 M atm-1 at 280 K).

Direct Probes of π-Delocalization in Prototypical Resonance-Stabilized Radicals: Hyperfine-Resolved Microwave Spectroscopy of Isotopic Propargyl and Cyanomethyl.

Delocalization of the unpaired electron in π-conjugated radicals has profound implications for their chemistry, but direct and quantitative characterization of this electronic structure in isolated molecules remains challenging. We apply hyperfine-resolved microwave rotational spectroscopy to rigorously probe π-delocalization in propargyl, CH2CCH, a prototypical resonance-stabilized radical and key reactive intermediate. Using the spectroscopic constants derived from the high-resolution cavity Fourier transform microwave measurements of an exhaustive set of 13C- and 2H-substituted isotopologues, together with high-level ab initio calculations of zero-point vibrational effects, we derive its precise semiexperimental equilibrium geometry and quantitatively characterize the spatial distribution of its unpaired electron. Our results highlight the importance of considering both spin-polarization and orbital-following contributions when interpreting the isotropic hyperfine coupling constants of π radicals. These physical insights are strengthened by a parallel analysis of the isoelectronic species cyanomethyl, CH2CN, using new 13C measurements also reported in this work. A detailed comparison of the structure and electronic properties of propargyl, cyanomethyl, and other closely related species allows us to correlate trends in their chemical bonding and electronic structure with critical changes in their reactivity and thermochemistry.

On the performance of composite schemes in determining equilibrium molecular structures.

Determination of equilibrium molecular structures is an essential ingredient in predicting spectroscopic parameters that help in identifying molecular carriers of microwave transitions. Here, the performance of two different ab initio composite approaches for obtaining equilibrium structures, "energy scheme" and "geometry scheme," is explored and compared to semi-experimental equilibrium structures. This study is performed for a set of 11 molecules which includes diatomics, linear triatomics, and a few non-linear molecules. The ab initio calculations were performed using three tiers of composite chemical recipes. The current results show that as the overall rigor of calculation is increased, the semi-experimental and the ab initio numbers agree to within 0.0003 Å for all molecules in the test set. The composite approach based on correcting the potential energy surface (energy scheme) and the one based on correcting the geometry directly (geometry scheme) show excellent agreement with each other. This work represents a step toward development of efficient and highly accurate procedures for computing ab initio equilibrium structures.

Influence of fourth-order vibrational corrections on semi-experimental (reSE) structures of linear molecules.

Semi-experimental structures (reSE) are derived from experimental ground state rotational constants combined with theoretical vibrational corrections. They permit a meaningful comparison with equilibrium structures based on high-level ab initio calculations. Typically, the vibrational corrections are evaluated with second-order vibrational perturbation theory (VPT2). The amount of error introduced by this approximation is generally thought to be small; however, it has not been thoroughly quantified. Herein, we assess the accuracy of theoretical vibrational corrections by extending the treatment to fourth order (VPT4) for a series of small linear molecules. Typical corrections to bond distances are on the order of 10-5 Å. Larger corrections, nearly 0.0002 Å, are obtained to the bond lengths of NCCN and CNCN. A borderline case is CCCO, which will likely require variational computations for a satisfactory answer. Treatment of vibrational effects beyond VPT2 will thus be important when one wishes to know bond distances confidently to four decimal places (10-4 Å). Certain molecules with shallow bending potentials, e.g., HOC+, are not amenable to a VPT2 description and are not improved by VPT4.

Photodissociation of dicarbon: How nature breaks an unusual multiple bond.

The dicarbon molecule (C2) is found in flames, comets, stars, and the diffuse interstellar medium. In comets, it is responsible for the green color of the coma, but it is not found in the tail. It has long been held to photodissociate in sunlight with a lifetime precluding observation in the tail, but the mechanism was not known. Here we directly observe photodissociation of C2 From the speed of the recoiling carbon atoms, a bond dissociation energy of 602.804(29) kJ·mol[Formula: see text] is determined, with an uncertainty comparable to its more experimentally accessible N2 and O2 counterparts. The value is within 0.03 kJ·mol-1 of high-level quantum theory. This work shows that, to break the quadruple bond of C2 using sunlight, the molecule must absorb two photons and undergo two "forbidden" transitions.

Precise Equilibrium Structure of Benzene.

Recent advances in gas-phase structure determination afford outstanding agreement between the CCSD(T)/cc-pCVTZ-corrected semi-experimental (reSE) equilibrium structures and their corresponding best theoretical estimates (BTEs) of the equilibrium structures (re) based upon corrections to the CCSD(T)/cc-pCV5Z geometries for the aromatic heterocycles pyrimidine and pyridazine. Herein, that same analysis is extended to the fundamental aromatic molecule benzene, using published experimental spectroscopic data for a total of 11 available isotopologues. The incorporation of rotational constants from all of these isotopologues and CCSD(T) corrections to address the impacts of both the vibration-rotation interaction and electron-mass distribution results in a highly precise and accurate reSE structure. The CCSD(T)/cc-pCV5Z optimized geometry has been further corrected to address a finite basis set, untreated electron correlation, relativistic effects, and a breakdown of the Born-Oppenheimer approximation. This analysis achieves outstanding agreement between the re (BTE) and reSE structural parameters of benzene to a highly satisfying level (0.0001 Å), an agreement that surpasses our recently published structures of the aforementioned nitrogen-substituted benzene analogues. The D6h geometry of benzene is now known to an unprecedented precision: RC-C = 1.3913 (1) Å and RC-H = 1.0809 (1) Å. The mutual agreement between theory and experiment presented in this work validates both, substantially resolving all discrepancies between the reSE and theoretical re structures available in the literature.

Connections between the accuracy of rotational constants and equilibrium molecular structures.

Rotational spectroscopy is the technique of choice for investigating molecular structures in the gas phase. Indeed, rotational constants are strongly connected to the geometry of the molecular system under consideration. Therefore, they are powerful tools for assessing the accuracy that quantum chemical approaches can reach in structural determinations. In this review article, it is shown how it is possible to measure the accuracy of a computed equilibrium geometry based on the comparison of rotational constants. But, it is also addressed what accuracy is required by computations for providing molecular structures and thus rotational constants that are useful to experiment. Quantum chemical methodologies for obtaining the "0.1% accuracy" for rotational constants are reviewed for systems ranging in size from small molecules to small polycyclic aromatic hydrocarbons. This accuracy for systems containing two dozen or so atoms opens the way towards future applications such as the accurate characterization of non-covalent interactions, which play a key role in several biological and technological processes.

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).

Rotamers of Methanediol: Composite Ab Initio Predictions of Structures, Frequencies, and Rovibrational Constants.

Geminal diols are known to be important intermediates in atmospheric ozonolysis and the aerosol cycle. Recently, the simplest member of this class, methanediol, was interrogated in the gas phase with infrared spectroscopy. To aid in future spectroscopic investigations of methanediol, including in the interstellar medium, we report fundamental frequencies and rovibrational constants for the two rotamers of this molecule using ab initio composite methods along with vibrational perturbation theory. Sensitivity of the predictions to the level of theory and the treatment of anharmonic resonances are carefully assessed. The OH stretching harmonic frequencies of both rotamers are particularly sensitive to the level of theory. The CH stretches of the Cs rotamer are sensitive to the treatment of anharmonic resonances with VPT2-based effective Hamiltonian models. Equilibrium bond distances and harmonic frequencies are converged conservatively to within 0.0005 Å and 3 cm-1, respectively. The effect of tunneling on the rotational constants is investigated with a 2D variational calculation, based on a relaxed hydroxyl torsional potential energy surface. Tunneling is found to be negligible in the lower energy C2 rotamer but should modify the rotational constants of the Cs rotamer on the order of MHz, giving rise to rotational line splittings of the same order. The rovibrational constants of the Cs rotamer are dominated by hydroxyl torsional effects, and here we see evidence for the breakdown of vibrational perturbation theory.

Relation Between Bond Angle and Carbon-Oxygen Stretching Frequencies in CO2-Containing Compounds.

The symmetric (νs) and antisymmetric (νas) O-C-O stretching modes of CO2-containing compounds encode structural information that can be difficult to decipher, due to the sensitivity of these spectral features to small shifts in charge distribution and structure, as well as the anharmonicities of these two vibrational modes. In this work, we discuss the relation between the frequency of these modes and the geometry of the O-C-O group, showing that the splitting between νs and νas (Δνas-s = νas - νs) can be predicted based only on the O-C-O bond angle obtained from quantum chemical calculations with reasonable accuracy (±46 cm-1, R2 = 0.994). The relationship is shown to hold for the infrared spectra of a variety of CO2-containing molecules measured in vacuo. The origins of this model are discussed in the framework of elementary mode analysis.

Calculated and Empirical Values of Vibronic Transition Dipole Moments of Reactive Chemical Intermediates for Determination of Concentrations.

Absorption spectroscopy has long been known as a technique for making molecular concentration measurements and has received enhanced visibility in recent years with the advent of new techniques, like cavity ring-down spectroscopy, that have increased its sensitivity. To apply the method, it is necessary to have a known molecular absorption cross section for the species of interest, which typically is obtained by measurements of a standard sample of known concentration. However, this method fails if the species is highly reactive, and indirect means for attaining the cross section must be employed. The HO2 and alkyl peroxy radicals are examples of reactive species for which absorption cross sections have been reported. This work explores and describes for these peroxy radicals the details of an alternative approach for obtaining these cross sections using quantum chemistry methods for the calculation of the transition dipole moment upon whose square the cross section depends. Likewise, details are given for obtaining the transition moment from the experimentally measured cross sections of individual rovibronic lines in the near-IR Ã-X̃ electronic spectrum of HO2 and the peaks of the rotational contours in the corresponding electronic transitions for the alkyl (methyl, ethyl, and acetyl) peroxy radicals. In the case of the alkyl peroxy radicals, good agreement for the transition moments, ≈20%, is found between the two methods. However, rather surprisingly, the agreement is significantly poorer, ≈40%, for the HO2 radical. Possible reasons for this disagreement are discussed.

Vibrationally excited states of 1H- and 2H-1,2,3-triazole isotopologues analyzed by millimeter-wave and high-resolution infrared spectroscopy with approximate state-specific quartic distortion constants.

In this work, we present the spectral analysis of 1H- and 2H-1,2,3-triazole vibrationally excited states alongside provisional and practical computational predictions of the excited-state quartic centrifugal distortion constants. The low-energy fundamental vibrational states of 1H-1,2,3-triazole and five of its deuteriated isotopologues ([1-2H]-, [4-2H]-, [5-2H]-, [4,5-2H]-, and [1,4,5-2H]-1H-1,2,3-triazole), as well as those of 2H-1,2,3-triazole and five of its deuteriated isotopologues ([2-2H]-, [4-2H]-, [2,4-2H]-, [4,5-2H]-, and [2,4,5-2H]-2H-1,2,3-triazole), are studied using millimeter-wave spectroscopy in the 130-375 GHz frequency region. The normal and [2-2H]-isotopologues of 2H-1,2,3-triazole are also analyzed using high-resolution infrared spectroscopy, determining the precise energies of three of their low-energy fundamental states. The resulting spectroscopic constants for each of the vibrationally excited states are reported for the first time. Coupled-cluster vibration-rotation interaction constants are compared with each of their experimentally determined values, often showing agreement within 500 kHz. Newly available coupled-cluster predictions of the excited-state quartic centrifugal distortion constants based on fourth-order vibrational perturbation theory are benchmarked using a large number of the 1,2,3-triazole tautomer isotopologues and vibrationally excited states studied.

Improved semi-experimental equilibrium structure and high-level theoretical structures of ketene.

The millimeter-wave rotational spectrum of ketene (H2C=C=O) has been collected and analyzed from 130 to 750 GHz, providing highly precise spectroscopic constants from a sextic, S-reduced Hamiltonian in the Ir representation. The chemical synthesis of deuteriated samples allowed spectroscopic measurements of five previously unstudied ketene isotopologues. Combined with previous work, these data provide a new, highly precise, and accurate semi-experimental (reSE) structure for ketene from 32 independent moments of inertia. This reSE structure was determined with the experimental rotational constants of each available isotopologue, together with computed vibration-rotation interaction and electron-mass distribution corrections from coupled-cluster calculations with single, double, and perturbative triple excitations [CCSD(T)/cc-pCVTZ]. The 2σ uncertainties of the reSE parameters are ≤0.0007 Å and 0.014° for the bond distances and angle, respectively. Only S-reduced spectroscopic constants were used in the structure determination due to a breakdown in the A-reduction of the Hamiltonian for the highly prolate ketene species. All four reSE structural parameters agree with the "best theoretical estimate" (BTE) values, which are derived from a high-level computed re structure [CCSD(T)/cc-pCV6Z] with corrections for the use of a finite basis set, the incomplete treatment of electron correlation, relativistic effects, and the diagonal Born-Oppenheimer breakdown. In each case, the computed value of the geometric parameter lies within the statistical experimental uncertainty (2σ) of the corresponding semi-experimental coordinate. The discrepancies between the BTE structure and the reSE structure are 0.0003, 0.0000, and 0.0004 Å for rC-C, rC-H, and rC-O, respectively, and 0.009° for θC-C-H.

Ring-Opening Dynamics of the Cyclopropyl Radical and Cation: the Transition State Nature of the Cyclopropyl Cation.

We provide compelling experimental and theoretical evidence for the transition state nature of the cyclopropyl cation. Synchrotron photoionization spectroscopy employing coincidence techniques together with a novel simulation based on high-accuracy ab initio calculations reveal that the cation is unstable via its allowed disrotatory ring-opening path. The ring strains of the cation and the radical are similar, but both ring opening paths for the radical are forbidden when the full electronic symmetries are considered. These findings are discussed in light of the early predictions by Longuet-Higgins alongside Woodward and Hoffman; we also propose a simple phase space explanation for the appearance of the cyclopropyl photoionization spectrum. The results of this work allow the refinement of the cyclopropane C-H bond dissociation energy, in addition to the cyclopropyl radical and cation cyclization energies, via the Active Thermochemical Tables approach.

Mechanism, thermochemistry, and kinetics of the reversible reactions: C2H3 + H2 ⇌ C2H4 + H ⇌ C2H5.

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations, was used in a study of the title reactions, which play an important role in the combustion of hydrocarbons. In the set of radical/radical reactions leading to soot formation in flames, the addition of H-atoms to alkenes is likely a common reaction, triggering the isomerization of complex hydrocarbons to aromatics. The heats of formation of C2H3, C2H4, and C2H5 are established to be 301.26 ± 0.30 at 0 K (297.22 ± 0.30 at 298 K), 60.89 ± 0.11 (52.38 ± 0.11), and 131.38 ± 0.22 (120.63 ± 0.22) kJ mol-1, respectively. The calculated rate constants from first principles agree well with experiments where they are available. Under conditions typical of high temperature combustion - where experimental work is very challenging with a consequent dearth of accurate data - we provide high-level theoretical results for kinetic modeling.

Theoretical Rovibrational Spectroscopy of Magnesium Tricarbide-Multireference Character Thwarts a Full Analysis of All Isomers.

Magnesium tricarbide isomers are studied herein with coupled cluster theory and multireference configuration interaction to support their possible detection in astrochemical environments such as the circumstellar envelope surrounding the star IRC +10216 or in terrestrial laboratories. Magnesium-bearing species may abound in the interstellar medium (ISM), but only eight (MgNC, MgCN, HMgNC, MgC2H, MgC3N, MgC4H, MgC5N, and MgC6H) have been directly identified thus far. Several possible isomers for the related MgC3 system are explored in their singlet and triplet spin multiplicities. Overall, this work offers quantum chemical insight of rovibrational spectroscopic data for MgC3 using quartic force fields (QFFs) based on the CCSD(T) and CCSD(T)-F12 levels of theory at the complete basis set (CBS) limit. Additional corrections with small basis set CCSDT(Q) and scalar relativistic effects are also included in the analysis. Salient multireference character is found in the singlet diamond electronic state, which makes a definitive assignment of the ground state challenging. Nevertheless, coupled cluster-based composite energies and multireference configuration interaction both predict that the 1A1 diamond isomer is 1.6-2.2 kcal mol-1 lower in energy than the 3A1 diamond isomer. Furthermore, highly accurate binding energies of various isomers MgC3 are provided for comparison to photodetachment experiments. Dipole moments along with harmonic infrared intensities will guide efforts for astronomical and spectroscopic characterization.

PAH Growth in Flames and Space: Formation of the Phenalenyl Radical.

Polycyclic aromatic hydrocarbons (PAHs) are intermediates in the formation of soot particles and interstellar grains. However, their formation mechanisms in combustion and interstellar environments are not fully understood. The production of tricyclic PAHs and, in particular, the conversion of a PAH containing a five-membered ring to one with a six-membered ring are of interest to explain PAH abundances in combustion processes. In the present work, resonant ionization mass spectrometry in conjunction with isotopic labeling is used to investigate the formation of the phenalenyl radical from acenaphthylene and methane in an electrical discharge. We show that in this environment the CH cycloaddition mechanism converts a five-membered ring to a six-membered ring. This mechanism can occur in tandem with other PAH formation mechanisms such as hydrogen abstraction/acetylene addition (HACA) to produce larger PAHs in flames and the interstellar medium.

Pure Rotational Spectroscopy of the CH2CN Radical Extended to the Sub-Millimeter Wave Spectral Region.

We present a thorough pure rotational investigation of the CH2CN radical in its ground vibrational state. Our measurements cover the millimeter and sub-millimeter wave spectral regions (79-860 GHz) using a W-band chirped-pulse instrument and a frequency multiplication chain-based spectrometer. The radical was produced in a flow cell at room temperature by H abstraction from acetonitrile using atomic fluorine. The newly recorded transitions of CH2CN (involving N″ and Ka″ up to 42 and 8, respectively) were combined with the literature data, leading to a refinement of the spectroscopic parameters of the species using a Watson S-reduced Hamiltonian. In particular, the A rotational constant and K-dependent parameters are significantly better determined than in previous studies. The present model, which reproduces all experimental transitions to their experimental accuracy, allows for confident searches for the radical in cold to warm environments of the interstellar medium.

Millimeter/Submillimeter-wave Spectroscopy and the Semi-experimental Equilibrium (reSE) Structure of 1H-1,2,4-Triazole (c-C2H3N3).

The millimeter/submillimeter spectrum of 1H-1,2,4-triazole is reported from 70 to 700 GHz, providing spectral frequencies directly comparable to radio telescopes and enabling an astronomical search. Using four deuteriated samples of 1,2,4-triazole, we measured, assigned, and least-squares fit transitions for 26 isotopologues to sextic A- and S-reduced Hamiltonians. An accurate and precise semi-experimental (reSE) structure from 50 independent moments of inertia has been obtained. Structural parameters are provided with 2σ uncertainties within 0.0009 Å for bond distances and 0.09° for bond angles. The structural parameters are in quite good agreement with the best theoretical estimate (BTE) obtained using CCSD(T)/cc-pCV5Z, where an agreement within the 2σ uncertainty is observed for all but one case. Despite the large number of isotopologues already included in this structure, more may be useful. One isotopologue, [1,3-2H]-1H-1,2,4-triazole, is observed to closely approach the oblate asymmetric-top limit, resulting in a clear breakdown of the A-reduction Hamiltonian. The highly accurate reSE structure and subsequent analysis demonstrates that the S-reduction is also unable to adequately model the spectrum of this isotopologue.

Atomic isotropic hyperfine properties for first row elements (B-F) revisited.

Benchmark quality isotropic hyperfine properties have been obtained for first row elements (B-F) using a systematic composite approach consisting of a sequence of core/valence correlation consistent basis sets, up through aug-cc-pCV8Z, along with configuration interaction and coupled cluster theory methods. The best nonrelativistic final values (in MHz) are 10.64 (B), 20.22 (C), 10.59 (N), -31.74 (O), and 318.30 (F) and are in very good agreement with available experimental values for these difficult-to-describe properties. Agreement is especially close in the case of N, which has the most accurate experimental value. The spin densities derived from the best composite level of theory were found to closely follow a simple quadratic scaling with the atomic number, Z. Observed convergence rates in the 1-particle and n-particle expansions obtained here may be useful in judging likely accuracy that can be expected in studies of molecular systems.