The multichannel i-propyl + O2 reaction system: A model of secondary alkyl radical oxidation.

The i-propyl + O2 reaction mechanism has been investigated by definitive quantum chemical methods to establish this system as a benchmark for the combustion of secondary alkyl radicals. Focal point analyses extrapolating to the ab initio limit were performed based on explicit computations with electron correlation treatments through coupled cluster single, double, triple, and quadruple excitations and basis sets up to cc-pV5Z. The rigorous coupled cluster single, double, and triple excitations/cc-pVTZ level of theory was used to fully optimize all reaction species and transition states, thus, removing some substantial flaws in reference geometries existing in the literature. The vital i-propylperoxy radical (MIN1) and its concerted elimination transition state (TS1) were found 34.8 and 4.4 kcal mol-1 below the reactants, respectively. Two ß-hydrogen transfer transition states (TS2, TS2') lie above the reactants by (1.4, 2.5) kcal mol-1 and display large Born-Oppenheimer diagonal corrections indicative of nearby surface crossings. An α-hydrogen transfer transition state (TS5) is discovered 5.7 kcal mol-1 above the reactants that bifurcates into equivalent α-peroxy radical hanging wells (MIN3) prior to a highly exothermic dissociation into acetone + OH. The reverse TS5 → MIN1 intrinsic reaction path also displays fascinating features, including another bifurcation and a conical intersection of potential energy surfaces. An exhaustive conformational search of two hydroperoxypropyl (QOOH) intermediates (MIN2 and MIN3) of the i-propyl + O2 system located nine rotamers within 0.9 kcal mol-1 of the corresponding lowest-energy minima.

Identification and Reactivity of s-cis,s-cis-Dihydroxycarbene, a New [CH2O2] Intermediate.

We report the first preparation of the s-cis,s-cis conformer of dihydroxycarbene (1cc) by means of pyrolysis of oxalic acid, isolation of the lower-energy s-trans,s-trans (1tt) and s-cis,s-trans (1ct) product conformers at cryogenic temperatures in a N2 matrix, and subsequent narrow-band near-infrared (NIR) laser excitation to give 1cc. Carbene 1cc converts quickly to 1ct via quantum-mechanical tunneling with an effective half-life of 22 min at 3 K. The potential energy surface features around 1 were pinpointed by convergent focal point analysis targeting the AE-CCSDT(Q)/CBS level of electronic structure theory. Computations of the tunneling kinetics confirm the time scale of the 1cc → 1ct rotamerization and suggest that direct 1cc → H2 + CO2 decomposition may also be a minor pathway. The intriguing latter possibility cannot be confirmed spectroscopically, but hints of it may be present in the measured kinetic profiles.

Riddles of the structure and vibrational dynamics of HO3 resolved near the ab initio limit.

The hydridotrioxygen (HO3) radical has been investigated in many previous theoretical and experimental studies over several decades, originally because of its possible relevance to the tropospheric HOx cycle but more recently because of its fascinating chemical bonding, geometric structure, and vibrational dynamics. We have executed new, comprehensive research on this vexing molecule via focal point analyses (FPA) to approach the ab initio limit of optimized geometric structures, relative energies, complete quartic force fields, and the entire reaction path for cis-trans isomerization. High-order coupled cluster theory was applied through the CCSDT(Q) and even CCSDTQ(P) levels, and CBS extrapolations were performed using cc-pVXZ (X = 2-6) basis sets. The cis isomer proves to be higher than trans by 0.52 kcal mol-1, but this energetic ordering is achieved only after the CCSDT(Q) milestone is reached; the barrier for cis → trans isomerization is a minute 0.27 kcal mol-1. The FPA central re(O-O) bond length of trans-HO3 is astonishingly long (1.670 Å), consistent with the semiexperimental re distance we extracted from microwave rotational constants of 10 isotopologues using FPA vibration-rotation interaction constants (αi). The D0(HO-O2) dissociation energy converges to a mere 2.80 ± 0.25 kcal mol-1. Contrary to expectation for such a weakly bound system, vibrational perturbation theory performs remarkably well with the FPA anharmonic force fields, even for the torsional fundamental near 130 cm-1. Exact numerical procedures are applied to the potential energy function for the torsional reaction path to obtain energy levels, tunneling rates, and radiative lifetimes. The cis → trans isomerization occurs via tunneling with an inherent half-life of 1.4 × 10-11 s and 8.6 × 10-10 s for HO3 and DO3, respectively, thus resolving the mystery of why the cis species has not been observed in previous experiments executed in dissipative environments that allow collisional cooling of the trans-HO3 product. In contrast, the pure ground eigenstate of the cis species in a vacuum is predicted to have a spontaneous radiative lifetime of about 1 h and 5 days for HO3 and DO3, respectively.

Intricate Conformational Tunneling in Carbonic Acid Monomethyl Ester.

Disentangling internal and external effects is a key requirement for understanding conformational tunneling processes. Here we report the s- trans/ s- cis tunneling rotamerization of carbonic acid monomethyl ester (1) under matrix isolation conditions and make comparisons to its parent carbonic acid (3). The observed tunneling rate of 1 is temperature-independent in the 3-20 K range and accelerates when using argon instead of neon as the matrix material. The methyl group increases the effective half life (τeff) of the energetically disfavored s- trans-conformer from 3-5 h for 3 to 11-13 h for 1. Methyl group deuteration slows the rotamerization further (τeff ≈ 35 h). CCSD(T)/cc-pVQZ//MP2/aug-cc-pVTZ computations of the tunneling probability suggest that the rate should be almost unaffected by methyl substitution or its deuteration. Thus the observed relative rates are puzzling, and they disagree with previous explanations involving fast vibrational relaxation after the tunneling event facilitated by the alkyl rotor.

Ring-Walking of Zerovalent Nickel on Aryl Halides.

While ring-walking is a critical step in transition metal catalyzed cross-coupling reactions, the associated metastable intermediates are often difficult to isolate and characterize. In this work, theoretical structures and energetics for ring-walking and oxidative addition of zerovalent nickel with 1-bromo-2-methylbenzene, 2-bromopyridine, 2-bromo-3-methyl-thiophene, and 2-bromopyrrole were computed at the B3LYP-D3/TZ2P-LANL2TZ(f)-LANL08d level. The mechanisms vary qualitatively with substrate ring size and type-the catalyst weaves along the edges of the benzene and pyridine rings, cuts through the interior of the thiophene ring, and arcs along the bond opposite the nitrogen atom in the pyrrole ring. Analogous computations on the ring-walking and oxidative addition of zerovalent palladium with 1-bromo-2-methylbenzene reveal an energetic profile similar to that of Ni but with much weaker overall binding to the arene. In all cases, dispersion corrections are found to be very important for computing accurate metal-substrate binding energies.