Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transition state”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The P( 4 S) + NH( 3 Σ – ) and N( 4 S) + PH( 3 Σ – )reactions as sources of interstellar phosphorus nitride

Phosphorus nitride (PN) is believed to be one of the major reservoirs of phosphorus in the interstellar medium (ISM). For this reason, understanding which reactions produce PN in space and predicting their rate coefficients is important for modelling the relative abundances of P-bearing species and clarifying the role of phosphorus in astrochemistry. Here, in this work, we explore the potential energy surfaces of the $\textrm{P}(^4\textrm{S}) + \textrm{NH}(^3\Sigma^-)$ and $\textrm{N}(^4\textrm{S}) + \textrm{PH}(^3\Sigma^-)$ reactions and the formation of $\textrm{H}(^2\textrm{S}) + \textrm{PN}(^1\Sigma^+)$ through high accuracy ab initio calculations and the variable reaction coordinate transition state theory (VRC-TST). We found that both reactions proceed without an activation barrier and with similar rate coefficients that can be described by a modified Arrhenius equation ( $k(T)=\alpha\!\left( T/300 \right)^{\beta} \exp\!{(\!-\!\gamma/T)})$ with $\alpha=0.93\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=0.24\, \rm K$ for the $\textrm{P} + \textrm{NH} \longrightarrow \textrm{H} + \textrm{PN}$ reaction and $\alpha=0.88\times 10^{-10}\rm cm^3\,s^{-1}$ , $\beta=-0.18$ and $\gamma=1.01\, \rm K$ for the $\textrm{N} + \textrm{PH} \longrightarrow \textrm{H} + \textrm{PN}$ one. Both reactions are expected to be relevant for modelling PN abundances even in the cold environments of the ISM. Given the abundance of hydrogen in space, we have also predicted rate coefficients for the destruction of PN via H + PN collisions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

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

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗
DOE OSTI2021

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

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

Excited-State Dynamics and Nonlinear Optical Properties of Hyperpolarizable Chromophores Based on Conjugated Bis(terpyridyl)Ru(II) and Palladium and Platinum Porphyrinic Components: Impact of Heavy Metals upon Supermolecular Electro-Optic Properties

A new series of strongly coupled oscillators based upon (porphinato)Pd, (porphinato)Pt, and bis(terpyridyl)-ruthenium(II) building blocks is described. These RuPPd, RuPPt, RuPPdRu, and RuPPtRu chromophores feature bis- (terpyridyl)Ru(II) moieties connected to the (porphinato)metal unit via an ethyne linker that bridges the 4!-terpyridyl and porphyrin macrocycle meso-carbon positions. Pump–probe transient optical data demonstrate sub-picosecond excited singlet-to-triplet-state relaxation. The relaxed lowest-energy triplet (T 1 ) excited states of these chromophores feature absorption manifolds that span the 800–1200 nm spectral region, microsecond triplet-state lifetimes, and large absorptive extinction coefficients [ε(T 1 → T n ) > 4 × 10 4 M –1 cm –1 ]. Dynamic hyperpolarizability (β λ ) values were determined from hyper-Rayleigh light scattering (HRS) measurements carried out at several incident irradiation wavelengths over the 800–1500 nm spectral region. Relative to benchmark RuPZn and RuPZnRu chromophores which showed large βHRS values over the 1200–1600 nm range, RuPPd, RuPPt, RuPPdRu, and RuPPtRu displayed large βHRS values over the 850–1200 nm region. Generalized Thomas–Kuhn sum (TKS) rules and experimental hyperpolarizability values were utilized to determine excited state-to-excited state transition dipole terms from experimental electronic absorption data and thus assessed frequency-dependent β λ values, including two- and three-level contributions for both β zzz and β xzx tensor components to the RuPPd, RuPPt, RuPPdRu, and RuPPtRu hyperpolarizability spectra. Furthermore, these analyses qualitatively rationalize how the βzzz and βxzx tensor elements influence the observed irradiation wavelength-dependent hyperpolarizability magnitudes. The TKS analysis suggests that supermolecules related to RuPPd, RuPPt, RuPPdRu, and RuPPtRu will likely feature intricate dependences of experimentally determined βHRS values as a function of irradiation wavelength that derive from substantial singlet–triplet mixing, and complex interactions among multiple different β tensor components that modulate the long wavelength regime of the nonlinear optical response.

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

Alkali Metal Control of Triplet-Mediated C–H Activation in Iron-Mediated Coupling of Dinitrogen and Benzene

