Search NASA⌕ Search

SEARCH · Search NASA

Results for “Chemical kinetics and 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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids

We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.

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

Chemical Thermodynamics and the Mathematical Integration of Reaction Kinetics

Key advances in the development of numerical methods for non-reacting compressible flows have been enabled by translating physical requirements into concrete numerical guidelines, such as the satisfaction of entropy inequalities for shock-capturing techniques [Lax, Contributions to Nonlinear Functional Analysis (1971) 603-634]. In the present work, we present nonlinear numerical analysis tools that draw from Chemical Thermodynamics , the branch of Nonequilibrium Thermodynamics that deals with chemical reactions. Through Gibbs formalism, chemical thermodynamics provides a well-known theoretical expression for the chemical equilibrium constant of a reaction in terms of reduced chemical potentials. A less-known, yet extremely valuable result, due to [Krambeck, Arch. Ration. Mech. Anal. , 38 (1970) 317], states that when this expression is implemented, mass-action kinetic models are consistent with the dynamical prescriptions of the 2nd law of thermodynamics. For fixed-temperature ordinary differential equations modeling constant-volume reacting gas mixtures, this leads to a decreasing Helmholtz free energy. If the temperature is allowed to vary in accordance with conservation of energy (1st law), this leads to the statement of increasing entropy. These nonlinear prescriptions can, and should be, used to further develop temporal integration techniques for reaction kinetics. We demonstrate that Krambeck's result holds even when the equilibrium constants are approximated from data. We prove this result by constructing the implicit free energy and the implicit entropy inherent to a given approximation. This is first done for a 5-species, 17-reaction model problem for air. With this structure established, elements of discrete entropy-stability theory [Tadmor, Acta Numer. , 12 (2003) 451] are leveraged to examine the consistency of time-integration schemes with these prescriptions. Using chemical potentials, one can compute the respective contributions of the kinetics model and of the temporal scheme to free energy/entropy variations. We introduce a nonlinear-stable version of the Discontinuous-Galerkin (DG) scheme in time which shows robustness improvements over the standard linearly-stable version. Most notably, the maximum timestep that can be resolved with the nonlinearly-stable variant tends to grow with polynomial order, in contrast to the linearly-stable variant. We generalize our constructions to arbitrary systems of reversible chemical reactions, ultimately showing that the compressible reacting Euler system admits the opposite of the implicitly constructed thermodynamic entropy as a mathematical entropy . This lays important theoretical foundations towards robust scheme development [Harten, J. Comput. Phys. 49 (1983) 151-164].

STMD↗
DOE OSTI2023

Ultrafast dynamics of CN radical reactions with chloroform solvent under vibrational strong coupling

Polariton chemistry may provide a new means to control molecular reactivity, permitting remote, reversible modification of reaction energetics, kinetics, and product yields. A considerable body of experimental and theoretical work has already demonstrated that strong coupling between a molecular vibrational mode and the confined electromagnetic field of an optical cavity can alter chemical reactivity without external illumination. However, the mechanisms underlying cavity-altered chemistry remain unclear in large part because the experimental systems examined previously are too complex for detailed analysis of their reaction dynamics. Here, we experimentally investigate photolysis-induced reactions of cyanide radicals with strongly-coupled chloroform (CHCl 3 ) solvent molecules and examine the intracavity rates of photofragment recombination, solvent complexation, and hydrogen abstraction. Here we use a microfluidic optical cavity fitted with dichroic mirrors to facilitate vibrational strong coupling (VSC) of the C–H stretching mode of CHCl 3 while simultaneously permitting optical access at visible wavelengths. Ultrafast transient absorption experiments performed with cavities tuned on- and off-resonance reveal that VSC of the CHCl 3 C–H stretching transition does not significantly modify any measured rate constants, including those associated with the hydrogen abstraction reaction. This work represents, to the best of our knowledge, the first experimental study of an elementary bimolecular reaction under VSC. We discuss how the conspicuous absence of cavity-altered effects in this system may provide insights into the mechanisms of modified ground state reactivity under VSC and help bridge the divide between experimental results and theoretical predictions in vibrational polariton chemistry.

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

