Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Point Kinetics”

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 775 records · Page 43

Hydroxylation and Cation Segregation in (La 0.5 Sr 0.5 )FeO 3-$δ$ Electrodes

Simultaneously achieving high activity and stability is the primary challenge when engineering (electro)catalysts. Transition metal perovskite oxides are employed as air electrodes for solid-oxide fuel cells and electrolyzers. However, degradation of oxygen exchange kinetics at the solid–gas interface, often linked to alkaline-earth cation segregation and precipitation, limits widespread commercialization. Here, we systematically investigated the surface degradation mechanism induced by gas-phase impurities in La 0.5 Sr 0.5 )FeO 3-$δ$ (LSF55) thin-film electrodes by varying the concentration of H 2 O, SO 2 , and CO 2 . Degradation of the area-specific resistance in ambient and humidified synthetic air is significantly greater than in dry ambient and dry synthetic air, pointing to the importance of water vapor. Time-resolved, in situ ambient pressure X-ray photoelectron spectroscopy performed in O 2 showed that nonbulk Sr is present on the surface before the exposure to water vapor. Upon introduction of water vapor, neither additional Sr segregation nor precipitation driven by water vapor is a necessary condition for degradation. Rather, hydroxylation of the surface induces irreversible and significant degradation. At the same time, we show that Sr migration driven by water vapor is partially reversible. These fundamental insights can be used for the rational design of electrodes with improved catalytic stability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Delving into the depths of NGC 3783 with XRISM

We present our study of the X-Ray Imaging and Spectroscopy Mission (XRISM) observation of the Seyfert-1 galaxy NGC 3783. XRISM’s Resolve microcalorimeter has enabled, for the first time, a detailed characterization of the highly ionized outflows in this active galactic nucleus. Our analysis constrains their outflow and turbulent velocities, along with their ionization parameter (ξ) and column density (N H ). The high-resolution Resolve spectrum reveals a distinct series of Fe absorption lines between 6.4 and 7.8 keV, ranging from Fe XVIII to Fe XXVI. At lower energies (1.8−3.3 keV), absorption features from Si, S, and Ar are also detected. Our spectroscopy and photoionization modeling of the time-averaged Resolve spectrum uncovers six outflow components, five of which exhibit relatively narrow absorption lines with outflow velocities ranging from 560 to 1170 km s −1 . In addition, a broad absorption feature is detected, which is consistent with Fe XXVI outflowing at 14 300 km s −1 (0.05 c). The kinetic luminosity of this component is 0.8−3% of the bolometric luminosity. Our analysis of the Resolve spectrum shows that more highly ionized absorption lines are intrinsically broader than those of lower-ionization species, indicating that the turbulent velocity of the six outflow components (ranging from 0 to 3500 km s −1 ) increases with ξ. Furthermore, we find that the column density (N H ) of the outflows generally declines with the ionization parameter up to log ξ = 3.2 but rises beyond this point, suggesting a complex ionization structure. The absorption profile of the Fe XXV resonance line is intriguingly similar to UV absorption lines (Lyα and C IV) observed by the Hubble Space Telescope, from which we infer that the outflows are clumpy in nature. Our XRISM/Resolve results from lower- and higher-ionization regimes support a “hybrid wind” scenario in which the observed outflows have multiple origins and driving mechanisms. We explore various interpretations of our findings within active galactic nucleus wind models.

X-rays: galaxies↗

The mechanism driving a solid–solid phase transition in a biomacromolecular crystal

A solid-solid phase transition (SSPT) occurs between distinguishable crystalline forms. Because of its importance in application and theory in material science and condensed matter physics, SSPT has been studied most extensively in metallic alloys, inorganic salt or small organic molecular crystals, but much less so in biomacromolecular crystals. In general, the mechanism of SSPT at atomic and molecular levels is not well understood. Here, we describe the ordered molecular rearrangements in biomacromolecular crystals of the adenine riboswitch (riboA) aptamer using real-time serial crystallography and solution atomic force microscopy (AFM). The large, ligand-induced conformational changes drive the initial phase transition from the apo unit cell (AUC) to the trans unit cell 1 (TUC1). During this transition, coaxial stacking of P1 duplexes becomes the dominant packing interface, whereas P2-P2 interactions are almost completely disrupted, resulting in “floating” layers of molecules. The coupling points in TUC1 and their local conformational flexibility allow the molecules to reorganize to achieve the more densely packed and energetically favorable bound unit cell (BUC). Our study thus reveals the interplay between the conformational changes and the crystal phases—the underlying mechanism that drives the phase transition. Using polarized video microscopy (PVM) to monitor the SSPT in small crystals at high ligand concentration, we have identified the time window during which the major conformational changes take place, and simulated the in crystallo kinetics. Together, these results provide the spatiotemporal information necessary for informing time-resolved crystallography (TRX) experiments. Moreover, this study illustrates a practical approach to characterize SSPT in transparent crystals.

