Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Complex chemical systems”

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

A Visual Understanding of Circular Dichroism Spectroscopy

Mapping chemical and structural properties to electronic and magnetic responses is critical to many applications such as quantum information science, where the precise storage and transmission of unique information is paramount. Specifically, constructing molecules and materials that provide strong polarized responses at tunable frequencies and with large anisotropies is key to optical processing of quantum information. Chiral molecules provide chiroptical response to circularly polarized light, making them attractive for quantum information science and other applications related to sensing, polarized photodetectors, and spintronics. Predicting a molecular design, a priori, with large anisotropies to circularly polarized light is challenging due to the complex interplay between electric and magnetic components of the optical response. In this work, we explore a visual representation of the electronic chiroptical response by decomposing the rotary strength into its constituent components. Here, we make use of the intuitive electronic oscillator framework to develop classical intuition regarding the rotary strength and its constituents. We explore three model chemical systems that exhibit local and global chirality. Our analysis reveals that local chirality necessarily exhibits competition between the local chiral center and chirality induced in other fragments of the molecule, resulting in both unexpected nonmonotonic trends and sign flips in chemically adjacent geometries. Furthermore, we can visually distinguish between local and global chirality via examination of the transition chiral tensor. Interestingly, we make strong connections to ferromagnetic and antiferromagnetic spin systems in that chiroptically inactive transitions exhibit antiferromagnetic-like alternating orbital patterns while active transitions show domain formation in an ferromagnetic-like alignment that produces a net chiroptical response.

36 MATERIALS SCIENCE↗
DOE OSTI2021

In situ ambient pressure XPS study of Pt/Cu(111) single-atom alloy in catalytically relevant reaction conditions

The presence of multiple reactant gases as well as reaction intermediates in a heterogeneous catalytic reaction results in a complex interaction between different components of the catalyst with each gas, which can alter the surface and chemical state of the catalyst differently than in the presence of an individual gas alone. In this study, we used in situ ambient pressure x-ray photoelectron spectroscopy to study the surface state of Pt/Cu(111) single-atom alloy model system in two catalytically relevant reaction conditions: CO 2 hydrogenation and CO oxidation. We found that the activation of CO 2 results in the formation of CO, which adsorbs on Pt sites at up to 400 K. In the presence of CO 2 and H 2 , Pt catalyzes the reverse water–gas shift reaction, which produces more CO and further stabilizes surface Pt atoms at 450 K. On the other hand, in CO oxidation condition, the presence of O 2 results in the formation of a thick Cu 2 O layer at higher temperatures, and Pt atoms are no longer detected in the surface and subsurface layers. When O 2 is introduced to the sample before CO, the formation of a complete Cu 2 O layer that covers all Pt atoms occurs immediately at room temperature. Finally, however, when CO is introduced at room temperature before O 2 , the presence of adsorbed CO on Pt sites stabilizes the surface Pt atoms and prevents the formation of a complete Cu 2 O layer, thus exposing the Pt atoms in 'holes' in the Cu 2 O layer.

36 MATERIALS SCIENCE↗
DOE OSTI2019

A two-phase substrate model for enzymatic hydrolysis of lignocellulose: application to batch and continuous reactors

Abstract Background Enzymatic hydrolysis continues to have a significant projected production cost for the biological conversion of biomass to fuels and chemicals, motivating research into improved enzyme and reactor technologies in order to reduce enzyme usage and equipment costs. However, technology development is stymied by a lack of accurate and computationally accessible enzymatic-hydrolysis reaction models. Enzymatic deconstruction of cellulosic materials is an exceedingly complex physico-chemical process. Models which elucidate specific mechanisms of deconstruction are often too computationally intensive to be accessible in process or multi-physics simulations, and empirical models are often too inflexible to be effectively applied outside of their batch contexts. In this paper, we employ a phenomenological modeling approach to represent rate slowdown due to substrate structure (implemented as two substrate phases) and feedback inhibition, and apply the model to a continuous reactor system. Results A phenomenological model was developed in order to predict glucose and solids concentrations in batch and continuous enzymatic-hydrolysis reactors from which liquor is continuously removed by ultrafiltration. A series of batch experiments were performed, varying initial conditions (solids, enzyme, and sugar concentrations), and best-fit model parameters were determined using constrained nonlinear least-squares methods. The model achieved a good fit for overall sugar yield and insoluble solids concentration, as well as for the reduced rate of sugar production over time. Additionally, without refitting model coefficients, good quantitative agreement was observed between results from continuous enzymatic-hydrolysis experiments and model predictions. Finally, the sensitivity of the model to its parameters is explored and discussed. Conclusions Although the phenomena represented by the model correspond to behaviors that emerge from clusters of mechanisms, and hence a set of model coefficients are unique to the substrate and the enzyme system, the model is efficient to solve and may be applied to novel reactor schema and implemented in computational fluid dynamics (CFD) simulations. Hence, this modeling approach finds the right balance between model complexity and computational efficiency. These capabilities have broad application to reactor design, scale-up, and process optimization.

