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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Automated theoretical chemical kinetics: Predicting the kinetics for the initial stages of pyrolysis

Large scale implementation of high level computational theoretical chemical kinetics offers the prospect for dramatically improving the fidelity of combustion chemical modeling. To facilitate such efforts, we have developed a suite of codes, collectively referred to as AutoMech, that allow for the automatic prediction of the kinetics for large sets of reactions via ab initio transition-state-theory based master-equation calculations. The primary input is simply the mechanism, a dictionary relating chemically identifiable species descriptors (e.g., SMILES or InChIs) to species labels in the mechanism, and a specification of the electronic structure and transition state theory models to be implemented. Here we illustrate the current utility of AutoMech through a study of the initial stages of pyrolysis for 3 sets of fuels: sequences of alkanes, alcohols, and aldehydes. For simplicity, the analysis focuses on abstractions from the fuel by H, CH 3 , and OH, and the decomposition of the resulting radicals. Altogether, there are a total of 166 input channels in these sets (more than 363 forward reactions when expanded to the full set of elementary reactions). The code successfully produces high quality rate estimates (with apparent uncertainties less than a factor of two in limited comparisons with experiment) for > 95% of these. For the radical decomposition reactions, the analysis includes predictions for the pressure dependence of the kinetics. This wide-ranging exploration illustrates (i) the effect of different levels of prediction on the expected accuracy, (ii) the branching between abstractions at different sites for different abstractors, (iii) the dependence of the rates on the chemical structure, and (iv) the variation in radical stabilities across chemical families. These results, as well as the demonstrated feasibility of the methodology, should find further utility in the development of accurate rate expressions for arbitrary fuels. (c) 2020 The Combustion Institute. Published by Elsevier Inc. All rights reserved.

Elliott, Sarah N.↗
DOE OSTI2021

Theoretical kinetics predictions for NH 2 + HO 2

Recent modeling studies of NH 3 oxidation, which are motivated by the prospective role of ammonia as a zero-carbon fuel, have indicated significant discrepancies between existing literature mechanisms. In this study high level theoretical kinetics predictions have been obtained for the reaction of NH 2 with HO 2 , which has previously been highlighted as an important reaction with high sensitivity and high uncertainty. The potential energy surface is explored with coupled cluster calculations including large basis sets and high-level corrections to yield high accuracy (~0.2 kcal/mol) estimates of the stationary point energies. Variational transition state theory is used to predict the microcanonical rate constants, which are then incorporated in master equation treatments of the temperature and pressure dependent kinetics. For the radical-radical channels, the microcanonical rates are obtained from variable reaction coordinate transition state theory implementing directly evaluated multireference electronic energies. The analysis yields predictions for the total rate constant as well as the branching to the NH 3 + O 2 , H 2 NO + OH, and HNO + H 2 O channels. Rate constants are also reported for the H 2 NO + OH reaction as they arise naturally from the analysis. The rate constant and branching fraction determined in this work for the NH 2 + HO 2 reaction deviate significantly from values used in most previous modeling studies. The fact that the main product channel is chain terminating, rather than propagating, has strong implications for modeling NH 3 ignition and oxidation, in particular at intermediate temperatures and elevated pressure.

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

Transition Path Sampling Study of the Feruloyl Esterase Mechanism

Serine hydrolases cleave peptide and ester bonds and are ubiquitous in nature, with applications in biotechnology, in materials, and as drug targets. The serine hydrolase two-step mechanism employs a serine–histidine–aspartate/glutamate catalytic triad, where the histidine residue acts as a base to activate poor nucleophiles (a serine residue or a water molecule) and as an acid to allow the dissociation of poor leaving groups. This mechanism has been the subject of debate regarding how histidine shuttles the proton from the nucleophile to the leaving group. To elucidate the reaction mechanism of serine hydrolases, we employ quantum mechanics/molecular mechanics-based transition path sampling to obtain the reaction coordinate using the Aspergillus niger feruloyl esterase A (AnFaeA) as a model enzyme. The optimal reaction coordinates include terms involving nucleophilic attack on the carbonyl carbon and proton transfer to, and dissociation of, the leaving group. During the reaction, the histidine residue undergoes a reorientation on the time scale of hundreds of femtoseconds that supports the “moving histidine” mechanism, thus calling into question the “ring flip” mechanism. We find a concerted mechanism, where the transition state coincides with the tetrahedral intermediate with the histidine residue pointed between the nucleophile and the leaving group. Moreover, motions of the catalytic aspartate toward the histidine occur concertedly with proton abstraction by the catalytic histidine and help stabilize the transition state, thus partially explaining how serine hydrolases enable poor nucleophiles to attack the substrate carbonyl carbon. Rate calculations indicate that the second step (deacylation) is rate-determining, with a calculated rate constant of 66 s –1 . Altogether, these results reveal the pivotal role of active-site dynamics in the catalytic mechanism of AnFaeA, which is likely similar in other serine hydrolases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1980

