SEARCH · Engineering Papers

Results for “density functional theory development”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,297 records · Page 72

All sources NASA NTRS DOE OSTI

DOE OSTI2022

13 C NMR study of amino acid salts in facilitated transport membranes for post-combustion carbon capture

Some of amino acid salts (AASs) have been demonstrated to be effective mobile carriers in facilitated transport membranes (FTMs) to achieve superior CO 2 permeance and CO 2 /N 2 selectivity for CO 2 capture from flue gas. Understandings of how the structures of different AASs affect the chemistry of the amine–CO 2 reaction are essential for the future development of more efficient AAS mobile carriers. In this study, the reaction chemistry of selected AASs with CO 2 was investigated by 13 C nuclear magnetic resonance (NMR) spectroscopy. The CO 2 loading and the distribution of major reaction products, including carbamate and bicarbonate products of each studied AAS were quantitatively analyzed. The positive correlation between the CO 2 loading of AAS (mol CO 2 /g AAS) and the CO 2 permeance of FTMs suggested that AAS with a higher CO 2 loading may improve the performance of FTMs. Our results also showed that increasing the steric hindrance of AAS could be a practical way to promote the bicarbonate reaction pathway and thus potentially increase the CO 2 loading (mol CO 2 /mol AAS). Moreover, 2-(1-piperazinyl)ethylamine (PZEA) was found to be a more effective multi-amine than piperazine (PZ) for synthesizing AASs due to the presence of more effective nitrogen sites per molecule. The superior CO 2 permeance of the PZEA-Sar membrane is also attributed to its more uniform membrane formation. As a result, the knowledge gained from this study will inform the rational design of more effective AAS carriers for CO 2 capture.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Role of the isovector spin-orbit potential in mitigating the CREX-PREX dilemma

Pioneering electroweak measurements of the neutron skin thickness in lead-208 and calcium-48 are challenging our understanding of nuclear dynamics. Many theoretical models suggest that the slope of the symmetry energy controls the development of a neutron skin in neutron-rich nuclei. This led to the expectation that if lead-208 exhibits a large neutron skin, calcium-48 should as well. Given that the PREX Collaboration reported a relatively thick neutron skin in lead, we anticipated that calcium would also have a significant neutron skin. Instead, the CREX Collaboration reported a thin neutron skin in calcium. Although many suggestions have been proposed, the “CREX-PREX dilemma” remains unsolved. Recently, an intriguing scenario has emerged, suggesting that an enhanced isovector spin-orbit interaction could simultaneously account for both results. Following this approach, we performed relativistic mean-field calculations with an increased isovector spin-orbit potential. Our findings indicate that, while this modification significantly affects the structure of calcium-48, it has only a marginal impact on lead-208, thereby bringing the results into better agreement with experiment. Furthermore, the strong enhancement required to mitigate the CREX-PREX dilemma destroys the agreement with a successful spin-orbit phenomenology, primarily by modifying the well-known ordering of spin-orbit partners.

Neutron skin thickness↗

DOE OSTI2024

How well do one-electron self-interaction-correction methods perform for systems with fractional electrons?

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z σ ) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew–Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of $H$$^{+}_{2}$ and $H$$^{+}_{2}$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy vs number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li + to F – at the expected distance from the experimental data and accurate ab initio data. Overall, both the PZSIC and LSIC methods reduce the delocalization errors substantially.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Machine learning Frenkel Hamiltonian parameters to accelerate simulations of exciton dynamics

In this manuscript, we develop multiple machine learning (ML) models to accelerate a scheme for parameterizing site-based models of exciton dynamics from all-atom configurations of condensed phase sexithiophene systems. This scheme encodes the details of a system’s specific molecular morphology in the correlated distributions of model parameters through the analysis of many single-molecule excited-state electronic-structure calculations. These calculations yield excitation energies for each molecule in the system and the network of pair-wise intermolecular electronic couplings. In this work, we demonstrate that the excitation energies can be accurately predicted using a kernel ridge regression (KRR) model with Coulomb matrix featurization. We present two ML models for predicting intermolecular couplings. The first one utilizes a deep neural network and bi-molecular featurization to predict the coupling directly, which we find to perform poorly. The second one utilizes a KRR model to predict unimolecular transition densities, which can subsequently be analyzed to compute the coupling. We find that the latter approach performs excellently, indicating that an effective, generalizable strategy for predicting simple bimolecular properties is through the indirect application of ML to predict higher-order unimolecular properties. Such an approach necessitates a much smaller feature space and can incorporate the insight of well-established molecular physics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Charged Defects in UO 2 Bulk and Surface: A First-Principles Study

