Search NASA⌕ Search

SEARCH · Search NASA

Results for “rate constants”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Electrochemistry at Back-Gated, Ultrathin ZnO Electrodes: Field-Effect Modulation of Heterogeneous Electron Transfer Rate Constants by 30× with Enhanced Gate Capacitance

We report steady-state voltammetry of outer-sphere redox species at back-gated ultrathin ZnO working electrodes in order to determine electron transfer rate constants k ET as a function of independently-controlled gate bias, V G . We demonstrate that k ET can be modulated as much as 30-fold by application of V G ≤ 8 V. Key to this demonstration was integrating the ultrathin (5 nm) ZnO on a high dielectric constant (k) insulator, HfO 2 (30 nm), which was grown on a Pd metal gate. The high-k HfO 2 dramatically decreased the required V G values and increased the gate-induced charge in ZnO compared to previous studies. Importantly, the enhanced gating power of the Pd/HfO 2 /ZnO stack meant it was possible to observe a non-monotonic dependence of k ET on V G , which reflects the inherent density of redox acceptor states in solution. Furthermore, this work adds to the growing body of literature demonstrating that electrochemical kinetics (i.e., rate constants and overpotentials) at ultrathin working electrodes can be tuned by V G , independent of the conventional electrochemical working electrode potential.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Re-evaluation of rate constants for the reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M)

Rate constants for the dissociation/recombination reaction N 2 H 4 (+ M) $^→_←$ NH 2 + NH 2 (+ M) are determined by a combination of quantum-chemical calculations and statistical unimolecular rate theory. Between 1100 and 2500 K, limiting low-pressure rate constants for hydrazine dissociation in the bath gas Ar of k 0 = [Ar] 6.1 × 10 20 (T/1000 K) –7.3 exp (-34490 K/T) cm 3 mol –1 s –1 , limiting high-pressure rate constants of k ∞ = 7.6 × 10 16 (T/1000 K) –1.0 exp(-33600 K/T) s –1 , and center broadening factors of the falloff curves (between 1100 and 1600 K in Ar) of F cent = 0.71 exp(- T/1460 K) + 0.29 exp(- T/21 K) + exp(-13400 K/T) were calculated. Using equilibrium constants K C = 1.7 × 10 3 (T/1000 K) –1.5 exp(-33,460 K/T) mol cm –3 , between 300 and 600 K limiting low-pressure rate constants for the reverse recombination of NH 2 radicals in the bath gas N 2 of k rec,0 = [N 2 ] 4.4 × 10 20 (T/400 K) –6.9 exp (-1630 K/T) cm 6 mol –2 s –1 , limiting high-pressure rate constants of k rec,∞ = 4.4 × 10 13 (T/1000 K) 0.44 exp(-140 K/T) cm 3 mol -1 s –1 , and center broadening factors (between 300 and 600 K in N 2 ) of F cent = exp(- T/1130 K) + exp (-10340 K/T) were obtained. A comparison with experimental results for hydrazine dissociation from the literature suggests incomplete falloff extrapolations toward k ∞ and experimental problems in the determination of k 0 at temperatures above about 1600 K. Finally, implications of the present re-evaluated rate constants for the modeling of high temperature ammonia oxidation kinetics are discussed, showing an only small influence of their precise values on the overall properties of the process.

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

Dissociation-induced depletion of high-energy reactant molecules as a mechanism for pressure-dependent rate constants for bimolecular reactions

In 1922, Lindemann proposed the now-well-known mechanism for pressure-dependent rate constants for unimolecular reactions: reactant molecules with sufficiently high energies dissociate more quickly than collisions can reestablish the Boltzmann distribution of the internal energies of the molecule during its dissociation at low pressures – yielding pressure-dependent rate constants for unimolecular reactions due to the preferential depletion of the high energy states capable of dissociation. In the last century, incredible progress has been made in achieving a far greater understanding of and quantitative predictions for unimolecular and association reactions. In the modern era, pressure-dependent phenomenological rate constants are now nearly universally used to describe the rates of unimolecular and associative reactions in phenomenological kinetic modeling. However, there is a second, more indirect, implication of Lindemann’s mechanism that relates to how these dissociation-induced non-equilibrium distributions impact bimolecular reactions, including non-associative bimolecular reactions – which are generally not considered to have pressure-dependent rate constants. Yet, as we show herein, the same high energy states depleted due to dissociation would otherwise react most rapidly in high-activation-energy bimolecular reactions – yielding a mechanism for pressure-dependent rate constants for bimolecular reactions (including non-associative reactions). Here, we present results from a case study for CH 2 O dissociation, isomerization, and bimolecular reaction with O2 to explore this question. Results from our master equation calculations indicate that the effect of dissociation-induced non-equilibrium distributions on bimolecular reactions can be substantial – even when chemical timescales are well separated from internal energy relaxational timescales (i.e. when the traditional rate constant description would be thought to apply). This effect is found to be more pronounced – and more complex – for bimolecular reactions involving molecular entities whose chemical timescales are merged with the internal energy relaxational timescales. Lastly, we present some ideas for discussion regarding what should be considered as “chemical species” in phenomenological kinetic models.

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

