Search NASA⌕ Search

SEARCH · Search NASA

Results for “reactive force fields”

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

Bond order predictions using deep neural networks

Machine learning is an extremely powerful tool for the modern theoretical chemist since it provides a method for bypassing costly algorithms for solving the Schrödinger equation. Already, it has proven able to infer molecular and atomic properties such as charges, enthalpies, dipoles, excited state energies, and others. Most of these machine learning algorithms proceed by inferring properties of individual atoms, even breaking down total molecular energy into individual atomic contributions. In this paper, we introduce a modified version of the Hierarchically Interacting Particle Neural Network (HIP-NN) capable of making predictions on the bonds between atoms rather than on the atoms themselves. We train the modified HIP-NN to infer bond orders for a large number of small organic molecules as computed via the Natural Bond Orbital package. We demonstrate that the trained model is extensible to molecules much larger than those in the training set by studying its performance on the COMP6 dataset. This method has applications in cheminformatics and force field parameterization and opens a promising future for machine learning models to predict other quantities that are defined between atoms such as density matrix elements, Hamiltonian parameters, and molecular reactivities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2002

Evaluation of Transport in the Lower Tropical Stratosphere in a Global Chemistry and Transport Model

A general circulation model (GCM) relies on various physical parameterizations and provides a solution to the atmospheric equations of motion. A data assimilation system (DAS) combines information from observations with a GCM forecast and produces analyzed meteorological fields that represent the observed atmospheric state. An off-line chemistry and transport model (CTM) can use winds and temperatures from a either a GCM or a DAS. The latter application is in common usage for interpretation of observations from various platforms under the assumption that the DAS transport represents the actual atmospheric transport. Here we compare the transport produced by a DAS with that produced by the particular GCM that is combined with observations to produce the analyzed fields. We focus on transport in the tropics and middle latitudes by comparing the age-of-air inferred from observations of SF6 and CO2 with the age-of-air calculated using GCM fields and DAS fields. We also compare observations of ozone, total reactive nitrogen, and methane with results from the two simulations. These comparisons show that DAS fields produce rapid upward tropical transport and excessive mixing between the tropics and middle latitudes. The unrealistic transport produced by the DAS fields may be due to implicit forcing that is required by the assimilation process when there is bias between the GCM forecast and observations that are combined to produce the analyzed fields. For example, the GCM does not produce a quasi-biennial oscillation (QBO). The QBO is present in the analyzed fields because it is present in the observations, and systematic implicit forcing is required by the DAS. Any systematic bias between observations and the GCM forecast used to produce the DAS analysis is likely to corrupt the transport produced by the analyzed fields. Evaluation of transport in the lower tropical stratosphere in a global chemistry and transport model.

Douglass, Anne R.↗
DOE OSTI2017

Crack propagation in silica from reactive classical molecular dynamics simulations

Mechanistic insight into the process of crack growth can be obtained through molecular dynamics (MD) simulations. In this investigation of fracture propagation, a slit crack was introduced into an atomistic amorphous silica model and mode I stress was applied through far–field loading until the crack propagates. Atomic displacements and forces and an Irving–Kirkwood method with a Lagrangian kernel estimator were used to calculate the J –integral of classical fracture mechanics around the crack tip. The resulting fracture toughness (K IC ), 0.76 ± 0.16 MPa√m, agrees with experimental values. In addition, the stress fields and dissipation energies around the slit crack indicate the development of an inelastic region ~30Å in diameter. This is one of the first reports of KIC values obtained from up–scaled atomic–level energies and stresses through the J –integral. In conclusion, the application of the ReaxFF classical MD force field in this study provides the basis for future research into crack growth in multicomponent oxides in a variety of environmental conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Effect of Domain Size, Boundary, and Loading Conditions on Mechanical Properties of Amorphous Silica: A Reactive Molecular Dynamics Study

Mechanical properties are very important when choosing a material for a specific application. They help to determine the range of usefulness of a material, establish the service life, and classify and identify materials. The size effect on mechanical properties has been well established numerically and experimentally. However, the role of the size effect combined with boundary and loading conditions on mechanical properties remains unknown. In this paper, by using molecular dynamics (MD) simulations with the state-of-the-art ReaxFF force field, we study mechanical properties of amorphous silica (e.g., Young’s modulus, Poisson’s ratio) as a function of domain size, full-/semi-periodic boundary condition, and tensile/compressive loading. We found that the domain-size effect on Young’s modulus and Poisson’s ratio is much more significant in semi-periodic domains compared to full-periodic domains. The results, for the first time, revealed the bimodular and anisotropic nature of amorphous silica at the atomic level. We also defined a “safe zone” regarding the domain size, where the bulk properties of amorphous silica can be reproducible, while the computational cost and accuracy are in balance.