Chemistry of Sugar Formation in the Gas Phase: Following the Activated Aldehyde

Sugars are produced by living organisms, and are required building blocks for life as we know it, which raises the foundational question of how sugars formed in a prebiotic environment. The abiotic formose reaction produces sugars from formaldehyde, but our understanding of its initiation step remains murky, with chemists invoking the concept of an “activated aldehyde” to seed this reaction. Singlet hydroxycarbenes, high-energy isomers of aldehydes, were recently reported to facilitate sugar formation under cold, nonaqueous conditions relevant to interstellar environments. Here, we generate singlet methylhydroxycarbene ( 1 CH 3 –C̈–OH) from the photodissociation of pyruvic acid and experimentally measure its gas-phase reaction with d 4 -acetaldehyde using multiplexed photoionization mass spectrometry. The C 4 H 4 D 4 O 2 isomer d 4 -acetoin is the sole product, which we kinetically link to the reactant CH 3 –C̈–OH, and attribute to a carbonyl-ene formation mechanism. We see no evidence of 3-hydroxybutanal, the C–H insertion product expected in carbene chemistry. Using automated exploration we calculate stationary points on the potential energy surface and report master equation rate coefficients from T = 20–600 K, providing quantitative kinetics of this fast reaction for use in chemical models. The prereactive complex in this reaction is stabilized by both hydrogen bonding and electrophilic carbene C═O interactions. These effects create a short-range dynamical bottleneck for the reaction besides the long- and midrange barrierless bottlenecks. Combined with recent reports of 1 HC̈OH production from methanol photodissociation and pyruvic acid production in cold irradiated ices, this work provides evidence that singlet hydroxycarbene + aldehyde chemistry is a feasible path to prebiotic sugar formation.

aldehydes↗
NASA NTRS2003

Dynamical Study of Prebiotic Processing by Comet Impact

Comets and meteoroids that bombarded the Earth, other planets and moons are considered possible deliverers of prebiotic materials manufactured in space. Simultaneously, chemical processing can be initiated by the large kinetic energy imparted to the planetary atmosphere during comet/meteoroid entry. The efficiency of organic synthesis and the diversity of products by impact shock are known to decrease as the reducing power of the atmosphere decreases. It is generally accepted that rich organic products are produced in a methane atmosphere whereas a carbon dioxide atmosphere is reported to yield a dearth of products In order to understand the details of impact chemistry and how it depends on the atmospheric composition, we carried out simulations of the chemistry initiated by comet/meteoroid impact upon a planetary atmosphere using different atmospheric compositions. The simulations were done by solving the set of coupled equations for mass, momentum, and energy conservations, chemical kinetics, and transport, that describe a high-energy impact shock, subsequent expansion and cooling of the hot shocked gas by mixing with the ambient gas, and the eventual steady state composition.

Huo, Winifred M.↗
DOE OSTI2023

Chemical looping conversion of CH4/CO2 to syngas on 5wt.%Ni@Ce0.6Zr0.4O2 catalyst: Impact of dynamic accumulation of surface carbon and oxygen vacancies