09 BIOMASS FUELS↗
DOE OSTI2020

Light yield quenching and quenching remediation in liquid scintillator detectors

Quenching of light emission from an LAB based scintillator by the addition of organic amines and carboxylic acids is examined. Chemical functional groups of the quenching agents play an important role in this reduction. It is shown that "salt" formation at a 1:1 mole ratio in a mixed amine-acid system, reduces quenching by a factor of ~ 2. Supporting NMR spectra are presented. Finally, this "quenching neutralization" has the potential to reduce the light loss incurred when metals complexed with quenching agents are loaded into organic scintillators.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Microscale Mechanical‐Chemical Modeling of Granular Salt: Insights for Creep

Abstract Microscale numerical modeling is potentially an effective approach for understanding salt creep, but it is quite challenging because natural salt grain boundaries can have complex and evolving geometry, and because contacts between these deformable salt mineral grains are dynamic as a result of compaction, chemical reaction, fluid flow, and heat transfer. In this study, we have overcome these challenges and developed a new microscale mechanical‐chemical (MC) model to analyze creep of salt at the microscale, accounting for coupled deformation, dynamic contacts, and chemical reaction in granular systems with realistic geometric representations. The MC model was realized by linking a new microscale mechanical code based on the numerical manifold method (NMM) to a reactive transport code named Crunch. We simulated the processes of reorganization of the salt grains, microfracturing, and pressure solution that contribute to creep at larger scales. Based on this first quantitative microscale model, we found that sharp corners of mineral grains can dominate the contact dynamics, microfracturing, and pressure solution when salt is compacted, thus governing the structural changes and porosity loss of the system. We found that pressure solution, which preferentially dissolves sharp corners and edges, can lead to relatively high porosity loss in the system, thus playing an important role in the creep of salt. Our analysis shows that the dynamic changes of salt granular systems involving grain relocation and pressure solution can occur repeatedly and continuously, thus contributing to longer‐term creep of salt at larger scales.

58 GEOSCIENCES↗
DOE OSTI2026

Co-Design of Charge Transport Superhighways to Connect Catalytic Sites in Soft Photoelectrochemical Systems

Efficient photon-to-electron-to-molecule conversion requires multi-length scale control over charge transport pathways, where electronic charges are delivered to catalytic sites under high mass transport flux. A fundamental question is how can we co-design charge transport pathways to promote efficient charge transfer to/from catalytic sites in complex three-dimensional architectures? Soft conducting polymer systems offer exceptional promise to provide three-dimensional charge transport networks, where electrolyte (ion and solvent) can interdiffuse to promote long-lived charge carriers and the molecular nature allows for strategic synthetic design of catalytic sites. Herein we combine theoretical and experimental approaches to investigate the earliest stages of photoelectrochemical deposition of near-surface catalytic sites (Pt) on soft bulk heterojunction polymeric semiconductors composed of a prototype donor (PTB7-Th) and a prototype acceptor (N2200) as a model system towards better understanding molecular catalyst-polymer site interactions. We focus initially on photoelectrochemical deposition of low Pt loadings, nanoparticle sizes (formed by progressive nucleation) below 20 nm, for both density functional theory (DFT) modeling studies and for spectroscopic characterization using surface-sensitive X-ray and UV-photoemission (XPS/UPS). DFT modeling of “n-type” N2200 slabs reveal for the first time that sulfur atoms in the thiophene units serve as the lowest-energy adsorption sites for single Pt atoms, while larger Pt clusters engage more complexly with both thiophene and naphthalene diimide (NDI) core sites. Changes in chemical composition observed by X-ray photoelectron spectroscopy (XPS) support the DFT predictions, and the angle-resolved measurements reveal that Pt nucleation initiates at subsurface sites which appear to be localized active domains that promote charge transport/transfer and enable vertical growth toward the surface. These results suggest that light-activated Pt nanoparticle deposition decorates energetically distinct sites, where photoactivity is dictated by the local energetics of those sites, and the fact that they represent the termini of charge transport “super-highways” – a small percentage of the total volume of the donor/acceptor polymeric active layer which carries most of the photocurrent generated during both Pt deposition and photoelectrochemical HER. We posit that these initial studies provide a foundational strategy for design of catalytic sites in the near surface regions of complex polymeric materials and advancing soft semiconductor-based photoelectrochemical systems. Achieving a nanometer-scale understanding of catalyst deposition and the impact of local composition and energetics on that placement, should ultimately provide the design guidelines (co-design) for a broad array of catalysts at sites that optimize that efficiency and maximize platform durability.

