Search NASA⌕ Search

SEARCH · Search NASA

Results for “atomic 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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Resolving the Heat Generated from ZrO 2 Atomic Layer Deposition Surface Reactions

Abstract In situ pyroelectric calorimetry and spectroscopic ellipsometry were used to investigate surface reactions in atomic layer deposition (ALD) of zirconium oxide (ZrO 2 ). Calibrated and time‐resolved in situ ALD calorimetry provides new insights into the thermodynamics and kinetics of saturating surface reactions for tetrakis(dimethylamino)zirconium(IV) (TDMAZr) and water. The net ALD reaction heat ranged from 0.197 mJ cm −2 at 76 °C to 0.155 mJ cm −2 at 158 °C, corresponding to an average of 4.0 eV/Zr at all temperatures. A temperature dependence for reaction kinetics was not resolved over the range investigated. The temperature dependence of net reaction heat and distribution among metalorganic and oxygen source exposure is attributed to factors including growth rate, equilibrium surface hydroxylation, and the extent of the reaction. ZrO 2 ‐forming surface reactions were investigated computationally using DFT methods to better understand the influence of surface hydration on reaction thermodynamics.

Bielinski, Ashley R.↗
DOE OSTI2025

Design Strategies Based on Electronic Interactions for Effective Catalysts in Lithium–Sulfur Batteries

Abstract Lithium–sulfur batteries (LSBs) are considered promising next‐generation batteries due to their high energy density (>500 W h kg −1 ). However, LSBs exhibit an unsatisfactory energy density (<400 W h kg −1 ) and cycle life (<300 cycles) because of the shuttle effect caused by soluble lithium polysulfide (LiPS) intermediates and the sluggish conversion reaction kinetics caused by insulating sulfur (S 8 ) and lithium sulfide (Li 2 S). Although various types of catalysts, including metal‐based compounds to single‐atom catalysts, have been reported to address these issues, most catalysts exhibited limited catalytic activity under practical lean electrolyte conditions (<5 µL mg −1 ). A comprehensive understanding of the synthetic strategy and catalytic mechanism of catalysts is essential for their design, but understanding the electronic effects of the catalysts and LiPS is more important. Furthermore, the electronic design of these catalysts is not well understood. In this review, we introduce the catalytic mechanisms in LSBs and discuss catalyst design strategies in terms of electronic effects on the interactions between reactants and catalysts, with a primary focus on heterogeneous catalytic systems. We additionally consider how the electronic property of homogeneous systems, particularly redox mediators, affects catalytic behavior under lean electrolyte conditions and propose future research directions for catalyst development in LSBs.

Chemistry↗
DOE OSTI2020

Nonequilibrium Solvent Effects during Photodissociation in Liquids: Dynamical Energy Surfaces, Caging and Chemical Identity

In the gas phase, potential energy surfaces can be used to provide insight into the details of photochemical reaction dynamics. In solution, however, it is unclear what potential energy surfaces, if any, can be used to describe even simple chemical reactions such as the photodissociation of a diatomic solute. In this paper, we use mixed quantum/classical (MQC) molecular dynamics (MD) to study the photodissociation of $Na$ $^+_2$ in both liquid Ar and liquid tetrahydrofuran (THF). We examine both the gas-phase potential surfaces and potentials of mean force (PMF), which assume that the solvent remains at equilibrium with the solute throughout the photodissociation process and show that neither resemble a nonequilibrium dynamical energy surface that is generated by taking the time integral of work. For the photodissociation of $Na$ $^+_2$ in liquid Ar, the dynamical energy surface shows clear signatures of solvent caging, and the degree of caging is directly related to the mass of the solvent atoms. For $Na$ $^+_2$ in liquid THF, local specific interactions between the solute and solvent lead to changes in chemical identity that create a kinetic trap that effectively prevents the molecule from dissociating. Here, the results show that nonequilibrium effects play an important role even in simple solution-phase reactions, requiring the use of dynamical energy surface to understand such chemical events.

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

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Catalytic reforming of methane with H 2 S via dynamically stabilized sulfur on transition metal oxides and sulfides