Syngas, a combination of carbon monoxide (CO) and hydrogen (H2), is a precursor to many chemicals and fuels, contributing to the billion-dollar global hydrocarbon industry. Chemical looping reforming (CLR) of the greenhouse gases methane (CH4) and carbon dioxide (CO2) allows for energy-efficient production of syngas. 5wt.% nickel (Ni) on ceria-zirconia (5wt%Ni@Ce0.6Zr0.4O2) mixed metal oxide catalyst was investigated here to explore pathways for enhanced syngas production on sustainable earth-abundant transition metal supported catalysts. The role of reduction-oxidation (redox) state of the catalyst, and carbon formation on the catalyst surface during chemical looping is explored to drive superior reaction kinetics, conversions, selectivity, and syngas ratios (H2/CO). The bulk and surface structure of the catalyst, along with carbon deposition features were characterized by electron microscopy, X-Ray diffraction, and ex situ Raman spectroscopy. The dynamic evolution of catalysts under CLR reaction conditions and the intrinsic reaction mechanisms were probed with in situ Raman spectroscopy, in situ Fourier transform spectroscopy, dynamic oxygen storage capacity (DOSC), and the Temporal Analysis of Products (TAP) reactor studies. Intrinsic kinetics of syngas production via CLR was correlated to the redox state of catalyst and the participation of nickel-catalyzed multiwalled carbon nanotube (CNT) growth, allowing enhanced CLR reaction performance.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Single-shot ultrafast dynamics of nanosecond pulsed plasmas: Transition from ps–ns nonequilibrium to near-full ionization

The ultrafast multi-stage evolution of state-defining properties in atmospheric-pressure nanosecond pulsed plasmas is quantified using single-discharge, jitter-free, continuous streak-sweep spectroscopy of N 2 (C → B) molecular and N + /O + ionic species with a single-shot time resolution as short as ∼160 ps. An early period of extreme nonequilibrium is observed with vibrational temperatures [T V (C)] dropping from ∼8000 K at <200 ps after plasma breakdown to <4000 K within ∼1 ns, with near-ambient rotational temperatures [T R (C)] ∼ 300 K due to limited collisional energy transfer. This early-time T V (C) trend is representative of a direct and unquenched look into the high-energy tail of the electron energy distribution function; thus, it tracks real-time changes in the mean electron energy via the N2(C) emission signatures. This is followed by the rapid onset of N + /O + ionic emission after a distinct time delay of ∼14.9 ns due to multistep chemical kinetics. The ionic emission enables determination of electron densities (n e ) >2 × 10 19 cm −3 and electron temperatures (T e ) >36 000 K, indicating the transition to a nearly fully ionized regime. The ps–ns temporal dynamics are also compared between air and N 2 plasmas to assess the influence of collisional partners, as well as across the anode, cathode, and central gap regions to identify spatial variations in plasma behavior. Finally, this work demonstrates versatile continuously probing approach capable of spectro-temporal and spatially resolved single-shot measurements of key state variables in the ps–ns evolution of atmospheric-pressure nanosecond pulsed plasmas.

Electric discharges↗
DOE OSTI2023

On the Interfacial Assembly of Anisotropic Amphiphilic Janus Particles

It has been shown both theoretically and experimentally that amphiphilic Janus particles are the most effective solid surfactants to stabilize interfaces. In most cases, the Janus particles investigated have uniform morphologies with Janus boundaries dividing the particle into halves. However, there are many examples of Janus particles where the hydrophilic and hydrophobic domains are not equally distributed. The effects of this uneven domain distribution on the mechanism and kinetics of Janus particle assembly, and final equilibrium state are not well-understood. Dynamic pendant drop tensiometry offers a means to probe both the equilibrium assembly and the kinetics and mechanism of assembly. Here, in this work, the interfacial kinetics and assembly of spherical anisotropic Janus particles are investigated using dynamic pendant drop tensiometry. Systematic studies quantifying the time-dependent interfacial behavior as a function of Janus particle morphology, chemical composition, particle concentration, and NaOH and HCl concentration are performed. These studies shed light on the assembly mechanism of more complex Janus particle morphologies and highlight their effectiveness as interface stabilizers.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Thermal Oxidation on the Structure, Surface Texturing, and Microstructure Evolution in Nanocrystalline Ga-O-N Films

