Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab Initio Kinetics”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Radical–Radical Reactions in Molecular Weight Growth: The Phenyl + Propargyl Reaction

The mechanism for hydrocarbon ring growth in sooting environments is still the subject of considerable debate. The reaction of phenyl radical (C 6 H 5 ) with propargyl radical (H 2 CCCH) provides an important prototype for radical–radical ring-growth pathways. We studied this reaction experimentally over the temperature range of 300–1000 K and pressure range of 4–10 Torr using time-resolved multiplexed photoionization mass spectrometry. We detect both the C 9 H 8 and C 9 H 7 + H product channels and report experimental isomer-resolved product branching fractions for the C 9 H 8 product. We compare these experiments to theoretical kinetics predictions from a recently published study augmented by new calculations. Here, these ab initio transition state theory-based master equation calculations employ high-quality potential energy surfaces, conventional transition state theory for the tight transition states, and direct CASPT2-based variable reaction coordinate transition state theory (VRC-TST) for the barrierless channels. At 300 K only the direct adducts from radical–radical addition are observed, with good agreement between experimental and theoretical branching fractions, supporting the VRC-TST calculations of the barrierless entrance channel. As the temperature is increased to 1000 K we observe two additional isomers, including indene, a two-ring polycyclic aromatic hydrocarbon, and a small amount of bimolecular products C 9 H 7 + H. Our calculated branching fractions for the phenyl + propargyl reaction predict significantly less indene than observed experimentally. We present further calculations and experimental evidence that the most likely cause of this discrepancy is the contribution of H atom reactions, both H + indenyl (C 9 H 7 ) recombination to indene and H-assisted isomerization that converts less stable C 9 H 8 isomers into indene. Especially at low pressures typical of laboratory investigations, H-atom-assisted isomerization needs to be considered. Regardless, the experimental observation of indene demonstrates that the title reaction leads, either directly or indirectly, to the formation of the second ring in polycyclic aromatic hydrocarbons.

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

Computing Reaction Kinetics with MC-PDFT–OPESf: Combining Multireference Electronic Structure Theory and Enhanced Sampling

Accurate rate constants are crucial for understanding and optimizing catalytic reactions mediated by enzymes, metalloproteins, and heterogeneous catalysts. These systems frequently present a dual computational challenge. Multiconfigurational reaction sites require multireference techniques for the accurate treatment of the electronic structure, and high activation barriers prevent efficient sampling of unbiased reactive transitions. In this work, we combine multiconfiguration pair-density functional theory (MC-PDFT) as an accurate and efficient multireference electronic structure method with on-the-fly probability-enhanced sampling flooding (OPESf) as an enhanced sampling method capable of accelerating reactive transitions. We demonstrate the approach on the Diels–Alder [4+2] cycloaddition between cis-butadiene and ethene as a reaction characterized by a large activation barrier and multireference character. MC-PDFT–OPESf provides reaction rates in agreement with experiments at a fraction of the computational cost required by conventional unbiased ab initio calculations. Here, we propose MC-PDFT–OPESf as an efficient approach for computing kinetics in strongly correlated molecular systems.

Chemical calculations↗
DOE OSTI2024

Quantum Molecular Charge-Transfer Model for Multistep Auger–Meitner Decay Cascade Dynamics

The fragmentation of molecular cations following inner-shell decay processes in molecules containing heavy elements underpins the X-ray damage effects observed in X-ray scattering measurements of biological and chemical materials, as well as in medical applications involving Auger electron-emitting radionuclides. Traditionally, these processes are modeled using simulations that describe the electronic structure at an atomic level, thereby omitting molecular bonding effects. This work addresses the gap by introducing a novel approach that couples Auger–Meitner decay to nuclear dynamics across multiple decay steps, by developing a decay spawning dynamics algorithm and applying it to potential energy surfaces characterized with ab initio molecular dynamics simulations. We showcase the approach on a model decay cascade following K-shell ionization of IBr and subsequent Kβ fluorescence decay. We examine two competing channels that undergo two decay steps, resulting in ion pairs with a total 3+ charge state. This approach provides a continuous description of the electron transfer dynamics occurring during the multistep decay cascade and molecular fragmentation, revealing the combined inner-shell decay and charge transfer time scale to be approximately 75 fs. In conclusion, our computed kinetic energies of ion fragments show good agreement with experimental data.

Ab initio molecular dynamics↗
DOE OSTI2021

The Reaction of o-Benzyne with Vinylacetylene: An Unexplored Way to Produce Naphthalene