Direct coupling of N 2 with abundant feedstocks like benzene to form N-containing organic compounds is a promising strategy for N 2 fixation pathways. The challenge of coupling N 2 activation and C–H bond oxidative addition was recently solved by introducing a reversible benzene C–H bond activation process mediated by a low-valent Fe(0) complex, which gave an organometallic product that could couple with partially reduced N 2 . Interestingly, the energetics of the C–H oxidative addition/reductive elimination step depends on the choice of alkali metal. However, the reason why the alkali metal influences the C–H bond activation remained elusive. In this work, we present a comprehensive study on this Fe(0)-mediated reversible C–H activation. Through density functional theory combined with high-level coupled cluster calculations, we discovered that the intermediate-spin triplet (S = 1) controls the energy of the transition state for C–H cleavage, while the high-spin quintet (S = 2) controls the position of the equilibrium. Na + drives the equilibrium toward oxidative addition due to an electrostatic effect, while K + and Rb + are dominated by a steric effect that favors the iron(0) species. The key role played by nonbonding interactions in the Fe-mediated C–H activation provides a conceptual model for alkali control over organometallic transformations.

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

Electronic Structure Theory Calculations Using Modern Architectures: KNL vs Haswell

The time to solution and parallel efficiency of several commonly used electronic structure methods (Hartree–Fock, density functional theory, second order perturbation theory, resolution of the identity second order perturbation theory, coupled cluster) are evaluated on both the Intel Xeon Haswell and the Intel Xeon Phi Knights Landing (KNL) architectures. The Haswell completes the benchmark calculations with a faster time to solution than the KNL for all molecules and methods tested. While the Haswell exhibits an average speedup of at least 3.5 relative to the KNL for all nonthreaded computations, the KNL has a better parallel efficiency than the Haswell with increasing core counts. The architectures are further tested using a more computationally costly coupled cluster method on a transition state reaction. The Haswell appears to be the best choice to minimize the time to solution, though for very large systems and high levels of theory that require memory intensive processes the superior memory hierarchy and larger on node memory of the KNL can make it a better choice. These results are used to showcase aspects of novel architectures that will increase efficiency for quantum chemistry applications.

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

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

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

Gas-Phase Complexes of Americium and Lanthanides with a Bis-triazinyl Pyridine: Reactivity and Bonding of Archetypes for F-Element Separations

Bis-triazinyl pyridines (BTPs) exhibit solution selectivity for trivalent americium over lanthanides (Ln), the origins of which remain uncertain. Here in this paper, electrospray ionization was used to generate gas-phase complexes [ML 3 ] 3+ , where M = La, Lu, or Am and L is EtBTP 2,6-bis(5,6-diethyl-1,2,4-triazin-3-yl)-pyridine. Collision-induced dissociation (CID) of [ML 3 ] 3+ in the presence of H 2 O yielded a protonated ligand [L(H)] + and hydroxide [ML 2 (OH)] 2+ or hydrate [ML(L–H)(H 2 O)] 2+ , where (L–H) - is a deprotonated ligand. Although solution affinities indicate stronger binding of BTPs toward Am 3+ versus Ln 3+ , the observed CID process is contrastingly more facile for M = Am versus Ln. To understand the disparity, density functional theory was employed to compute potential energy surfaces for two possible CID processes, for M = La and Am. In accordance with the CID results, both the rate determining transition state barrier and the net energy are lower for [AmL 3 ] 3+ versus [LaL 3 ] 3+ and for both product isomers, [ML 2 (OH)] 2+ and [ML(L–H)(H 2 O)] 2+ . More facile removal of a ligand from [AmL 3 ] 3+ by CID does not necessarily contradict stronger Am 3+ –L binding, as inferred from solution behavior. In particular, the formation of new bonds in the products can distort kinetics and thermodynamics expected for simple bond cleavage reactions. In addition to correctly predicting the seemingly anomalous CID behavior, the computational results indicate greater participation of Am 5f versus La 4f orbitals in metal–ligand bonding.

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

Domain Separated Density Functional Theory for Reaction Energy Barriers and Optical Excitations

We recently proposed Domain Separated Density Functional Theory (DS-DFT), a framework that allows for the combination of different levels of theory for the computation of the electronic structure of molecules. Here, this work discusses the application of DS-DFT to the computation of transition-state energy barriers and optical absorption spectra. We considered several hydrogen abstraction reactions and optical spectra of molecule/metal cluster systems, including the absorption of individual species such as carbon monoxide, methane, and molecular hydrogen to a Li6 cluster. We present and discuss two domain-separated methods: i), the screened-density approximation (SDA), and, ii), linearly weighted exchange (LWE). We find that SDA, which is applied as a hybridization based on atomic domains, could be useful to computing energy barriers, whereas LWE is suited for the analysis of electronic properties such as ground-state gaps, excitation energies, and oscillator strengths.

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

