Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “chemical rate theory”

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 OSTI2021

Radiation-Enhanced Anion Transport in Hematite

The influence of radiation-induced (1 MeV energy H+ to ~0.1 dpa at 450°C), non-equilibrium point defect populations on mass transport is studied with an integrated campaign of experimental and theoretical methods. Using epitaxial thin films of hematite with embedded 18O tracer layers and nanoscale atom probe tomography measurements, it is shown that anion self-diffusion is enhanced by at least 2 orders of magnitude under irradiation compared to thermal diffusion alone. Complementary scanning transmission electron microscopy of vacuum annealed specimens, reveals associated microstructural changes in the oxide films, including local phase transformation to Fe3O4 and the development of nanoscale voids from vacancy coalescence. Point defect formation and migration energies were computed from density functional theory and applied within the context of chemical rate theory to analyze contributions from both interstitial and vacancy mechanisms to self-diffusion in thermal and irradiation conditions. Comparisons are made between calculated, literature and newly measured self-diffusion values, revealing good agreement on the magnitude of radiation-enhanced anion diffusion. Further, the model suggests a transition from vacancy to interstitialcy mechanisms at low temperatures and high oxygen activity, consistent with the varied activation energies reported from prior studies.

iron oxides, isotopic tracers, atom probe tomograp↗
DOE OSTI2021

Radiation Enhanced Anion Diffusion in Chromia

Cr 2 O 3 is a ubiquitous passivating film on alloys that is vital to their stability and performance in hostile environments. In this study, we directly observe and measure how anion transport through the oxide is affected by 400 keV Ar 2+ ion irradiation as a function of temperature (30–500 °C) and dose (0.33–2.00 dpa) in single-crystal Cr 2 O 3 films using embedded 18 O isotopic tracers and atom probe tomography. Diffusion coefficients are experimentally determined and compared to a chemical rate theory model. Experimental and theoretical results are broadly in agreement, showing anion diffusion increases by at least 4–5 orders of magnitude upon irradiation with a weak dependence on temperature and a stronger dose-driven sink-strength dependence. These results reveal that radiation could significantly reduce the protectiveness of Cr 2 O 3 films when passivity relies upon limited anion diffusivity, even at relatively modest levels of irradiation and low temperatures.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

On the origin of ground-state vacuum-field catalysis: Equilibrium consideration

Recent experiments suggest that vibrational strong coupling (VSC) may significantly modify ground-state chemical reactions and their rates even without external pumping. The intrinsic mechanism of this “vacuum-field catalysis” remains largely unclear. Generally, modifications of thermal reactions in the ground electronic states can be caused by equilibrium or non-equilibrium effects. The former are associated with modifications of the reactant equilibrium distribution as expressed by the transition state theory of chemical reaction rates, while the latter stem from the dynamics of reaching and leaving transition state configurations. In this work, we examine how VSC can affect chemical reactions rates in a cavity environment according to transition state theory. Our approach is to examine the effect of coupling to cavity mode(s) on the potential of mean force (PMF) associated with the reaction coordinate. Within the context of classical nuclei and classical photons and also assuming no charge overlap between molecules, we find that while the PMF can be affected by the cavity environment, this effect is negligible for the usual micron-length cavities used to examine VSC situations.

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

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↗
DOE OSTI2024

Phenalenyl growth reactions and implications for prenucleation chemistry of aromatics in flames

The energetics and kinetics of phenalene and phenalenyl growth reactions were studied theoretically. Rate constants of phenalene and phenalenyl H-abstraction and C 2 H 2 addition to the formed radicals were evaluated through quantum-chemical and rate-theory calculations. The obtained values, assigned to all π radicals, were tested in deterministic and kinetic Monte Carlo simulations of aromatics growth under conditions of laminar premixed flames. Kekulé and non-Kekulé structures of the polycyclic aromatic hydrocarbons (PAHs) evolving in the stochastic simulations were identified by on-the-fly constrained optimization. The numerical results demonstrated an increased PAH growth and qualitatively reproduced experimental observations of Homann and co-workers of non-decaying PAH concentrations with nearly equal abundances of even and odd carbon-atom PAHs. The analysis revealed that the PAH growth proceeds via alternating and sterically diverse acetylene and methyl HACA additions. Furthermore, the rapid and diverse spreading in the PAH population supports a nucleation model as PAH dimerization, assisted by the non-equilibrium phenomena, forming planar aromatics first and then transitioning to the PAH–PAH stacking with size.

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

Comprehensive Kinetics on the C 7 H 7 Potential Energy Surface under Combustion Conditions