Charge transfer rate constants for the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran: The spin-boson model vs the linearized semiclassical approximation

Charge transfer rate constants were calculated for the carotenoid-porphyrin-C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran. The calculation was based on mapping the all-atom anharmonic Hamiltonian of this system onto the spin-boson Hamiltonian. The mapping was based on discretizing the spectral density from the time correlation function of the donor–acceptor potential energy gap, as obtained from all-atom molecular dynamics simulations. Different spin-boson Hamiltonians were constructed for each of the possible transitions between the three excited electronic states in two different triad conformations. The rate constants of three possible transitions were calculated via the quantum-mechanically exact Fermi’s golden rule (FGR), as well as a progression of more approximate expressions that lead to the classical Marcus expression. The advantage of the spin-boson approach is that once the mapping is established, the quantum-mechanically exact FGR and the hierarchy of approximations are known in closed form. The classical Marcus charge transfer rate constants obtained with the spin-boson Hamiltonians were found to reproduce those obtained from all-atom simulations with the linearized semiclassical approximation, thereby confirming the equivalence of the two approaches for this system. Furthermore, within the spin-boson Hamiltonian, we also found that the quantum-mechanically exact FGR rate constants were significantly enhanced compared to the classical Marcus theory rate constants for two out of three transitions in one of the two conformations under consideration. The results confirm that mapping to the spin-boson model can yield accurate predictions for charge transfer rate constants in a system as complex as CPC 60 dissolved in tetrahydrofuran.

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

Inefficient intramolecular vibrational energy redistribution for the H + HO2 reaction and negative internal energy dependence for its rate constant

Quasiclassical trajectories (QCT) and newly constructed global potential energy surfaces are used to compute thermal and nonthermal rate constants for the H + HO 2 reaction. The thermal QCTs rate constants are up to 50% smaller than transition state theory (TST) rate constants based on the same level of electronic structure theory. This reduction is demonstrated to result from inefficient intramolecular vibrational energy redistribution (IVR) in the transient H 2 O 2 well, with a significant fraction of trajectories that reach the H 2 O 2 well promptly dissociating back to reactants instead of via the heavily statistically favored 2OH channel. The nonstatistical reduction factor, κ IVR , that quantifies this effect is shown to increase in importance with temperature, with κ IVR = 0.81 at 300 K and 0.47 at 2500 K. Finally, we show that inefficient IVR causes H + HO 2 rate constants mediated by H 2 O 2 to depend inversely on the initial vibrational excitation of HO 2 .

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

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

How to calculate the rate constants for nonradiative transitions between the M S components of spin multiplets?

Predicting the rates of spin-dependent processes characterized by nonradiative transitions between electronic states with different spin multiplicities is important for understanding the mechanisms of many photochemical and catalytic reactions. To calculate these rates, it is necessary to define the spin state representation and the couplings between these states that drives the interstate transitions. In this work, we describe three different approaches to calculating the spin-orbit coupling (SOC), transition probabilities, and rate constants between the MS components of the electronic states with different spin multiplicities. We implemented these approaches in our nonadiabatic statistical theory (NAST) software package, which predicts the transition probabilities and rate constants of spin-dependent processes using information obtained from electronic structure calculations. Here, we discuss the advantage and drawbacks of each approach and, as an example, calculate the rate constants for transitions between the spin states in the active site model of the protein rubredoxin.

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

Shock tube/laser absorption measurement of the rate constant of the reaction: H 2 O 2 + CO 2 2OH + CO 2