Theoretical Study of the Phenoxy Radical Recombination with the O( 3 P) Atom, Phenyl plus Molecular Oxygen Revisited

Quantum chemical calculations of the C 6 H 5 O 2 potential energy surface (PES) were carried out to study the mechanism of the phenoxy + O( 3 P) and phenyl + O 2 reactions. CASPT2(15e,13o)/CBS//CASSCF(15e,13o)/DZP multi-reference calculations were utilized to map out the minimum energy path for the entrance channels of the phenoxy + O( 3 P) reaction. Stationary points on the C 6 H 5 O 2 PES were explored at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311++G** level for the species with a single-reference character of the wave function and at the CASPT2(15e,13o)/CBS//B3LYP/6-311++G** level of theory for the species with a multi-reference character of the wave function. Conventional, variational, and variable reaction coordinate transition state theories were employed in Rice-Ramsperger-Kassel-Marcus Master Equation calculations to assess temperature- and pressure-dependent phenomenological rate constants and product branching ratios. Here, the main bimolecular product channels of the phenoxy + O( 3 P) reaction are concluded to be para/ortho-benzoquinone + H, 2,4-cyclopentadienone + HCO, and, at high temperatures, also phenyl + O 2 . The main bimolecular product channels of the phenyl + O 2 reaction include 2,4-cyclopentadienone + HCO at lower temperatures, and phenoxy + O( 3 P) at higher temperatures. Both for the phenoxy + O( 3 P) and phenyl + O 2 reactions, the collisional stabilization of peroxybenzene at low temperatures and high pressures competes with the bimolecular product channels.

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

Mechanistic Investigation of Dimethylmercury Formation Mediated by a Sulfide Mineral Surface

Mercury (Hg) pollution is a global environmental problem, but many key transformations among 29 various Hg species are not well understood. Although anaerobic microorganisms are the primary 30 producers of monomethylmercury (MMeHg), the abiotic formation of dimethylmercury (DMeHg) 31 from MMeHg have been suggested to account for a large portion of DMeHg in the ocean. It has 32 been shown that abiotic formation of DMeHg from MMeHg can be facilitated by reduced sulfur 33 groups on mineral surfaces. However, a detailed mechanistic understanding of this reaction is 34 lacking. Here we perform density functional theory (DFT) calculations to explore mechanisms of 35 DMeHg formation on the 110 surface of a CdS(s) (hawleyite) nanoparticle. We show that 36 coordination of MMeHg substituents to adjacent reduced sulfur groups protruding from the surface 37 facilitate DMeHg formation and that the reaction proceeds through direct transmethylation from 38 one MMeHg substituent to another. Solvation effects were found to play an important role in the 39 surface reconstruction of the nanoparticle and in decreasing the energetic barrier for DMeHg 40 formation. We also investigate DMeHg formation from MMeHg mediated by organic dithiols and 41 compare the results to those mediated by the mineral surface. In contrast to the surface reaction, 42 the calculations suggest that dithiol-mediated formation of DMeHg instead proceeds through a 43 mechanism in which Hg–Hg interactions activate a methyl group from one MMeHg substituent 44 for intramolecular transfer to the adjacent MMeHg substituent. These findings refine previously 45 proposed mechanisms of abiotic DMeHg formation and fill gaps in our understanding of global 46 Hg cycling in the environment.

Lian, Peng↗
DOE OSTI2022

