Structural and Thermochemical Properties of Hydroxymethyl (CH2OH) Radicals and Cations Derived from Observations of B̃2A‘(3p) ← X̃2A‘‘ Electronic Spectra and fromab InitioCalculations

Stereodynamical control of product branching in multi-channel barrierless hydrogen abstraction of CH3OH by F

Comprehensive dynamical simulations of a prototypical multi-channel reaction on a globally accurate potential energy surface show that the non-statistical product branching is dictated by unique stereodynamics in the entrance channels.

Hydrogen abstraction from methanol (CH₃OH) by F atoms presents an ideal proving ground to investigate dynamics of multi-channel reactions, because two types of hydrogen can be abstracted from the methanol molecule leading to the HF + CH₃O and HF + CH₂OH products. Using the quasi-classical trajectory approach on a globally accurate potential energy surface based on high-level ab initio calculations, this work reports a comprehensive dynamical investigation of this multi-channel reaction, yielding measurable attributes including integral and differential cross sections, as well as branching ratios. It is shown that while complex-forming and direct mechanisms coexist at low collision energies, these barrierless reaction channels are dominated at high energies by the direct mechanism, in which the reaction is only possible for trajectories entering into the respective dynamical cones of acceptance. Perhaps more importantly, the non-statistical product branching is found to be dictated by unique stereodynamics in the entrance channels.

Active vs. spectator modes in nonadiabatic photodissociation dynamics of the hydroxymethyl radical via the 22A(3s) Rydberg state

The choice of the active degrees of freedom (DOFs) is a pivotal issue in a reduced-dimensional model of quantum dynamics when a full-dimensional one is not feasible. Here, several five-dimensional (5D) models are used to investigate the nonadiabatic photodissociation dynamics of the hydroxymethyl (CH₂OH) radical, which possesses nine internal DOFs, in its lowest absorption band. A normal-mode based scheme is used to identify the active and spectator modes, and its predictions are confirmed by 5D quantum dynamical calculations. Our results underscore the important role of the CO stretching mode in the photodissociation dynamics of CH₂OH, originating from the photo-induced promotion of an electron from the half-occupied π*_CO antibonding orbital to a carbon Rydberg orbital.

Sub-Doppler infrared spectroscopy of CH2OH radical in a slit supersonic jet: Vibration-rotation-tunneling dynamics in the symmetric CH stretch manifold

The sub-Doppler CH-symmetric stretch (ν₃) infrared absorption spectrum of a hydroxymethyl (CH₂OH) radical is observed and analyzed with the radical formed in a slit-jet supersonic discharge expansion (T_rot = 18 K) via Cl atom mediated H atom abstraction from methanol. The high sensitivity of the spectrometer and reduced spectral congestion associated with the cooled expansion enable first infrared spectroscopic observation of hydroxymethyl transitions from both ± symmetry tunneling states resulting from large amplitude COH torsional motion. Nuclear spin statistics due to exchange of the two methyl H-atoms aid in unambiguous rovibrational assignment of two A-type K_a = 0 ← 0 and K_a = 1 ← 1 bands out of each ± tunneling state, with additional spectral information obtained from spin-rotation splittings in P, Q, and R branch K_a = 1 ← 1 transitions that become resolved at low N. A high level ab initio potential surface (CCSD(T)-f12b/cc-pvnzf12 (n = 2,3)/CBS) is calculated in the large amplitude COH torsional and CH₂ wag coordinates, which in the adiabatic approximation and with zero point correction predicts ground state tunneling splittings in good qualitative agreement with experiment. Of particular astrochemical interest, a combined fit of the present infrared ground state combination differences with recently reported millimeter-wave frequencies permits the determination of improved accuracy rotational constants for the ground vibrational state, which will facilitate ongoing millimeter/microwave searches for a hydroxymethyl radical in the interstellar medium.

Active Thermochemical Tables: Sequential Bond Dissociation Enthalpies of Methane, Ethane, and Methanol and the Related Thermochemistry

Active Thermochemical Tables (ATcT) thermochemistry for the sequential bond dissociations of methane, ethane, and methanol systems were obtained by analyzing and solving a very large thermochemical network (TN). Values for all possible C-H, C-C, C-O, and O-H bond dissociation enthalpies at 298.15 K (BDE298) and bond dissociation energies at 0 K (D0) are presented. The corresponding ATcT standard gas-phase enthalpies of formation of the resulting CHn, n = 4-0 species (methane, methyl, methylene, methylidyne, and carbon atom), C2Hn, n = 6-0 species (ethane, ethyl, ethylene, ethylidene, vinyl, ethylidyne, acetylene, vinylidene, ethynyl, and ethynylene), and COHn, n = 4-0 species (methanol, hydroxymethyl, methoxy, formaldehyde, hydroxymethylene, formyl, isoformyl, and carbon monoxide) are also presented. The ATcT thermochemistry of carbon dioxide, water, hydroxyl, and carbon, oxygen, and hydrogen atoms is also included, together with the sequential BDEs of CO2 and H2O. The provenances of the ATcT enthalpies of formation, which are quite distributed and involve a large number of relevant determinations, are analyzed by variance decomposition and discussed in terms of principal contributions. The underlying reasons for periodic appearances of remarkably low and/or unusually high BDEs, alternating along the dissociation sequences, are analyzed and quantitatively rationalized. The present ATcT results are the most accurate thermochemical values currently available for these species.