Herein, the mechanisms and kinetics of the reaction of ortho-benzyne with vinylactylene have been studied by ab initio and density functional CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G(d,p) calculations of the pertinent potential energy surface combined with Rice-Ramsperger-Kassel-Marcus - Master Equation calculations of reaction rate constants at various temperatures and pressures. Under prevailing combustion conditions, the reaction has been shown to predominantly proceed by the biradical acetylenic mechanism initiated by the addition of C4H4 to one of the C atoms of the triple bond in ortho-benzyne by the acetylenic end, with a significant contribution of the concerted addition mechanism. Following the initial reaction steps, an extra six-membered ring is produced and the rearrangement of H atoms in this new ring leads to the formation of naphthalene, which can further dissociate to 1- or 2-naphthyl radicals. The o-C 6 H 4 + C 4 H 4 reaction is highly exothermic, by ~143 kcal/mol to form naphthalene and by 31-32 kcal/mol to produce naphthyl radicals + H, but features relatively high entrance barriers of 9-11 kcal/mol. Although the reaction is rather slow, much slower than the reaction of phenyl radical with vinylacetylene, it forms naphthalene and 1- and 2-naphthyl radicals directly, with their relative yields controlled by the temperature and pressure, and thus represents a viable source of the naphthalene core under conditions where ortho-benzyne and vinylacetylene are available.

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

Enhanced Menshutkin S N 2 Reactivity in Mesoporous Silica: The Influence of Surface Catalysis and Confinement

A significant enhancement in the Menshutkin S N 2 reaction between 1-methylimidazole (MeIm) and methyl thiocyanate (MeSCN) is observed when the reaction is confined in the nanoscale silica pores of MCM41 and SBA15. The experiments in the silica pores are conducted without surrounding bulk reaction mixture. The influences of temperature, pore radius, and surface chemistry on the kinetics of the confined reaction are analyzed with time-dependent infrared spectroscopy, molecular dynamics simulations, and ab initio calculations. The rate constant of the pseudo-first order reaction increases with decreasing pore size, and the activation energy is found to decrease by 5.6 kJ/mol in the smallest pore studied (2.8 nm) relative to the bulk reaction. The rate constant dependence on pore size is accurately described by a two-state model in which molecules within the 4.6 Å interfacial layer experience a 2.4-fold rate constant increase relative to those reacting at the bulk rate further away from the interface. The removal of polar silanol groups from the silica surface via passivation with trimethylsilyl chloride results in bulk-like kinetics despite a reduction in the pore diameter, demonstrating the role of silanols as catalytic sites. Here, electronic structure calculations of the energy profile on a model silica surface confirm that silanol groups, particularly those of the vicinal type, can reduce the activation energy and reaction endothermicity through the donation of hydrogen bonds to the reactant, transition state, and product complexes.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

First-principle investigation on catalytic hydrogenation of benzaldehyde over Pt-group metals

We report that understanding the hydrogenation of organic compounds in the aqueous phase has always been fundamentally important for improving carbon neutral pathways to fuels and value-added chemicals. In this study, we investigated both thermodynamic and kinetic profiles of benzaldehyde hydrogenation over the Pd(111) and Pt(111) metal surfaces using density functional theory (DFT) and ab initio molecular dynamic (AIMD) simulations. The adsorption of H 2 shows the mixed preference of H adsorption sites on the Pt(111), while the fcc adsorption site is dominant for H on the Pd(111). When benzaldehyde is added to the systems, we observe a strong reduction of benzaldehyde on charged Pd (111) surface compared with that on neutral surface. In contrast, charged state of the Pt(111) surface does not change their interaction. Subsequent hydrogenation reaction of benzaldehyde over Pd(111), proceeding via Langmuir-Hinshelwood mechanism, is affected by two major factors: the presence of H 2 O solvent and surface charge. The presence of H 2 O solvent greatly reduces the activation energy of C-H and O-H bond formation during the hydrogenation process. Furthermore, the hydrogenation step via C-H bond formation is preferred thermodynamically and kinetically over O-H bond formation during thermocatalytic hydrogenation, while the opposite trend holds true during electrocatalytic hydrogenation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Nuclear quantum effects of metal surface-mediated C–H activation

The nuclear quantum effects of surface-mediated C–H activation of surface CH 3 are considered for the pristine Pt(111) and Au(111) surfaces at 300 K. The kinetic barriers without nuclear quantum effects are calculated using both static density functional theory calculations and ab initio molecular dynamics. Static calculations are performed using the harmonic approximation while the free energy pathway is calculated using enhanced sampling molecular dynamics. Machine learning potentials are trained using generated datasets and validated against the ab initio molecular dynamics generated free energy pathways. The machine learning potentials are used to perform centroid molecular dynamics to consider the nuclear quantum effects of C–H activation. Nuclear quantum effects are found to have a very significant effect on the free energy pathway, with reduced importance at higher temperatures and in the CD 3 case.