Reforming of methane with H 2 S is a promising path to directly utilize sour natural gas reserves, although some aspects of the mechanism and structure–function relations remain elusive. In this report we show that metal oxides of group 4–6 elements, which are inert for steam and dry methane reforming reactions, are active and stable (pre)catalysts for the H 2 S reforming of methane. The key active sites are sulfur species (S*) that are dynamically bound to metal cations during catalysis. Similar H/D isotope exchange patterns and universal rate inhibition by H 2 on representative catalysts indicate that H 2 S decomposition and recombination of surface hydrogen atoms are quasi-equilibrated, whereas CH 4 dissociation steps are reversible, yet not quasi-equilibrated. An in-depth analysis of the kinetic data and isotopic substitution effects identifies S*-mediated C–H bond cleavage as the most plausible rate-limiting step common for all catalysts, with subtle yet essential differences between 3d and 4d/5d catalysts in the thermodynamic stability of S*.

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

Ultrafast conformational dynamics of Rydberg-excited N -methyl piperidine

We have observed the ultrafast conformational dynamics of electronically excited N-methyl piperidine (NMP) using time-resolved Rydberg fingerprint spectroscopy. Optical excitation at various wavelengths ranging from 212 nm to 229 nm leads to the 3s or 3p Rydberg states and induces coherent oscillatory motions with periods of about 700 fs. Furthermore, these coherent motions survive the internal conversion from 3p to 3s but then dephase on a time scale of a few oscillations. Intramolecular vibrational energy redistribution on a picosecond time scale leads to an equilibrium between two conformeric structures that are separated in binding energy by 0.09 eV. Model calculations using the DFT-SIC method are in excellent agreement with the experiments and identify the conformers as the chair and twist structures of NMP. The analysis of the equilibrium parameters at long time delays as a function of excitation wavelength allows for the extraction of thermodynamic parameters for the conformeric transformation. We derive an enthalpy of the chair to twist reaction in the 3s excited state of 62 meV with an entropy of 19.70 J mol –1 K –1 . An activation energy of 276 meV is also obtained with a kinetic model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

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

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2023

Bioinspired stability enhancement in deuterium-substituted organic–inorganic hybrid perovskite solar cells

Abstract In hybrid perovskite solar cells (PSCs), the reaction of hydrogens (H) located in the amino group of the organic A-site cations with their neighboring halides plays a central role in degradation. Inspired by the retarded biological activities of cells in heavy water, we replaced the light H atom with its abundant, twice-as-heavy, nonradioactive isotope, deuterium (D) to hamper the motion of H. This D substitution retarded the formation kinetics of the detrimental H halides in Pb-based PSCs, as well as the H bond-mediated oxidation of Sn2+ in Sn–Pb-based narrow-bandgap PSCs, evidenced by accelerated stability studies. A computational study indicated that the zero point energy of D-based formamidinium (FA) is lower than that of pristine FA. In addition, the smaller increase in entropy in D-based FA than in pristine FA accounts for the increased formation free energy of the Sn2+ vacancies, which leads to the retarded oxidation kinetics of Sn2+. In this study, we show that substituting active H with D in organic cations is an effective way to enhance the stability of PSCs without sacrificing photovoltaic (PV) performance. This approach is also adaptable to other stabilizing methods.

14 SOLAR ENERGY↗
DOE OSTI2022

Nuclear cusps and singularities in the nonadditive kinetic potential bifunctional from analytical inversion

The non-additive kinetic potential $v^{\text{NAD}}$ is a key quantity in density-functional theory (DFT) embedding methods, such as frozen density embedding theory and partition DFT. $v^{\text{NAD}}$ is a bi-functional of electron densities $\rho_{\rm B}$ and $\rho_{\rm tot} = \rho_{\rm A} + \rho_{\rm B}$. It can be evaluated using approximate kinetic-energy functionals, but accurate approximations are challenging. The behavior of $v^{\text{NAD}}$ in the vicinity of the nuclei has long been questioned, and singularities were seen in some approximate calculations. In this article, the existence of singularities in $v^{\text{NAD}}$ is analyzed analytically for various choices of $\rho_{\rm B}$ and $\rho_{\rm tot}$, using the nuclear cusp conditions for the density and Kohn-Sham potential. It is shown that no singularities arise from smoothly partitioned ground-state Kohn-Sham densities. We confirm this result by numerical calculations on diatomic test systems HeHe, HeLi$^+$, and H$_2$, using analytical inversion to obtain a numerically exact $v^{\rm NAD}$ for the local density approximation. We examine features of $v^{\rm NAD}$ which can be used for development and testing of approximations to $v^{\rm NAD}[\rho_{\rm B},\rho_{\rm tot}]$ and kinetic-energy functionals.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Single Atomic Iron Site Catalysts via Benign Aqueous Synthesis for Durability Improvement in Proton Exchange Membrane Fuel Cells