bimodular materials↗
DOE OSTI2022

The influence of particle size on the fluid dynamics of a laser-induced plasma

The interaction of a laser-induced shock wave with nanoparticles and microparticles of aluminum oxide is investigated through experiments and modeling. The chemistry and physics of the interaction between the particles and plasma generated from laser ablation shows similarities and discrete differences for the two particle sizes. For both particle sizes, early stage (< 10 μs) ionization was dominant and evidenced by higher concentrations of Al II. While both sizes exhibit ionization over the same duration, the intensity of emission was greater for nanoparticles indicating greater concentrations of ionized species. Moreover, the dispersion of species was notably more elongated for microparticles while radial dispersion was more pronounced for nanoparticles with elevated drag forces. At later stages (i.e., > 10 μs), oxidation reactions were dominant for both particle sizes, but the same distinctions in flow field were observed and attributed to particle drag. In all stages of interaction, microparticles expand axially with less drag that suppresses their radial expansion. As a result, the dispersion of reactive species was mapped over an up to 80% larger area for nanoparticles relative to microparticles. Results shown here can be applied toward advancing experimental diagnostics and particle-shock wave modeling and simulation efforts for energetic materials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nanoscale wetting controls reactive Pd ensembles in synthesis of dilute PdAu alloy catalysts

The performance of bimetallic dilute alloy catalysts is largely determined by the size of minority metal ensembles on the nanoparticle surface. By analyzing the synthesis of catalysts comprising Pd 8 Au 92 nanoparticles supported on silica using surface-sensitive techniques, we report that whether Pd overgrowth occurs before or after Au nanoparticle deposition onto the support controls the surface Pd ensemble size and abundance. These differences in Pd ensembles influence catalytic reactivity in H 2 –D 2 isotope exchange and benzaldehyde hydrogenation, which, in correlation with theoretical calculations, is used to elucidate the active site(s) in each reaction. To clarify how the synthetic sequence controls the formation of Pd ensembles, we combine numerical wetting calculations and molecular dynamics simulations (with a machine-learned force field) to visualize Pd deposition and migration on the nanoparticle surface, respectively. Our results suggest that the nanoparticle–support interface restricts nanoparticle accessibility to Pd deposition, which consequently controls the Pd ensemble size, illustrating the critical role of nanoscale wetting phenomena during bimetallic catalyst preparation.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Reactive optical matter: light-induced motility in electrodynamically asymmetric nanoscale scatterers

From Newton’s third law, which is known as the principle of actio et reactio, we expect the forces between interacting particles to be equal and opposite for closed systems. Otherwise, “nonreciprocal” forces can arise. This has been shown theoretically in the interaction between dissimilar optically trapped particles that are mediated by an external field. As a result, despite the incident external field not having a transverse component of momentum, the particle pair experiences a force in a direction that is transverse to the light propagation direction. In this letter, we directly measure the net nonreciprocal forces in electrodynamically interacting asymmetric nanoparticle dimers and nanoparticle structures that are illuminated by plane waves and confined to pseudo one-dimensional geometries. We show via electrodynamic theory and simulations that interparticle interactions cause asymmetric scattering from heterodimers. Therefore, the nonreciprocal forces are a consequence of momentum conservation. Our study demonstrates that asymmetric scatterers exhibit directed motion due to the breakdown of mirror symmetry in open systems that involve incident external fields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1990

Exploring the notion of space coupling propulsion

All existing methods of space propulsion are based on expelling a reaction mass (propellant) to induce motion. Alternatively, 'space coupling propulsion' refers to speculations about reacting with space-time itself to generate propulsive forces. Conceivably, the resulting increases in payload, range, and velocity would constitute a breakthrough in space propulsion. Such speculations are still considered science fiction for a number of reasons: (1) it appears to violate conservation of momentum; (2) no reactive media appear to exist in space; (3) no 'Grand Uniform Theories' exist to link gravity, an acceleration field, to other phenomena of nature such as electrodynamics. The rationale behind these objectives is the focus of interest. Various methods to either satisfy or explore these issues are presented along with secondary considerations. It is found that it may be useful to consider alternative conventions of science to further explore speculations of space coupling propulsion.

Millis, Marc G.↗
DOE OSTI2022

Modular Approach for the Synthesis of Bottlebrush Diblock Copolymers from Poly(Glycidyl Methacrylate)-block-Poly(Vinyldimethylazlactone) Backbones

