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

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2024

Mechanistic study of direct coupling of CO 2 and C 2 H 4 over atomically dispersed metal at graphene edges

Direct coupling of CO 2 and ethylene (hereinafter DCCE) to acrylic acid is valuable for valorizing CO 2 to manufacture acrylate-derived products. However, previous studies in DCCE have been limited on molecular catalysts with challenges in improving catalytic performance. In this work, we employed density functional theory calculations and ab initio molecular dynamics simulations to investigate the heterogeneous catalysis of DCCE over atomically dispersed metal centers at nitrogen-doped zigzag edge of graphene. Based on competitive adsorption and structural stability, Mo, Cr, V, Ru, and Ni active sites are chosen to explore the reaction kinetics. Here, we find that the activation barriers are determined by the charge redistribution at transition states, which explains the trend of activity for the C-C coupling and the hydrogen transfer, two key steps in DCCE. Furthermore, we show that the intramolecular hydrogen transfer (rate-limiting step) is hindered due to the lack of local coordinate at the active sites. We thus propose to use co-adsorbed water as a “proton-exchanger” following a water-assisted route, and show that the activation barriers are reduced over all metal centers. Particularly, water promotes the hydrogen transfer over metals with strong CO 2 -ethylene co-activation and facile C-C coupling kinetics, which could be considered promising for DCCE. In both mechanisms, the stability of metallactone intermediate can be used to predict the catalytic activity. It is anticipated that the insights from this work can provide guidelines for mimicking well-defined multifunctional active sites in molecular catalysts to design heterogeneous catalysts for such C-C coupling, which advances catalytic utilization of CO 2 .

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

Towards accurate orbital-free simulations: A generalized gradient approximation for the noninteracting free energy density functional

For orbital-free ab initio molecular dynamics, especially on systems in extreme thermodynamic conditions, we provide the first pseudo-potential-adapted generalized gradient approximation (GGA) functional for the non-interacting free energy. This is achieved by systematic finite-temperature extension of our recent LKT ground state non-interacting kinetic energy GGA functional (Phys. Rev. B 98, 041111(R) (2018)). We test the performance of the new functional first via static lattice calculations on crystalline aluminum and silicon. Then we compare deuterium equation of state results against both path-integral Monte Carlo and conventional (orbital-dependent) Kohn-Sham results. The new functional, denoted LKTF, outperforms the previous best semi-local free energy functional, VT84F (Phys. Rev. B 88, 161108(R) (2013)), and provides modestly faster simulations. Finally, we discuss subtleties of identification of kinetic and entropic contributions to non-interacting free-energy functionals obtained by extension from ground state orbital-free kinetic energy functionals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗
DOE OSTI2023

Case Study: Calorimetry

This chapter presents the use of adsorption microcalorimetry to study acid-base catalysis over mixed metal oxides. Understanding the acid-base properties of catalysts is of paramount importance for steering product selectivity in a wide variety of catalytic industrial processes. In particular, catalytic conversion of biomass into biofuels and chemicals demands comprehensive understanding of the acid-base properties of the catalyst to selectively convert functional groups. Herein, the use of probe molecules to measure the strength and density of acid and basic sites on mixed metal oxide surfaces is introduced. Further, the challenges when characterizing mixed metal oxides, presenting surface reconstruction, are explained. An example displays how combining adsorption calorimetry with ab initio calculations explains the activation of a reactant on the surface of a mixed metal oxide. Acid-base catalysis over nanoshapes is studied to unveil that different nanoshapes lead to different extents of surface reconstruction, which ultimately affect acid-base catalysis. The use of steady-state isotopic transient kinetic analysis (SSITKA), in combination with adsorption microcalorimetry, unveils that surface reconstruction of a set of perovskites impacts the density of active sites to a greater extent than the intrinsic reactivity of the active site. At the end of this chapter, fundamental insights into catalyst transient behavior during isomerization of n-butane, and the reaction mechanism for conversion of 4-methylpentan-2-ol, are presented.

Polo Garzon, Felipe↗
DOE OSTI2023

Nickel Isonicotinate Framework with Optimal Pore Structure for Complete Discrimination of Hexane Isomers

Efficient separation of physicochemically similar alkanes is of vital importance. Adsorptive separation utilizing porous materials such as metal–organic frameworks with tunable pore structure and surface functionality represents an energy-efficient technology. In this study, we demonstrate successful separation of alkanes with varying degree of branching using a microporous nickel isonicotinate framework, Ni(4-PyC) 2 . Its 2-fold interpenetrated diamondoid structure with well-suited pore size enables selective adsorption of linear and monobranched hexane isomers, while excluding dibranched isomer. Breakthrough experiments validated its capability to completely discriminate all three hexane isomers. Ab initio calculations combined with in situ infrared spectroscopic analysis unveiled the nature of host–guest interactions and differences in the binding energies and diffusion barriers among the isomers. Furthermore, having well-balanced adsorption uptakes (146 and 79 mg g –1 of nHEX and 3MP), high nHEX/DMB uptake ratio (12.2) and fast kinetics, Ni(4-PyC) 2 stands out as a promising adsorbent for complete separation of hexane isomers under ambient conditions.

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