An extensive examination of the nanoscale, crystallographic growth dynamics of the system, which is impacted by the thermal energy given to the GaN, is carried out to derive a deeper understanding of the growth kinetics, morphology and microstructure evolution, chemical bonding, and optical properties of Ga-O-N films. Thermal annealing of GaN films is performed in the temperature range of 900–1200 °C. Crystal structure, phase formation, chemical composition, surface morphology, and microstructure evolution of Ga-O-N films are investigated as a function of temperature. Increasing temperature induces surface oxidation, which results in the formation of stable β-Ga2O3 phase in the GaN matrix, where the overall film composition evolves from nitride (GaN) to oxynitride (Ga-O-N). While GaN surfaces are smooth, planar, and featureless, oxidation induced granular-to-rod shaped morphology evolution is seen with increasing temperature to 1200 °C. The considerable texturing and stability of the nanocrystalline Ga-O-N on Si substrates can be attributed to the surface and interface driven modification because of thermal treatment. Corroborating with structure and chemical changes, Raman spectroscopic analyses also indicate that the chemical bonding evolution progresses from fully Ga-N bonds to Ga-O-N. While the GaN oxidation process starts with the formation of β-Ga 2 O 3 at an annealing temperature of 1000 °C, higher annealing temperatures induce structural distortion with the potential formation of Ga-O-N bonds. The structure-phase-chemical composition correlation, which will be useful for nanocrystalline materials for selective optoelectronic applications, is established in Ga-O-N films made by thermal treatment of GaN.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Design of light- and chemically responsive protein assemblies through host-guest interactions

Host-guest (HG) interactions have been widely used to build responsive materials and molecular machines owing to their inherently dynamic nature, interaction specificity, and responsiveness to diverse stimuli. Here, in this work, we have set out to exploit these advantages of HG chemistry in the design of dynamic protein assemblies, using a C 4 symmetric protein, C98 RhuA, as a building block. We show that a C98 RhuA variant individually modified with β-cyclodextrin (βCD) (host) or azobenzene (guest) functionalities can specifically pair with each other to form highly ordered 1D and 2D assemblies. Association and dissociation of βCD RhuA- azo RhuA assemblies can be controlled by UV and visible light as well as by small-molecule modulators of βCD-azobenzene interactions. Kinetics analyses reveal that βCD RhuA- azo RhuA nanotubes assemble without a nucleation barrier, a highly unusual occurrence for helical supramolecular systems. Taken together, our findings provide a compelling example for achieving complex structural and dynamic outcomes in protein assembly through simple chemical design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗
NASA NTRS1972

The Aladdin experiment. I - Dynamics.

A coordinated rocket experiment was carried out on Nov. 20, 1970 at Eglin AFB, Fla. (30 deg N) to determine atmospheric layering and density distributions of ions and neutrals. The program measured vertical profiles between 70 and 150 km altitude of total neutrals, ions and electrons, individual species, temperature, density, diffusion coefficients, horizontal winds, and turbulent structure. A dynamic computer code model using these parameters was compared with the observed profiles. Wind fields were measured from seven chemical trails, separated up to 60 km horizontally and 600 sec in time. Kinetic temperatures were measured from fluorescence of chemically-released AlO, ranging from 440 to 680 K between 126 and 161 km. Densities derived from diffusion of the AlO ranged from 10 to 5 times 10 to the -12th power per cu cm between 126 and 139 km. A rocket-borne pitot probe showed wavelike density deviations from current standard atmosphere models, amounting to -25% at 80 and 110 km and to +25% at 100 and 125 km. The derived temperature profile from the pitot probe showed a double minimum at 74 and 102 km.

Rosenberg, N. W.↗
DOE OSTI2022

Kinetic model development and Bayesian uncertainty quantification for the complete reduction of Fe-based oxygen carriers with CH 4 , CO, and H 2 for chemical looping combustion

In this work, three kinetic models are developed and calibrated for the complete multi-step reduction of an Fe-based oxygen carrier (OC) particle with CH 4 , CO, and H 2 , using data from thermogravimetric analysis. The complete reduction rate profiles exhibit complex dynamics whose trajectory is significantly different depending on the reducing gas. A Bayesian model building and parameter estimation framework is applied for simultaneous parameter and model structure uncertainty quantification. The final models show excellent agreement between model predictions and calibration data, as well as new data not used for calibration (for the reduction of the OC with CH 4 ). Parameter uncertainty is quantified by determining joint posterior distribution, and model structure uncertainty is addressed by incorporating Gaussian process stochastic functions (represented by Bayesian smoothing splines) into the kinetic models. The final kinetic models with discrepancy functions are readily employable in equation-oriented simulation and optimization platforms.

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