A strategy that enables the facile synthesis of bottlebrush block copolymers with flexible backbones was developed. A demonstration of the strategy’s utility was carried out by grafting end-functionalized polymethylmethacrylate (PMMA) and polystyrene (PS) to the dually reactive block copolymer, poly(glycidyl methacrylate)-block-poly(vinyldimethylazlactone) (PGMA-b-PVDMA). Five different bottlebrush diblock copolymers were investigated by size-exclusion chromatography (SEC), 1H NMR, Fourier transform infrared (FT-IR), differential scanning calorimetry (DSC), X-ray scattering methods, atomic force microscopy (AFM), rheology and computational simulations using molecular dynamics (MD), and self-consistent field theory (SCFT). A relationship between the glass transition temperature and the fraction of chain ends was demonstrated by DSC and highlights the potential of this synthetic method to tailor polymer properties. All five samples were found to be in a disordered phase exhibiting multiscale structures revealed by two broad peaks in small-angle X-ray scattering (SAXS) that can be attributed to graft-to-graft and backbone-to-backbone density correlations using MD simulations. The SCFT-based simulations justify the observation of a disordered phase due to its stabilization by the grafts. Additionally, this modular approach can be easily extended to other grafts, including responsive, conducting, and charged polymers with the prerequisite end groups. The versatility and ease of assembling these functional bottlebrushes constitute a powerful “toolbox” method for the rapid and scalable synthesis of novel bottlebrush block copolymers with desired properties.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Li + solvation and kinetics of Li + –BF 4 - /PF 6 - ion pairs in ethylene carbonate. A molecular dynamics study with classical rate theories

Using our polarizable force-field models and employing classical rate theories of chemical reactions, we examine the ethylene carbonate (EC) exchange process between the first and second solvation shells around Li+ and the dissociation kinetics of ion pairs Li+-[BF4] and Li+-[PF6] in this solvent. We calculate the exchange rates using transition state theory and correct them with transmission coefficients computed by the reactive flux; Impey, Madden, and McDonald approaches; and Grote-Hynes theory. We found the residence times of EC around Li+ ions varied from 70 to 450 ps, depending on the correction method used. We found the relaxation times changed significantly from Li+-[BF4] to Li+-[PF6] ion pairs in EC. Our results also show that, in addition to affecting the free energy of dissociation in EC, the anion type also significantly influence the dissociation kinetics of ion pairing. This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, and Biosciences. The calculations were carried out using computer resources provided by the Office of Basic Energy Sciences.

Chang, Tsun-Mei↗
DOE OSTI2021

Accurate simulation of surfaces and interfaces of ten FCC metals and steel using Lennard–Jones potentials

Abstract The earlier integration of validated Lennard–Jones (LJ) potentials for 8 fcc metals into materials and biomolecular force fields has advanced multiple research fields, for example, metal–electrolyte interfaces, recognition of biomolecules, colloidal assembly of metal nanostructures, alloys, and catalysis. Here we introduce 12-6 and 9-6 LJ parameters for classical all-atom simulations of 10 further fcc metals (Ac, Ca (α), Ce (γ), Es (β), Fe (γ), Ir, Rh, Sr (α), Th (α), Yb (β)) and stainless steel. The parameters reproduce lattice constants, surface energies, water interfacial energies, and interactions with (bio)organic molecules in 0.1 to 5% agreement with experiment, as well as qualitative mechanical properties under standard conditions. Deviations are reduced up to a factor of one hundred in comparison to earlier Lennard–Jones parameters, embedded atom models, and density functional theory. We also explain a quantitative correlation between atomization energies from experiments and surface energies that supports parameter development. The models are computationally very efficient and applicable to an exponential space of alloys. Compatibility with a wide range of force fields such as the Interface force field (IFF), AMBER, CHARMM, COMPASS, CVFF, DREIDING, OPLS-AA, and PCFF enables reliable simulations of nanostructures up to millions of atoms and microsecond time scales. User-friendly model building and input generation are available in the CHARMM-GUI Nanomaterial Modeler. As a limitation, deviations in mechanical properties vary and are comparable to DFT methods. We discuss the incorporation of reactivity and features of the electronic structure to expand the range of applications and further increase the accuracy.

Kanhaiya, Krishan (ORCID:0000000236222655)↗
DOE OSTI2019

Sodium Halide Adsorption and Water Structure at the α-Alumina(0001)/Water Interface