Uranium dioxide (UO 2 ) is the primary fuel used in nuclear reactors. Under the extreme heat and radiation inside a reactor, this material inevitably develops defects in its crystal structure. To investigate the nature and behavior of these defects, DFT+U calculations were employed to investigate charged point defects in both bulk UO 2 and its most stable surface, the (111) plane. The formation of defects and their impact on the electronic structure were systematically examined. The results reveal that these defects introduce localized electronic states, alter magnetic behavior, and modify the structural properties. In general, such defects act as deep traps capable of capturing and retaining charge carriers. The stability of these defects depends strongly on the chemical environment and the position of the Fermi level. Surface defect calculations reveal that oxygen vacancies form more readily at the surface than in the bulk over a wide range of electron chemical potential, with subsurface oxygen vacancies being more stable than those in the top layer. Overall, the findings demonstrate how charged defects influence magnetism, transport, and stability in UO 2 , providing insights that may guide improvements in the safety and efficiency of nuclear fuel.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Magnetoelastic coupling and spin contributions to entropy and thermal transport in biferroic yttrium orthochromite

Direct engineering of material properties through exploitation of spin, phonon, and charge-coupled degrees of freedom is an active area of development in materials science. However, the relative contribution of the competing orders to controlling the desired behavior is challenging to decipher. In particular, the independent role of phonons, magnons, and electrons, quasiparticle coupling, and relative contributions to the phase transition free energy largely remain unexplored, especially for magnetic phase transitions. Here, we study the lattice and magnetic dynamics of biferroic yttrium orthochromite using Raman, infrared, and inelastic neutron spectroscopy techniques, supporting our experimental results with first-principles lattice dynamics and spin-wave simulations across the antiferromagnetic transition at T N ~ 138 K. Spectroscopy data and simulations together with the heat capacity (C p ) measurements, allow us to quantify individual entropic contributions from phonons (0.01 ± 0.01k B atom -1 ), dilational (0.03 ± 0.01k B atom -1 ), and magnons (0.11 ± 0.01k B atom -1 ) across TN. High-resolution phonon measurements conducted in a magnetic field show that anomalous T-dependence of phonon energies across TN originates from magnetoelastic coupling. Additionally, phonon scattering is primarily governed by the phonon–phonon coupling, with little contribution from magnon–phonon coupling, short-range spin correlations, or magnetostriction effects; a conclusion further supported by our thermal conductivity measurements conducted up to 14 T, and phenomenological modeling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗

DOE OSTI2025

Facet-dependent structure and dissociation of water at pristine IrO 2 /water interfaces

Understanding the microscopic structure of water at metal oxide interfaces is crucial for advancing electrocatalysis. IrO 2 , specifically, has shown exceptional activity for electrochemical water oxidation, but we currently lack a fundamental understanding of how the surface structure of IrO 2 impacts water reactivity. In this work, we developed a machine learning potential trained to first-principles accuracy for modeling IrO 2 /water interfaces across different facets: (110), (100), (101), and (001). Using extensive machine learning molecular dynamics simulations, we investigated the spontaneous dissociation of water molecules at these interfaces. Our results reveal a distinct dissociation probability trend: (110) > (100) ≈ (101) > (001), which we attribute primarily to the reaction thermodynamics of surface water dissociation. A strong correlation is observed between the surface Ir–O bond distances and the dissociation probabilities, highlighting the role of surface geometry in modulating reactivity. As a consequence, the interfacial solvation structures and hydrogen bonding environments are dynamically tuned by the varying water dissociation capabilities across facets. This work elucidates how water dissociation energetics depend on surface orientation and interfacial structure, offering atomistic insights into manipulating reaction chemistry at electrocatalytic interfaces.

organic↗

DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Boosting CO 2 reduction on Fe-N-C with sulfur incorporation: Synergistic electronic and structural engineering

