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 55 records · Page 3

Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories

Here, the accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH 3 ) molecule and the nucleophilic substitution (S N 2) reaction F – + CH 3 F → FCH 3 + F – . We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.

Iyer, Gopal R.↗

Restricted Partition Function Semiclassical Transition State Theory (RPF-SCTST): Applications to Reactions Involving Hydrogen Cyanide, Hydrogen Peroxide, and Formaldehyde Oxide

We demonstrate the efficiency of the numerical calculation of thermal semiclassical transition state theory (SCTST) rates across several representative chemical reactions using the restricted partition function (RPF)─viz a function that depends only on the imaginary action associated with a reaction. Here, we treat the potential energy surfaces (PESs) through fourth order expansions around the saddle point and integrate the Hamiltonian up to second order in employing vibrational perturbation theory. We apply this formalism to uni- and bimolecular reactions with different─though relatively high─barrier heights to uncover the influence of the barrier properties, minimal energies, and separability of the rotational motion on rate constants and tunneling corrections. Although all modes are coupled within the RPF, we found in our numerical examples that the rotational component in the absence of significant rotational distortions can be separated from the vibrational contributions. Moreover, a classical treatment of the rotational contribution is adequate over the temperature range from approximately 100 to 1000 K. We also found that the choice of DFT basis set can lead to variations in rate constants of up to an order of magnitude. Lastly, different schemes for counting eligible energy levels result in rate constants that differ by no more than a factor of 2, with the remaining discrepancies attributed to the treatment of near-convergent levels in perturbation theory.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Numerical Method for Computing the State Transition Matrix Using Poincare Integral Invariants

The Poincare integral invariants describe the volumes of sets in Hamiltonian phase space. We use these invariants to derive a new numerical procedure for obtaining the state transition matrix (STM), which can be applied to both conservative and nonconservative systems. The method is analogous to a finite difference approximation of the STM, where perturbed states are numerically propagated along with the reference trajectory. We discuss the mathematical similarities between this new STM and existing methods, show numerical results for orbital motion and uncertainty propagation, and discuss new insights afforded by the Hamiltonian properties of phase flow.

state transition matrix↗

Approximate Cartesian state transition matrix

Based on a Taylor series expansion, an easily-computed approximation to the Cartesian state transition matrix is presented for a general velocity-independent force field. Suitable for the short time intervals encountered in onboard navigation applications of the extended Kalman filter, it provides approximately five decimal digits of accuracy for earth orbiting spacecraft with update intervals of one minute, and better accuracy for shorter intervals.

Markley, F. L.↗

Using integrals of the state transition matrix for efficient transient-response computations

An iterative technique for the computation of lower-order state systems with poles at the origin (such as aircraft lateral dynamics) or with time-polynomial inputs is developed and demonstrated. The integrated input or output signals are computed using the integrals of the state transition matrix, thus reducing the system order and computation time. In an application to the lateral response of a DC-8 aircraft (Grunwald, 1981), this method required only 7824 multiplication or addition operations (at fourth order) to achieve the accuracy obtained in 15,012 operations by a conventional sixth-order computation.

Grunwald, A. J.↗

A Numerical Scheme for Ordinary Differential Equations Having Time Varying and Nonlinear Coefficients Based on the State Transition Matrix

A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.

Bartels, Robert E.↗

Modeling Magnetic Disk Wind State Transitions in Black Hole X-Ray Binaries

We analyze three prototypical black hole X-ray binaries, 4U 1630–472, GRO J1655–40, and H1743–322, in an effort to systematically understand the intrinsic state transition of the observed accretion disk winds between wind-on and wind-off states by utilizing state-of-the-art Chandra/HETGS archival data from multi-epoch observations. We apply our magnetically driven wind models in the context of magnetohydrodynamic(MHD)calculations to constrain(1)their global density slope(p),(2)their density (n17) at the foot point of the inner most launching radius, and(3)the abundances of heavier elements(AFe,S,Si). Incorporating the MHD winds into xstar photoionization calculations in a self-consistent manner, we create a library of synthetic absorption spectra given the observed X-ray continua. Our analysis clearly indicates a characteristic bimodal transition of multi-ion X-ray winds; i.e., the wind density gradient is found to steepen(fro mp∼1.2–1.4 to∼1.4–1.5)while its density normalization declines as the source transitions from the wind-onto the wind-off state. The model implies that the ionized wind remains physically present even in the wind-off state, despite its apparent absence in the observed spectra. Super solar abundances for heavier elements are also favored. Our global multi-ion wind models, taking into account soft X-ray ions as well as Fe K absorbers, show that the internal wind condition plays an important role in wind transitions besides photoionization changes. Simulated XRISM/Resolve and Athena/X-IFU spectra are presented to demonstrate a high fidelity of the multi-ion wind model for a better understanding of these powerful ionized winds in the coming decades