The automated kinetics workflow code, KinBot, was used to explore and characterize the regions of the C 7 H 7 potential energy surface that are relevant to combustion environments and especially soot inception. We first explored the lowest-energy region, which includes the benzyl, fulvenallene + H, and cyclopentadienyl + acetylene entry points. We then expanded the model to include two higher-energy entry points, vinylpropargyl + acetylene and vinylacetylene + propargyl. The automated search was able to uncover the pathways from the literature. In addition, three important new routes were discovered: a lower-energy pathway connecting benzyl with vinylcyclopentadienyl, a decomposition mechanism from benzyl that results in side-chain hydrogen atom loss to produce fulvenallene + H, and shorter and lower energy routes to the dimethylene-cyclopentenyl intermediates. We systematically reduced the extended model to a chemically relevant domain composed of 63 wells, 10 bimolecular products, 87 barriers, and 1 barrierless channel and constructed a master equation using the CCSD(T)-F12a/cc-pVTZ//ωB97X-D/6-311++G(d,p) level of theory to provide rate coefficients for chemical modeling. Our calculated rate coefficients show excellent agreement with measured ones. We also simulated concentration profiles and calculated branching fractions from the important entry points to provide an interpretation of this important chemical landscape.

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

Concepts Relevant for the Kinetic Analysis of Reversible Reaction Systems

The net rate of a reversible chemical reaction is the difference between unidirectional rates of traversal along forward and reverse reaction paths. In a multistep reaction sequence, the forward and reverse trajectories, in general, are not the microscopic reverse of one another; rather, each unidirectional route is comprised of distinct rate-controlling steps, intermediates, and transition states. Consequently, traditional descriptors of rate (e.g., reaction orders) do not reflect intrinsic kinetic information but instead conflate unidirectional contributions determined by (i) the microscopic occurrence of forward/reverse reactions (i.e., unidirectional kinetics) and (ii) the reversibility of reaction (i.e., nonequilibrium thermodynamics). This review aims to provide a comprehensive resource of analytical and conceptual tools which deconvolute the contributions of reaction kinetics and thermodynamics to disambiguate unidirectional reaction trajectories and precisely identify rate- and reversibility-controlling molecular species and steps in reversible reaction systems. The extrication of mechanistic and kinetic information from bidirectional reactions is accomplished through equation-based formalisms (e.g., De Donder relations) grounded in principles of thermodynamics and interpreted in the context of theories of chemical kinetics developed in the past 25 years. Here, the aggregate of mathematical formalisms detailed herein is general to thermochemical and electrochemical reactions and encapsulates a diverse body of scientific literature encompassing chemical physics, thermodynamics, chemical kinetics, catalysis, and kinetic modeling.

42 ENGINEERING↗
DOE OSTI2021

Microstructure-Dependent Rate Theory Model of Radiation-Induced Segregation in Binary Alloys

Conventional rate theory often uses the mean field concept to describe the effect of inhomogeneous microstructures on the evolution of radiation induced defect and solute/fission product segregation. However, the spatial and temporal evolution of defects and solutes determines the formation and spatial distribution of radiation-induced second phase such as precipitates and gas bubbles/voids, especially in materials with complicated microstructures and subject to high dose radiation. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed to investigate the effect of inhomogeneous thermodynamic and kinetics properties of defects on diffusion and accumulations of solute A in AB binary alloys. Four independent concentrations: atom A, interstitial A, interstitial B, and vacancy on [A, B] sublattice are used as field variables to describe temporal and spatial distribution and evolution of defects and solute A. The independent concentrations of interstitial A and interstitial B allow to describe their different generation rates, thermodynamic and kinetic properties, and release the assumptions of interstitial generation and sink strength used in the conventional rate theory. Microstructure and concentration dependent chemical potentials of defects are used to calculate the driving forces of defect diffusions. With the model, the effects of defect chemical potentials and mobilities on the RIS in polycrystalline AB model alloys have been simulated. The results demonstrated the model capability in predicting defect evolution in materials with inhomogeneous thermodynamic and kinetic properties of defects. The model can be extended to materials with complicated microstructures such as a wide range of grain size distribution, coating structure and multiphases as well as radiation-induced precipitation subject to severe radiation damage.

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 OSTI2020

Resolving Heterogeneous Dynamics of Excess Protons in Aqueous Solution with Rate Theory