Developing earth-abundant efficient catalysts for CO 2 reduction reaction (CO 2 RR) is of paramount importance for electrochemical conversion of CO 2 into value-added products. Despite numerous studies on iron and nitrogen codoped carbon (Fe-N-C) catalysts, grand challenges exist due to limited performance and understanding of catalytic mechanisms. This study reports a general strategy to boost electrocatalytic CO 2 RR activity of Fe-N-C with the incorporation of S atoms to engineer carbon support structure and electronic properties of active Fe-N sites simultaneously via a copolymer-assisted synthetic approach. The employment of N,S comonomers significantly increases the numbers of micropores and surface area, enabling dense atomic Fe-N and enhanced utilization efficiency. We report the first-principles calculations reveal that S modulation upraises the Fermi energy of Fe 3d and increases charge density on Fe atoms of Fe-N 4 , thereby enhancing intrinsic catalytic reactivity and selectivity for CO 2 reduction by strengthening the binding interaction between the Fe site and key COOH* intermediate. These integrated structural and electronic merits endow Fe-NS-C with outstanding activity (e.g., CO Faradaic efficiency of 98% at an overpotential of 490 mV) and stability (without deactivation in 30 h), ranking it one of the most active Fe-N-C reported to date. The finding offers an innovative design strategy to enable the design of advanced catalysts for CO 2 conversion.

25 ENERGY STORAGE↗

DOE OSTI2025

Defect Complexes in CrSBr Revealed Through Electron Microscopy and Deep Learning

Atomic defects underpin the properties of van der Waals materials, and their understanding is essential for advancing quantum and energy technologies. Scanning transmission electron microscopy is a powerful tool for defect identification in atomically thin materials, and extending it to multilayer and beam-sensitive materials would accelerate their exploration. Here, we establish a comprehensive defect library in a bilayer of the magnetic quasi-1D semiconductor CrSBr by combining atomic-resolution imaging, deep learning, and calculations. We apply a custom-developed machine learning work flow to detect, classify, and average point vacancy defects. This classification enables us to uncover several distinct Cr interstitial defect complexes, combined Cr and Br vacancy defect complexes, and lines of vacancy defects that extend over many unit cells. We show that their occurrence is in agreement with our computed structures and binding energy densities, reflecting the intriguing layer interlocked crystal structure of CrSBr. Our ab initio calculations show that the interstitial defect complexes give rise to highly localized electronic states. These states are of particular interest due to the reduced electronic dimensionality and magnetic properties of CrSBr and are, furthermore, predicted to be optically active. Our results broaden the scope of defect studies in challenging materials and reveal new defect types in bilayer CrSBr that can be extrapolated to the bulk and to over 20 materials belonging to the same FeOCl structural family.

deep learning↗

DOE OSTI2025

A New Family of Ternary Intermetallic Compounds with Dualistic Atomic Ordering – The ZIP Phases

A new family of nanostructured ternary intermetallic compounds − named the ZIP phases − is introduced in this work. The ZIP phases exhibit dualistic atomic ordering, i.e., they form two structural variants: one with the fcc diamond cubic structure (space group Fd$\bar{3}$m) and one with the hexagonal structure (space group P6 3 /mmc). They are also characterized by metallic behavior, ionic bonding, and atomic zigzagging. Powder metallurgical routes involving pressure-assisted densification are adopted to demonstrate ZIP phase synthesis in the Nb-Si-Ni, Nb-Si-Co, Ta-Si-Ni, V-Si-Ni, and Nb-Si-Fe ternary systems. Crucially, reactive hot pressing is capable of producing high-purity ZIP phase materials after the judicious, elemental system-specific optimization of the processing route. Synthesis of phase-pure materials – demonstrated in the Nb-Si-Ni ternary system by the synthesis of quasi phase-pure Nb 3 SiNi 2 and Ni 3 SiNb 2 ZIP phase-based materials – is a steppingstone to the prospective exploitation of the ZIP phases. Characterization of Nb 3 SiNi 2 and Ni 3 SiNb 2 involves crystal structure determination, spatially resolved chemical analysis, and determination of select thermal, electrical, magnetic, mechanical, and physical properties. Density functional theory is used to assess the stability of Nb 3 SiNi 2 & Ni 3 SiNb 2 and derivative binary compounds at different temperatures, also exploring the exfoliation of these two ZIP phases along specific surfaces to produce 2D derivatives.

