Search

A new strategy is presented for computing anharmonic partition functions for the motion of adsorbates relative to a catalytic surface. Importance sampling is compared with conventional Monte Carlo. The importance sampling is significantly more efficient. This new approach is applied to CH3* on Ni(111) as a test case. The motion of methyl relative to the nickel surface is found to be anharmonic, with significantly higher entropy compared to the standard harmonic oscillator model. The new method is freely available as part of the Minima-Preserving Neural Network within the A D T HERM package. [ABSTRACT FROM AUTHOR]

Fluorinated species have a pivotal role in semiconductor material chemistry and some of them have been detected beyond the Earth's atmosphere. Achieving good energy accuracy on fluorinated species using quantum chemical calculations has long been a challenge. In addition, obtaining direct experimental thermochemical quantities has also proved difficult. Here, we report the threshold photoelectron and photoion yield spectra of SiF and CF radicals generated with a fluorine reactor. The spectra were analysed with the support of ab initio calculations, resulting in new experimental values for the adiabatic ionisation energies of both CF (9.128 ± 0.006 eV) and SiF (7.379 ± 0.009 eV). Using these values, the underlying thermochemical network of Active Thermochemical Tables was updated, providing further refined enthalpies of formation and dissociation energies of CF, SiF, and their cationic counterparts. [ABSTRACT FROM AUTHOR]

Empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N2, and H2O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm−1 (0.12 kJ mol−1), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches. [ABSTRACT FROM AUTHOR]

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 cm³ 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 HCHO + H2O - HOCH2OH (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). [ABSTRACT FROM AUTHOR]

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). [ABSTRACT FROM AUTHOR]

A number of economical modifications to the high-accuracy extrapolated ab initio thermochemistry (HEAT) model chemistry are evaluated. The two resulting schemes, designated as mHEAT and mHEAT+, are designed for efficient and pragmatic evaluation of molecular energies in systems somewhat larger than can be practically studied by the unapproximated HEAT scheme. It is found that mHEAT+ produces heats of formation with nearly subchemical (±1 kJ/mol) accuracy at a substantially reduced cost relative to the full scheme. Total atomization energies calculated using the new thermochemical recipes are compared to the results of the HEAT-345(Q) model chemistry, and enthalpies of formation for the three protocols are also compared to Active Thermochemical Tables. Finally, a small selection of transition states is studied using mHEAT and mHEAT+, which illuminates some interesting features of reaction barriers and serves as an initial benchmark of the performance of these model chemistries for chemical kinetics applications. [ABSTRACT FROM AUTHOR]

A combined (fixed-J) two-dimensional master-equation/semi-classical transition state theory/variational Rice-Ramsperger-Kassel-Marcus approach has been used to compute reaction rate coefficients of •OH with CH3OH over a wide range of temperatures (10–2500 K) and pressures (10−1–104 Torr) based on a potential energy surface that has been constructed using a modification of the high accuracy extrapolated ab initio thermochemistry (HEAT) protocol. The calculated results show that the title reaction is nearly pressure-independent when T > 250 K but depends strongly on pressure at lower temperatures. In addition, the preferred mechanism and rate constants are found to be very sensitive to temperature. The reaction pathway CH3OH + •OH → CH3O• + H2O proceeds exclusively through tunneling at exceedingly low temperatures (T ≤ 50 K), typical of those established in interstellar environments. In this regime, the rate constant is found to increase with decreasing temperature, which agrees with low-temperature experimental results. The thermodynamically favored reaction pathway CH3OH + •OH → •CH2OH + H2O becomes dominant at higher temperatures (T ≥ 200 K), such as those found in Earth's atmosphere as well as combustion environments. By adjusting the ab initio barrier heights slightly, experimental rate constants from 200 to 1250 K can be satisfactorily reproduced. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The best available enthalpy of formation of H(O.sub.2) is obtained through the use of the Active Thermochemical Table. The new thermochemistry of H(O.sub.2), together with other arguments suggests that the former equilibrium constant for NO + H(O.sub.2) -- OH + N(O.sub.2), which is of major importance in stratospheric ozone depletion, was underestimated by the factor of ~2.

The best currently available values for the enthalpy of formation of H[O.sub.2], and new Active Thermochemical Table (ATcT)values for several other species, such as NO, N[O.sub.2], and OH were used on the basis of ATcT approach. The experiments were performed in [[C.sub.2][H.sub.5]I/N[O.sub.2]] mixtures, using the reflected shock tube technique and OH-radical electronic absorption detection and using a multipass optical system.

Born-Oppenheimer direct dynamics simulations were performed to study atomistic details of the F + CH3CN→HF + CH2CN H-atom abstraction reaction. The simulation trajectories were calculated with a combined M06-2X/MP2 algorithm utilizing the 6-311++G** basis set. The experiments were performed at 300 K, and assuming the accuracy of transition state theory (TST), the trajectories were initiated at the F… HCH2CN abstraction TS with a 300 K Boltzmann distribution of energy and directed towards products. Recrossing of the TS was negligible, confirming the accuracy of TST. HF formation was rapid, occurring within 0.014 ps of the trajectory initiation. The intrinsic reaction coordinate (IRC) for reaction involves rotation of HF about CH2CN and then trapping in the CH2CN… HF post-reaction potential energy well of ~10 kcal/mol with respect to the HF + CH2CN products. In contrast to this IRC, five different trajectory types were observed: the majority proceeded by direct H-atom transfer and only 11% approximately following the IRC. The HF vibrational and rotational quantum numbers, n and J, were calculated when HF was initially formed and they increase as potential energy is released in forming the HF +CH2CNproducts. The population of the HF product vibrational states is only in qualitative agreement with experiment, with the simulations showing depressed and enhanced populations of the n = 1 and 2 states as compared to experiment. Simulations with an anharmonic zero-point energy constraint gave product distributions for relative translation, HF rotation, HF vibration, CH2CN rotation, and CH2CN vibration as 5%, 11%, 60%, 7%, and 16%, respectively. In contrast, the experimental energy partitioning percentages to HF rotation and vibration are6%and 41%. Comparisons are made between the current simulation and those for other F + H-atom abstraction reactions. The simulation product energy partitioning and HF vibrational population for F +CH3CN→HF+CH2CNresemble those for other reactions. A detailed discussion is given of possible origins of the difference between the simulation and experimental energy partitioning dynamics for F +CH3CN→HF +CH2CN. The F +CH3CNreaction also forms the CH3C(F)N intermediate, in which the F-atom adds to the C#8801;N bond. However, this intermediate and F… CH3CN and CH3CN… F van der Waals complexes are not expected to affect the F + CH3CN→HF + CH2CN product energy partitioning. [ABSTRACT FROM AUTHOR]

The concept behind active thermochemical tables(ATcT), a paradigm of how to derive accurate, reliable, and internally consistent thermochemistry, is presented. As opposed to the traditional sequential thermochemistry, ATcT uses the thermochemical network approach.