Fine-tuning catalytic selectivity by modulating catalyst-environment interactions: CO 2 hydrogenation over Pd-based catalysts

Capturing catalytic behaviors under operational conditions is pivotal to gaining a mechanistic understanding and promoting the design of robust catalysts. The challenge lies in the difficulty of monitoring real-time surface dynamics driven by catalyst-environment interactions. Here, in this work, we introduce a framework based on density functional calculations and kinetic modeling. This framework significantly improves the accuracy of theoretical models’ descriptions of experimental observations by quantifying environmental impacts on surface phases and active sites. CO 2 hydrogenation over Pd-based catalysts is taken as a showcase. The observed selectivity variations of Pd and Pd-M bimetallic catalysts strongly correlate with hydrogen coverage maintained under typical CO 2 hydrogenation conditions. By reducing the amount of surface hydrogen, the selectivity tuned effectively from formic acid toward CO and methanol. This study not only deepens the comprehension of dynamics of active sites under active chemical conditions but also introduces an alternative opportunity for catalytic tuning by modulating catalyst-environment interactions.

03 NATURAL GAS↗
DOE OSTI2022

Structural, Dynamic, and Chemical Complexities in Zinc Anode of an Operating Aqueous Zn-Ion Battery

Aqueous Zn-ion battery is a promising technology for electrochemical energy storage. The formation of Zn dendrites, however, can jeopardize the cell cycle life and thus, hinders the industrial adoption of this technology. A fundamental understanding of the kinetic mechanisms is crucial for improving the Zn-ion battery. Here, in situ and operando X-ray microscopy methods are utilized to visualize the Zn plating and stripping behaviors under different electrochemical conditions. It is demonstrated that the substrate curvature, local morphology, electrochemical protocols, and the surface chemistry can collectively affect the Zn plating behavior. Furthermore, these results provide new insights for developing the next-generation dendrite-free and long-span aqueous Zn-ion battery.

25 ENERGY STORAGE↗
DOE OSTI2023

Development and evaluation of a skeletal mechanism for EHN additized gasoline mixtures in large Eddy simulations of HCCI combustion

Advanced Low Temperature Combustion modes, such as the Sandia proposed Additive-Mixing Fuel Injection (AMFI), can unlock significant potential to boost fuel conversion efficiency and ultimately improve the energy conversion of internal combustion engines. This is a novel improved combustion process that is enabled by supplying small (<5%) variable amounts of autoignition improver to the fuel to enhance the engine operation and control. Common, diesel-fuel ignition-quality enhancing additive, 2-ethylexyl nitrate (EHN), is doped into gasoline to enable Sandia LTGC + AMFI combustion. This manuscript focuses on the development of a reduced sub-mechanism for EHN chemical kinetics at engine relevant conditions that is implemented into a skeletal mechanism for chemical kinetic studies of gasoline surrogate fuels. The mechanism validation utilized zero-dimensional numerical simulations and comparison to shock tube ignition-delay data of pure and EHN-doped n-heptane. Additional validation is presented with Homogeneous Charge Compression-Ignition (HCCI) engine data of pure and EHN-doped research-grade E10 gasoline. Then, the mechanism was deployed in a 3-D computational fluid dynamics (CFD) using Large Eddy Simulations (LES) to model the HCCI engine experiments of 0.4% vol EHN additized E10 gasoline at several equivalence ratios. Simulations showed a very good performance of the mechanism, and the model accurately reproduced (a) the ignition point, (b) combustion phasing, (c) combustion duration, and (d) the peak of the heat release rates of the engine experiments. The results show that EHN promotes Low-Temperature Heat Release, ultimately driving the gasoline to autoignite at thermodynamic conditions where the fuel would not otherwise ignite. Overall, this work demonstrates a viable reduced chemical-kinetic mechanism for EHN and shows that it can be combined with a skeletal gasoline mechanism for CFD-LES analysis of well-mixed LTGC that matches well with experimental results. The CFD-LES analysis also shows the spatial distribution of EHN-fuel interactions that control the autoignition throughout the combustion chamber.