Intermetallic compounds (IMCs)↗

DOE OSTI2025

Data-Driven Kinetic Reaction Networks for Separation Chemistry

Understanding complex, multistep chemical reactions at the molecular level is a major challenge whose solution would greatly benefit the design and optimization of numerous chemical processes. The separation of rare-earth (4f) and actinide (5f) elements is an example where improving our chemical understanding is important for designing and optimizing new chemistries, even with a limited number of observations. Here, in this work, we leverage data-driven artificial intelligence and machine-learning approaches to develop kinetic reaction networks that describe the liquid–liquid extraction mechanism of uranium using N,N-di-2-ethylhexyl-isobutyramide (DEHiBA). Specifically, we compare and contrast the properties of two classes of models: (1) purely data-driven models that are regularized using chemistry-agnostic, L1 regression and (2) chemistry-informed models that are regularized using relative reaction energies provided by quantum mechanical calculations. We observe that purely data-driven models are unbiased, simple, and accurate in their predictions of experimental measurements when provided with sufficient data but are difficult to fully constrain and interpret. In contrast, chemistry-informed models exhibit significantly improved chemical interpretability and consistency, providing a detailed description of the separation process while achieving high accuracy through ensemble averaging. Overall, the dominant species predicted to be extracted into the organic phase is UO 2 (NO 3 ) 2 (DEHiBA) 2 , agreeing with experimental slope analysis, thermodynamic modeling, EXAFS, and crystal structures. This work demonstrates that leveraging the fundamental structure of the problem can lead to efficient learning schemes that provide both accurate predictions and chemical insights at a low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Promoting the cleavage of C-O bonds at the interface between a metal oxide cluster and a Co(0001) support

As a first step toward the rational design of Co-based catalysts with a higher activity and selectivity, we determine how one can activate a C-O bond at the interface between a metal oxide cluster and a Co(0001) support. The hypothesis here is that the metal ions in metal oxide clusters on a Co(0001) support enhance the adsorption of CO and weaken the C-O bond strength, which can then facilitate the dissociation of the CO reactant. To test this hypothesis, we developed three computational models of Ti4O8/Co(0001), Zr4O8/Co(0001), and Mn8O8/Co(0001). We quantify the CO adsorption behavior at the interface sites between an oxide cluster and the Co(0001) support as well as the corresponding IR spectra. We correlate the computed CO stretch frequencies with their CO adsorption energies, as well as the CO stretch frequency with the C-O bond length, and related these findings to the changes in the chemical bonding in the bound CO. The interface is the most favorable site for CO adsorption. Adsorption results in an increase of the C-O bond length and a decrease in its vibrational frequency. From a chemical bonding analysis, the bond order in CO at this site drops from 3 (in the gas phase) to 1. This decrease in bond order is a necessary precursor stage for CO dissociation. The experimental measurements of the corresponding FTIR spectra support this point. The favorability of CO adsorption at the interface sites is due to an electron transfer from the metal ion in the metal oxide cluster to the O atom in CO. We establish a linear relationship between C-O bond length and CO frequency and this relationship is found to be independent of the support, type of metal oxide cluster, or the adsorption site.

Fischer-Tropsch Synthesis, Heterogeneous Catalysis↗

DOE OSTI2020

High Temperature Erosion Modeling in Particle Based CSP Systems