Alumina is one of the most abundant minerals and has a wide range of industrial applications, with catalysis as one of the most important. Of particular relevance for catalysis is the structure of the mineral/water interface. In this work, water structure and sodium halide adsorption at the neutral α-alumina(0001)/water interface are investigated using molecular dynamics simulations. This work demonstrates the accuracy of the chosen model of the alumina/water interface and shows that high charge density monovalent ions, such as Na+ and F–, have a strong affinity for the interface due to the specific pattern of alumina surface OH groups, such that the adsorbed ions displace waters that are hydrogen-bonded to the surface in their absence. A significant portion of the driving force for anion adsorption arises from surface bound Na+, which reverse the intrinsic surface dipole field and drive the accumulation of halides at the interface. The resulting electrolytic interfacial structure reorients water molecules as far as 1 nm from the surface. Although ion adsorption does not alter the global orientation of surface OH groups, it significantly affects their local geometry. This in turn may affect the reactivity of surface groups and thus play a role in chemical processes occurring at the interface.

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

Recent Advances in the Application of Functionalized Lignin in Value-Added Polymeric Materials

The quest for converting lignin into high-value products has been continuously pursued in the past few decades. In its native form, lignin is a group of heterogeneous polymers comprised of phenylpropanoids. The major commercial lignin streams, including Kraft lignin, lignosulfonates, soda lignin and organosolv lignin, are produced from industrial processes including the paper and pulping industry and emerging lignocellulosic biorefineries. Although lignin has been viewed as a low-cost and renewable feedstock to replace petroleum-based materials, its utilization in polymeric materials has been suppressed due to the low reactivity and inherent physicochemical properties of lignin. Hence, various lignin modification strategies have been developed to overcome these problems. Herein, we review recent progress made in the utilization of functionalized lignins in commodity polymers including thermoset resins, blends/composites, grafted functionalized copolymers and carbon fiber precursors. In the synthesis of thermoset resins such as polyurethane, phenol-formaldehyde and epoxy, they are covalently incorporated into the polymer matrix, and the discussion is focused on chemical modifications improving the reactivity of technical lignins. In blends/composites, functionalization of technical lignins is based upon tuning the intermolecular forces between polymer components. In addition, grafted functional polymers have expanded the utilization of lignin-based copolymers to biomedical materials and value-added additives. Different modification approaches have also been applied to facilitate the application of lignin as carbon fiber precursors, heavy metal adsorbents and nanoparticles. These emerging fields will create new opportunities in cost-effectively integrating the lignin valorization into lignocellulosic biorefineries.

36 MATERIALS SCIENCE↗
DOE OSTI2018

The Unified Electrochemical Band Diagram Framework: Understanding the Driving Forces of Materials Electrochemistry

Recent work that establishes a picture of the driving forces that govern material transformations and degradation in electrochemical environments to enable the ab initio design of electrochemical materials is highlighted. Select prototype systems are used to describe how the interplay between materials properties such as crystal field splitting, band edge energies, surface termination, material length scale, dielectric constant, and isoelectric point, and electrolyte properties such as pH and ion type, impacts electrochemical behavior - i.e., redox potentials, reaction enthalpies, reactivity, and decoupled ionic/electronic processes. Ab initio modeling of charged defects and intercalants within the grand canonical unified electrochemical band-diagram (UEB) framework is shown to enable the quantitative prediction of electrochemical materials behavior. UEB combines electrochemical theory, charged defect theory, and band diagram descriptions and can be used both for materials discovery and development. First, a pedagogical description of the UEB framework is presented, and then the application of this framework to reveal mechanisms for high rate electronic charge storage in cation incorporated α.-MnO 2 and λ-MnO 2 , high desalination efficiency of thin-film NaMn 4 O 8 , and the flat charge/discharge profile of FePO 4 is reviewed. Finally, new prospects for the application of the UEB framework to electrolyte design, interfacial engineering, and catalysis are suggested.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electric-Field-Induced Assembly of an Ionic Liquid–Water Interphase Enables Efficient Heavy Metal Electrosorption

Controlling ion desolvation, transport, and charge transfer at the electrode-electrolyte interface (EEI) is critical to enabling the rational design of efficient and selective separation of targeted heavy metals and decontamination of industrial waste water. One challenge is to sufficiently resolve and interrogate the intermediate transformation steps between solvated metal cations and bare cations for electroreduction at the EEI and establish pathways to modulate these steps to achieve efficient energy transfer for targeted reactive separations. Herein, we obtained a predictive understanding of the effect of adding a hydrophobic ionic liquid, 1-ethyl-3-methylimidazolium chloride (EMIMCl), to an aqueous electrolyte on modulating the desolvation and electrosorption of Pb 2+ cations using a combination of experimental and theoretical techniques. Here, we reveal formation of a compact interphase layer consisting of EMIMCl-Pb complex clusters under an applied electric field using operando electrochemical Raman spectroscopy, atomic force microscopy, and electrochemical impedance spectroscopy measurements combined with classical molecular dynamics simulations. The application of a larger negative potential is shown to result in formation of a well oriented layer with the positive imidazolium ring of the EMIMCl lying on the electrode and the long hydrophobic alkyl chain extending into the bulk electrolyte. This oriented layer is demonstrated to facilitate desolvation of incoming solvated Pb 2+ cations and decrease the charge transfer resistance for Pb electrodeposition, which has important implications for selective removal of Pb from contaminated mixtures. Overall, our studies open up new opportunities to modulate ion desolvation using hydrophobic ionic liquids in aqueous electrolytes for efficient heavy metal separation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Summary Report of the Reactive CO 2 Capture: Process Integration for the New Carbon Economy Workshop