A solution of the variational equations for elliptic orbits in rotating coordinates

For elliptic reference orbits, formulas are given for the perturbation state transition matrix of the two-body problem. The formulas relate perturbations expressed in a local vertical rotating coordinate system and are valid for motion in the linear neighborhood of reference orbits with e in the range of 0 to 1. The elements of the state transition matrix are expressed in terms of natural parameters (horizontal and radial velocity, radius, eccentricity, true anomaly, etc.) at the initial and final points. In addition to the general form, a simplified version, valid for small eccentricity orbits, is given.

Jones, J. B.↗
NASA NTRS2006

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'prop-burn-prop' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons are made with three methods using modified linearized covariance propagation. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. The second method does not model the maneuver but includes a process noise matrix to account for the uncertainty in its magnitude. The third method, which is essentially a hybrid of the first two, includes the nominal portion of the maneuver via the state transition matrix and uses a process noise matrix to account for the magnitude uncertainty. The first method is unable to produce the final orbit covariance except in the case of zero maneuver uncertainty. The second method yields good accuracy for the final covariance matrix but fails to model the final orbital state accurately. Agreement between the simulated covariance data produced by this method and the Monte Carlo truth data fell within 0.5-2.5 percent over a range of maneuver sizes that span two orders of magnitude (0.1-20 m/s). The third method, which yields a combination of good accuracy in the computation of the final covariance matrix and correct accounting for the presence of the maneuver in the nominal orbit, is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, PC. However, applications for the two other methods exist and are briefly discussed. Although the process model ("prop-burn-prop") that was studied is very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗
DOE OSTI2022

A Conical Intersection Influences the Ground State Rearrangement of Fulvene to Benzene

The rearrangement of fulvene to benzene is believed to play an important role in the formation of soot during hydrocarbon combustion. Previous work has identified two possible mechanisms for the rearrangement-a unimolecular path and a hydrogen-atom-assisted, bimolecular path. Computational results to date have suggested that the unimolecular mechanism faces a barrier of about 74 kcal/mol, which makes it unable to compete with the bimolecular mechanism under typical combustion conditions. This computed barrier is about 10 kcal/mol higher than the experimental value, which is an unusually large discrepancy for modern electronic structure theory. In the present work, we have reinvestigated the unimolecular mechanism computationally, and we have found a second transition state that is approximately 10 kcal/mol lower in energy than the previously identified one and, therefore, in excellent agreement with the experimental value. The existence of two transition states for the same rearrangement arises because there is a conical intersection between the two lowest singlet states which occurs in the vicinity of the reaction coordinates. The two possible paths around the cone on the lower adiabatic surface give rise to the two distinct saddle points. The lower barrier for the unimolecular mechanism now makes it competitive with the bimolecular one, according to our calculations. In support of this conclusion, we have reanalyzed some previous experimental results on anisole pyrolysis, which leads to benzene as a significant product and have shown that the unimolecular and bimolecular mechanisms for fulvene ? benzene must be occurring competitively in that system. Finally, we have identified that similar conical intersections arise during the isomerizations of benzofulvene and isobenzofulvene to naphthalene.

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

Kinetic analyses for solid-state phase transition of metastable amorphous-AlO x (2.5 < x ≤ 3.0) nanostructures into crystalline alumina polymorphs

