Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Science Opportunities and Capabilities Enabled by LCLS-II and LCLS-II HE

X-Ray Science with LCLS-II and LCLS-II HE: Chemical dynamics: Reaction dynamics, charge transfer, molecular photocatalysts, natural & artificial photosynthesis; Catalysis: Homogeneous and heterogeneous catalysis, interfacial & geo/environmental chemistry; Materials Physics: Heterogeneity, spontaneous fluctuations, nonequilibrium dynamics, extreme environments; Quantum Materials: Emergent phenomena & collective excitations; Biological Function & Structural Dynamics: Dynamics in physiological environments.

43 PARTICLE ACCELERATORS↗
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 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 OSTI2019

Stochastic self-tuning hybrid algorithm for reaction-diffusion systems

Many biochemical phenomena involve reactants with vastly different concentrations, some of which are amenable to continuum-level descriptions, while the others are not. We present a hybrid self-tuning algorithm to model such systems. The method combines microscopic (Brownian) dynamics for diffusion with mesoscopic (Gillespie-type) methods for reactions and remains efficient in a wide range of regimes and scenarios with large variations of concentrations. Furthermore, its accuracy, robustness, and versatility are balanced by redefining propensities and optimizing the mesh size and time step. We use a bimolecular reaction to demonstrate the potential of our method in a broad spectrum of scenarios: from almost completely reaction-dominated systems to cases where reactions rarely occur or take place very slowly. The simulation results show that the number of particles present in the system does not degrade the performance of our method. This makes it an accurate and computationally efficient tool to model complex multireaction systems.

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

Physics-based, neural network force fields for reactive molecular dynamics: Investigation of carbene formation from [EMIM + ][OAc - ]

Reactive molecular dynamics simulations enable detailed understanding of solvent effects on chemical reaction mechanisms and reaction rates. While classical molecular dynamics using reactive force fields allows significantly longer simulation time scales and larger system sizes compared with ab initio molecular dynamics, constructing reactive force fields is a difficult and complex task. In this work, we describe a general approach following the Empirical Valence Bond (EVB) framework for constructing ab initio reactive force fields for condensed phase simulations by combining physics-based methods with neural networks (PB/NN). The physics-based terms ensure correct asymptotic behavior of electrostatic, polarization, and dispersion interactions, and are compatible with existing solvent force fields. Neural networks are utilized for versatile description of short-range orbital interactions within the transition state region, and accurate rendering of vibrational motion of the reacting complex. Herein, we demonstrate our methodology for a simple deprotonation reaction of the 1-ethyl-3-methylimidazolium (EMIM+) cation with acetate to form 1-ethyl- 3-methylimidazol-2-ylidene and acetic acid. Our PB/NN force field exhibits ~ 1 kJ/mol MAE accuracy within the transition state region for the gas-phase complex. To characterize solvent modulation of the reaction profile, we compute potentials of mean force (PMFs) for the gas-phase reaction as well as the reaction within a four ion cluster, and benchmark against ab initio molecular dynamics simulations. We find that the surrounding ionic environment significantly destabilizes formation of the carbene product, and we show that this effect is accurately captured by the reactive force field. By construction, the PB/NN potential may be directly employed for simulations of other solvents/chemical environments without additional parameterization.

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

Systematic development of reduced reaction mechanisms for dynamic modeling

A method for systematically developing a reduced chemical reaction mechanism for dynamic modeling of chemically reactive flows is presented. The method is based on the postulate that if a reduced reaction mechanism faithfully describes the time evolution of both thermal and chain reaction processes characteristic of a more complete mechanism, then the reduced mechanism will describe the chemical processes in a chemically reacting flow with approximately the same degree of accuracy. Here this postulate is tested by producing a series of mechanisms of reduced accuracy, which are derived from a full detailed mechanism for methane-oxygen combustion. These mechanisms were then tested in a series of reactive flow calculations in which a large-amplitude sinusoidal perturbation is applied to a system that is initially quiescent and whose temperature is high enough to start ignition processes. Comparison of the results for systems with and without convective flow show that this approach produces reduced mechanisms that are useful for calculations of explosions and detonations. Extensions and applicability to flames are discussed.

Frenklach, M.↗
DOE OSTI2021

Using Dissipative Particle Dynamics to Model Effects of Chemical Reactions Occurring within Hydrogels

Computational models that reveal the structural response of polymer gels to changing, dissolved reactive chemical species would provide useful information about dynamically evolving environments. However, it remains challenging to devise one computational approach that can capture all the interconnected chemical events and responsive structural changes involved in this multi-stage, multi-component process. Here, we augment the dissipative particle dynamics (DPD) method to simulate the reaction of a gel with diffusing, dissolved chemicals to form kinetically stable complexes, which in turn cause concentration-dependent deformation of the gel. Using this model, we also examine how the addition of new chemical stimuli and subsequent reactions cause the gel to exhibit additional concentration-dependent structural changes. Through these DPD simulations, we show that the gel forms multiple latent states (not just the “on/off”) that indicate changes in the chemical composition of the fluidic environment. Hence, the gel can actuate a range of motion within the system, not just movements corresponding to the equilibrated swollen or collapsed states. Moreover, the system can be used as a sensor, since the structure of the layer effectively indicates the presence of chemical stimuli.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

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

Real-space observation of the dissociation of a transition metal complex and its concurrent energy redistribution

