Engineering Papers⌕ Search

Engineering topics

Ruscic, Branko

Publications and source records attributed to Ruscic, Branko.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Dynamics of HCN, HNC, and HNCO Formation in the 193 nm Photodissociation of Formamide

Formamide (NH 2 CHO) is the simplest molecule containing a peptide linkage [-NH-C($=$O)-], and it plays an essential role in the study of prebiotic chemistry. Exposure to UV irradiation allows formamide to decompose and act as a prebiotic feedstock in the formation of nucleobases and other necessary starting materials. The photodissociation mechanism of gaseous formamide at 193 nm is studied using (a) chirped-pulse Fourier transform millimeter-wave spectroscopy in the 260-290 GHz spectral region in a room-temperature flow-tube reactor at 1 mu bar pressure, (b) a combination of electronic structure theory, transition state theory, and quasiclassical trajectories, and (c) the Active Thermochemical Tables. The HCN and HNC photoproducts of hydrogenated (NH 2 CHO) and deuterated (NH 2 CDO and ND 2 CHO) formamide precursors are examined to gain insight into the photodissociation mechanism. Here, the theoretical investigation has characterized the main pathway leading to each of the HCN/HNC isomers from the precursor isotopologues. The theoretical branching ratio [HNC]/[HCN] = 2.1 for nascent photofragments agrees with the experiment. The effect of the postphotolysis HNC ↔ HCN isomerization on the [HNC]/[HCN] ratio is predicted. We report the experimental branching ratio [HNCO]: ([HNC] + [HCN]) = 12 ± 3 and propose that most of HNCO originates from dissociation on the S-1 electronic state of formamide.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

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, k 1 = 7.5 × 10 −12 cm 3 s −1 . High-level accuracy coupled cluster calculations in combination with E,J -resolved two-dimensional master equation analyses yield k 1 = (2.4 ± 0.5) × 10 −12 cm 3 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 + H 2 O ⇌ HOCH 2 OH (2) with a HOCH 2 OH to HCHO ratio of only ca . 2%. Although HOCH 2 OH 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 HOCH 2 OH reacting with OH. The best estimate, 3.3 Tg/y, is about 30 times less than recently reported. The theoretical equilibrium constant K eq (2) determined in this work also allows us to estimate the Henry’s law constant of methanediol (8.1 × 10 5 M atm −1 at 280 K).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

The C 2 H 4 O isomers in the oxidation of ethylene

A combined rotational spectroscopy, thermochemistry, and kinetic modeling study explores the mechanism of acetaldehyde (ethanal), vinyl alcohol (ethenol), and ethylene oxide (oxirane) formation in the oxidation of ethylene (ethene). Multiplexed quantitative detection of oxidation intermediates and products exiting a SiO 2 / SiC microreactor at 1700 K is demonstrated with the BrightSpec W-band chirped-pulse Fourier transform millimeter-wave spectrometer. The broadband rotational spectrum contains transitions of formaldehyde (CH 2 O), methoxy (CH 3 O), methanol (CH 3 OH), ketene (CH 2 CO), acetaldehyde (CH 3 CHO), syn-vinyl alcohol (syn- CH 2 CHOH), anti-vinyl alcohol (anti-CH 2 CHOH), oxirane (c-CH 2 OCH 2 ), propyne (CH 3 CCH), and syn-propanal (syn-CH 3 CH 2 CHO). We focus on the three C 2 H 4 O species and deduce their concentration ratio [CH 3 CHO]:[CH 2 CHOH]:[c-CH 2 OCH 2 ] = 1:0.7(2):0.06(2) by comparing the observed line intensities to those simulated with the PGOPHER software. Detailed thermochemistry of the C 2 H 4 O isomers and conformers is provided by Active Thermochemical Tables. The observed excess concentrations of vinyl alcohol and oxirane relative to the more stable acetaldehyde compared to the equilibrium ratio at 1700 K (1:0.087:0.000024), point to direct chemical pathways to these higher energy isomers. The mechanism for the formation of the three C 2 H 4 O isomers is analyzed using a 0-D homogenous reactor kinetics simulation for ethylene oxidation. The ratios of the C 2 H 4 O isomers concentrations predicted by the kinetic model are compared to the experimental values.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Refining the thermochemical properties of CF, SiF, and their cations by combining photoelectron spectroscopy, quantum chemical calculations, and the Active Thermochemical Tables approach

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Active Thermochemical Tables: Enthalpies of Formation of Bromo- and Iodo-Methanes, Ethenes and Ethynes

Here, the thermochemistry of halocarbon species containing iodine and bromine is examined through an extensive interplay between new Feller–Peterson–Dixon (FPD) style composite methods and a detailed analysis of all available experimental and theoretical determinations using the thermochemical network that underlies the Active Thermochemical Tables (ATcT). From the computational viewpoint, a slower convergence of the components of composite thermochemistry methods is observed relative to species that solely contain first row elements, leading to a higher computational expense for achieving comparable levels of accuracy. Potential systematic sources of computational uncertainty are investigated, and, not surprisingly, spin-orbit coupling is found to be a critical component, particularly for iodine containing molecular species. The ATcT analysis of available experimental and theoretical determinations indicates that prior theoretical determinations have significantly larger uncertainties than originally reported, particularly in cases where molecular spin-orbit effects were ignored. Accurate and reliable heats of formation are reported for 38 halogen containing systems, based on combining the current computations with previous experimental and theoretical work via the ATcT approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Configuration Space Integration for Adsorbate Partition Functions: The Effect of Anharmonicity on the Thermophysical Properties of CO–Pt(111) and CH 3 OH–Cu(111)

A method for computing anharmonic thermophysical properties for adsorbates on metal surfaces has been extended to include libration, or frustrated rotation. Classical phase space integration is used with Monte Carlo sampling of the configuration space to obtain the partition function of CO on Pt(111) and CH 3 OH on Cu(111). A minima-preserving neural network potential energy surrogate is used within the integration routines. Direct state counting using discrete variable representation is used to benchmark the results. We find that the phase space integration approach is in excellent agreement with the direct state counting results. Comparison with standard models such as the harmonic oscillator indicates that anharmonicity contributes significantly to the thermodynamic properties of CH 3 OH on Cu(111). We find that there is also a considerable difference between the harmonic oscillator and phase space integration for CO on Pt(111), although the discrepancy can largely be attributed to the presence of multiple binding sites within the unit cell. We demonstrate that a multisite harmonic oscillator model might be sufficient for CO-Pt(111). A more thorough description of the potential energy surface, which can be achieved with phase space integration, is necessary for weakly bound adsorbates such as CH 3 OH. In conclusion, the thermophysical properties were used to calculate free energies of adsorption on the respective metals, and subsequently the equilibrium constants and Langmuir isotherms in relevant temperature ranges. The results show that the choice of model to obtain partition functions greatly affects the resulting surface coverages in kinetic models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Sub 20 cm -1 computational prediction of the CH bond energy – a case of systematic error in computational thermochemistry

Here, the bond dissociation energy of methylidyne, D 0 (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 D 0 (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 D 0 (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 D 0 (CH + ), and 85 831.0 ± 6.0 cm -1 for IE(CH).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

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

Here, 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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Mechanism, thermochemistry, and kinetics of the reversible reactions: C 2 H 3 + H 2 ⇌ C 2 H 4 + H ⇌ C 2 H 5

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 C 2 H 3 , C 2 H 4 , and C 2 H 5 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 –here we provide high-level theoretical results for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