Solid-solid phase change materials (SS-PCMs) hold promise for energy storage/dissipation in batteries and energetic materials. Yet, phase change kinetics for SS-PCMs undergoing metastable to semi-stable/stable phase transformations remain relatively ill-studied because trapping metastable phases remain challenging. Recently, we demonstrated the kinetic entrapment and stabilization of a highly disordered and amorphous Al-oxide phase m-AlO x @C (x~2.5-3.0) via laser ablation synthesis in solution (LASiS). We report here, to our knowledge, the first chemical kinetics analysis for S-S phase transition of the m-AlO 3 @C nanocomposites (< 5–8 nm sizes) into semi-stable equilibrium alumina phases (θ/γ-Al 2 O 3 ) via disproportionation reaction, while releasing excess trapped gases. Our results indicate the atomic density of the AlO 3 structures to be ~5–10 times less than that of the final Al 2 O 3 phases, which led to the hypothesis of a volume shrinkage process during their phase transition. Temperature-dependent X-ray diffraction studies reveal the high-temperature phase transition for m-AlO 3 → θ/γ-Al 2 O 3 to follow contracting volume kinetics model, thereby validating our earlier hypothesis. Using the geometric volume contraction model, reaction kinetics analyses from Arrhenius plots reveal the activation energy barrier for the phase transition to be ~270±11 kJ/mol. This makes the activation energy barrier nearly identical to the oxidation of micron-sized Al particles.

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

Facet-Dependent Water Inhibition of Alkanol Dehydration on TiO2 via Distinct Water–Alkanol Complexes

Water is ubiquitous in biomass-derived feeds, yet its molecular impact on oxygen-elimination reactions remains poorly understood, particularly for catalysts exposing different facets. Here, we utilize well-defined TiO2 nanocrystals with dominant (101) and (001) facets to reveal a pronounced facet-dependent effect of water, where inhibition for dehydration of isopropanol (IPA) on the TiO2(001) surface is about four times more severe than TiO2(101). Through a combination of in situ solid state NMR, in situ infrared spectroscopy, kinetics studies, and theoretical calculations, we demonstrate that this disparity arises from the formation of distinct alkanol-water complex intermediates. On TiO2(001), IPA undergoes dissociative adsorption to form an isopropoxide-H2O complex that readily drives the surface into a complex-dominated regime. This pathway increases the activation barrier for C–H cleavage by 40 kJ mol-1 by inducing a disordered transition state. In contrast, TiO2(101) favors molecular IPA adsorption with weak hydrogen bonding to water, resulting in a smaller complex formation constant and a much smaller activation barrier increase (25 kJ mol-1). By quantitatively linking facet-dependent complex coverage to transition-state destabilization, this work moves beyond simple site-blocking models and provides a conceptual framework for designing catalysts that remain active in water-containing environments

Hu, Wenda↗
DOE OSTI2024

Demystifying group-4 polyolefin hydrogenolysis catalysis: Gaseous propane hydrogenolysis mechanism over the same catalysts

A kinetic/mechanistic investigation of gaseous propane hydrogenolysis over the single-site heterogeneous polyolefin depolymerization catalysts AlS/ZrNp 2 and AlS/HfNp 2 (AlS = sulfated alumina, Np = neopentyl), is use to probe intrinsic catalyst properties without the complexities introduced by time- and viscosity-dependent polymer medium effects. In a polymer-free automated plug-flow catalytic reactor, propane hydrogenolysis turnover frequencies approach 3,000 h −1 at 150 °C. Both catalysts exhibit approximately linear relationships between rate and [H 2 ] at substoichiometric [H 2 ] with rate law orders of 0.66 ± 0.09 and 0.48 ± 0.07 for Hf and Zr, respectively; at higher [H 2 ], the rates approach zero-order in [H 2 ]. Reaction orders in [C 3 H 8 ] and [catalyst] are essentially zero-order under all conditions, with the former implying rapid, irreversible alkane binding/activation. This rate law, activation parameter, and DFT energy span analysis support a scenario in which [H 2 ] is pivotal in one of two plausible and competing rate-determining transition states—bimolecular metal-alkyl bond hydrogenolysis vs. unimolecular β-alkyl elimination. The Zr and Hf catalyst activation parameters, ΔH ‡ = 16.8 ± 0.2 kcal mol −1 and 18.2 ± 0.6 kcal mol −1 , respectively, track the relative turnover frequencies, while ΔS ‡ = −19.1 ± 0.8 and −16.7 ± 1.4 cal mol −1 K −1 , respectively, imply highly organized transition states. These catalysts maintain activity up to 200 °C, while time-on-stream data indicate multiday activities with an extrapolated turnover number ~92,000 at 150 °C for the Zr catalyst. This methodology is attractive for depolymerization catalyst discovery and process optimization.