Atomically-dispersed iron-nitrogen-carbon (Fe–N–C) catalysts have arisen as promising candidates for replacing the costly precious metal catalysts in fuel cells but still face some grand challenges, such as insufficient site density and durability. Herein, we report a self-assembly method in an aqueous solution to develop an atomically-dispersed iron catalyst with high oxygen reduction reaction (ORR) activity and stability in acidic electrolytes. As determined by high-resolution transmission electron microscopy (HR-TEM), X-ray absorption spectroscopy (XAS), and high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM), this benign aqueous synthesis strategy facilitates the formation of homogeneous atomic nitrogen-coordinated iron sites embedded in a popcorn-like porous graphitic carbon matrix. These catalyst properties contribute to the improved ORR kinetic current density and mass transport. By controlling synthesis chemistry, the correlation between structure and property is systematically investigated. In this work, the iron content is the most critical material property and can regulate site density and graphitic carbon structures in the catalyst, impacting catalytic activity and stability. The enhanced performance and durability were examined in both acidic aqueous electrolytes and membrane electrode assemblies.

25 ENERGY STORAGE↗
NASA NTRS1972

Elasticity of some mantle crystal structures. I - Pleonaste and hercynite spinel.

The elasticity of high-pressure mantle phases can be characterized by using data for chemically similar crystal compounds. The single-crystal elastic constants are determined as a function of pressure and temperature for pleonaste spinel and at room conditions for hercynite spinel. The bulk modulus increases from 1.95 Mb for pleonaste spinel to 2.10 Mb for hercynite spinel. Low or negative values of the pressure derivatives of shear constants are characteristic of the spinel structure and imply a low kinetic barrier to phase transformations and diffusion. Compressional and shear velocities of the spinel phase of olivine are estimated as a function of mean atomic weight by using the pleonaste and hercynite data.

Wang, H.↗
NASA NTRS1975

Experimental study of the C/3P/ + N/4S/ + M yields CN/B 2Sigma+/ + M recombination

The kinetics of the CN violet band emission produced by the reaction of nitrogen with cynogen is investigated as a function of the densities of atomic carbon and atomic nitrogen and of total pressure. It is shown that CN is formed in a three-body recombination reaction from C and N atoms. The reactions studied have relevance to investigations of the chemistry of interstellar molecules.

Washida, N.↗
NASA NTRS1976

The broadening of the calcium resonance line in a high-temperature helium atmosphere

The paper describes experimental measurements of the shape of the Ca I resonance line at 4227 A in a high-temperature helium atmosphere. A ballistic piston compressor was used to produce hot helium in the temperature range from 3000 to 7000 K and the number-density range from 2 to 4 by 10 to the 20th power per cu cm, which conditions approximate those postulated for the atmospheres of certain cool white-dwarf stars. The Boltzmann temperature of the calcium atoms was measured by the brightness-emissivity method, the absorption line shape was measured with a twelve-channel polychromator, and the helium kinetic temperature was calculated from the equation of an ideal gas. The observed deviation from thermodynamic equilibrium is illustrated. The value of the ratio of the damping constant to the helium number density at 5000 K is found to be 45 (+ or - 4) by 10 to the -22nd power A cu cm, or 4.7 (+ or - 0.4) by 10 to the -9th power rad/s cu cm. It is noted that no existing line-broadening calculation for the calcium-helium system is consistent with the present results. Recent studies of the 4227-A line in spectra of the sun and the white dwarf van Maanen 2 are discussed.

Driver, R. D.↗
NASA NTRS1985

Rapid evaluation of ion thruster lifetime using optical emission spectroscopy

A major life-limiting phenomenon of electric thrusters is the sputter erosion of discharge chamber components. Thrusters for space propulsion are required to operate for extended periods of time, usually in excess of 10,000 hr. Lengthy and very costly life-tests in high-vacuum facilities have been required in the past to determine the erosion rates of thruster components. Alternative methods for determining erosion rates which can be performed in relatively short periods of time at considerably lower costs are studied. An attempt to relate optical emission intensity from an ion bombarded surface (screen grid) to the sputtering rate of that surface is made. The model used a kinetic steady-state (KSS) approach, balancing the rates of population and depopulation of ten low-lying excited states of the sputtered molybdenum atom (MoI) with those of the ground state to relate the spectral intensities of the various transitions of the MoI to the population densities. Once this is accomplished, the population density can be related to the sputtering rate of the target. Radiative and collisional modes of excitation and decay are considered. Since actual data has not been published for MoI excitation rate and decay constants, semiempirical equations are used. The calculated sputtering rate and intensity is compared to the measured intensity and sputtering rates of the 8 and 30 cm ion thrusters.