14 SOLAR ENERGY↗
DOE OSTI2025

Chemically and mechanically recyclable polyester-based multilayer plastics

Approximately 100 million tons of multilayered plastics (MLPs) are produced each year worldwide but are not recycled due to their complex structure. Here, this work aims to design polyester-based multilayer plastics (80–100 % polyester) that provide barrier performance comparable to typical 9–12-layer commercial MLPs, while also enabling both chemical recycling (back to parent monomer) and mechanical recycling (grind-and-melt reprocessing). Such dual recyclability is not achievable with non-polyester multilayers, such as all-polyolefin systems. Furthermore, we emphasize how the multilayer architecture was tailored to balance barrier properties, mechanical integrity, and end-of-life recyclability for both flexible and rigid packaging applications. Two main categories of polyester-based MLPs are reported; in the first type, poly(butylene adipate-co-terephthalate) (PBAT)-70 % polyglycolic acid (PGA) is used as middle barrier layer, while in second type, middle barrier layer is Ethylene-vinyl alcohol (EVOH) copolymers. Polyethylene terephthalate (PET) was used as a structural layer, while either PBAT or poly(butylene succinate) (PBS) was used to enable thermal sealing and serve as the product contact layer. These MLPs are recycled by both chemical and mechanical recycling processes. Techno-economic analysis (TEA) shows that MLPs incorporating EVOH as barrier layer have similar or lower selling costs (0.32 $\$$/m 2 ) than commercial MLPs. Life cycle assessment (LCA) indicates EVOH-based MLPs have a lower carbon footprint and lower energy consumption relative to commercial MLP benchmarks. This work offers simplified MLPs that are easy to manufacture and ready to recycle, which will significantly reduce environmental impact of MLP packaging while also providing a cost-effective and practical solution for industry.

Barrier properties↗
DOE OSTI2025

Non-traditional stable isotope measurements using laser desorption Orbitrap mass spectrometry: Implications for planetary missions.

Isotopic fractionation recorded in planetary materials provides insights into physical, chemical, and/or potential biological processes occurring on Solar System bodies. As we enter into the next decades of planetary exploration, the crucial information revealed by isotopic compositions of rocky and icy samples mandates that next generation spaceflight instrumentation possess the capability to measure isotope ratios in situ with sufficient precision/accuracy to distinguish between such processes. Here, in addition to identifying and fingerprinting complex organic materials with high accuracy and ultrahigh mass resolutions, laser desorption Orbitrap™ mass spectrometry (LD-O-MS) has the capability to determine the elemental and isotopic composition of solid planetary materials such as rock, regolith, organics, ice, etc. Here, we use a space-qualified LD-O-MS instrument comprising a 266 nm ultraviolet (UV) laser and Orbitrap mass analyzer ruggedized for planetary applications to investigate the stable isotopic composition of Ti and Zn metal plates. Based on the isotopic analyses, we constrain the performance of the CORALS instrument to sub per mille (‰) level accuracy and precision for Ti and at the per mille (‰) level for Zn. An LD-O-MS instrument is a versatile instrument capable of measuring isotopic composition of a variety of planetary samples and would constitute a critical instrument in the exploration of various planetary bodies including but not limited to Moon, Mars, Enceladus and other ocean or icy worlds, Ceres and other asteroids, and comets, thereby answering several high-priority questions pertaining to the formation and evolution of our Solar System.

58 GEOSCIENCES↗
DOE OSTI2022

Multitask Machine Learning of Collective Variables for Enhanced Sampling of Rare Events