03 NATURAL GAS↗
DOE OSTI2026

Transient hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation. II. Isomer-resolved unimolecular dynamics

Transient carbon-centered hydroperoxyalkyl intermediates (•QOOH) in isopentane oxidation are characterized by their time- and energy-resolved unimolecular dissociation dynamics to hydroxyl (OH) and cyclic ether products. Two distinct •QOOH isomers are examined with radical sites at a primary carbon of one of the methyl groups (β-Me) or a secondary carbon (β-Et) of the ethyl group. Energy-dependent unimolecular rates are obtained from the time-dependent appearance of OH products for the two isomers and compared with statistical microcanonical rates computed using RRKM theory, including heavy-atom tunneling, based on high-level theoretical calculations. A benchmark-corrected approach is utilized to compute high-accuracy stationary-point energies, most importantly, transition-state barriers, for the •QOOH_Me and •QOOH_Et isomers in isopentane oxidation, building on higher-level reference calculations for the oxidation of ethane (C 2 H 5 O 2 ) and propane (C 3 H 7 O 2 ), respectively. Here, the measured rates are compared with RRKM calculations incorporating the benchmark-corrected transition-state parameters, a vibrationally adiabatic multidimensional hindered-rotor treatment of key torsions, and quantum tunneling. Agreement between experiment and theory validates the statistical description and shows faster decay for •QOOH_Et due to its lower barrier. Both β-QOOH isomers decay almost exclusively to OH + cyclic ether products under the conditions studied.

Oxidation processes↗
DOE OSTI2026

Nuclear–Electronic Orbital General Rate Theory: Predicting Hydrogen Kinetic Isotope Effects in the Deep Tunneling Regime

Hydrogen transfer is a critical component of many chemical and biological processes. The ratio of rate constants for hydrogen and deuterium transfer defines the H/D kinetic isotope effect (KIE), which is a powerful tool for elucidating hydrogen transfer mechanisms. Interpretation of experimental H/D KIEs relies on accurate and affordable computational methods. However, due to their light mass, hydrogen and deuterium can undergo tunneling, which is challenging to describe in multidimensional molecular systems. Herein, we introduce the nuclear–electronic orbital general rate theory (NEO-GRT), which enables the efficient prediction of H/D KIEs based on full-dimensional molecular quantum chemistry calculations. The NEO-GRT approach describes the hydrogen transfer rate constant with a general expression that spans the vibrationally adiabatic and nonadiabatic hydrogen tunneling regimes. The input quantities are computed using NEO density functional theory, which treats the transferring hydrogen or deuterium nucleus quantum mechanically on the same level as the electrons. We investigate two intramolecular proton transfer reactions in organic molecules at temperatures down to 50 K to evaluate the performance of NEO-GRT by comparison to transition state theory and ring-polymer instanton theory. The KIEs computed with NEO-GRT agree with those calculated using ring-polymer instanton theory for the full-dimensional molecular systems at the same level of electronic structure theory. This agreement indicates that NEO-GRT captures the deep hydrogen tunneling effects, in contrast to transition state theory, which neglects such effects. Given its relatively low computational cost, NEO-GRT is a promising approach for predicting H/D KIEs in large organic and organometallic systems.

Hydrogen↗
DOE OSTI2021

Role of Hydrogen-bonded Bimolecular Formic Acid–Formate Complexes for Formic Acid Decomposition on Copper: A Combined First-Principles and Microkinetic Modeling Study

Hydrogen bonding interactions alter the nanoscale reaction mechanisms of many chemistries. Yet, it remains unclear how they affect heterogeneously catalyzed decomposition of formic acid (FA), a reaction of intense interest since FA is a promising hydrogen carrier. In this work, we elucidate how hydrogen bonding affects the reaction mechanisms for FA decomposition on Cu(111) by combining first-principles density functional theory calculations to calculate reaction energetics, Latin-hypercube sampling to elucidate stable high-coverage adsorbate configurations, and coverage self-consistent mean-field microkinetic models to predict reaction kinetics. We demonstrate that hydrogen-bonded complexes of FA with formate (bimolecular FA–HCOO complexes) can play a dominant role in FA decomposition. Specifically, our first-principles calculations show that hydrogen bonding of FA with HCOO may stabilize the crucial monodentate HCOO intermediate and the transition states for HCOO decomposition, especially at low coverages. We predict that, depending on the reaction conditions, 40–80% of the reaction flux goes through pathways involving the bimolecular FA–HCOO complexes. Additionally, the active site for FA decomposition on Cu(111) involves a high coverage (~0.4 monolayers (ML)) of these complexes, which unexpectedly stabilize intermediates and transition states via van der Waals interactions. Our work provides molecular insights consistent with previous experimental observations on supported Cu/Al 2 O 3 catalysts. This paves the way toward the development of novel catalysts for FA decomposition as well as for other industrially important chemistries with intermediates capable of hydrogen bonding, such as ammonia electrooxidation and CO 2 hydrogenation.

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