We address the role of the linear mixing rule in the kinetics of the H 2 O 2 decomposition system by reporting the rate constant for H 2 O 2 + M = 2OH + M (M = Ar and CO 2 ) in the temperature range of 1087-1234 K at low pressures in a mixture of 20% CO 2 in Argon. The reaction rate constant was inferred from H 2 O concentrations monitored by using a laser-absorption spectroscopy-based water diagnostic. To the best of our knowledge, this is the first measurement of the rate constant of this reaction in a mixture to be reported in literature. A significant discrepancy was found between the rate constants derived using the traditional linear mixing rule and the reduced pressure linear mixing rule. Furthermore, this discrepancy can have serious implications on the predictive accuracy of these kinetic models, especially under conditions relevant to the operation of supercritical CO 2 (sCO 2 ) power cycles that rely on oxy-fuel combustion in a working fluid comprised almost entirely of CO 2 .

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2021

CTRAMER: An open-source software package for correlating interfacial charge transfer rate constants with donor/acceptor geometries in organic photovoltaic materials

In this paper, we present CTRAMER (Charge-Transfer RAtes from Molecular dynamics, Electronic structure, and Rate theory)—an opensource software package for calculating interfacial charge-transfer (CT) rate constants in organic photovoltaic (OPV) materials based on ab initio calculations and molecular dynamics simulations. The software is based on identifying representative donor/acceptor geometries within interfacial structures obtained from molecular dynamics simulation of donor/acceptor blends and calculating the corresponding Fermi’s golden rule CT rate constants within the framework of the linearized-semiclassical approximation. While the methods used are well established, the integration of these state-of-the-art tools originating from different disciplines to study photoinduced CT processes with explicit treatment of the environment, in our opinion, makes this package unique and innovative. The software also provides tools for investigating other observables of interest. After outlining the features and implementation details, the usage and performance of the software are demonstrated with results from an example OPV system.

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

Neutral gas pressure dependence of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations

In this computational study, we describe a self-consistent trajectory simulation approach to capture the effect of neutral gas pressure on ion–ion mutual neutralization (MN) reactions. The electron transfer probability estimated using Landau–Zener (LZ) transition state theory is incorporated into classical trajectory simulations to elicit predictions of MN cross sections in vacuum and rate constants at finite neutral gas pressures. Electronic structure calculations with multireference configuration interaction and large correlation consistent basis sets are used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of background gas on ion transport. For H + – H – and Li + – H(D) – , our approach quantitatively agrees with measured speed-dependent cross sections for up to ~10 5 m/s. For the ion pair Ne + – Cl – , our predictions of the MN rate constant at ~1 Torr are a factor of ~2 to 3 higher than the experimentally measured value. Similarly, for Xe + – F – in the pressure range of ~20 000–80 000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. Here, the paradigm of using trajectory simulations to self-consistently capture the effect of gas pressure on MN reactions advanced here provides avenues for the inclusion of additional nonclassical effects in future work.

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

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2021

Tutorial: Langevin Dynamics methods for aerosol particle trajectory simulations and collision rate constant modeling

The Langevin Dynamics (LD) method (also known in the literature as Brownian Dynamics) is routinely used to simulate aerosol particle trajectories for transport rate constant calculations as well as to understand aerosol particle transport in internal and external fluid flows. This tutorial intends to explain the methodological details of setting up a LD simulation of a population of aerosol particles and to deduce rate constants from an ensemble of classical trajectories. We discuss the applicability and limitations of the translational Langevin equation to model the combined stochastic and deterministic motion of particles in fields of force or fluid flow. The drag force and stochastic “diffusion” force terms that appear in the Langevin equation are discussed elaborately, along with a summary of common forces relevant to aerosol systems (electrostatic, gravity, van der Waals, …); a commonly used first order and a fourth order Runge-Kutta time stepping schemes for linear stochastic ordinary differential equations are presented. A MATLAB® implementation of a LD code for simulating particle settling under gravity using the first order scheme is included for illustration. Scaling analysis of aerosol transport processes and the selection of timestep and domain size for trajectory simulations are demonstrated through two specific aerosol processes: particle diffusion charging and coagulation. Fortran® implementations of the first order and fourth order time-stepping schemes are included for simulating the 3D motion of a particle in a periodic domain. Lastly, potential applications and caveats to the usage of LD are included as a summary.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Measurements of Radical Reactivity with an Imine, (CF 3 ) 2 CNH: Rate Constants for Chlorine Atoms and Hydroxyl Radicals and the Global Warming Potential

