Arbitrary-Order Derivatives of Quantum Chemical Methods via Automatic Differentiation
Not provided.
Chemistry↗
Engineering topics
Schaefer, Henry F.
Publications and source records attributed to Schaefer, Henry F..
Energetics and kinetics of various cyano radical hydrogen abstractions
The cyano radical (CN) is an abundant, open-shell molecule found in a variety of environments, including the atmosphere, the interstellar medium and combustion processes. In these environments, it often reacts with small, closed-shell molecules via hydrogen abstraction. Both carbon and nitrogen atoms of the cyano radical are reactive sites, however the carbon is more reactive with reaction barrier heights generally between 2–15 kcal mol -1 lower than those of the analogous nitrogen. The CN + HX → HCN/HNC + X, with X = H, CH 3 , NH 2 , OH, F, SiH 3 , PH 2 , SH, Cl, C 2 H, CN reactions have been studied at a high-level of theory, including CCSD(T)-F12a. Finally, kinetics were obtained over the 100–1000 K temperature range, showing excellent agreement with those rate constants that have been determined experimentally.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Isomer-Dependent Reaction Mechanisms of Cyclic Ether Intermediates: cis -23-Dimethyloxirane and trans -23-Dimethyloxirane.
Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β-hydroperoxyalkyl radicals (β-˙QOOH) or from reactions of HÒO with alkenes. The cis- and trans-isomers of 2,3-dimethyloxirane are intermediates of n-butane oxidation, and while rate coefficients for β-˙QOOH → 2,3-dimethyloxirane + OH are reported extensively, subsequent reaction mechanisms of the cyclic ethers are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe the consumption of 2,3-dimethyloxirane by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. The present research examines the isomerdependence of 2,3-dimethyloxirane reaction mechanisms in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred via the detection of products from Cl-initiated oxidation of both cis-2,3-dimethyloxirane and trans-2,3-dimethyloxirane using multiplexed photoionization mass spectrometry (MPIMS). The experiments were conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, the enthalpies of stationary points on the ˙R + O 2 surfaces were computed at the ccCA-PS3 level of theory. In total, 28 barrier heights were computed on the 2,3-dimethyloxiranylperoxy surfaces. Two notable aspects are low-lying pathways that form resonance-stabilized ketohydroperoxide-type radicals caused by ˙QOOH ring-opening when the unpaired electron is localized adjacent to the ether group, and cis-trans isomerization of ˙R and ˙QOOH radicals, via inversion, which enable reaction pathways otherwise restricted by stereochemistry. Several species were identified in the MPIMS experiments from ring opening of 2,3-dimethyloxiranyl radicals. Neither of the two conjugate alkene isomers prototypical of ˙R + O 2 reactions were detected. Products were also identified from decomposition of ketohydroperoxide-type radicals. The present work provides the first analysis of 2,3-dimethyloxirane oxidation chemistry and reveals that consumption pathways are complex and require the expansion of submechanisms in chemical kinetics mechanisms.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The atmospheric importance of methylamine additions to Criegee intermediates
We report that Criegee intermediates are important targets for study in atmospheric chemistry because of their capacity to oxidize airborne species. Among these species, ammonia has received critical attention for its presence in polluted agricultural or industrial areas and its role in forming particulate matter and condensation nuclei. Although methylamine has been given less attention than ammonia, both theoretical and experimental studies have demonstrated that the additional methyl substitution on the ammonia derivatives increases the rate constants for some systems. This suggests that the methylamine addition to Criegee intermediates could be more significant to atmospheric processes. In this work, geometries are optimized at the DF-CCSD(T)/ANO1 level for the methylamine addition reactions to the simplest Criegee intermediate and the anti- and syn-methylated Criegee intermediates. Energies for each stationary point were computed at the CCSD(T)/CBS level with corrections from the CCSDT(Q) method. Rate constants are obtained for each reaction using canonical transition state theory. Although methylamine addition proved to be a more favorable reaction relative to ammonia addition, the significantly lower concentration of atmospheric methylamine limits the prevalence of these reactions, even in the most optimal conditions. It is unlikely that the methylamine addition to Criegee intermediates will contribute significantly to the consumption of Criegee intermediates in the atmosphere.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Reaction mechanisms of a cyclic ether intermediate: Ethyloxirane
Oxiranes are a class of cyclic ethers formed in abundance during low-temperature combustion of hydrocarbons and biofuels, either via chain-propagating steps that occur from unimolecular decomposition of β -hydroperoxyalkyl radicals ( β -˙QOOH) or from reactions of HOÓ with alkenes. Ethyloxirane is one of four alkyl-substituted cyclic ether isomers produced as an intermediate from n -butane oxidation. While rate coefficients for β -˙QOOH → ethyloxirane + ÓH are reported extensively, subsequent reaction mechanisms of the cyclic ether are not. As a result, chemical kinetics mechanisms commonly adopt simplified chemistry to describe ethyloxirane consumption by convoluting several elementary reactions into a single step, which may introduce mechanism truncation error—uncertainty derived from missing or incomplete chemistry. Here we provide fundamental insight on reaction mechanisms of ethyloxirane in support of ongoing efforts to minimize mechanism truncation error. Reaction mechanisms are inferred from the detection of products during chlorine atom-initiated oxidation experiments using multiplexed photoionization mass spectrometry conducted at 10 Torr and temperatures of 650 K and 800 K. To complement the experiments, calculations of stationary point energies were conducted using the ccCA-PS3 composite method on ˙R + O 2 potential energy surfaces for the four ethyloxiranyl radical isomers, which produced barrier heights for 24 reaction pathways. In addition to products from ˙QOOH → cyclic ether + ÓH and ˙R + O 2 → conjugate alkene + HOÓ, both of which were significant pathways and are prototypical to alkane oxidation, other species were identified from ring-opening of both ethyloxiranyl and ˙QOOH radicals. The latter occurs when the unpaired electron is localized on the ether group, causing the initial ˙QOOH structure to ring-open and form a resonance-stabilized ketohydroperoxide-type radical. The present work provides the first analysis of ethyloxirane oxidation chemistry, which reveals that consumption pathways are complex and may require an expansion of submechanisms to increase the fidelity of chemical kinetics mechanisms.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The water trimer reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O
All important stationary points on the potential energy surface (PES) for the reaction OH + (H 2 O) 3 → (H 2 O) 2 OH + H 2 O have been fully optimized using the “gold standard” CCSD(T) method with the large Dunning correlation-consistent cc-pVQZ basis sets. Three types of pathways were found. For the pathway without hydrogen abstraction, the barrier height of the transition state (TS1) is predicted to lie 5.9 kcal mol -1 below the reactants. The two major complexes (H 2 O) 3 ∙∙∙OH (CP1 and CP2a) are found to lie 6.3 and 11.0 kcal mol -1 , respectively, below the reactants [OH + (H 2 O) 3 ]. For one of the H-abstraction pathways the lowest classical barrier height is predicted to be much higher, 6.1 kcal mol -1 (TS2a) above the reactants. For the other H-abstraction pathway the barrier height is even higher, 15.0 (TS3) kcal mol -1 . Vibrational frequencies and the zero-point vibrational energies connected to the PES are also reported. Fianlly, the energy barriers for the H-abstraction pathways are compared with those for the OH + (H 2 O) 2 and OH + H 2 O reactions, and the effects of the third water on the energetics are usually minor (0.2 kcal mol -1 ).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Energetics and mechanisms for the acetonyl radical + O 2 reaction: An important system for atmospheric and combustion chemistry
The acetonyl radical (•CH 2 COCH 3 ) is relevant to atmospheric and combustion chemistry due to its prevalence in many important reaction mechanisms. One such reaction mechanism is the decomposition of Criegee intermediates in the atmosphere that can produce acetonyl radical and OH. In order to understand the fate of the acetonyl radical in these environments and to create more accurate kinetics models, we have examined the reaction system of the acetonyl radical with O 2 using highly reliable theoretical methods. Structures were optimized using coupled cluster theory with singles, doubles, and perturbative triples [CCSD(T)] with an atomic natural orbital (ANO0) basis set. Energetics were computed to chemical accuracy using the focal point approach involving perturbative treatment of quadruple excitations [CCSDT(Q)] and basis sets as large as cc-pV5Z. The addition of O 2 to the acetonyl radical produces the acetonylperoxy radical, and multireference computations on this reaction suggest it to be barrierless. Additionally, no submerged pathways were found for the unimolecular isomerization of the acetonylperoxy radical. Besides dissociation to reactants, the lowest energy pathway available for the acetonylperoxy radical is a 1-5 H shift from the methyl group to the peroxy group through a transition state that is 3.3 kcal mol -1 higher in energy than acetonyl radical + O 2 . The ultimate products from this pathway are the enol tautomer of the acetonyl radical along with O 2 . Multiple pathways that lead to OH formation are considered; however, all of these pathways are predicted to be energetically inaccessible, except at high temperatures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Relatives of cyanomethylene: replacement of the divalent carbon by B − , N + , Al − , Si, P + , Ga − , Ge, and As +
The lowest lying singlet and triplet states of nine relatives of cyanomethylene are studied with highly rigorous ab initio methods, and periodic trends in their electronic structures are analyzed.
Abbott, Boyi Z.↗
Important features of the potential energy surface of the methylamine plus O( 1 D) reaction
This research presents an ab initio characterization of the potential energy surface for the methylamine plus 1D oxygen atom reaction, which may be relevant to interstellar chemistry. Geometries and harmonic vibrational frequencies were determined for all stationary points at the CCSD(T)/aug-cc-pVTZ level of theory. The focal point method along with several additive corrections was used to obtain reliable CCSDT(Q)/CBS potential energy surface features. Extensive conformational analysis and intrinsic reaction coordinate computations were performed to ensure accurate chemical connectivity of the stationary points. Additionally, five minima were determined to be possible products of this reaction and three novel transition states were found that were previously unreported or mislabeled in the literature. The pathways we present can be used to guide further searches for NH 2 containing species in the interstellar medium.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Characterization of the 2-methylvinoxy radical + O 2 reaction: A focal point analysis and composite multireference study
Vinoxy radicals are involved in numerous atmospheric and combustion mechanisms. High-level theoretical methods have recently shed new light on the reaction of the unsubstituted vinoxy radical with O 2 . The reactions of 1-methylvinoxy radical and 2-methylvinoxy radical with molecular oxygen have experimental high pressure limiting rate constants, k ∞ , 5–7 times higher than that of the vinoxy plus O 2 reaction. In this work, high-level ab initio quantum chemical computations are applied to the 2-methylvinoxy radical plus O 2 system, namely, the formation and isomerization of the 1-oxo-2-propylperoxy radical, the immediate product of O 2 addition to the 2-methylvinoxy radical. Multireference methods were applied to the entrance channel. No barrier to O 2 addition could be located, and more sophisticated treatment of dynamic electron correlation shows that the principal difference between O 2 addition to the vinoxy and 2-methylvinoxy radicals is a larger steric factor for 2-methylvinoxy + O 2 . This is attributed to the favorable interaction between the incoming O 2 molecule and the methyl group of the 2-methylvinoxy radical. Via the focal point approach, energetics for this reaction were determined, in most cases, to chemical accuracy. The coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] correlation energy and Hartree–Fock energies were independently extrapolated to the complete basis set limit. A correction for the effect of higher excitations was computed at the CCSDT(Q)/6-31G level. Corrections for the frozen-core approximation, the Born–Oppenheimer approximation, the nonrelativistic approximation, and the zero-point vibrational energy were included. From the 1-oxo-2-propylperoxy radical, dissociation to reactants is competitive with the lowest energy isomerization pathway. The lowest energy isomerization pathway ultimately forms acetaldehyde, CO, and ·OH as the final products.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Riddles of the structure and vibrational dynamics of HO 3 resolved near the ab initio limit
The hydridotrioxygen (HO 3 ) radical has been explored in many previous theoretical and experimental studies over several decades, originally because of its possible relevance to the tropospheric HO x 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 D 0 (HO–O 2 ) 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 HO 3 and DO 3 , 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-HO 3 product. However, 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 HO 3 and DO 3 , respectively.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The addition of methanol to Criegee intermediates
Bimolecular reactions involving stabilized Criegee intermediates (SCI) have been the target of many studies due to the role these molecules play in atmospheric chemistry. Recently, kinetic rates for the addition reaction of the simplest SCI (formaldehyde oxide) and its methylated analogue (acetone oxide) with methanol were reported both experimentally and theoretically. We re-examine the energy profile of these reactions by employing rigorous ab initio methods. Optimized CCSD(T)/ANO1 geometries are reported for the stationary points along the reaction path. Energies are obtained at the CCSD(T)/CBS level of theory. Contributions of full triple and quadruple excitations are computed to assess the convergence of this method. Rate constants are obtained using conventional canonical transition state theory under the rigid rotor harmonic oscillator approximation and with the inclusion of a one-dimensional hindered rotor treatment. These corrections for internal rotations have a significant impact on computed kinetic rate constants. With this approach, we compute rate constants for the addition of methanol to formaldehyde oxide (H 2 COO) and acetone oxide [(CH 3 ) 2 COO] at 298.15 K as (1.2 ± 0.8) × 10 -13 and (2.8 ± 1.3) × 10 -15 cm 3 s -1 , respectively. Additionally, we investigate the temperature dependence of the rate constant, concluding that the transition state barrier height and tunneling contributions shape the qualitative behaviour of these reactions.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Convergent energies and anharmonic vibrational spectra of Ca2H2 and Ca2H4 constitutional isomers
Three constitutional isomers of both Ca 2 H 2 and Ca 2 H 4 have been characterized with molecular electronic structure theory. Correlation methods as complete as CCSDT(Q) and basis sets as large as cc-pwCV5Z have been used to converge the relative energies within chemical accuracy (≤1 kcal mol -1 ). Anharmonic vibrational frequencies were computed using second-order vibrational perturbation theory employing CCSD(T)/cc-pwCVTZ cubic and quartic force-fields and a CCSD(T)/cc-pwCVQZ quadratic force field. The monobridged [Ca(μ 2 -H)CaH] and dibridged [Ca(μ 2 -H) 2 Ca] isomers of Ca 2 H 2 were predicted to lie 6.5 and 12.9 kcal mol -1 below the energy of the classical HCaCaH linear isomer, respectively. Despite the energetic favorability of the bridged Ca 2 H 2 isomers, we conclude (surprisingly) that only the higher energy linear structure has been observed in the laboratory. At 0 K, the tribridged [Ca(μ 2 -H) 3 CaH] isomer of Ca 2 H 4 is predicted to be enthalpically favored by 0.9 kcal mol -1 in comparison to the enthalpy of the dibridged [HCa(μ 2 -H) 2 CaH] structure. Comparison of experiment with our computed frequencies suggests that the observed vibrational features arise from both the dibridged and the tribridged Ca 2 H 4 structures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
tert -Butyl peroxy radical: ground and first excited state energetics and fundamental frequencies
Alkylperoxy radicals (RO 2 ˙) are key intermediates in combustion and atmospheric oxidation processes. As such, reliable detection and monitoring of these radicals can provide a wealth of information about the underlying chemistry. The tert-butyl peroxy radical is the archetypal tertiary peroxy radical, yet its vibrational spectroscopy is largely unexplored. To aid in future experimental investigations, we have performed high-level theoretical studies of the fundamental vibrational frequencies of the ground- and first excited states. A conformer search on both electronic surfaces reveals single minimum-energy structures. We predict an à 2 A' ← X 2 A" adiabatic excitation energy of 7738 cm -1 via focal point analysis, approximating the CCSDT(Q)/CBS level of theory. This excitation energy agrees to within 17 cm -1 of the most accurate experimental measurement. We compute CCSD(T) fundamental vibrational frequencies via second-order vibrational perturbation theory (VPT2), using a hybrid force field in which the quadratic (cubic/quartic) force constants are evaluated with the ANO1 (ANO0) basis set. Anharmonic resonance polyads are treated with the VPT2 + K effective Hamiltonian approach. Among the predicted fundamental frequencies, the ground state O–O stretch, excited state O–O stretch, and excited state C–O–O bend fundamentals are predicted at 1138, 959, and 490 cm -1 , respectively. Basis set sensitivity is found to be particularly great for the O–O stretches, similar to what has already been noted in smaller, unbranched peroxy radicals. Exempting these O–O stretches, agreement with the available experimental fundamentals is generally good (±10 cm -1 ).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