1. Introduction. Wear and erosion damage of materials from solid particle and surface interactions is a major issue in various industries. Although more common in mining and oil and gas production, erosion is becoming a critical issue in renewable technologies as well such as particle based concentrated solar thermal power (CSP) systems. In particle based CSP systems, solid particles are used to absorb solar energy and as thermal storage. However, these particles may cause significant amount of wear to system components while through the system at high temperatures. This damage can be costly, and therefore, requires a greater understanding of solid particle erosion in CSP systems. Throughout the years, models and tools have been developed to predict and control erosion in industries such as oil and gas production. However, these erosion models and erosion prediction tools have been mainly developed based on erosion data for much higher velocities and lower temperatures, i.e. the operating conditions for which the existing models have been built are not pertinent to those expected in CSP system. It is known that erosion depends on many factors such as material properties, erodent particle properties, and particle impact speed and angle. However, the effect of temperature on erosion is not vastly investigated. In this work, a temperature-based correlation is introduced that will be used to modify the existing erosion models to predict erosion rates at velocities and temperatures relevant to Gen3 CSP systems. The existing models are first validated against erosion experiments run at low temperatures and low velocities. Subsequently, a temperature correction term is developed that can extend the existing models to high temperatures, based on the available experimental data showing the effect of temperature on impact erosion. 2. Erosion Models. Continuing sequence of impacts from solid particles on surfaces would result in loss of material due to mechanical interaction between solid surface and particles. Erosion equations have been developed to predict erosion under different conditions and erosion mechanism, including cutting and deformation erosion. These models are mainly either mechanistic, empirical, and semi-mechanistic models, when the latter combine the theory of the erosion mechanism and particle motion with the available experimental data. One of the first empirical correlation was introduced by American Petroleum Institute (API) Recommended Practice (RP) 14E [1]. This correlation was very conservative, and several improved empirical and semi-mechanistic equations were developed in the years following. More recently, the models introduced by Erosion/Corrosion Research Center (E/CRC) at the University of Tulsa and Arabnejad. et al. [2] are frequently used in the literature and industry, as they account for many parameters affecting erosion including particle impact speed and angle, material density and hardness, and particle size and shape. However, both these models were developed based on data at relatively high velocities and also do not take into account any mechanical changes in material associated with thermal cycling at high temperatures. 3. Erosion Prediction at High Temperatures. In this work, the accuracy of Arabnejad et al. model and E/CRC model is investigated to predict impact erosion at conditions relevant to Gen3 CSP systems. Two significant deviations expected in Gen3 systems compared to operating conditions used to build these models are particle impact velocities and system temperature. As a first step, the performance of these models was validated at low velocities and temperatures. Impact erosion experiments were conducted on SS316 coupons using HSP 40/70 ceramic particles. Assuming particle rate of 1 kg/s/m, 0.0254 m of particle curtain thickness, and particle velocity of 1.5 m/s in the system, an overall erosion of 0.59 mm/year was calculated. The experimental results were subsequently compared to those from computational simulations and erosion of 0.47 mm/year was obtained using the Arabnejad et al. model. The results indicated that the models work well at low velocities and low temperature conditions. To account for temperature effect, a mathematical correlation was developed using data published by DUCOM [3] for Inconel 600 eroded by alumina particles at high velocities. The correlation was applied to both Arabnejad et. al model and the E/CRC model. The prediction results from these modified models were within 20% of this experimental data. Predictions of of erosion by the modified E/CRC model at three different temperatures are made. Similar to the calculation at low temperature, a typical CSP system with 1 kg/s/m of particle flow rate per unit length of the particle curtain and a curtain thickness of 0.0254 m is assumed. The annual thickness loss calculations were run for different particle impact velocities, assuming a uniform particle impact area equal to the cross-sectional area of the curtain (particle-particle interactions and dispersion of particles are not considered in obtaining the results). It is observed that erosion increases exponentially as temperature increases. Furthermore, it is also noted that, erosion changes non-linearly with impact velocity. Currently, experiments are also being conducted to measure erosion of SS316 at 800 ºC with HSP 40/70 particles at low impact velocities. We expect to use the results from high temperature testing to further improve the temperature correlation function. Similar models are also being developed for abrasion erosion resulting from particle sliding along the surfaces as well as attrition from particle to particle and particle to surface interactions. References. [1] Institute, A. P. (1991). API Recommended Practice for Design and Installation of Offshore Production Platform Piping System, API RP 14E. [2] Arabnejad, H., Mansouri, A., Shirazi, S. A., and McLaury, B. S. (2015a). Development of mechanistic erosion equation for solid particles. Wear, 332–333, 1044–1050. http://doi.org/10.1016/j.wear.2015.01.031. [3] https://ducom.com/high-temperature-erosion-evaluating-sample-wear/

14 SOLAR ENERGY↗

DOE OSTI2023

Synthetic Access to a Framework-Stabilized and Fully Sulfided Analogue of an Anderson Polyoxometalate that is Catalytically Competent for Reduction Reactions

