Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction dynamics”

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
NASA NTRS1971

The molecular dynamics of atmospheric reaction

Detailed information about the chemistry of the upper atmosphere took the form of quantitative data concerning the rate of reaction into specified states of product vibration, rotation and translation for exothermic reaction, as well as concerning the rate of reaction from specified states of reagent vibration, rotation and translation for endothermic reaction. The techniques used were variants on the infrared chemiluminescence method. Emphasis was placed on reactions that formed, and that removed, vibrationally-excited hydroxyl radicals. Fundamental studies were also performed on exothermic reactions involving hydrogen halides.

Polanyi, J. C.↗
DOE OSTI2021

A Chemical Dynamics Study of the Reaction of the Methylidyne Radical (CH, X 2 Π) with Dimethylacetylene (CH 3 CCCH 3 , X 1 A 1g )

The gas-phase reaction of the methylidyne (CH; X 2 Π) radical with dimethylacetylene (CH 3 CCCH 3 ; X 1 A 1g ) was studied at a collision energy of 20.6 kJ mol -1 under single collision conditions with experimental results merged with ab initio calculations of the potential energy surface (PES) and ab initio molecule dynamics (AIMD) simulations. Here, the crossed molecular beams experiment reveals that the reaction proceeds barrierless via indirect scattering dynamics through long-lived C 5 H 7 reaction intermediate(s) ultimately dissociating to C 5 H 6 isomers along with atomic hydrogen with atomic hydrogen predominantly released from the methyl groups as verified by replacing the methylidyne with the D1-methylidyne reactant. AIMD simulations reveal that the reaction dynamics are statistical leading predominantly to p28 (1-methyl-3-methylene-cyclopropene, 13 %) and p8 (1-penten-3-yne, 81 %) plus atomic hydrogen with a significant amount of available energy being channeled into the internal excitation of the polyatomic reaction products.

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

Signatures of Non-universal Quantum Dynamics of Ultracold Chemical Reactions of Polar Alkali Dimer Molecules with Alkali Metal Atoms: Li( 2 S) + NaLi( a 3 Σ + ) → Na( 2 S) + Li 2 ( a 3 Σ u + )

Ultracold chemical reactions of weakly bound triplet-state alkali metal dimer molecules have recently attracted much experimental interest. Here we perform rigorous quantum scattering calculations with a new ab initio potential energy surface to explore the chemical reaction of spin-polarized NaLi(a 3 Σ + ) and Li( 2 S) to form Li 2 (a 3 Σ u + ) and Na( 2 S). The reaction is exothermic and proceeds readily at ultralow temperatures. Significantly, we observe strong sensitivity of the total reaction rate to small variations of the three-body part of the Li 2 Na interaction at short range, which we attribute to a relatively small number of open Li 2 (a 3 Σ u + ) product channels populated in the reaction. This provides the first signature of highly non-universal dynamics seen in rigorous quantum reactive scattering calculations of an ultracold exothermic insertion reaction involving a polar alkali dimer molecule, opening up the possibility of probing microscopic interactions in atom+molecule collision complexes via ultracold reactive scattering experiments.

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

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

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

Fast predictions of liquid-phase acid-catalyzed reaction rates using molecular dynamics simulations and convolutional neural networks

The rates of liquid-phase, acid-catalyzed reactions relevant to the upgrading of biomass into high-value chemicals are highly sensitive to solvent composition and identifying suitable solvent mixtures is theoretically and experimentally challenging. We show that the complex atomistic configurations of reactant–solvent environments generated by classical molecular dynamics simulations can be exploited by 3D convolutional neural networks to enable accurate predictions of Brønsted acid-catalyzed reaction rates for model biomass compounds. We develop a 3D convolutional neural network, which we call SolventNet, and train it to predict acid-catalyzed reaction rates using experimental reaction data and corresponding molecular dynamics simulation data for seven biomass-derived oxygenates in water–cosolvent mixtures. We show that SolventNet can predict reaction rates for additional reactants and solvent systems an order of magnitude faster than prior simulation methods. This combination of machine learning with molecular dynamics enables the rapid, high-throughput screening of solvent systems and identification of improved biomass conversion conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1983

Dynamic fracture toughnesses of reaction-bonded silicon nitride

The room-temperature dynamic fracture response of reaction-bonded silicon nitride is investigated using a hybrid experimental-numerical procedure. In this procedure, experimentally determined crack velocities are utilized to drive a dynamic finite-element code or dynamic finite-difference code in its generation mode in order to extract numerically the dynamic stress intensity factor of the fracturing specimen. Results show that the dynamic fracture toughness vs crack velocity relations of the two reaction-bonded silicon nitrides do not follow the general trend in those relations of brittle polymers and steel. A definite slow crack velocity during the initial phase of dynamic crack propagation is observed in reaction-bonded silicon nitride, which results in a nonunique dynamic fracture toughness vs crack velocity relation. In addition, it is found that a propagating crack will continue to propagate under a static stress intensity factor substantially lower than K(IC).