Integrated Physics-based Modeling and Experiments for Improved Prediction of Combustion Dynamics in Low-Emission Systems

Concurrent simulation and experiment was undertaken to assess the ability of a hybrid RANS-LES model to predict combustion dynamics in a single-element lean direct-inject (LDI) combustor showing self-excited instabilities. High frequency pressure modes produced by Fourier and modal decomposition analysis were compared quantitatively, and trends with equivalence ratio and inlet temperature were compared qualitatively. High frequency OH PLIF and PIV measurements were also taken. Submodels for chemical kinetics and primary and secondary atomization were also tested against the measured behavior. For a point-wise comparison, the amplitudes matched within a factor of two. The dependence on equivalence ratio was matched. Preliminary results from simulation using an 18-reaction kinetics model indicated instability amplitudes closer to measurement. Analysis of the simulations suggested a band of modes around 1400 Hz were due to a vortex bubble breakdown and a band of modes around 6 kHz were due to a precessing vortex core hydrodynamic instability. The primary needs are directly coupled and validated ab initio models of the atomizer free surface flow and the primary atomization processes, and more detailed study of the coupling between the 3D swirling flow and the local thermoacoustics in the diverging venturi section.

combustion↗
DOE OSTI2023

Understanding palladium–tellurium cluster formation on WTe2: From a kinetically hindered distribution to thermodynamically controlled monodispersity

Abstract A fundamental understanding of the transition metal dichalcogenide (TMDC)–metal interface is critical for their utilization in a broad range of applications. We investigate how the deposition of palladium (Pd), as a model metal, on WTe2(001), leads to the assembly of Pd into clusters and nanoparticles. Using X-ray photoemission spectroscopy, scanning tunneling microscopy imaging, and ab initio simulations, we find that Pd nucleation is driven by the interaction with and the availability of mobile excess tellurium (Te) leading to the formation of Pd-Te clusters at room temperature. Surprisingly, the nucleation of Pd-Te clusters is not affected by intrinsic surface defects, even at elevated temperatures. Upon annealing, the Pd-Te nanoclusters adopt an identical nanostructure and are stable up to ∼523 K. Density functional theory calculations provide a foundation for our understanding of the mobility of Pd and Te atoms, preferential nucleation of Pd-Te clusters, and the origin of their annealing-induced monodispersity. These results highlight the role the excess chalcogenide atoms may play in the metal deposition process. More broadly, the discoveries of synthetic pathways yielding thermally robust monodispersed nanostructures on TMDCs are critical to the manufacturing of novel quantum and microelectronics devices and catalytically active nano-alloy centers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Solvent structure and dynamics over Brønsted acid MWW zeolite nanosheets

In the liquid phase of heterogeneous catalysis, solvent plays an important role and governs the kinetics and thermodynamics of a reaction. Although it is often difficult to quantify the role of the solvent, it becomes particularly challenging when a zeolite is used as the catalyst. This difficulty arises from the complex nature of the liquid/zeolite interface and the different solvation environments around catalytically active sites. Here, we use ab initio molecular dynamics simulations to probe the local solvation structure and dynamics of methanol and water over MWW zeolite nanosheets with varying Brønsted acidity. We find that the zeolite framework and the number and location of the acid sites in the zeolite influence the structure and dynamics of the solvent. In particular, methanol is more likely to be in the vicinity of the aluminum (Al 3+ ) at the T4 site than at T1 due to easy accessibility. The methanol oxygen binds strongly to the Al at the T4 site, weakening the Al–O for the bridging acid site, which results in the formation of the silanol group, significantly reducing the acidity of the site. The behavior of methanol is in direct contrast to that of water, where protons can easily propagate from the zeolite to the solvent molecules regardless of the acid site location. Furthermore, our work provides molecular-level insights into how solvent interacts with zeolite surfaces, leading to an improved understanding of the catalytic site in the MWW zeolite nanosheet.

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

An Accurate Potential Energy Surface for H2O