The rate constant kOH for the reaction of 1,1,1,3,3,3-hexafluoroprop-2-imine with OH radicals was measured relative to two reference compounds, CH 3 F and CH 3 CHF 2 , to be k OH = (4.2 ± 1.1) × 10 −14 cm 3 molecule −1 s −1 at 295 K. This implies an atmospheric lifetime with respect to consumption by OH of 0.75 years. Reaction with Cl atoms yielded k Cl = (7.9 ± 1.7) × 10 −16 cm 3 molecule −1 s −1 at 295 K, and reaction with O 3 has an upper limit of k O3 < 4 × 10 −23 cm 3 molecule −1 s −1 , so that the atmospheric consumption by Cl and O 3 is negligibly slow. Absolute infrared cross sections of the imine yield a radiative efficiency of 0.34 W m −2 ppb −1 , which is corrected to 0.23 W m −2 ppb −1 for the effects of atmospheric lifetime. The imine’s corresponding 100-year global warming potential is 64 ± 19. This value is an upper limit, given that heterogenous atmospheric removal paths, such as hydrolysis in water droplets, are not included.

climate metric↗
DOE OSTI2021

Correlations of the Steam Oxidation Rate Constant of BWR Alloy Zircaloy-2 at 800–1400 °C

Steam oxidation experiments were conducted at 800–1400 °C with boiling water reactor alloy Zircaloy-2 strip specimens. Sample weight gain measurements were performed on the oxidized specimens before and after the test and were compared to oxygen pickup calculations using the Cathcart–Pawel correlation. The results showed that Zircaloy-2 follows the parabolic law at temperatures above 1000 °C. At or below 1000 °C, the oxidation rate was very low when compared to Cathcart–Pawel correlation and can be represented by a cubic expression. Arrhenius expressions are given to describe the parabolic rate constants at temperatures above 1000 °C and cubic rate constants are provided for temperatures at or below 1000 °C. Here, the weight gains calculated by our Arrhenius correlations are in excellent agreement with the measured sample weight gains at all test temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Progress in modeling hydrogen assisted ammonia oxidation with new experiments and a further reconciliation of the NH 3 + OH rate constant

Hydrogen-assisted oxidation of ammonia in a premixed, laminar flow tubular reactor under reducing conditions was investigated experimentally and through chemical kinetic modeling. Due to its impact on the competition for OH among ammonia and hydrogen, the rate constant for NH 3 + OH (R1) was determined through state-of-the-art theoretical kinetics calculations, employing composite energies that include the effects of higher order electronic excitations on the electronic energies along the variational reaction path and treating the limitations in the kinetics posed by the passage through a hydrogen-bonded complex. The resulting rate constant was in close agreement with the recent experimental value from Zaczek et al. (2025), settling a long-term dispute about the high-temperature value of k 1 and confirming within 20% the value previously used in modeling. The chemical kinetic model, with no other changes, captured well measured concentrations of NH 3 , H 2 , NO, and N 2 O from flow reactor oxidation of NH 3 /H 2 at slightly reducing conditions over a range of temperature (900-1350 K) and NH 3 /H 2 ratios (0.5-2.0). Comparison of the present results with reported data from a non-premixed setup indicates that for laminar flow tubular reactors, the reactor configuration may have implications for the observed H 2 consumption due to the possibility of preferential oxidation during mixing.

Ab initio theory↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2022

Determining the kinetic rate constants of Fe 3 O 4 -to-Fe and FeO-to-Fe reduction by H 2

Steel production using coal accounts for ~ 8% of global carbon emissions. “Green Steel” is a new grand concept proposed recently to make steel from iron ores using renewable derived “Green Hydrogen” to achieve zero carbon emission. The kinetics and rate-limiting steps of iron ore reduction into iron with H 2 as a reducing agent is critically important to the success of this new technology. While reduction of Fe 2 O 3 into Fe by H 2 follows multiple steps, the past research on this topic mainly deals with the overall averaged kinetics, giving little information on the elemental and rate-limiting steps. In this work, we report a kinetic study specifically design to attain kinetic rate constants of one-step reduction of Fe 3 O 4 -to-Fe and FeO-to-Fe. Guided by thermodynamics, we show first how to create in situ the desirable starting oxide phases, i.e., Fe 3 O 4 and FeO, with precisely controlled the ratio of partial pressures of H 2 O and H 2 . We then show time-dependent raw H 2 O content data collected by a mass spectrometer and the processed reduction data to extract kinetic rate constants. We found that the kinetics of the two one-step reduction reactions follows nicely the Johnson-Mehl-Avrami (JMA) phase transformation model. The one-step reduction mechanisms and activation energy are also discussed.

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

Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. Further, the accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