Rate theories have found great utility across the chemical sciences by providing a physically transparent way to analyze dynamical processes. Here we demonstrate the benefits of using transition state theory and Marcus theory to study the rate of proton transfer in HCl solutions. By using long ab initio molecular dynamics simulations we show that good agreement is obtained between these two different formulations of rate theory and how they can be used to study the pathways and life-time of proton transfer in aqueous solution. Since both rate theory formulations utilize identical sets of molecular data, this provides a self-consistent theoretical picture of the rates of aqueous phase proton transfer. Specifically, we isolate and quantify the rates of proton transfer, ion-pair dissociation, and solvent exchange in concentrated HCl solutions. Our analysis predicts a concentration dependence to both proton transfer and ion-pairing. Moreover, our estimate of the life-time for the Zundel species is 0.8 ps and 1.3 ps for 2M and 8M HCl, respectively. In this paper, we demonstrate that concentration effects, can indeed be quantified through the combination of state-of-the-art simulation and theory and provides a picture of the important correlations between the cation (hydronium) and the counter-ion in acid solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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

Substitution Reactions in the Pyrolysis of Acetone Revealed through a Modeling, Experiment, Theory Paradigm

The development of high-fidelity mechanisms for chemically reactive systems is a challenging process that requires the compilation of rate descriptions for a large and somewhat ill-defined set of reactions. The present unified combination of modeling, experiment, and theory provides a paradigm for improving such mechanism development efforts. Here we combine broadband rotational spectroscopy with detailed chemical modeling based on rate constants obtained from automated ab initio transition state theory-based master equation calculations and high-level thermochemical parametrizations. Broadband rotational spectroscopy offers quantitative and isomer-specific detection by which branching ratios of polar reaction products may be obtained. Using this technique, we observe and characterize products arising from H atom substitution reactions in the flash pyrolysis of acetone (CH 3 C(O)CH 3 ) at a nominal temperature of 1800 K. The major product observed is ketene (CH 2 CO). Minor products identified include acetaldehyde (CH 3 CHO), propyne (CH 3 CCH), propene (CH 2 CHCH 3 ), and water (HDO). Literature mechanisms for the pyrolysis of acetone do not adequately describe the minor products. The inclusion of a variety of substitution reactions, with rate constants and thermochemistry obtained from automated ab initio kinetics predictions and Active Thermochemical Tables analyses, demonstrates an important role for such processes. Furthermore, the pathway to acetaldehyde is shown to be a direct result of substitution of acetone's methyl group by a free H atom, while propene formation arises from OH substitution in the enol form of acetone by a free H atom.

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 OSTI2024

Unconventional Error Cancellation Explains the Success of Hartree–Fock Density Functional Theory for Barrier Heights

Energy barriers, which control the rates of chemical reactions, are seriously underestimated by computationally efficient semilocal approximations for the exchange-correlation energy. The accuracy of a semilocal density functional approximation is strongly boosted for reaction barrier heights by evaluating that approximation non-self-consistently on Hartree–Fock electron densities, which has been known for ~30 years. Here, the conventional explanation is that the Hartree–Fock theory yields the more accurate density. This work presents a benchmark Kohn–Sham inversion of accurate coupled-cluster densities for the reaction H 2 + F → HHF → H + HF and finds a strong, understandable cancellation between positive (excessively overcorrected) density-driven and large negative functional-driven errors (expected from stretched radical bonds in the transition state) within this Hartree–Fock density functional theory. This confirms earlier conclusions (Kaplan, A. D., et al. J. Chem. Theory Comput. 2023, 19, 532–543) based on 76 barrier heights and three less reliable, but less expensive, fully nonlocal density functional proxies for the exact density.

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

Microstructure modeling of nuclear structural materials: Recent progress and future directions

Modeling and simulation of microstructures are essential to understand the complex responses and behaviors of nuclear materials in extreme environments. The needs to assess the extended life operation as well as the growing interest in accelerating nuclear materials development and qualification have stimulated the use of high-fidelity multiscale models aided by empirical and ab initio data. This paper reviews the role of various models across different length and time scales in investigating irradiation effects on microstructure evolution and degradation, in particular the embrittlement caused by radiation induced or enhanced formation of nanoscale chemical heterogeneities. The strength and limitations of these models, including classical rate theories, cluster dynamics, phase-field methods, and atomistic models informed by ab initio energies, are discussed with seminal examples. Challenges regarding the lack of thermo-kinetic data and theoretical treatments considering chemical complexities and magnetic excitations, as well as the stabilizing effect by excess point defects in nuclear structural materials are presented, along with potential solutions based on ab initio informed surrogate energy models and statistical sampling by Monte Carlo simulations. Further, the review then highlights the opportunities to leverage the advantages of different methods by establishing hybrid models by shared variables or coupled codes and applications. Finally, the review concludes with forward-looking remarks on how the use of physics-based models can aid the improvement of machine-learning models of property degradation and vice versa.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Quantifying the Decomposition Kinetics of Linear C 2 –C 4 Perfluoroalkanes