Rock, B. A.↗
NASA NTRS1985

Rapid evaluation of ion thruster lifetime using optical emission spectroscopy

A major life-limiting phenomenon of electric thrusters is the sputter erosion of discharge chamber components. Thrusters for space propulsion are required to operate for extended periods of time, usually in excess of 10,000 hr. Lengthy and very costly life-tests in high-vacuum facilities have been required in the past to determine the erosion rates of thruster components. Alternative methods for determining erosion rates which can be performed in relatively short periods of time at considerably lower costs are studied. An attempt to relate optical emission intensity from an ion bombarded surface (screen grid) to the sputtering rate of that surface is made. The model used a kinetic steady-state (KSS) approach, balancing the rates of population and depopulation of ten low-lying excited states of the sputtered molybdenum atom (MoI) with those of the ground state to relate the spectral intensities of the various transitions of the MoI to the population densities. Once this is accomplished, the population density can be related to the sputting rate of the target. Radiative and collisional modes of excitation and decay are considered. Since actual data has not been published for MoI excitation rate and decay constants, semiempirical equations are used. The calculated sputtering rate and intensity is compared to the measured intensity and sputtering rates of the 8 and 30 cm ion thrusters.

Rock, B. A.↗
NASA NTRS1990

A unifying picture of gas-phase formation and growth of PAH (Polycyclic Aromatic Hydrocarbons), soot, diamond and graphite

A variety of seemingly different carbon formation processes -- polycyclic aromatic hydrocarbons and diamond in the interstellar medium, soot in hydrocarbon flames, graphite and diamond in plasma-assisted-chemical vapor deposition reactors -- may all have closely related underlying chemical reaction mechanisms. Two distinct mechanisms for gas-phase carbon growth are discussed. At high temperatures it proceeds via the formation of carbon clusters. At lower temperatures it follows a polymerization-type kinetic sequence of chemical reactions of acetylene addition to a radical, and reactivation of the resultant species through H-abstraction by a hydrogen atom.

Frenklach, Michael↗
DOE OSTI2023

OH Roaming during the Ozonolysis of α-Pinene: A New Route to Highly Oxygenated Molecules?

The formation of low-volatility organic compounds in the ozonolysis of α-pinene, the dominant atmospheric monoterpene, provides an important route to aerosol formation. Here, in this work, we consider a previously unexplored set of pathways for the formation of highly oxygenated molecules in α-pinene ozonolysis. Pioneering, direct experimental observations of Lester and co-workers have demonstrated a significant production of hydroxycarbonyl products in the dissociation of Criegee intermediates. Theoretical analyses indicate that this production arises from OH roaming-induced pathways during the OO fission of the vinylhydroperoxides (VHPs), which in turn come from internal H transfers in the Criegee intermediates. Ab initio kinetics computations are used here to explore the OH roaming-induced channels that arise from the ozonolysis of α-pinene. For computational reasons, the calculations consider a surrogate for α-pinene, where two spectator methyl groups are replaced with H atoms. Multireference electronic structure calculations are used to illustrate a variety of energetically accessible OH roaming pathways for the four VHPs arising from the ozonolysis of this α-pinene surrogate. Ab initio transition-state theory-based master equation calculations indicate that for the dissociation of stabilized VHPs, these OH roaming pathways are kinetically significant with a branching that generally increases from ~20% at room temperature up to ~70% at lower temperatures representative of the troposphere. For one of the VHPs, this branching already exceeds 60% at room temperature. For the overall ozonolysis process, these branching ratios would be greatly reduced by a limited branching to the stabilized VHP, although there would also be some modest roaming fraction for the nonthermal VHP dissociation process. The strong exothermicities of the roaming-induced isomerizations/additions and abstractions suggest new routes to fission of the cyclobutane rings. Such ring fissions would facilitate further autoxidation reactions, thereby providing a new route for producing highly oxygenated nonvolatile precursors to aerosol formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