Computational Mechanistic Study of the Gas Phase Oxidation of Methanol with Ozone

The reaction of methanol and ozone on the singlet potential energy surface has been carried out using the RMP2 theoretical approach in connection with the 6–311++G(d, p) basis set. To obtain a more reliable energy, single point calculations were performed for all species at the CCSD(T) level. One pre-reactive complex C1 is formed between the reactants. From a variety of the complexes, five types of product are obtained, of which four types have enough thermodynamic stability at the RMP2 level. In the thermodynamic approach, the most favoured route begins with the formation of a pre-reactive complex C1 and gives H2CO+H2O+3O2 as the final adducts. This process is exothermic by – 73.243 kcal mol–1 in standard enthalpy and exergonic by – 83.382 kcal mol–1.

Size effects on cation heats of formation. II. Methyl substitutions in oxygen compounds

A relation between the heat of formation of molecular ions and cation size is used to study the effects of methyl substitution on a series of species containing oxygen. These include methanol, the hydroxymethyl radical, formaldehyde and related isomers, acetaldehyde and related isomers, ketene, formic acid, and acetic acid. This size-dependent relation is found to be valid in most cases, enabling a choice to be made among conflicting reported ion and neutral heat of formation values, but it holds much less well in methyl-substituted formaldehyde and its carbene isomers. The cations formed in methanol, hydroxymethyl, hydroxycarbene, and formic acid by methyl substitution at carbon sites were found to be more stable than those substituted at oxygen sites. Suggestions are made for investigating or reinvestigating the ionization energies and the heats of formation of several of the molecules studied in cases where multiple choices are available in the literature or where our analysis suggests that more reliable values are needed.

Theoretical Prediction of the Heats of Formation of C2H5O• Radicals Derived from Ethanol and of the Kinetics of β-C−C Scission in the Ethoxy Radical

Thermochemical parameters of three C(2)H(5)O* radicals derived from ethanol were reevaluated using coupled-cluster theory CCSD(T) calculations, with the aug-cc-pVnZ (n = D, T, Q) basis sets, that allow the CC energies to be extrapolated at the CBS limit. Theoretical results obtained for methanol and two CH(3)O* radicals were found to agree within +/-0.5 kcal/mol with the experiment values. A set of consistent values was determined for ethanol and its radicals: (a) heats of formation (298 K) DeltaHf(C(2)H(5)OH) = -56.4 +/- 0.8 kcal/mol (exptl: -56.21 +/- 0.12 kcal/mol), DeltaHf(CH(3)C*HOH) = -13.1 +/- 0.8 kcal/mol, DeltaHf(C*H(2)CH(2)OH) = -6.2 +/- 0.8 kcal/mol, and DeltaHf(CH(3)CH(2)O*) = -2.7 +/- 0.8 kcal/mol; (b) bond dissociation energies (BDEs) of ethanol (0 K) BDE(CH(3)CHOH-H) = 93.9 +/- 0.8 kcal/mol, BDE(CH(2)CH(2)OH-H) = 100.6 +/- 0.8 kcal/mol, and BDE(CH(3)CH(2)O-H) = 104.5 +/- 0.8 kcal/mol. The present results support the experimental ionization energies and electron affinities of the radicals, and appearance energy of (CH(3)CHOH+) cation. Beta-C-C bond scission in the ethoxy radical, CH(3)CH2O*, leading to the formation of C*H3 and CH(2)=O, is characterized by a C-C bond energy of 9.6 kcal/mol at 0 K, a zero-point-corrected energy barrier of E0++ = 17.2 kcal/mol, an activation energy of Ea = 18.0 kcal/mol and a high-pressure thermal rate coefficient of k(infinity)(298 K) = 3.9 s(-1), including a tunneling correction. The latter value is in excellent agreement with the value of 5.2 s(-1) from the most recent experimental kinetic data. Using RRKM theory, we obtain a general rate expression of k(T,p) = 1.26 x 10(9)p(0.793) exp(-15.5/RT) s(-1) in the temperature range (T) from 198 to 1998 K and pressure range (p) from 0.1 to 8360.1 Torr with N2 as the collision partners, where k(298 K, 760 Torr) = 2.7 s(-1), without tunneling and k = 3.2 s(-1) with the tunneling correction. Evidence is provided that heavy atom tunneling can play a role in the rate constant for beta-C-C bond scission in alkoxy radicals.

Unimolecular processes in CH2OH below the dissociation barrier: O–H stretch overtone excitation and dissociation