Kobayashi, A. S.↗
DOE OSTI2026

Covalent Adaptable Networks: Reprocessable Cross-Linked Polymers

Thermoset polymers have desirable properties, such as excellent thermal and mechanical stability, but their covalent cross-links typically prevent repair or recycling. By enabling and controlling dynamic exchange reactions within polymer networks, their covalent bonds rearrange and allow the polymer to be reshaped. These viscoelastic polymer networks, now known as covalent adaptable networks (CANs), are an important frontier for improving plastic circularity, as well as for designing valuable stimuli-responsive materials. This Review describes the history of CANs, dating back to the early days of polymer science, and the evolution of their classification and nomenclature. A comprehensive survey of dynamic reactions and linkage chemistries is provided, as well as methods to characterize and reprocess CANs. Beyond straightforward reprocessing, many advanced applications of CANs and their composites are now emerging. Lastly, we provide perspective on how the development of new chemistries, strategies to control stimuli-responsive bond exchange and mechanical properties, and a deep understanding of exchange reactions will advance this field toward scalable, sustainable, and high-value materials.

Sun, Molly [Northwestern University, Evanston, IL ↗
DOE OSTI2020

Revealing nanoscale dynamics during an epoxy curing reaction with x-ray photon correlation spectroscopy

In this work the evolution of nanoscale properties is measured during the thermally triggered curing of an industrial epoxy adhesive. We use x-ray photon correlation spectroscopy (XPCS) to track the progression of the curing reaction through the local dynamics of filler particles that reflect the formation of a thermoset network. Out-of-equilibrium dynamics are resolved through identification and analysis of the intensity– intensity autocorrelation functions obtained from XPCS. The characteristic time scale and local velocity of the filler is calculated as functions of time and temperature. We find that the dynamics speed up when approaching the curing temperature (T cure ), and decay rapidly once T cure is reached. We compare the results from XPCS to conventional macroscale characterization by differential scanning calorimetry (DSC). The demonstration and implementation of nanoscale characterization of curing reactions by XPCS proves useful for future development and optimization of epoxy thermoset materials and other industrial adhesive systems.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Elucidating the chemical dynamics of the elementary reactions of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with 2-methylpropene ((CH 3 ) 2 CCH 2 ; X 1 A 1 )

Exploiting the crossed molecular beam technique, we studied the reaction of the 1-propynyl radical (CH 3 CC; X 2 A 1 ) with 2-methylpropene (isobutylene; (CH 3 ) 2 CCH 2 ; X 1 A 1 ) at a collision energy of 38 ± 3 kJ mol –1 . The experimental results along with ab initio and statistical calculations revealed that the reaction has no entrance barrier and proceeds via indirect scattering dynamics involving C 7 H 11 intermediates with lifetimes longer than their rotation period(s). The reaction is initiated by the addition of the 1-propynyl radical with its radical center to the π-electron density at the C1 and/or C2 position in 2-methylpropene. Further, the C 7 H 11 intermediate formed from the C1 addition either emits atomic hydrogen or undergoes isomerization via [1,2-H] shift from the CH 3 or CH 2 group prior to atomic hydrogen loss preferentially leading to 1,2,4-trimethylvinylacetylene (2-methylhex-2-en-4-yne) as the dominant product. The molecular structures of the collisional complexes promote hydrogen atom loss channels. RRKM results show that hydrogen elimination channels dominate in this reaction, with a branching ratio exceeding 70%. Since the reaction of the 1-propynyl radical with 2-methylpropene has no entrance barrier, is exoergic, and all transition states involved are located below the energy of the separated reactants, bimolecular collisions are feasible to form trimethylsubstituted 1,3-enyne (p1) via a single collision event even at temperatures as low as 10 K prevailing in cold molecular clouds such as G+0.693. As a result, the formation of trimethylsubstituted vinylacetylene could serve as the starting point of fundamental molecular mass growth processes leading to di- and trimethylsubstituted naphthalenes via the HAVA mechanism.

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

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Insight into SEI Growth in Li-Ion Batteries using Molecular Dynamics and Accelerated Chemical Reactions

The performance of lithium-ion batteries (LIB) using organic electrolytes strongly depends on the formation of a stable solid electrolyte interphase (SEI) film. Elucidating the dynamic evolution and spatial composition of the SEI can be very useful to study the stability of SEI components and help optimize the formation cycles of LIB. We propose a classical molecular dynamics simulation protocol for predicting the first stages of SEI formation using a reaction method involving the decomposition of EC and LiPF 6 molecules in the electrolyte. We accelerate the formation of SEI components near the anode surface by increasing the probability of reactions, implemented through a geometry matching scheme, followed by a force-field reconfiguration. We observe the formation of gases (C 2 H 4 ), inorganic (Li 2 CO 3 and LiF) and organic (LEDC) components. Finally, this protocol shows promise to be able to evaluate the effects of varying electrolyte compositions and additives on SEI layer structure and composition.

25 ENERGY STORAGE↗
DOE OSTI2020

Reactive Molecular Dynamics Simulations and Quantum Chemistry Calculations To Investigate Soot-Relevant Reaction Pathways for Hexylamine Isomers

Sooting tendencies of a series of nitrogen-containing hydrocarbons (NHCs) have been recently characterized experimentally using the yield sooting index (YSI) methodology. This work aims to identify soot-relevant reaction pathways for three selected C6H15N amines, namely, dipropylamine (DPA), diisopropylamine (DIPA), and 3,3-dimethylbutylamine (DMBA) using ReaxFF molecular dynamics (MD) simulations and quantum mechanical (QM) calculations and to interpret the experimentally observed trends. ReaxFF MD simulations are performed to determine the important intermediate species and radicals involved in the fuel decomposition and soot formation processes. QM calculations are employed to extensively search for chemical reactions involving these species and radicals based on the ReaxFF MD results and also to quantitatively characterize the potential energy surfaces. Specifically, ReaxFF simulations are carried out in the NVT ensemble at 1400, 1600, and 1800 K, where soot has been identified to form in the YSI experiment. These simulations account for the interactions among test fuel molecules and pre-existing radicals and intermediate species generated from rich methane combustion, using a recently proposed simulation framework. ReaxFF simulations predict that the reactivity of the amines decrease in the order DIPA > DPA > DMBA, independent of temperature. Both QM calculations and ReaxFF simulations predict that C2H4, C3H6, and C4H8 are the main nonaromatic soot precursors formed during the decomposition of DPA, DIPA, and DMBA, respectively, and the associated reaction pathways are identified for each amine. Both theoretical methods predict that sooting tendency increases in the order DPA, DIPA, and DMBA, consistent with the experimentally measured trend in YSI. This work demonstrates that sooting tendencies and soot-relevant reaction pathways of fuels with unknown chemical kinetics can be identified efficiently through combined ReaxFF and QM simulations. Overall, predictions from ReaxFF simulations and QM calculations are consistent, in terms of fuel reactivity, major intermediates, and major nonaromatic soot precursors.

fuel decomposition↗
DOE OSTI2023

Gas phase synthesis of the C40 nano bowl C 40 H 10

Nanobowls represent vital molecular building blocks of end-capped nanotubes and fullerenes detected in combustion systems and in deep space such as toward the planetary nebula TC-1, but their fundamental formation mechanisms have remained elusive. By merging molecular beam experiments with electronic structure calculations, we reveal a complex chain of reactions initiated through the gas-phase preparation of benzocorannulene (C 24 H 12 ) via ring annulation of the corannulenyl radical (C 20 H 9 • ) by vinylacetylene (C 4 H 4 ) as identified isomer-selectively in situ via photoionization efficiency curves and photoion mass-selected threshold photoelectron spectra. In silico studies provided compelling evidence that the benzannulation mechanism can be expanded to pentabenzocorannulene (C 40 H 20 ) followed by successive cyclodehydrogenation to the C40 nanobowl (C 40 H 10 ) – a fundamental building block of buckminsterfullerene (C 60 ). This high-temperature pathway opens up isomer-selective routes to nanobowls via resonantly stabilized free-radical intermediates and ring annulation in circumstellar envelopes of carbon stars and planetary nebulae as their descendants eventually altering our insights of the complex chemistry of carbon in our Galaxy.

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

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2005

Optimal Sampling of a Reaction Coordinate in Molecular Dynamics

Estimating how free energy changes with the state of a system is a central goal in applications of statistical mechanics to problems of chemical or biological interest. From these free energy changes it is possible, for example, to establish which states of the system are stable, what are their probabilities and how the equilibria between these states are influenced by external conditions. Free energies are also of great utility in determining kinetics of transitions between different states. A variety of methods have been developed to compute free energies of condensed phase systems. Here, I will focus on one class of methods - those that allow for calculating free energy changes along one or several generalized coordinates in the system, often called reaction coordinates or order parameters . Considering that in almost all cases of practical interest a significant computational effort is required to determine free energy changes along such coordinates it is hardly surprising that efficiencies of different methods are of great concern. In most cases, the main difficulty is associated with its shape along the reaction coordinate. If the free energy changes markedly along this coordinate Boltzmann sampling of its different values becomes highly non-uniform. This, in turn, may have considerable, detrimental effect on the performance of many methods for calculating free energies.

Pohorille, Andrew↗