Machine Learning Accelerated First-Principles Study of the Hydrodeoxygenation of Propanoic Acid

The complex reaction network of catalytic biomass conversions often involves hundreds of surface intermediates and thousands of reaction steps, greatly hindering the rational design of metal catalysts for these conversions. Here, we present a framework of machine learning (ML)-accelerated first-principles studies for the hydrodeoxygenation (HDO) of propanoic acid over transition metal surfaces. The microkinetic model (MKM) is initially parametrized by ML-predicted energies and iteratively improved by identifying the rate-determining species and steps (RDS), computing their energies by density functional theory (DFT), and reparameterizing the MKM until all the RDS are computed by DFT. The Gaussian process (GP) model performs significantly better than the linear ridge regression model for predicting both the adsorption free energies and transition state free energies. Parameterized with energies from the GP model, only 5–20% of the full reaction network has to be computed by DFT for the MKM to possess DFT-level accuracy for the TOF and dominant reaction pathway. While the linear ridge regression model performs worse than the GP model, its performance is greatly improved when only transition states are predicted by the regression model and adsorption energies are computed by DFT. Overall, we find that a high accuracy in adsorption free energies is more important for a reliable MKM than a high accuracy in TS free energies. Lastly, based on the GP model with GOH and GCHCHCO as catalyst descriptors, we build two-dimensional volcano plots in activity and selectivity that can help design promising alloy catalysts for HDO reactions of organic acids.

adsorption↗
NASA NTRS1993

Near-Threshold Rovibrational Excitation of HF by Electron Impact

The pronounced threshold peaks found in cross sections for electron-impact excitation of the vibrational levels of HF have been studied using an R-matrix method, in which electron exchange, correlation effects and polarization are included. Vibration is modelled inside the R-matrix sphere using the energy-modified adiabatic phase matrix method (EMAP). Rotation of the target molecule is treated outside the R-matrix sphere, using the frame transformation method. The external rovibronic close-coupling equations are solved by analytic methods. Analysis of the multichannel threshold structures indicates that the peaks found for rovibrational excitation of nu = 1 can be interpreted as virtual states while those for nu = 2 correspond to nuclear-excited Feshbach resonances. These two different results correspond simply to a change of sign of a small parameter that is determined by the threshold analysis. Whereas the overall shapes of the excitation peaks agree with earlier theoretical results, we have found detailed rotational fine structures for scattering energies on the order of the rotational constant above and below the opening of the vibrational channels. The partial and total cross sections for rovibrational state-to-state transitions in the energy region of the peaks are analysed. Characteristic cusp structures and excitation peaks are found at the rovibrational thresholds for transitions with small angular momentum transfer (Delta)j =0, 1 and 2. Differential cross sections are similar to recent measurements. These experimental data do not show the relatively sharp rise at the vibrational thresholds indicated by theory. Theory and experiment are in good agreement for the shape and magnitude of the angular dependence of cross sections for rovibrational state-to-state transitions, especially at larger energies.

Thummel, Helmar T.↗
DOE OSTI2026

Selective Chemical Looping Combustion of Terminal Alkynes in Mixtures with Alkenes