We have carried out extensive high quality ab initio electronic structure calculations of the ground state potential energy surface (PES) and dipole moment function (DMF) for H2O. A small adjustment is made to the PES to improve the agreement of line positions from theory and experiment. The theoretical line positions are obtained from variational ro-vibrational calculations using the exact kinetic energy operator. For the lines being fitted, the root-mean-square error was reduced from 6.9 to 0.08 /cm. We were then able to match 30,092 of the 30,117 lines from the HITRAN 96 data base to theoretical lines, and 80% of the line positions differed less than 0.1 /cm. About 3% of the line positions in the experimental data base appear to be incorrect. Theory predicts the existence of many additional weak lines with intensities above the cutoff used in the data base. To obtain results of similar accuracy for HDO, a mass dependent correction to the PH is introduced and is parameterized by simultaneously fitting line positions for HDO and D2O. The mass dependent PH has good predictive value for T2O and HTO. Nonadiabatic effects are not explicitly included. Line strengths for vibrational bands summed over rotational levels usually agree well between theory and experiment, but individual line strengths can differ greatly. A high temperature line list containing about 380 million lines has been generated using the present PES and DMF

Schwenke, David W.↗
DOE OSTI2023

Energetics of high temperature degradation of fentanyl into primary and secondary products

Fentanyl is a synthetic opioid used for managing chronic pain. Due to its higher potency (50–100×) than morphine, fentanyl is also an abused drug. A sensor that could detect illicit fentanyl by identifying its thermally degraded fragments would be helpful to law enforcement. While experimental studies have probed the thermal degradation of fentanyl, little theoretical work has been done to understand the mechanism. Here, we studied the thermal degradation pathways of fentanyl using extensive ab initio molecular dynamics simulations combined with enhanced sampling via multiple-walker metadynamics. We calculated the free energy profile for each bond suggested earlier as a potential degradation point to map the thermodynamic driving forces. In conclusion, we also estimated the forward attempt rate of each bond degradation reaction to gain information about degradation kinetics.

Poudel, Bharat↗
DOE OSTI2022

Understanding Hydrogenation Chemistry at MgB 2 Reactive Edges from Ab Initio Molecular Dynamics

Solid-state hydrogen storage materials often operate via transient, multistep chemical reactions at complex interfaces that are difficult to capture. Here, we use direct ab initio molecular dynamics simulations at accelerated temperatures and hydrogen pressures to probe the hydrogenation chemistry of the candidate material MgB 2 without a priori assumption of reaction pathways. Focusing on highly reactive ($10\overline{1}0$) edge planes where initial hydrogen attack is likely to occur, we track mechanistic steps toward the formation of hydrogen-saturated BH 4 – units and key chemical intermediates, involving H 2 dissociation, generation of functionalities and molecular complexes containing BH 2 and BH 3 motifs, and B–B bond breaking. Further, the genesis of higher-order boron clustering is also observed. Different charge states and chemical environments at the B-rich and Mg-rich edge planes are found to produce different chemical pathways and preferred speciation, with implications for overall hydrogenation kinetics. The reaction processes rely on B–H bond polarization and fluctuations between ionic and covalent character, which are critically enabled by the presence of Mg 2+ cations in the nearby interphase region. Our results provide guidance for devising kinetic improvement strategies for MgB 2 -based hydrogen storage materials, while also providing a template for exploring chemical pathways in other solid-state energy storage reactions.

08 HYDROGEN↗
DOE OSTI2020

Conformational Panorama and Chirality Controlled Structure–Energy Relationship in a Chiral Carboxylic Acid Dimer

Chirality recognition in dimers of tetrahydro-2-furoic acid (THFA) was studied in a conformer-specific manner using rotational spectroscopy and theoretical approaches. THFA shows a strong preference for the trans- over the cis-COOH configuration. Two drastically different scenarios are possible for the detectable (THFA) 2 : a kinetically preferred dimer bound by feeble interactions between two trans-THFAs or a thermodynamically favored dimer with a double hydrogen-bonded ring structure between two cis subunits. To identify the conformers responsible for the extremely dense rotational spectra observed, it was essential not only to locate several hundred homo/heterochiral (THFA) 2 minima in ab initio calculations but also to evaluate the energetic connectivities among the minima. The study further reveals an interesting chirality dependent structure-energy ordering relationship. A method for enantiomeric excess (ee) determination of THFA is presented using a recently proposed chiral self-tag approach.

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

Oxidation of methylamine

A detailed chemical kinetic model for oxidation of methylamine has been developed, based on theoretical work and a critical evaluation of data from the literature. In this work, the rate coefficients for the reactions of CH 3 NH 2 + O 2 → CH 2 NH 2 / CH 3 NH + HO 2 , CH 3 NH 2 + H → CH 3 + NH 3 , CH 3 NH → CH 2 NH 2 , and CH 3 NH + O 2 → CH 2 NH + HO 2 were calculated from ab initio theory. The mechanism was validated against experimental results from batch reactors, flow reactors, shock tubes, and premixed flames. The model predicts satisfactorily explosion limits for CH 3 NH 2 and its oxidation in a flow reactor. However, oxidation in the presence of nitric oxide, which strongly promotes reaction at lower temperatures, is only described qualitatively. Furthermore, calculated flame speeds are higher than reported experimental values; the model does not capture the inhibiting effect of the NH 2 group in CH 3 NH 2 compared to CH 4 . More work is desirable to confirm the products of the CH 3 NH + NO reaction and to look into possible pathways to NH 3 in methylamine oxidation.