Bunting, Rhys J. [Lawrence Livermore National Labo↗
DOE OSTI2020

Electronic correlations and flattened band in magnetic Weyl semimetal candidate Co 3 Sn 2 S 2

The interplay between electronic correlations and topological protection may offer a rich avenue for discovering emergent quantum phenomena in condensed matter. However, electronic correlations have so far been little investigated in Weyl semimetals (WSMs) by experiments. Here, we report a combined optical spectroscopy and theoretical calculation study on the strength and effect of electronic correlations in a magnet Co 3 Sn 2 S 2 . The electronic kinetic energy estimated from our optical data is about half of that obtained from single-particle ab initio calculations in the ferromagnetic ground state, which indicates intermediate-strength electronic correlations in this system. Furthermore, comparing the energy and side-slope ratios between the interband-transition peaks at high energies in the experimental and single-particle-calculation-derived optical conductivity spectra with the bandwidth-renormalization factors obtained by many-body calculations enables us to estimate the Coulomb-interaction strength (U ~ 4 eV) in Co 3 Sn 2 S 2 . Besides, a sharp experimental optical conductivity peak at low energy, which is absent in the single-particle-calculation-derived spectrum but is consistent with the optical conductivity peaks obtained by many-body calculations with U ~ 4 eV, indicates that an electronic band connecting the two Weyl cones is flattened by electronic correlations and emerges near the Fermi energy in Co 3 Sn 2 S 2 . Our work paves the way for exploring flat-band-generated quantum phenomena in WSMs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

SAXS analysis of the intrinsic tenase complex bound to a lipid nanodisc highlights intermolecular contacts between factors VIIIa/IXa

Abstract The intrinsic tenase (Xase) complex, formed by factors (f) VIIIa and fIXa, forms on activated platelet surfaces and catalyzes the activation of factor X to Xa, stimulating thrombin production in the blood coagulation cascade. The structural organization of the membrane-bound Xase complex remains largely unknown, hindering our understanding of the structural underpinnings that guide Xase complex assembly. Here, we aimed to characterize the Xase complex bound to a lipid nanodisc with biolayer interferometry (BLI), Michaelis–Menten kinetics, and small-angle X-ray scattering (SAXS). Using immobilized lipid nanodiscs, we measured binding rates and nanomolar affinities for fVIIIa, fIXa, and the Xase complex. Enzyme kinetic measurements demonstrated the assembly of an active enzyme complex in the presence of lipid nanodiscs. An ab initio molecular envelope of the nanodisc-bound Xase complex allowed us to computationally model fVIIIa and fIXa docked onto a flexible lipid membrane and identify protein–protein interactions. Our results highlight multiple points of contact between fVIIIa and fIXa, including a novel interaction with fIXa at the fVIIIa A1–A3 domain interface. Lastly, we identified hemophilia A/B-related mutations with varying severities at the fVIIIa/fIXa interface that may regulate Xase complex assembly. Together, our results support the use of SAXS as an emergent tool to investigate the membrane-bound Xase complex and illustrate how mutations at the fVIIIa/fIXa dimer interface may disrupt or stabilize the activated enzyme complex.

Hematology↗
DOE OSTI2019

A Review on Lignin Liquefaction: Advanced Characterization of Structure and Microkinetic Modeling

Lignin liquefaction microkinetics is a move toward a more first principles (i.e., ab initio) based understanding at the molecular level in reaction engineering. While the microkinetic modelling of reactions to obtain kinetic rate parameters of chemical reactions have been widely used in the field of gas phase combustion and heterogeneous catalysis, this approach has not been as thoroughly developed in the area of biomass thermochemical reactions (e.g., lignin pyrolysis, hydrothermal liquefaction). The difficulties in establishing the structure of complex heterogeneous materials, like lignin, is perhaps the main challenge in developing rational microkinetic descriptions of biomass thermochemical reactions. In this manuscript we review the current state of the art and the challenges to develop microkinetic models for lignin liquefaction technologies (e.g., pyrolysis, hydrothermal liquefaction, solvolysis). A general strategy for the development of microkinetic models for lignin liquefaction technologies is discussed. Here, the first hurdle is to obtain sufficiently rich experimental data of lignin underlying polymeric structure and methodologies to use this data to build realistic lignin structural representations. Some analytical techniques for lignin structural characterization and their associated data, as well as a correlation for calculating the degree of macromolecular lignin branching are discussed. The presence of small lignin oligomeric structures and the role of these structures in lignin pyrolysis is also addressed. The ways in which elementary deconstruction and repolymerization reactions occur within this structure to form a liquid intermediate, and how these deconstruction products continue to interact with each other until they are removed from the liquid intermediate is thoroughly discussed. Further, experimental work with model compounds and the effect of reaction parameters (e.g., temperature, pressure, vapor residence time) are reviewed.

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

Rehybridization dynamics into the pericyclic minimum of an electrocyclic reaction imaged in real-time

Abstract Electrocyclic reactions are characterized by the concerted formation and cleavage of both σ and π bonds through a cyclic structure. This structure is known as a pericyclic transition state for thermal reactions and a pericyclic minimum in the excited state for photochemical reactions. However, the structure of the pericyclic geometry has yet to be observed experimentally. We use a combination of ultrafast electron diffraction and excited state wavepacket simulations to image structural dynamics through the pericyclic minimum of a photochemical electrocyclic ring-opening reaction in the molecule α-terpinene. The structural motion into the pericyclic minimum is dominated by rehybridization of two carbon atoms, which is required for the transformation from two to three conjugated π bonds. The σ bond dissociation largely happens after internal conversion from the pericyclic minimum to the electronic ground state. These findings may be transferrable to electrocyclic reactions in general.

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

Theoretical and kinetic modeling study of hydrazine oxidation

The present work constitutes the first theoretical and kinetic modeling study of hydrazine oxidation, which may be important for burnout in ammonia-fueled combustion. The kinetics of the oxidation of N 2 H 4 , N 2 H 3 and tHNNH by molecular oxygen were investigated via a quantum chemistry/canonical transition state theory approach. Geometries and anharmonic frequencies were obtained with density functional theory, and energies from coupled cluster calculations (CCSD(T)) extrapolated to the infinite basis set limit, with corrections for core-valence electron correlation, scalar relativistic effects, and higher level correlation up to lambda coupled cluster, CCSDT(Q) Λ . The key reactions occurred on the N 2 H 4 O 2 potential energy surface, where the results indicated a fast reaction of N 2 H 3 with HO 2 via singlet adducts to yield tHNNH + H 2 O 2 and HNN(H)O + H 2 O, while reaction on the triplet surface proceeds via a bound complex followed by a tight, submerged barrier to yield N 2 H 4 + O 2 . The results were incorporated in a detailed reaction mechanism, which was used to interpret the shock tube results from Michel and Wagner (1965) on the effect of O 2 on hydrazine conversion at 1100–1400 K. The kinetic model captured qualitatively the observed behavior, but underestimated the reaction rate under oxidizing conditions. The hydrazine pyrolysis chemistry dominated conversion at reducing conditions and/or high temperature. At oxidizing conditions and intermediate temperatures (≲ 1400 K), reactions of N 2 -amines with HO 2 and O 2 were important for the oxidation rate.

Ab initio calculations↗
DOE OSTI2025

Theoretical and kinetic modeling study of H 2 S pyrolysis

Hydrogen sulfide pyrolysis was investigated theoretically and through chemical kinetic modeling. Reactions on the SHH potential energy surface, primarily S + H 2 (+Ar) ⇌ H 2 S (+ Ar) (R1) and S + H 2 ⇌ SH + H (R6b) were characterized by ab initio calculations. Results for k 1 were in good agreement with experiment, but deviated strongly below 2000 K from values previously used in modeling. Collider efficiencies for H 2 S, S 2 , and N 2 compared to Ar were calculated for R1. Hydrogen sulfide decomposition experiments reported in literature were re-examined in terms of an updated detailed chemical kinetic model. Concentration profiles for the atomic S at high temperature in shock tubes supported the present value of k 1 and served to constrain the rate constants for reaction of S with SH and H 2 S. To explain results from batch and flow reactors, conducted at high H 2 S concentrations in the 900–1400 K range, a very fast rate constant was required for HSS + H ⇌ SH + SH. Under dilute conditions, the gas-phase chemistry was too slow to compete and the decomposition of H 2 S was controlled by loss on the reactor surface.

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

Toward Higher Voltage Solid-State Batteries by Metastability and Kinetic Stability Design

The energy density of battery systems is limited largely by the electrochemical window of the electrolyte. Here, the combined thermodynamic and kinetic effects of mechanically induced metastability are shown to greatly widen the operational voltage window of solid-state batteries based on ceramic-sulfide electrolytes. Solid electrolyte voltage stability up to 10 V is achieved with minimal degradation, far beyond the capability of organic liquid electrolytes. Furthermore, combined experiment, ab initio computation, and theoretical modeling identify the nature of mechanically constrained Li 10 GeP 2 S 12 decomposition both within the bulk and at interfaces with cathode materials at very high voltages. Previously unclear kinetic processes are identified that, when properly implemented, can potentially allow solid-state full cells with remarkably high operational voltages.

25 ENERGY STORAGE↗