The selective combustion of terminal alkynes in mixtures with alkenes is demonstrated during anaerobic reduction half-cycles on bulk bismuth oxide (Bi 2 O 3 ) as an approach to remove alkynes, which act as inhibitors in olefin polymerization. Bi 2 O 3 combusts phenylacetylene in styrene, 3-methylphenylacetylene in 3-methylstyrene, propyne in propylene, 1-hexyne in 1-hexene, and 1-octyne in 1-octene, with alkyne combustion selectivities exceeding 96%. Near unity reaction orders for hydrocarbon consumption during reduction half-cycles are consistent with combustion pathways initiated by rate-determining initial C–H activation, which drive selective alkyne combustion through intrinsic differences in the first-order rate constants for alkyne and alkene combustion rather than preferential adsorption of alkynes on Bi 2 O 3 surfaces. Computational assessments of initial C–H activation pathways for alkynes and alkenes on (010) α-Bi 2 O 3 surfaces using density functional theory illustrate that heterolytic transition states which form proton-carbanion pairs on Bi–O sites kinetically favor the activation of alkynes rather than alkenes due to differences in C–H bond acidity, and the barrier for heterolytic C–H activation is dictated in part by the sum of the molecular deprotonation energy and the energy to bind an R – carbanion to a Bi site in its transition-state geometry. Finally, these heterolytic reactivity channels during selective chemical looping combustion present novel routes for purifying olefin gas streams containing alkyne impurities.

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

Probing the dynamics and bottleneck of the key atmospheric SO 2 oxidation reaction by the hydroxyl radical

SO 2 (Sulfur dioxide) is the major precursor to the production of sulfuric acid (H 2 SO 4 ), contributing to acid rain and atmospheric aerosols. Sulfuric acid formed from SO 2 generates light-reflecting sulfate aerosol particles in the atmosphere. This property has prompted recent geoengineering proposals to inject sulfuric acid or its precursors into the Earth’s atmosphere to increase the planetary albedo to counteract global warming. SO 2 oxidation in the atmosphere by the hydroxyl radical HO to form HOSO 2 is a key rate-limiting step in the mechanism for forming acid rain. However, the dynamics of the HO + SO 2 → HOSO 2 reaction and its slow rate in the atmosphere are poorly understood to date. Herein, we use photoelectron spectroscopy of cryogenically cooled HOSO 2 – anion to access the neutral HOSO 2 radical near the transition state of the HO + SO 2 reaction. Spectroscopic and dynamic calculations are conducted on the first ab initio-based full-dimensional potential energy surface to interpret the photoelectron spectra of HOSO 2 – and to probe the dynamics of the HO + SO 2 reaction. In addition to the finding of a unique pre-reaction complex (HO⋯SO 2 ) directly connected to the transition state, dynamic calculations reveal that the accessible phase space for the HO + SO 2 → HOSO 2 reaction is extremely narrow, forming a key reaction bottleneck and slowing the reaction rate in the atmosphere, despite the low reaction barrier. This study underlines the importance of understanding the full multidimensional potential energy surface to elucidate the dynamics of complex bimolecular reactions involving polyatomic reactants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2006

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'propagate-burn-propagate' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons between it and three methods using modified linearized covariance propagation are made. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. In the second method, the maneuver is not modeled but the uncertainty in its magnitude is accounted for by the inclusion of a process noise matrix. In the third method, which is essentially a hybrid of the first two, the nominal portion of the maneuver is included via the state transition matrix while a process noise matrix is used to account for the magnitude uncertainty. Since this method also correctly accounts for the presence of the maneuver in the nominal orbit, it is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, Pc. However, applications for the two other methods exist and are briefly discussed. Despite the fact that the process model ('propagate-burn-propagate') that was studied was very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗
NASA NTRS1984

Quantum chemical studies of a model for peptide bond formation. 3. Role of magnesium cation in formation of amide and water from ammonia and glycine

The SN2 reaction between glycine and ammonia molecules with magnesium cation Mg2+ as a catalyst has been studied as a model reaction for Mg(2+)-catalyzed peptide bond formation using the ab initio Hartree-Fock molecular orbital method. As in previous studies of the uncatalyzed and amine-catalyzed reactions between glycine and ammonia, two reaction mechanisms have been examined, i.e., a two-step and a concerted reaction. The stationary points of each reaction including intermediate and transition states have been identified and free energies calculated for all geometry-optimized reaction species to determine the thermodynamics and kinetics of each reaction. Substantial decreases in free energies of activation were found for both reaction mechanisms in the Mg(2+)-catalyzed amide bond formation compared with those in the uncatalyzed and amine-catalyzed amide bond formation. The catalytic effect of the Mg2+ cation is to stabilize both the transition states and intermediate, and it is attributed to the neutralization of the developing negative charge on the electrophile and formation of a conformationally flexible nonplanar five-membered chelate ring structure.

Magnesium/chemistry↗