Keigo Fukumura↗

Survival of itinerant excitations and quantum spin state transitions in YbMgGaO 4 with chemical disorder

A recent focus of quantum spin liquid (QSL) studies is how disorder/randomness in a QSL candidate affects its true magnetic ground state. The ultimate question is whether the QSL survives disorder or the disorder leads to a “spin-liquid-like” state, such as the proposed random-singlet (RS) state. Since disorder is a standard feature of most QSL candidates, this question represents a major challenge for QSL candidates. YbMgGaO 4 , a triangular lattice antiferromagnet with effective spin-1/2 Yb 3+ ions, is an ideal system to address this question, since it shows no long-range magnetic ordering with Mg/Ga site disorder. Despite the intensive study, it remains unresolved as to whether YbMgGaO 4 is a QSL or in the RS state. Here, through ultralow-temperature thermal conductivity and magnetic torque measurements, plus specific heat and DC magnetization data, we observed a residual κ 0 /T term and series of quantum spin state transitions in the zero temperature limit for YbMgGaO 4 . These observations strongly suggest that a QSL state with itinerant excitations and quantum spin fluctuations survives disorder in YbMgGaO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Generator coordinate method for transition-state dynamics in nuclear fission

Since its beginnings, fission theory has assumed that low-energy induced fission takes place through transition-state channels at the barrier tops. Nevertheless, up to now there is no microscopic theory applicable to those conditions. We suggest that modern reaction theory is suitable for this purpose, and propose a methodology based on a configuration-interaction framework using the generator coordinate method (GCM). Simple reaction-theoretic models are constructed with the Gaussian overlap approximation to parametrize both the dynamics within the channels and their incoherent couplings to states outside the barrier. The physical characteristics of the channels examined here are their effective bandwidths and the quality of the coupling to compound-nucleus states as measured by the transmission factor T. We also investigate the spacing of GCM states with respect to their degree of overlap. We find that a rather coarse mesh provides an acceptable accuracy for estimating the bandwidths and transmission factors. The common numerical stability problem in using the GCM is avoided due to the choice of meshes and the finite bandwidths of the channels. Here, the bandwidths of the channels are largely controlled by the zero-point energy with respect to the collective coordinate in the GCM configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Rotational state-to-state transition rate coefficients for H 2 O + H 2 O collisions at nonequilibrium conditions

Aims.The goal is to develop a database of rate coefficients for rotational state-to-state transitions in H 2 O + H 2 O collisions that is suitable for the modeling of energy transfer in nonequilibrium conditions, in which the distribution of rotational states of H 2 O deviates from local thermodynamic equilibrium. Methods.A two-temperature model was employed that assumed that although there is no equilibrium between all possible degrees of freedom in the system, the translational and rotational degrees of freedom can be expected to achieve their own equilibria independently, and that they can be approximately characterized by Boltzmann distributions at two different temperatures,T kin andT rot . Results.Upon introducing our new parameterization of the collisional rates, taking into account their dependence on bothT kin andT rot , we find a change of up to 20% in the H 2 O rotational level populations for both ortho and para-H 2 O for the part of the cometary coma where the nonequilibrium regime occurs.

Astronomy & Astrophysics↗

Radio astronomical determination of ground state transition frequencies of CH

From a comparison of the spectra of CH and other molecules toward the continuum source Cas A and four dark clouds, the ground state transition frequencies of CH have been determined and are reported. The relative errors in these frequencies are about twice as small as those obtained earlier by radio astronomical methods for the two main lines of ground state OH.

Rydbeck, O. E. H.↗

A Time Integration Algorithm Based on the State Transition Matrix for Structures with Time Varying and Nonlinear Properties

A variable order method of integrating the structural dynamics equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. When the time variation of the system can be modeled exactly by a polynomial it produces nearly exact solutions for a wide range of time step sizes. Solutions of a model nonlinear dynamic response exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with solutions obtained by established methods.