The pyrolysis of C 2 –C 4 linear perfluoroalkanes was investigated using laser schlieren densitometry (LS) in a diaphragmless shock tube. Mixtures of 1, 2, and 4% perfluoroethane (PFE), perfluoropropane (PFP), and perfluorobutane (PFB) in excess krypton were shock-heated to 1400–2500 K at pressures of 60, 120, and 240 Torr. The initiation step for all perfluoroalkanes involved carbon–carbon bond fission. The rate constants for unimolecular decomposition were computed using microcanonical rate theory, and these rate constants were combined into a detailed chemical kinetic mechanism. For PFE and PFP, there was one dissociation channel, while PFB dissociated via two competing bond-fission reactions. The perfluoroalkyl radicals subsequently recombined or underwent further unimolecular decomposition. In conclusion, the kinetic model accurately simulated the density gradient profiles of all perfluoroalkanes.

Guzman, Eduardo H. [Brown Univ., Providence, RI (U↗
DOE OSTI2024

Quantum Simulations of Radiation Damage in a Molecular Polyethylene Analog

Abstract An atomic‐level understanding of radiation‐induced damage in simple polymers like polyethylene is essential for determining how these chemical changes can alter the physical and mechanical properties of important technological materials such as plastics. Ensembles of quantum simulations of radiation damage in a polyethylene analog are performed using the Density Functional Tight Binding method to help bind its radiolysis and subsequent degradation as a function of radiation dose. Chemical degradation products are categorized with a graph theory approach, and occurrence rates of unsaturated carbon bond formation, crosslinking, cycle formation, chain scission reactions, and out‐gassing products are computed. Statistical correlations between product pairs show significant correlations between chain scission reactions, unsaturated carbon bond formation, and out‐gassing products, though these correlations decrease with increasing atom recoil energy. The results present relatively simple chemical descriptors as possible indications of network rearrangements in the middle range of excitation energies. Ultimately, the work provides a computational framework for determining the coupling between nonequilibrium chemistry in polymers and potential changes to macro‐scale properties that can aid in the interpretation of future radiation damage experiments on plastic materials.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Physiochemical Machine Learning Models Predict Operational Lifetimes of CH3NH3PbI3 Perovskite Solar Cells

Halide perovskites are promising photovoltaic (PV) materials with the potential to lower the cost of electricity and greatly expand the penetration of PV if they can demonstrate long-term stability under illumination in the presence of moisture and oxygen. The solar cell service lifetime as quantified by the T80 (the time required for the power conversion efficiency to drop to 80% of its starting value) is a useful metric to assess stability. The T80 for utility, commercial, or residential PV systems needs to be several decades in order to yield low-cost electricity, and thus it is not practical to directly measure the T80. It would be useful if T80 could be predicted from the initial dynamics of a solar cell’s performance, but until now no models have been developed to forecast T80. In this work, we report the development of machine learning models to predict T80 of ITO/NiOx/CH3NH3PbI3/C60/BCP/Ag solar cells operating at maximum power point under 1-sun equivalent photon flux in air at varying temperatures and relative humidities. Efficiency losses are driven by short-circuit current and fill factor, indicating that chemical decomposition of the perovskite is a major contributor to degradation. Spatial patterns evident from in situ dark field optical microscopy suggest that the electric field gradient at device edges plays a significant role in perovskite decomposition, along with photochemical reactions with O2 and H2O. Models are trained using a menu of features from three distinct categories: (i) features based on measurements of the initial rates of change of device parameters, (ii) features based on the ambient conditions during operation (temperature, & partial pressure of H2O), and (iii) features based on underlying physics and chemistry. We show that a theory-based physiochemical feature derived from a model of the chemical reaction kinetics of the rate of degradation of the CH3NH3PbI3 is particularly valuable for prediction. This physiochemical feature was selected as the first or second most dominant feature in the best performing models. With a dataset consisting of 45 accelerated degradation experiments with T80 that range over a factor of 30, the model predicts T80 with an accuracy of about 40% (|predicted T80 - observed T80| / observed T80) on samples not used in training. This hybrid ML approach should be effective when applied to other compositions, device architectures, and advanced packaging schemes.

14 SOLAR ENERGY↗