36 MATERIALS SCIENCE↗

Pressure-coupled vaporization and combustion responses of liquid-fuel droplets in high-pressure environments

The dynamic responses of liquid-fuel droplet vaporization and combustion to ambient pressure oscillations are examined. The analysis is based on the complete sets of conservation equations for both gas and liquid phases, and accommodates detailed treatments of finite-rate chemical kinetics and variable properties. With a full account of thermodynamic phase equilibrium at the droplet surface, the model enables a systematic examination of the effects of ambient flow conditions on the droplet behavior. The responses of hydrocarbon fuel droplets in both sub- and super-critical environments are investigated. Results indicate that the droplet gasification and burning mechanisms depend greatly on the ambient pressure. In particular, a rapid enlargement of the vaporization and combustion responses occurs when the droplet surface reaches its critical point, mainly due to the strong variations of latent heat of vaporization and thermophysical properties at the critical state.

Yang, Vigor↗

Iodine recombination in xenon solvent: Clusters in the gas to liquid-like state transition

Supercritical fluids (SCFs) have attracted significant attention as solvents for chemical reactions due to their unique properties, such as high diffusivity, low viscosity, and tunable solvation properties. These properties profoundly influence reaction kinetics and are often attributed to the formation of molecular clusters within SCFs. To study the effect of supercritical solvent on chemical reactivity and dynamics of reactions, one needs to understand the dynamics of clusters in supercritical fluid. Extensive experiments on the photodissociation and recombination of iodine in supercritical fluids served as a model system for understanding these effects. Experimental studies have been complemented by theoretical and computational investigations, which mostly employ Monte Carlo or empirical molecular dynamics simulations. However, computational studies using non-reactive force fields and ab initio approaches present challenges in capturing reactive processes at larger scales within supercritical fluids. Here, in this work, we developed the ReaxFF parameters by training against quantum mechanics data. ReaxFF reactive force field based molecular dynamics simulations were performed, studying the dynamics of a xenon solvent and cage effect at different thermodynamic conditions for the iodine recombination reaction. We show that the conditions near the critical point are the optimal conditions to study the cage effect. We show that the average lifetime of xenon clusters ranging between 5 and 11 ps is comparable to iodine geminate recombination. Our simulation results of iodine recombination in xenon solvent demonstrate the higher probability of iodine molecule formation in the presence of xenon clusters. Finally, we show that the supercritical condition exhibits the highest recombination rate for iodine atoms.

Cage effect↗

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

Demonstration of reconstruction-free static magnetic control of DIII-D plasma with deep reinforcement learning

This paper presents the development and experimental validation of a reinforcement learning (RL)-based magnetic controller on the DIII-D tokamak. The controller directly maps raw magnetic diagnostic signals to actuator commands, replacing the traditional isoflux control algorithm based on equilibrium reconstruction. Four RL controllers are trained using the Soft Actor–Critic algorithm with an asymmetric Actor–Critic architecture in the NSFsim simulator. All controllers are deployed in the DIII-D Plasma Control System and operated with a 4 kHz feedback loop. Two randomization strategies are evaluated during training: evolving kinetic profiles and fixed kinetic profiles within each episode. The latter approach is found to better capture experimental deviations in the current density profile and to provide overall improved control performance. Robust operation is demonstrated across heating power scans in both L- and H-mode plasmas, as well as during transient events such as L–H transitions and pellet injections. Control errors in plasma shape and radial position remained within 1.5–2.0 cm and 1 cm, respectively. A notable discrepancy was observed in the vertical X-point position, with errors of up to approximately 4 cm, attributed to the current density distribution mismatches between simulations and experiments.

DIII-D↗

Vanadyl Phosphates A x VOPO 4 (A = Li, Na, K) as Multielectron Cathodes for Alkali-Ion Batteries