The OH-stretch overtone spectroscopy and dynamics of the hydroxymethyl radical, CH(2)OH, are reported in the region of the second and third overtones, which is above the thermochemical threshold to dissociation to H+CH(2)O (D(0)=9600 cm(-1)). The second overtone spectrum at 10 484 cm(-1) is obtained by double resonance IR-UV resonance enhanced multiphoton ionization (REMPI) spectroscopy via the 3p(z) electronic state. It is rotationally resolved with a linewidth of 0.4 cm(-1) and displays properties of local-mode vibration. No dissociation products are observed. The third overtone spectra of CH(2)OH and CD(2)OH are observed at approximately 13 600 cm(-1) by monitoring H-atom photofragments while scanning the excitation laser frequency. No double resonance REMPI spectrum is detected, and no D fragments are produced. The spectra of both isotope analogs can be simulated with a linewidth of 1.3 cm(-1), indicating dissociation via tunneling. By treating the tunneling as one dimensional and using the calculated imaginary frequency, the barrier to dissociation is estimated at about 15 200 cm(-1), in good agreement with theoretical estimations. The Birge-Sponer plot is linear for OH-stretch vibrations 1nu(1)-4nu(1), demonstrating behavior of a one-dimensional Morse oscillator. The anharmonicity parameter derived from the plot is similar to the values obtained for other small OH containing molecules. Isomerization to methoxy does not contribute to the predissociation signal and the mechanism appears to be direct O-H fission via tunneling. CH(2)OH presents a unique example in which the reaction coordinate is excited directly and leads to predissociation via tunneling while preserving the local-mode character of the stretch vibration.

Statistical and nonstatistical nonadiabatic photodissociation from the first excited state of the hydroxymethyl radical

Photodissociation of hydroxymethyl (H(2)COD) from the first excited state, a Rydberg 3s state, can produce either H+HCOD or H(2)CO+D. These processes involve a seam of conical intersections. An analysis of the seam is reported based on a pointwise determination of the three gradient vectors that characterize a conical intersection: the energy difference gradient and the coupling gradient, which span the branching space, and the gradient of the average energy. These data are used to understand why H, but not D is produced at energies near threshold, the difference in the recoil anisotropy for H and D dissociation, and whether H or D will be produced impulsively or after formation of an intermediate complex.

Quantum Monte Carlo study of the reaction: Cl+CH3OH→CH2OH+HCl

A theoretical study is reported of the Cl+CH3OH-->CH2OH+HCl reaction based on the diffusion Monte Carlo (DMC) variant of the quantum Monte Carlo method. Using a DMC trial function constructed as a product of Hartree-Fock and correlation functions, we have computed the barrier height, heat of reaction, atomization energies, and heats of formation of reagents and products. The DMC heat of reaction, atomization energies, and heats of formation are found to agree with experiment to within the error bounds of computation and experiment. Møller-Plesset second order perturbation theory (MP2) and density functional theory, the latter in the B3LYP generalized gradient approximation, are found to overestimate the experimental heat of reaction. Intrinsic reaction coordinate calculations at the MP2 level of theory demonstrate that the reaction is predominantly direct, i.e., proceeds without formation of intermediates, which is consistent with a recent molecular beam experiment. The reaction barrier as determined from MP2 calculations is found to be 2.24 kcal/mol and by DMC it is computed to be 2.39(49) kcal/mol.

State-to-state dynamics of the Cl+CH3OH→HCl+CH2OH reaction

Molecular chlorine, methanol, and helium are co-expanded into a vacuum chamber using a custom designed "late-mixing" nozzle. The title reaction is initiated by photolysis of Cl2 at 355 nm, which generates monoenergetic Cl atoms that react with CH3OH at a collision energy of 1960 +/- 170 cm(-1) (0.24 +/- 0.02 eV). Rovibrational state distributions of the nascent HCl products are obtained via 2 + 1 resonance enhanced multiphoton ionization, center-of-mass scattering distributions are measured by the core-extraction technique, and the average internal energy of the CH3OH co-products is deduced by measuring the spatial anisotropy of the HCl products. The majority (84 +/- 7%) of the HCl reaction products are formed in HCl(v = 0) with an average rotational energy of [Erot] = 390 +/- 70 cm(-1). The remaining 16 +/- 7% are formed in HCl(v = 1) and have an average rotational energy of [Erot] = 190 +/- 30 cm(-1). The HCl(v = 1) products are primarily forward scattered, and they are formed in coincidence with CH2OH products that have little internal energy. In contrast, the HCl(v = 0) products are formed in coincidence with CH2OH products that have significant internal energy. These results indicate that two or more different mechanisms are responsible for the dynamics in the Cl + CH3OH reaction. We suggest that (1) the HCl(v = 1) products are formed primarily from collisions at high impact parameter via a stripping mechanism in which the CH2OH co-products act as spectators, and (2) the HCl(v = 0) products are formed from collisions over a wide range of impact parameters, resulting in both a stripping mechanism and a rebound mechanism in which the CH2OH co-products are active participants. In all cases, the reaction of fast Cl atoms with CH3OH is with the hydrogen atoms on the methyl group, not the hydrogen on the hydroxyl group.