Computing accurate reaction rates is a central challenge in computational chemistry and biology because of the high cost of free energy estimation with unbiased molecular dynamics. In this work, a data-driven machine learning algorithm is devised to learn collective variables with a multitask neural network, where a common upstream part reduces the high dimensionality of atomic configurations to a low dimensional latent space and separate downstream parts map the latent space to predictions of basin class labels and potential energies. Here, the resulting latent space is shown to be an effective low-dimensional representation, capturing the reaction progress and guiding effective umbrella sampling to obtain accurate free energy landscapes. This approach is successfully applied to model systems including a 5D Müller Brown model, a 5D three-well model, the alanine dipeptide in vacuum, and an Au(110) surface reconstruction unit reaction. It enables automated dimensionality reduction for energy controlled reactions in complex systems, offers a unified and data-efficient framework that can be trained with limited data, and outperforms single-task learning approaches, including autoencoders.

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

DECOVALEX-2023: An international collaboration for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems

The DECOVALEX initiative is an international research collaboration (www.decovalex.org), initiated in 1992, for advancing the understanding and modeling of coupled thermo-hydro-mechanical-chemical (THMC) processes in geological systems. DECOVALEX stands for “DEvelopment of COupled Models and VALidation against EXperiments”. The creation of this international initiative was motivated by the recognition that prediction of these coupled effects is an essential part of the performance and safety assessment of geologic disposal systems for radioactive waste and spent nuclear fuel. DECOVALEX emphasizes joint analysis and comparative modeling of the complex perturbations and coupled processes in geologic repositories and how these impact long-term performance predictions. The most recent phase of the DECOVALEX Project, here referred to as DECOVALEX-2023, started in early 2020 and ended in late 2023. More than fifty research teams associated with 17 international DECOVALEX partner organizations participated in the comparative evaluation of eight modeling tasks covering a wide range of spatial and temporal scales, geological formations, and coupled processes. This Virtual Special Issue on DECOVALEX-2023 provides an in-depth overview of these collaborative research efforts and how these have advanced the state-of-the-art of understanding and modeling coupled THMC processes. While primarily focused on radioactive waste, much of the work included here has wider application to many geoengineering topics.

Coupled processes↗
DOE OSTI2021

Defect-Driven Oxide Transformations and the Electrochemical Interphase

The redox reaction pathway is crucial to the sustainable production of the fuels and chemicals required for a carbon-neutral society. Our society is becoming increasingly dependent on devices using batteries and electrolyzers, all of which rely on a series of redox reactions. The overall properties of oxide materials make them very well suited for such electrochemical and catalytic applications due to their associated cationic redox properties and the static site-adsorbate interactions. As these technologies have matured, it has become apparent that defect-driven redox reactions, defect-coupled diffusion, and structural transformations that are both time- and rate-dependent are also critical materials processes. This change in focus, considering not only redox properties but also more complex, dynamic behaviors, represents a new research frontier in the molecular sciences as they are strongly linked to device operation and degradation and lie at the heart of various phenomena that take place at electrochemical interfaces. Fundamental studies of the structural, electronic, and chemical transformation mechanisms are key to the advancement of materials and technological innovations that could be implemented in various electrochemical systems. In this Account, we focus on recent studies and advances in characterizing and understanding the dynamic redox evolution and structural transformations that take place in model perovskites and layered oxides under reactive conditions and correlate them with degradation mechanisms and operations in electrolyzers and batteries. We show that the dynamic evolution of oxygen vacancies and cationic migration in the surface or bulk occurs at the solid-liquid interface, using a combination of different synchrotron-based Xray spectroscopies and scattering probes. Detailed redox-structure-reactivity correlation studies show how defects and diffusion processes can be tailored to drive various physical and chemical transformations in electrolyzers and batteries. We also highlight a strong correlation between oxygen redox reactivity and structural reorganization in both model thin films and particles, helping to bridge the gap between fundamental studies of the reaction mechanism and device applications. Finally, on the basis of these findings, we discuss strategies to probe and tune the redox reactivity and structural stability of the redox-active oxide interphase toward devising efficient pathways for energy and chemical harvesting.

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

Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software

The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).

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

Material analysis of metal hydrides for bulk hydrogen storage