Vanadyl phosphates comprise a class of multielectron cathode materials capable of cycling two Li + , about 1.66 Na + , and some K + ions per redox center. In this review, structures, thermodynamic stabilities, and ion diffusion kinetics of various A x VOPO 4 (A = Li, Na, K, NH 4 ) polymorphs are discussed. Both the experimental data and first-principle calculations indicate kinetic limitations for alkali metal ions cycling, especially between for 0 ≤ x ≤ 1, and metastability of phases with x > 1. This creates challenges for multiple-ion cycling, as the slow kinetics call for nanosized particles, which being metastable and reactive with organic electrolytes are prone to side reactions. Thus, various synthesis approaches, surface coating, and transition metal ion substitution strategies are discussed here as possible ways to stabilize A x VOPO 4 structures and improve alkali metal ion diffusion. Finally, the role of advanced characterization techniques, such as X-ray absorption spectroscopy, diffraction, pair distribution function analysis and 7 Li and 31 P NMR, in understanding the reaction mechanism from both structural and electronic points of view is emphasized.

25 ENERGY STORAGE↗

Prominence and Cool Loop Energetics Measured in the UV, EUV, and H-alpha

We investigate the thermal and dynamic properties of moving features in a prominence jet, a approx. 10(exp 5)K loop near a prominence channel, and an erupting prominence. In order to make measurements of the quickly moving features seen in prominences in the UV we use the SOHO/SUMER spectrograph to take a time series of exposures from a single pointing position, providing a measurement of spectral line properties as a function of time and position along the slit. The lines observed cover a broad range of temperatures from 80,000 - 1.6 million K. These measurements are combined with TRACE movies in transition region and coronal temperature bands and with ground based H-alpha data to obtain more complete information concerning prominence structure and motions. The resulting observations allow us perform DEM analysis and to calculate limits on densities, pressures, and the thermal and kinetic energies of the moving sources.

Kucera, Therese A.↗

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Supersonic shock wave/vortex interaction

Although shock wave/vortex interaction is a basic and important fluid dynamics problem, very little research has been conducted on this topic. Therefore, a detailed experimental study of the interaction between a supersonic streamwise turbulent vortex and a shock wave was carried out at the Penn State Gas Dynamics Laboratory. A vortex is produced by replaceable swirl vanes located upstream of the throat of various converging-diverging nozzles. The supersonic vortex is then injected into either a coflowing supersonic stream or ambient air. The structure of the isolated vortex is investigated in a supersonic wind tunnel using miniature, fast-response, five-hole and total temperature probes and in a free jet using laser Doppler velocimetry. The cases tested have unit Reynolds numbers in excess of 25 million per meter, axial Mach numbers ranging from 2.5 to 4.0, and peak tangential Mach numbers from 0 (i.e., a pure jet) to about 0.7. The results show that the typical supersonic wake-like vortex consists of a non-isentropic, rotational core, where the reduced circulation distribution is self similar, and an outer isentropic, irrotational region. The vortex core is also a region of significant turbulent fluctuations. Radial profiles of turbulent kinetic energy and axial-tangential Reynolds stress are presented. The interactions between the vortex and both oblique and normal shock waves are investigated using nonintrusive optical diagnostics (i.e. schlieren, planar laser scattering, and laser Doppler velocimetry). Of the various types, two Mach 2.5 overexpanded-nozzle Mach disc interactions are examined in detail. Below a certain vortex strength, a 'weak' interaction exists in which the normal shock is perturbed locally into an unsteady 'bubble' shock near the vortex axis, but vortex breakdown (i.e., a stagnation point) does not occur. For stronger vortices, a random unsteady 'strong' interaction results that causes vortex breakdown. The vortex core reforms downstream of the rear stagnation point, and the reduced circulation distribution once again becomes self-similar in this region. A-new model of this interaction is proposed. Finally, a curve defining the approximate limits of supersonic vortex breakdown is presented.

Settles, G. S.↗

Convective and radiative heating for vehicle return from the Moon and Mars