Gas-Phase Study of the Elementary Reaction of the D1-Ethynyl Radical (C 2 D; X 2 Σ + ) with Propylene (C 3 H 6 ; X 1 A') under Single-Collision Conditions

The bimolecular gas-phase reactions of the D1-ethynyl radical (C 2 D; X 2 Σ + ) with propylene (C 3 H 6 ; X 1 A’) and partially substituted D3-3,3,3-propylene (C 2 H 3 CD 3 ; X 1 A’) were studied under single collision conditions utilizing the crossed molecular beams technique. Combining our laboratory data with electronic structure and statistical calculations, the D1-ethynyl radical is found to add without barrier to the C1 and C2 carbons of the propylene reactant, resulting in doublet C 5 H 6 D intermediate(s) with lifetime(s) longer than their rotational period(s). These intermediates undergo isomerization and unimolecular decomposition via atomic hydrogen loss through tight exit transition states forming predominantly cis/trans-3-penten-1-yne ((HCC)CH=CH(CH 3 )) and to a minor amount 3-methyl-3-buten-1-yne ((HCC)C(CH 3 )=CH 2 ) via overall exoergic reactions. Although the title reaction does not lead to the cyclopentadiene molecule (c-C 5 H 6 , X 1 A 1 ), high temperature environments can convert the identified acyclic C 5 H 6 isomers through hydrogen atom assisted isomerization to cyclopentadiene (c-C 5 H 6 , X 1 A 1 ). Since both the ethynyl radical and propylene reactants have been observed in cold interstellar environments such as TMC-1 and the reaction is exoergic and all barriers lie below the energy of the separated reactants, the these C 5 H 6 product isomers are predicted to form in those low temperature regions.

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

Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase

Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Electron Binding Energy Spectra of Al n Pt – Clusters—A Combined Experimental and Computational Study

Results of size-selected electron photo-detachment experiments and density functional theory calculations on anionic Al n Pt – , n = 1–7, clusters are presented and analyzed. The measured and calculated spectra of electron binding energies are, overall, in excellent accord with each other. The analysis reveals the general importance of accounting for the multiplicity of structural forms of a given-size cluster that can contribute to its measured spectrum, especially when the clusters are fluxional and/or the conditions of the experiment allow for structural transitions. Here, we show that for the systems studied here, the size-specific peculiarities of the measured spectra can be understood in terms of the combined contributions of corresponding different accessible stable equilibrium conformations, bona-fide transition-state configurations, and electronic-crossing structures that may play the role of effective barriers in electronically nonadiabatic dynamics.

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

Probing High-Temperature Amine Chemistry: Is the Reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 Important?

The reaction NH 3 + NH 2 ⇌ N 2 H 3 + H 2 (R1) has been identified as a key step to explain experimental results for pyrolysis and oxidation of ammonia. However, no direct experimental or theoretical evidence for the reaction has been reported. In the present work, the reaction was studied by ab initio theory and by reinterpretation of experimental data. We could not locate a transition state for R1 occurring as a direct process, but alternative mechanisms yielded an upper bound to k 1 of 1.5 × 10 13 exp(–58.9 kcal mol –1 /RT) cm 3 mol –1 s –1 over 1000–2500 K, several orders of magnitude lower than values applied in modeling. Consistent with the theoretical work, re-evaluation of NH 3 pyrolysis data supported a very low value of k 1 . However, this finding opens up a novel unresolved issue. Current kinetic models cannot capture the NH 3 oxidation behavior in a number of laminar flow reactor and jet-stirred reactor experiments without adopting an improbably high value for k 1 . Furthermore, important oxidation steps might be underestimated or missing from mechanisms.

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

Conformer-Dependent Chemistry: Experimental Product Branching of the Vinyl Alcohol + OH + O 2 Reaction

The concentration of formic acid in Earth's tropo-sphere is underestimated by detailed chemical models compared to field observations. Phototautomerization of acetaldehyde to its less stable tautomer vinyl alcohol, followed by the OH-initiated oxidation of vinyl alcohol, has been proposed as a missing source of formic acid that improves the agreement between models and field measurements. Theoretical investigations of the OH + vinyl alcohol reaction in excess O 2 conclude that OH addition to the alpha carbon of vinyl alcohol produces formaldehyde + formic acid + OH, whereas OH addition to the β site leads to glycoaldehyde + HO 2 . Furthermore, these studies predict that the conformeric structure of vinyl alcohol controls the reaction pathway, with the anti-conformer of vinyl alcohol promoting α OH addition, whereas the syn-conformer promotes β addition. However, the two theoretical studies reach different conclusions regarding which set of products dominate. We studied this reaction using time-resolved multiplexed photoionization mass spectrometry to quantify the product branching fractions. Our results, supported by a detailed kinetic model, conclude that the glycoaldehyde product channel (arising mostly from syn-vinyl alcohol) dominates over formic acid production with a 3.6:1.0 branching ratio. Additionally, this result supports the conclusion of Lei et al. that conformer-dependent hydrogen bonding at the transition state for OH-addition controls the reaction outcome. As a result, tropospheric oxidation of vinyl alcohol creates less formic acid than recently thought, increasing again the discrepancy between models and field observations of Earth's formic acid budget.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