Guleria, Gaurav↗
DOE OSTI2026

Kinetic control of phase evolution and defect-mediated recombination in AACVD-grown copper antimony sulphide thin films

Copper antimony sulphide (CAS) is a multinary chalcogenide semiconductor in which small deviations from stoichiometry can drive phase competition and strong defect-mediated modulation of optoelectronic properties. However, the systematic roles of copper precursor fraction and deposition time in governing phase evolution, off-stoichiometry, and recombination dynamics in aerosol-assisted chemical vapour deposition (AACVD)-grown undoped CAS thin films remain insufficiently understood. In this work, undoped CAS thin films were deposited by AACVD using Cu(dedtc)2 and Sb(dedtc)3 single-source precursors at 550 °C and a carrier gas flow rate of 150 sccm, while the Cu(dedtc)2 mole fraction (x = 0.15 – 0.55) and deposition time (1 – 2 hours) were systematically varied to probe how growth kinetics influence phase composition, microstructure, and defect-mediated optical properties without post-deposition annealing or extrinsic doping. Increasing copper precursor content drives phase evolution toward tetrahedrite-dominant CAS films at intermediate Cu(dedtc)2 mole fractions, with the film deposited at x = 0.35 exhibiting the strongest tetrahedrite character within the parameter space examined. The films are also copper-rich, antimony-poor, and sulphur-deficient, consistent with off-stoichiometric growth and intrinsic defect formation, plausibly including copper interstitials, Cu-on-Sb antisites, and sulphur vacancies. These growth-dependent compositional deviations are accompanied by tunable indirect optical bandgaps of approximately 1.60 – 2.18 eV and weak visible photoluminescence governed by defect-mediated recombination. Time-resolved photoluminescence reveals bi-exponential decay behaviour with lifetimes of approximately 0.1 – 3.6 ns, with emission dominated by slower donor–acceptor pair recombination and a smaller contribution from faster trap-assisted pathways. Collectively, these results establish an explicit kinetic process–structure–defect–property relationship for AACVD-grown CAS thin films and provide a growth–structure–property framework relevant to future optimization of CAS-based optoelectronic and energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomic-Level Features for Kinetic Monte Carlo Models of Complex Chemistry from Molecular Dynamics Simulations

The high computational cost of evaluating atomic interactions recently motivated the development of computationally inexpensive kinetic models, which can be parametrized from MD simulations of complex chemistry of thousands of species or other processes and accelerate the prediction of the chemical evolution by up to four order of magnitude. Such models go beyond the commonly employed potential energy surface fitting methods in that they are aimed purely at describing kinetic effects. So far, such kinetic models utilize molecular descriptions of reactions and have been constrained to only reproduce molecules previously observed in MD simulations. Therefore, these descriptions fail to predict the reactivity of unobserved molecules, for example in the case of large molecules or solids. In this work, we propose a new approach for the extraction of reaction mechanisms and reaction rates from MD simulations, namely the use of atomic-level features. Using the complex chemical network of hydrocarbon pyrolysis as example, it is demonstrated that kinetic models built using atomic features are able to explore chemical reaction pathways never observed in the MD simulations used to parametrize them, a critical feature to describe rare events. Atomic-level features are shown to construct reaction mechanisms and estimate reaction rates of unknown molecular species from elementary atomic events. Through comparisons of the model ability to extrapolate to longer simulation timescales and different chemical compositions than the ones used for parameterization, it is demonstrated that kinetic models employing atomic features retain the same level of accuracy and transferability as the use of features based on molecular species, while being more compact and parametrized with less data. We also find that atomic features can better describe the formation of large molecules enabling the simultaneous description of small molecules and condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