The aerothermal environment is examined for two vehicle forebodies near the peak heating points of lunar and martian return-to-earth trajectories at several nominal entry velocities. The first vehicle forebody is that of a 70 deg aerobrake for entry into earth orbit; the second, a capsule of Apollo configuration for direct entry into the earth's atmosphere. The configurations and trajectories are considered likely candidates for such missions. Two-temperature, thermochemical nonequilibrium models are used in the flow field analyses. In addition to Park's empirical model for dissociation under conditions of thermal nonequilibrium, the Gordiets kinetic model for the homonuclear dissociation of N2 and O2 is also considered. Temperature and emission profiles indicate nonequilibrium effects in a 2 to 5 cm post shock region. Substantial portions of the shock layer flow appear to be in equilibrium. The shock layer over an aerobrake for return from the moon exhibits the largest extent of nonequilibrium effects of all considered missions. Differences between the Gordiets and Parks kinetic model were generally very small for the lunar return aerobrake case, the greatest difference of 6.1 percent occurring in the radiative heating levels.

Greendyke, Robert B.↗

A Multi-Phase Modeling Framework Suitable for Dynamic Applications

Under dynamic loading conditions and the associated extreme conditions many metals will undergo phase transformations. The change in crystal structure associated with solid–solid phase transformations can significantly alter the subsequent mechanical response of the material. For the interpretation of experiments involving dynamic loading it is beneficial to have a modeling framework that captures key features of the material response while remaining relatively simple. We introduce a candidate framework and apply it to the metal tin to highlight a range of behaviors that are captured by the model. We also discuss potential extensions to capture additional behaviors that could be important for certain materials and loading scenarios. The model is useful for analysis of results from dynamic experiments and offers a point of departure for more complex model formulations.

36 MATERIALS SCIENCE↗

Dose, LET, time and strain dependence of radiation-induced 53BP1 foci in 15 mouse strains ex vivo and associations to in vivo radiation susceptibility

We present a comparative analysis on the repair of radiation-induced DNA damage ex vivo in 15 strains of mice, including 5 inbred reference strains and 10 collaborative-cross strains, of both sexes. Non-immortalized primary skin fibroblasts derived from 76 mice were subjected to both low- and high-LET radiation (0.1, 1 and 4 Gy of X rays; 1.1 and 3 particles/100μm2 of 350 MeV/n 40Ar and 600 MeV/n 56Fe). Automated image quantification of 53BP1 radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET. Similarly to what we had previously reported for immortalized human cell lines [1], we observed a saturation of RIF number with dose at 4h post-irradiation, with more RIF/Gy for lower LET (X rays and 40Ar) compared to 56Fe. However at later time points (24h and above), the trend was inverted with more RIF/Gy for higher LET. Our data suggest that multiple DSBs cluster into RIF: as the linear density of DSBs increases with LET, so does the probability of having more DSBs per RIF, which makes it more difficult for cells to fully resolve high-LET-induced RIF, explaining the hypersensitivity to high-LET radiation despite a low number of RIF. Taking into account the amount of clustering at a given dose and LET, but also the kinetics of DNA damage repair, we introduced a novel mathematical formalism to evaluate the number of remaining RIF over time. We showed that the newly introduced kinetic metrics can be used as surrogate biomarkers for in vivo radiation toxicity, with potential applications in radiotherapy and human space exploration. In particular, we observed an association between the repairable fraction of RIF measured in vitro and survival levels of immune cells collected from irradiated mice. Moreover, the speed of DNA damage repair correlated with spontaneous cancer incidence data collected from the Mouse Tumor Biology database, suggesting a relationship between the efficiency of DSB repair after irradiation and cancer risk. In addition to the efficacy of repair and persistent RIF levels, even the amount of spontaneous foci without irradiation was shown to be strain dependent, indicating that these phenotypes are at least partially driven by genetics, and supporting their potential as indicators of individual radiation sensitivity. [1] Neumaier, T., et al., PNAS, 2012 (8) 109:443

Radiation, DNA damage, repair kinetics↗

Multiplicity scaling of light nuclei production in relativistic heavy-ion collisions

Using the nucleon coalescence model based on kinetic freeze-out nucleons from the 3D MUSIC+UrQMD and the 2D VISHNU hybrid model with a crossover equation of state, we study the multiplicity dependence of deuteron ( d ) and triton ( t ) production from central to peripheral Au+Au collisions at s NN = 7.7, 14.5, 19.6, 27, 39, 62.4 and 200 GeV and Pb+Pb at s NN = 2.76 TeV, respectively. It is found that the ratio N t N p / N d 2 of the proton yield N p , deuteron yield N d and triton yield N t exhibits a scaling behavior in its multiplicity dependence, i.e., decreasing monotonically with increasing charged-particle multiplicity. A similar multiplicity scaling of this ratio is also found in the nucleon coalescence calculation based on kinetic freeze-out nucleons from a multiphase transport (AMPT) model. The scaling behavior of N t N p / N d 2 can be naturally explained by the interplay between the sizes of light nuclei and the nucleon emission source. We further argue that the multiplicity scaling of N t N p / N d 2 can be used to validate the production mechanism of light nuclei, and the collision energy dependence of this yield ratio can further serve as a baseline in the search for the QCD critical point in relativistic heavy-ion collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