Decarbonization of our global economy is required to limit planetary warming to +1.5⁰C above pre-industrial levels, an ambitious goal set into motion by the Paris agreement. Given the scale and urgency, the solution demands international, cross-sector advancements spanning policy, social responsibility, and technology, with emphasis on step changes over incremental changes. To that end, technological revolutions that disrupt the status quo need to be envisioned and enacted. One potential technological revolution is the production of fuels, chemicals, and materials from carbon dioxide (CO 2 ) as the starting feedstock, leveraging renewable energy as the driving force. Over the past decades, significant research, development, and deployment has occurred on technologies for capturing CO 2 from point sources or the air and utilizing this CO 2 as a working fluid or as a chemical reactant; however, most of this work has been siloed in these two categories. Recently, an emerging field has started to explore the direct integration of CO 2 capture and conversion technologies as a means to reduce overall energy demand (i.e., avoid energy penalty of CO 2 desorption/regeneration of capture media) and capital expense through process intensification. This strategy represents an opportunity to leapfrog forward this technological revolution. However, the field is in its infancy and the technologies are at an early stage of development, thus it is critically important to define and assess the value proposition of this strategy relative to alternatives (e.g., separated capture and conversion technologies, fuels and chemicals derived from renewable feedstocks like biomass, and industrial electrification) to chart a path forward. To identify next steps, we organized a workshop titled “Reactive CO 2 Capture: Process Integration for the New Carbon Economy” which was held in Golden, Colorado, February 18–19, 2020. The focus of this workshop was to discuss approaches for merging CO 2 capture and CO 2 conversion/utilization systems into what we denoted as an integrated "reactive capture" strategy. By our definition, reactive capture of CO 2 is the coupled process of capturing CO 2 from a mixed gas stream and converting it into a valuable product without going through a purified CO 2 intermediate (see full definition in the Introduction section). This report seeks to summarize feedback from the approximately 125 participants and subject matter experts in attendance from academia, industry, U.S. Department of Energy (DOE), and DOE national laboratories. The workshop agenda is included in Appendix A and the full list of attendees can be found in Appendix B. To elucidate a path forward, we first asked the attendees to define what success would look like for reactive capture in the short term (0–5 years), midterm (5–10 years), and long term (10+ years) and then asked them to answer four questions related to how we could achieve that success: (1) What are the key barriers and challenges to success? (2) What are needed activities to overcome barriers and challenges? (3) What opportunities will arise from these activities? (4) What is a target outcome and what metrics need to be met?

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2019

Molecular Dynamics of Combustion Reactions in Supercritical Carbon Dioxide. Part 5: Computational Study of Ethane Dissociation and Recombination Reactions C 2 H 6 ⇌ CH 3 + CH 3

Fossil fuel oxy-combustion is an emerging technology where the habitual nitrogen diluent is replaced by high-pressure supercritical CO 2 (sCO 2 ), which increases the efficiency of energy conversion. In this study, the chemical kinetics of the combustion reaction C 2 H 6 ⇌ CH 3 + CH 3 in the sCO 2 environment is predicted at 30–1000 atm and 1000–2000 K. We adopt a multiscale approach, where the reactive complex is treated quantum mechanically in rigid rotor/harmonic oscillator approximation, while environment effects at different densities are taken into account by the potential of mean force, produced with classical molecular dynamics (MD). Here, we used boxed MD, where enhanced sampling of infrequent events of barrier crossing is accomplished without application of the bias potential. The multistate empirical valence bond model is applied to describe free radical formation accurately at the cost of the classical force field. Predicted rates at low densities agree well with the literature data. In conclusion, rate constants at 300 atm are 2.41 × 10 14 T –0.20 exp(–77.03 kcal/mol/RT) 1/s for ethane dissociation and 8.44 × 10 –19 T 1.42 exp(19.89 kcal/mol/RT) cm 3 /molecule/s for methyl–methyl recombination.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