Bartels, Robert E.↗

Quantum and anharmonic effects in non-adiabatic transition state theory

Quantitative descriptions of non-adiabatic transition rates at intermediate temperatures are challenging due to the simultaneous importance of quantum and anharmonic effects. Here, in this paper, the interplay between quantum effects—for motion across or along the seam of crossing—and anharmonicity in the seam potential is considered within the weak coupling limit. The well-known expression for quantized 1-D motion across the seam (i.e., tunneling) in the linear terms approximation is derived in the thermal domain using the Lagrangian formalism, which is then applied to the case when tunneling is distributed along the seam of crossing (treating motion along the seam classically). For high-frequency quantum modes, a vibrationally adiabatic (VA) approach is developed that introduces to the non-adiabatic rate constant a factor associated with high-frequency wavefunction overlap; this approach treats the high-frequency motion along the seam quantum mechanically. To test these methodologies, the reaction N 2 O ↔ N 2 + O( 3 P) was chosen. CCSD(T)-F12b/cc-pVTZ-F12 explorations of the 3 A'- 1 A' seam of N 2 O revealed that seam anharmonicity has a strong effect on the rate constant (a factor of ~20 at 2000 K). Several quantum effects were found to be significant at intermediate/lower temperatures, including the quantum N–N vibration that was coupled with seam anharmonicity using the VA approach. Finally, a 1-D approximation to non-adiabatic instanton theory is presented to estimate the validity limit of the linear terms model at low temperatures (⁠~ 250 K for N 2 O). We recommend that the assumptions built into many statistical theories for non-adiabatic reactions—harmonic behavior, classical motion, linear terms, and weak coupling—should be verified on a case-by-case basis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The CH(X 2 Π) + H 2 O reaction: two transition state kinetics

he reaction of ground state methylidyne (CH) with water vapor (H 2 O) is theoretically re-investigated using high-level coupled cluster computations in combination with semi-classical transition state theory (SCTST) and two-dimensional master equation simulations. Insertion of CH into a H–O bond of H 2 O over a submerged barrier via a well-skipping mechanism yielding solely H and CH 2 O is characterized. The reaction kinetics is effectively determined by the formation of a pre-reaction van der Waals complex (PRC, HC—OH 2 ) and its subsequent isomerization to activated CH 2 OH in competition with PRC re-dissociation. The tunneling effects are found to be minor, while variational effects in the PRC → CH 2 OH step are negligible. The calculated rate coefficient k ( T ) is nearly pressure-independent, but strongly depends on temperature with pronounced down-up behavior: a high value of 2 x 10 -10 cm 3 s -1 at 50 K, followed by a fairly steep decrease down to 8 × 10 -12 cm 3 s -1 at 900 K, but increasing again to 5 × 10 -11 cm 3 s -1 at 3500 K. Over the T-range of this work, k ( T ) can be expressed as: k ( T , P = 0) = 2.31 × 10 -11 ( T /300 K) -1.615 exp(-38.45/ T ) cm 3 s -1 for T = 50-400 K k ( T , P = 0) = 1.15 × 10 -12 ( T /300 K) 0.8637 exp(892.6/ T ) cm 3 s -1 for T = 400-1000 K k ( T , P = 0) = 4.57 × 10 -15 ( T /300 K) 3.375 exp(3477.4/ T ) cm 3 s -1 for T = 1000–3500 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Computation of a Continuous-thrust State Transition Matrix Incorporating Accurate Hardware and Ephemeris Models

A significant body of work exists showing that providing a nonlinear programming (NLP) solver with expressions for the problem constraint gradient substantially increases the speed of program execution and can also improve the robustness of convergence, especially for local optimizers. Calculation of these derivatives is often accomplished through the computation of spacecraft's state transition matrix (STM). If the two-body gravitational model is employed as is often done in the context of preliminary design, closed form expressions for these derivatives may be provided. If a high fidelity dynamics model, that might include perturbing forces such as the gravitational effect from multiple third bodies and solar radiation pressure is used then these STM's must be computed numerically. We present a method for the power hardward model and a full ephemeris model. An adaptive-step embedded eight order Dormand-Prince numerical integrator is discussed and a method for the computation of the time of flight derivatives in this framework is presented. The use of these numerically calculated derivatieves offer a substantial improvement over finite differencing in the context of a global optimizer. Specifically the inclusion of these STM's into the low thrust missiondesign tool chain in use at NASA Goddard Spaceflight Center allows for an increased preliminary mission design cadence.