Polyoxometalates (POMs) featuring 7, 12, 18, or more redox-accessible transition metal ions are ubiquitous as selective catalysts, especially for oxidation reactions. The corresponding synthetic and catalytic chemistry of stable, discrete, capping-ligand-free polythiometalates (PTMs), which could be especially attractive for reduction reactions, is much less well developed. Among the challenges are the propensity of PTMs to agglomerate and the tendency for agglomeration to block reactant access of catalyst active sites. Nevertheless, the pervasive presence of transition metal sulfur clusters metalloenzymes or cofactors that catalyze reduction reactions and the justifiable proliferation of studies of two-dimensional (2D) metal-chalcogenides as reduction catalysts point to the promise of well-defined and controllable PTMs as reduction catalysts. Here, we report the fabrication of agglomeration-immune, reactant-accessible, capping-ligand-free Co II Mo 6 IV S 24 n– clusters as periodic arrays in a water-stable, hierarchically porous Zr-metal–organic framework (MOF; NU1K) by first installing a disk-like Anderson polyoxometalate, Co III Mo 6 VI O 24 m– , in size-matched micropores where the siting is established via difference electron density (DED) X-ray diffraction (XRD) experiments. Flowing H 2 S, while heating, reduces molybdenum(VI) ions to Mo(IV) and quantitatively replaces oxygen anions with sulfur anions (S 2– , HS – , S 2 2– ). DED maps show that MOF-templated POM-to-PTM conversion leaves clusters individually isolated in open-channel-connected micropores. Importantly, the structure of the immobilized cluster as determined, in part, by X-ray photoelectron spectroscopy (XPS), X-ray absorption fine structure (XAFS) analysis, and pair distribution function (PDF) analysis of total X-ray scattering agrees well with the theoretically simulated structure. PTM@MOF displays both electrocatalytic and photocatalytic competency for hydrogen evolution. Nevertheless, the initially installed PTM appears to be a precatalyst, gaining competency only after the loss of ~3 to 6 sulfurs and exposure to hydride-forming metal ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Interpretation of Ion Irradiation and Neutron Irradiation Damage in Additively Manufactured 316 Stainless Steel using Multiscale Modeling

The accelerated adoption of nuclear energy necessitates advanced manufacturing technologies, such as additive manufacturing, to meet heightened supply chain requirements and support innovative reactor technologies. Due to the unique microstructural characteristics of additively manufactured materials under distinct solidification conditions, comprehensive evaluation of their performance in reactor environments is essential. The Advanced Materials and Manufacturing Technologies program under the Department of Energy's Office of Nuclear Energy focuses on understanding the irradiation performance and damage evolution of laser powder bed fusion 316 stainless steel, with an emphasis on integrating ion and neutron irradiation data to accelerate the development and qualification of materials for advanced nuclear reactor applications. While ion irradiation is a cost- and time-effective method, modeling and simulation are required to interpret the data for the broader range of irradiation conditions encountered in advanced reactors. In fiscal year 2025, integrated multiscale modeling and simulations were conducted to assess irradiation damage in additively manufactured 316 stainless steel. Key outcomes include predictions of chromium enrichment at grain boundaries, nickel enrichment at dislocation cell walls and void surfaces, and heterogeneous void evolution under ion and neutron irradiation conditions. Cluster dynamics simulations revealed the coarsening of voids at high irradiation temperatures and the suppression of void growth by high network dislocation density, while also demonstrating significant growth and coarsening of voids and self-interstitial atom loops at low dose rates. Machine learning-accelerated atomistic simulations highlighted the impact of the local environment and chromium concentration on vacancy diffusivity, providing key insights on the influence of composition on void swelling and radiation-induced segregation. Additionally, molecular dynamics simulations demonstrated the presence of defect production bias and a significant effect of carbon content on defect cluster behavior. These combined efforts aim to predict the performance of additively manufactured materials under various reactor conditions, supporting their qualification for nuclear reactor applications by interpreting ion irradiation data. This report underscores the potential of integrated multiscale modeling to analyze ion irradiation data in the effort to accelerate the qualification of additively manufactured materials for nuclear reactor components.

316 stainless steel↗