A kinetic evaluation on NO 2 formation in the post-flame region of pressurized oxy-combustion process

Pressurized oxy-combustion is a promising technology that can significantly re-duce the energy penalty associated with first generation oxy-combustion for CO 2 capture in coal-fired power plants. However, higher pressure enhances the production of strong acid gases, including NO 2 and SO 3 , aggravating the corrosion threat during flue gas re-circulation. In the flame region, high temperature NO x exists mainly as NO, while conversion from NO to NO 2 happened in post-flame region. In this study, the conversion of NO → NO 2 has been kinetically evaluated under representative post-flame conditions of pressurized oxy-combustion after validating the mechanism (80 species and 464 reactions), which includes nitrogen and sulfur chemistry based on GRI-MECH 3.0. The effects of residence time, temperature, pressure, major species (O 2 /H 2 O), and minor or trace species (CO/SO x ) on NO 2 formation are studied. The calculation results show that when pressure is increased from 1 to 15 bar, NO 2 is increased from 1 to 60 ppm, and the acid dew point increases by over 80°C. Higher pressure and temperature greatly reduce the time required to reach equilibrium. With increasing pressure and decreasing temperature, O plays a much more important role than HO 2 in the oxidation of NO. A higher water vapor content accelerates NO 2 formation in all cases by providing more O and HO 2 radicals. The addition of CO or SO 2 also promotes the formation of NO 2 . The NO 2 formation in a pressurized oxy-combustion furnace can be over 10 times that of an atmospheric air-combustion furnace.

01 COAL, LIGNITE, AND PEAT↗

Experimental evaluation of heat transfer on a 1030:1 area ratio rocket nozzle

A 1030:1 carbon steel, heat-sink nozzle was tested. The test conditions included a nominal chamber pressure of 2413 kN/sq m and a mixture ratio range of 2.78 to 5.49. The propellants were gaseous oxygen and gaseous hydrogen. Outer wall temperature measurements were used to calculate the inner wall temperature and the heat flux and heat rate to the nozzle at specified axial locations. The experimental heat fluxes were compared to those predicted by the Two-Dimensional Kinetics (TDK) computer model analysis program. When laminar boundary layer flow was assumed in the analysis, the predicted values were within 15% of the experimental values for the area ratios of 20 to 975. However, when turbulent boundary layer conditions were assumed, the predicted values were approximately 120% higher than the experimental values. A study was performed to determine if the conditions within the nozzle could sustain a laminar boundary layer. Using the flow properties predicted by TDK, the momentum-thickness Reynolds number was calculated, and the point of transition to turbulent flow was predicted. The predicted transition point was within 0.5 inches of the nozzle throat. Calculations of the acceleration parameter were then made to determine if the flow conditions could produce relaminarization of the boundary layer. It was determined that if the boundary layer flow was inclined to transition to turbulent, the acceleration conditions within the nozzle would tend to suppress turbulence and keep the flow laminar-like.

Kacynski, Kenneth J.↗

Theoretical H + O 3 rate coefficients from ring polymer molecular dynamics on an accurate global potential energy surface: assessing experimental uncertainties

Thermal rate coefficients and kinetic isotope effects have been calculated for an important atmospheric reaction H/D + O 3 → OH/OD + O2 based on an accurate permutation invariant polynomial-neural network potential energy surface, using ring polymer molecular dynamics (RPMD), quasi-classical trajectory (QCT) and variational transition-state theory (VTST) with multidimensional tunneling. The RPMD approach yielded results that are generally in better agreement with experimental rate coefficients than the VTST and QCT ones, especially at low temperatures, attributable to its capacity to capture quantum effects such as tunneling and zero-point energy. Additionally, the theoretical results support one group of existing experiments over the other. In addition, rate coefficients for the D + O 3 → OD + O 2 reaction are also reported using the same methods, which will allow a stringent assessment of future experimental measurements, thus helping to reduce the uncertainty in the recommended rate coefficients of this reaction.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