Low Thrust↗

The protonation of N2O reexamined - A case study on the reliability of various electron correlation methods for minima and transition states

The protonation of N2O and the intramolecular proton transfer in N2OH(+) are studied using various basis sets and a variety of methods, including second-order many-body perturbation theory (MP2), singles and doubles coupled cluster (CCSD), the augmented coupled cluster (CCSD/T/), and complete active space self-consistent field (CASSCF) methods. For geometries, MP2 leads to serious errors even for HNNO(+); for the transition state, only CCSD/T/ produces a reliable geometry due to serious nondynamical correlation effects. The proton affinity at 298.15 K is estimated at 137.6 kcal/mol, in close agreement with recent experimental determinations of 137.3 +/- 1 kcal/mol.

Martin, J. M. L.↗

Magnetic Field Perturbations to a Soft X-ray-Activated Fe (II) Molecular Spin State Transition

The X-ray-induced spin crossover transition of an Fe (II) molecular thin film in the presence and absence of a magnetic field has been investigated. The thermal activation energy barrier in the soft X-ray activation of the spin crossover transition for [Fe{H2B(pz)2}2(bipy)] molecular thin films is reduced in the presence of an applied magnetic field, as measured through X-ray absorption spectroscopy at various temperatures. The influence of a 1.8 T magnetic field is sufficient to cause deviations from the expected exponential spin state transition behavior which is measured in the field free case. We find that orbital moment diminishes with increasing temperature, relative to the spin moment in the vicinity of room temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Anaerobic benzene oxidation in Geotalea daltonii involves activation by methylation and is regulated by the transition state regulator AbrB

ABSTRACT Benzene is a widespread groundwater contaminant that persists under anoxic conditions. The aim of this study was to more accurately investigate anaerobic microbial degradation pathways to predict benzene fate and transport. Preliminary genomic analysis of Geotalea daltonii strain FRC-32, isolated from contaminated groundwater, revealed the presence of putative aromatic-degrading genes. G. daltonii was subsequently shown to conserve energy for growth on benzene as the sole electron donor and fumarate or nitrate as the electron acceptor. The hbs gene, encoding for 3-hydroxybenzylsuccinate synthase (Hbs), a homolog of the radical-forming, toluene-activating benzylsuccinate synthase (Bss), was upregulated during benzene oxidation in G. daltonii , while the bss gene was upregulated during toluene oxidation. Addition of benzene to the G. daltonii whole-cell lysate resulted in toluene formation, indicating that methylation of benzene was occurring. Complementation of σ 54 - (deficient) E. coli transformed with the bss operon restored its ability to grow in the presence of toluene, revealing bss to be regulated by σ 54 . Binding sites for σ 70 and the transition state regulator AbrB were identified in the promoter region of the σ 54 -encoding gene rpoN, and binding was confirmed. Induced expression of abrB during benzene and toluene degradation caused G. daltonii cultures to transition to the death phase. Our results suggested that G. daltonii can anaerobically oxidize benzene by methylation, which is regulated by σ 54 and AbrB. Our findings further indicated that the benzene, toluene, and benzoate degradation pathways converge into a single metabolic pathway, representing a uniquely efficient approach to anaerobic aromatic degradation in G. daltonii . IMPORTANCE The contamination of anaerobic subsurface environments including groundwater with toxic aromatic hydrocarbons, specifically benzene, toluene, ethylbenzene, and xylene, has become a global issue. Subsurface groundwater is largely anoxic, and further study is needed to understand the natural attenuation of these compounds. This study elucidated a metabolic pathway utilized by the bacterium Geotalea daltonii capable of anaerobically degrading the recalcitrant molecule benzene using a unique activation mechanism involving methylation. The identification of aromatic-degrading genes and AbrB as a regulator of the anaerobic benzene and toluene degradation pathways provides insights into the mechanisms employed by G. daltonii to modulate metabolic pathways as necessary to thrive in anoxic contaminated groundwater. Our findings contribute to the understanding of novel anaerobic benzene degradation pathways that could potentially be harnessed to develop improved strategies for bioremediation of groundwater contaminants.

Bullows, James E.↗