10 SYNTHETIC FUELS↗
DOE OSTI2023

Hydrogen in disordered titania: connecting local chemistry, structure, and stoichiometry through accelerated exploration

Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Separating Xylene Isomers with a Calcium Metal‐Organic Framework

Abstract The purification of p‐xylene (pX) from its xylene isomers represents a challenging but important industrial process. Herein, we report the efficient separation of pX from its ortho‐ and meta‐ isomers by a microporous calcium‐based metal–organic framework material (HIAM‐203) with a flexible skeleton. At 30 °C, all three isomers are accommodated but the adsorption kinetics of o‐xylene (oX) and m‐xylene (mX) are substantially slower than that of pX, and at an elevated temperature of 120 °C, oX and mX are fully excluded while pX can be adsorbed. Multicomponent column breakthrough measurements and vapor‐phase/liquid‐phase adsorption experiments have demonstrated the capability of HIAM‐203 for efficient separation of xylene isomers. Ab initio calculations have provided useful information for understanding the adsorption mechanism.

Yu, Liang↗
DOE OSTI2025

New Dimension in Ab Initio Electronic Structure Theory: Temperature, Pressure, and Chemical Potential

Ab initio electronic structure theory has transformed gas-phase molecular science with its predictive ability. In the attempt to bring such predictive ability to macroscopic systems and condensed matter, the theory must integrate quantum mechanics with statistical thermodynamics, so that thermodynamic functions such as free energy, internal energy, entropy, and chemical potentials are computed as functions of temperature in a systematically converging series of approximations. Here, a general, versatile strategy of elevating ab initio electronic structure theory to nonzero temperatures is introduced and discussed.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Data for "Plasmon-driven exciton formation in a non-equilibrium Fermi liquid"

This repository contains source data for key plots presented in the manuscript "Plasmon-driven exciton formation in a non-equilibrium Fermi liquid." Experimental data that was analyzed in Igor Pro 8 are presented as the .pxp files used to generate individual sub-plots. Electronic spectral function calculations are provided as .txt files, in which consecutive rows refer to the meshgrid x coordinate, y coordinate, spectral function (and, where relevant, axis-projected local angular momentum). We additionally include the Wannier model and DFT-obtained bulk band structure on which the Wannier model was based. Files are named as the number of the figure in the manuscript to which they correspond, with additional details included where necessary. Details of file names: 2a_DOS_Lxz_Ek_KGM_40layer_xnum_800kpt_tot.txt: Density of states, xz-axis projected local orbital angular momentum, for 800 points along the K-Gamma-M path, for a 40-layer model. 2c_composite_y.pxp: ARPES (angle-resolved photoemission spectroscopy) spectra along the ky axis, including both a scan near the Fermi level and a scan at high kinetic energies. 2d_LCP_RCP_diff_Sect_20K.pxp: difference between ARPES constant energy cuts at T=20 K at E0 + 0.23 eV taken with left- and right-circularly polarized photons. The polarization-integrated intensity at the constant energy cut is also included. 2e_DOS_L45_E11pt79_m0pt25to0pt25_xnum_800kpt_tot.txt: Density of states, xz-projected local orbital angular momentum, and corresponding k-points in two dimensions from ab-initio electronic structure calculations for a constant-energy cut. 3a_[x]_[y]ps: ARPES cut under excitation at a fluence of x uJ/cm2, measured y ps after photoexcitation. Measurements were performed at 9 K. 3b_[x]: Energy distribution curves under excitation at a fluence x uJ/cm2 at selected delay times after photoexcitation. 4a_ImSigma_vs_temperature.pxp: Imaginary self energy (extracted from ARPES linewidths) at different energies above E0 for selected lattice temperatures. 4b_EELS_lowE.pxp: Electron energy loss spectrum over a low energy range 5b_diff_55m15.pxp: Difference between momentum-integrated Tr-ARPES traces at 55 uJ/cm2 and 15 uJ/cm2 photoexcitation. Time-dependent intensity at each energy level has been normalized to a maximum of 1 for each individual fluence prior to subtraction. 5d_invtau_at_EX_vs_fluence.pxp: decay rate at a specified energy EX for different excitation fluences, from single exponential fits. NOTE: Analyses based on the Wannier model presented here should cite both the associated Article and this dataset. For all other files in the repository, citing the dataset alone is sufficient.

Acharya, Rishi [University of Illinois] (ORCID:000↗