The United States Department of Energy (DOE) Hydrogen Storage Engineering Center of Excellence (HSECoE), 2009–2014, was established to evaluate the viability of material-based hydrogen storage for light duty vehicles. The HSECoE concluded that there were no metal hydrides that simultaneously met the DOE targets for gravimetric and volumetric capacity, and rate of onboard refueling. However, it was noted by the DOE that potential applications existed for stationary hydrogen storage systems. An application of particular interest is using a metal-hydride to store a hydrogen supply for fuel cell backup power to a 5 MWe data center. A generalized numerical model that coupled mass and energy conservation, chemical kinetics and thermodynamics was developed to simulate the discharge process for a given metal-hydride. The model was applied to an interstitial and a complex metal-hydride. In conclusion, it was found that the interstitial metal-hydride, which is discharged by a pressure swing is preferrable in a practical sense to the complex metal-hydride, which requires a temperature swing discharge process.

08 HYDROGEN↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Electron–Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab-initio wavefunction calculations at the multi-reference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron – nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on non-relativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538–549] and therefore had only limited applicability. Furthermore, results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex, demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

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

An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast—Part 2: Circulation, Boundary Layer, and Clouds

The Western North Atlantic Ocean (WNAO) is a complex land-ocean-atmosphere system that experiences a wide range of atmospheric phenomena, which in turn drive unique aerosol transport pathways, cloud morphologies, and boundary layer variability. This work, Part 2 of a 2-part paper series, provides an overview of the atmospheric circulation, boundary layer variability, three-dimensional cloud structure, and precipitation over the WNAO; the companion paper (Part 1) focused on chemical characterization of aerosol, gas, and wet deposition. Seasonal changes in the circulation and sea surface temperature explain a clear transition in cloud morphologies from small shallow cumulus clouds, convective clouds, and tropical storms in summer, to stratus/stratocumulus and multi-layer cloud systems associated with winter storms. Synoptic variability in cloud fields is estimated using satellite-based weather states, and the role of postfrontal conditions (cold-air outbreaks) in the development of stratiform clouds is further analyzed. Precipitation is persistent over the ocean, with a regional peak over the Gulf Stream path, where offshore sea surface temperature gradients are large and surface fluxes reach a regional peak. Satellite data show a clear annual cycle in cloud droplet number concentration along the coast in winter, and minima in summer, suggesting a marked annual cycle in aerosol-cloud interactions. The performance of four climate models are assessed against satellite cloud retrievals. The models qualitatively reproduce the annual cycle in cloud cover and liquid water path, but with large discrepancies across models, especially in the extra-tropics. The paper concludes with a summary of outstanding issues and recommendations for future work.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2020

An extended methodology for automated calculations of non-Boltzmann kinetic sequences: H + C 2 H 2 + X and combustion impact

It is generally assumed in phenomenological kinetic models that bimolecular reactions only occur between species whose rovibrational energy follows a Boltzmann (thermal) distribution. That is, any complexes initially formed in non-Boltzmann distributions are assumed to be thermalized by energy-transferring collisions prior to bimolecular reactions. Given the high mole fractions of reactive species, X, in combustion environments, reactive collisions of the complexes with X often occur on the same timescale as energy-transferring collisions – yielding sequences proceeding through non-Boltzmann intermediates across multiple potential energy surfaces. Recent studies have shown that such non-Boltzmann kinetic sequences can have substantial impact on the global reactivity in combustion systems. Simulations of these non-Boltzmann reaction sequences, which can be described in phenomenological kinetic models via chemically termolecular reactions, require that rovibrational excitation from one potential energy surface be carried over to the next. Here this paper presents an extended theoretical and computational methodology that couples multiple master equations and derives rate constants for phenomenological reactions describing the conversion of thermal reactants to thermal products for use in phenomenological kinetic schemes. The methodology is then implemented using in-house scripts for non-Boltzmann sequences involving C2H$^*_3$ + X (with X = O 2 , H, and OH) where C 2 H$^*_3$ is formed via H + C 2 H 2 association – which were identified as having strong potential for influencing combustion predictions in a recent study. The results reveal that non-Boltzmann reaction sequences for X = O 2 (the primary focus of this paper) significantly alters the total conversion rate from H + C 2 H 2 to products and product branching fractions from those of thermal sequential pathways. Furthermore, the present results demonstrate that non-Boltzmann reaction sequences have significant impact – as high as an order of magnitude – on predicted ignition delay times. Similarly, they yield significantly different dependence of ignition delay times with temperature and O 2 mole fraction – yielding signatures that are likely observable experimentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