Mechanistic insights into photodissociation dynamics of transition metal carbonyls, like Fe(CO) 5 , are fundamental for understanding active catalytic intermediates. Although extensively studied, the structural dynamics of these systems remain elusive. Using ultrafast X-ray scattering, we uncover the photochemistry of Fe(CO) 5 in real space and time, observing synchronous oscillations in atomic pair distances, followed by a prompt rotating CO release preferentially in the axial direction. This behavior aligns with simulations, reflecting the interplay between the axial Fe-C distances’ potential energy landscape and non-adiabatic transitions between metal-to-ligand charge-transfer states. Additionally, we characterize a secondary delayed CO release associated with a reduction of Fe-C steady state distances and structural dynamics of the formed Fe(CO) 4 . Our results quantify energy redistribution across vibration, rotation, and translation degrees of freedom, offering a microscopic view of complex structural dynamics, enhancing our grasp on Fe(CO) 5 photodissociation, and advancing our understanding of transition metal catalytic systems.

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

The quantum dynamics of electronically nonadiabatic chemical reactions

Considerable progress was achieved on the quantum mechanical treatment of electronically nonadiabatic collisions involving energy transfer and chemical reaction in the collision of an electronically excited atom with a molecule. In the first step, a new diabatic representation for the coupled potential energy surfaces was created. A two-state diabatic representation was developed which was designed to realistically reproduce the two lowest adiabatic states of the valence bond model and also to have the following three desirable features: (1) it is more economical to evaluate; (2) it is more portable; and (3) all spline fits are replaced by analytic functions. The new representation consists of a set of two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for dynamics calculations on both the electronic quenching and reaction processes in collisions of Na(3p2p) with H2. The new two-state representation was obtained by a three-step process from a modified eight-state diatomics-in-molecules (DIM) representation of Blais. The second step required the development of new dynamical methods. A formalism was developed for treating reactions with very general basis functions including electronically excited states. Our formalism is based on the generalized Newton, scattered wave, and outgoing wave variational principles that were used previously for reactive collisions on a single potential energy surface, and it incorporates three new features: (1) the basis functions include electronic degrees of freedom, as required to treat reactions involving electronic excitation and two or more coupled potential energy surfaces; (2) the primitive electronic basis is assumed to be diabatic, and it is not assumed that it diagonalizes the electronic Hamiltonian even asymptotically; and (3) contracted basis functions for vibrational-rotational-orbital degrees of freedom are included in a very general way, similar to previous prescriptions for locally adiabatic functions in various quantum scattering algorithms.

Truhlar, Donald G.↗
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 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↗
NASA NTRS2017

Modeling Chemical Reactions in Classical Molecular Dynamics Simulations

An algorithm capable of incorporating multi-step reaction mechanisms into atomistic molecular dynamics (MD) simulations using traditional fixed valence force fields is proposed and implemented within the framework of LAMMPS (Large-scale Atomic Molecular Massively Parallel Simulator). This extension, referred to as fix bond/react, enables bonding topology modifications during a running MD simulation using pre- and post-reaction bonding templates to carry out a pre-specified reaction. Candidate reactants are first identified by interatomic separation, followed by the application of a generalized topology matching algorithm to confirm they match the pre-reaction template. This is followed by a topology conversion to match the post-reaction template and a dynamic relaxation to minimize high energy configurations. Two case studies, the condensation polymerization of nylon 6,6 and the formation of a highly-crosslinked epoxy, are simulated to demonstrate the robustness, stability, and speed of the algorithm. Improvements which could increase its utility are discussed.

Jacob R Gissinger↗
NASA NTRS1995

Advanced Multi-phase Flow CFD Model Development for Solid Rocket Motor Flowfield Analysis

A Navier-Stokes code, finite difference Navier-Stokes (FDNS), is used to analyze the complicated internal flowfield of the SRM (solid rocket motor) to explore the impacts due to the effects of chemical reaction, particle dynamics, and slag accumulation on the solid rocket motor (SRM). The particulate multi-phase flowfield with chemical reaction, particle evaporation, combustion, breakup, and agglomeration models are included in present study to obtain a better understanding of the SRM design. Finite rate chemistry model is applied to simulate the chemical reaction effects. Hermsen correlation model is used for the combustion simulation. The evaporation model introduced by Spalding is utilized to include the heat transfer from the particulate phase to the gase phase due to the evaporation of the particles. A correlation of the minimum particle size for breakup expressed in terms of the Al/Al2O3 surface tension and shear force was employed to simulate the breakup of particles. It is assumed that the breakup occurs when the Weber number exceeds 6. A simple L agglomeration model is used to investigate the particle agglomeration. However, due to the large computer memory requirements for the agglomeration model, only 2D cases are tested with the agglomeration model. The VOF (Volume of Fluid) method is employed to simulate the slag buildup in the aft-end cavity of the redesigned solid rocket motor (RSRM). Monte Carlo method is employed to calculate the turbulent dispersion effect of the particles. The flowfield analysis obtained using the FDNS code in the present research with finite rate chemical reaction, particle evaporation, combustion, breakup, agglomeration, and VOG models will provide a design guide for the potential improvement of the SRM including the use of materials and the shape of nozzle geometry such that a better performance of the SRM can be achieved. The simulation of the slag buildup in the aft-end cavity can assist the designer to improve the design of the RSRM geometry.

Liaw, Paul↗