Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Monte Carlo”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

CASM — A software package for first-principles based study of multicomponent crystalline solids

CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Monte Carlo simulations of Cu/Ni–Si–Mn co-precipitation in duplex stainless steels

First-Passage Kinetic Monte Carlo (FPKMC) simulations of species migration in duplex stainless steels were performed in order to establish relationships between alloy content, segregation behavior, and the formation of Ni-Si-Mn rich particles in cast duplex stainless steels during thermal aging. The Ni-Si-Mn-rich second phase forms after extended aging at reactor operating temperatures and degrades alloy properties. Simulations of Ni-Si-Mn cluster formation were validated through comparison with experimental results obtained through Atom Probe Tomography (APT) on similar alloy compositions, identifying several trends. First, Cu promotes the formation of Ni-Si-Mn clusters, but only when Ni or Mn prefer segregation to the surface of Cu particles; without the segregation of these species, the critical composition for the clusters to form was not achieved. Second, the width of precipitate-denuded zones near ?/d interfaces increases with decreasing Cu content. This finding was in strong agreement with APT results, further validating our model. Finally, our model predicts that Ni-Si-Mn cluster formation will be the most extensive when the Si:Mn ratio is approximately 1:1 and the least extensive when one of these key elements is less concentrated (Mn = 0.5 at.% or Si = 0.5 at.%). The implications of these results on how to improve alloy properties are discussed.

Cu co-precipitation, Atom probe tomography (APT), ↗
DOE OSTI2021

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

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

Graph-theoretical KMC Framework for Calculating Effective Diffusivity in TPBAR Components: Effective Diffusivity of Tritium in α-Zr Grain Boundaries

We report the development of a computer simulation tool based on the Graph-theoretical kinetic Monte Carlo (GT-KMC) approach that can be used to simulate diffusion and derive effective diffusivity of species in complex structures of TPBAR components. The GT-KMC framework was implemented in AKSOME, an on-lattice self-learning kinetic Monte Carlo tool developed at PNNL and was previously used to study solute diffusion in metal alloys. This newly developed tool will be helpful in simulating atom diffusion and extracting their effective diffusivities in grain boundaries and interfaces in various TPBAR components. Once benchmarked, the tool will be used to simulate the diffusion of tritium (hydrogen) along grain boundaries in α-Zr using the activation energy barrier data previously obtained from Density Functional Theory calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomistic modeling of meso -timescale processes with $\mathrm{SEAKMC}$: A perspective and recent developments

On-the-fly kinetic Monte Carlo (kMC) methods have recently garnered significant attentions after successful applications to various atomic-scale problems using a timescale outside the reach of classical molecular dynamics. These methods play a critical role in modeling atomistic meso-timescale processes, and it is therefore essential to further improve their capabilities. In this report we review one of the on-the-fly kMC methods, Self-Evolving Atomistic kinetic Monte Carlo (SEAKMC) and propose two schemes that considerably enhance the efficiency of saddle point searches (SPSs) during the simulations. The performance of these schemes is tested using the diffusion of point defects in bcc Fe. In addition, we discuss approaches to significantly mitigate limitations of these schemes, which further improves their efficiencies. Importantly, these schemes improve the SPS efficiency not only for SEAKMC but also for other on-the-fly kMC methods, broadening the applications of on-the-fly kMC simulations to complex meso-timescale problems.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Precipitation-site competition in duplex stainless steels: Cu clusters vs spinodal decomposition interfaces as nucleation sites during thermal aging

Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.

Duplex stainless steel, Thermal aging, Precipitati↗
DOE OSTI2021

Toward alcohol synthesis from CO hydrogenation on Cu(111)-supported MoS 2 – predictions from DFT+KMC

In the quest for cheap and efficient catalysts for alcohol synthesis from syngas, a material of interest is single-layer MoS2 owing to its low cost, abundancy, and flexible structure. Because of the inertness of its basal plane, however, it is essential to find ways that make it catalytically active. Herein, by means of density functional theory based calculations of reaction pathways and activation energy barriers and accompanying kinetic Monte Carlo simulations, we show that while S vacancy row structures activate the MoS 2 basal plane, further enhancement of chemical activity and selectivity can be achieved by interfacing the MoS 2 layer with a metallic support. When defect-laden MoS2 is grown on Cu(111), there is not only an increase in the active region (surface area of active sites) but also charge transfer from Cu to MoS 2 , resulting in a shift of the Fermi level such that the frontier states (d orbitals of the exposed Mo atoms) appear close to it, making the MoS 2 /Cu(111) system ready for catalytic activity. Finally, our calculated thermodynamics of reaction pathways lead to the conclusion that the Cu(111) substrate promotes both methanol and ethanol as the products, while kinetic Monte Carlo simulations suggest a high selectivity toward the formation of ethanol.

2D materials↗
DOE OSTI2024

Source of Processable Vitrimer Viscosities: Swap Frequencies and Steric Factors

Vitrimers exhibit high, processable viscosities, where other polymers do not, and are among the most promising polymers for closed-loop material circularity. We sought to investigate the underlying chemical kinetic factors that result in high viscosities for vitrimers, which are crucial to designing vitrimers with tunable viscosity. To interrogate these factors, we achieved the first simulated predictions of real vitrimer viscosities, using a novel kinetic Monte Carlo molecular dynamics method, overcoming the time and length scale gaps to predict experimental bulk viscosities. The vitrimer architecture investigated is based on poly(dimethylsiloxane) chains and vinylogous urethane bond swaps. We probed the effects of the extent of free swapping groups, %F, the activation energy, E A , and the steric factor, ρ. The steric factor is related to the intrinsic reaction probability for molecules with sufficient energy. All three factors were found to be significant, but the role of ρ was found to be the biggest and also the most underappreciated. The results show that the inclusion of accurate ρ is of critical importance for viscosity predictions, with the evidence suggesting that the typical assumption of ρ = 1 is not valid for vitrimers and that, indeed, very low steric factors are present in bond-swap vitrimers such that values of ρ < 10 –10 may be typical. This greatly influences the bond exchange rates and, ultimately, the viscosities. Recognition of this result is necessary for the prediction of vitrimer viscosities from molecular simulations and to make vitrimers by design from molecular dynamics. We also investigated the effects that E A , ρ, and the number of free swapping groups have upon vitreous range temperatures, TV, with respect to achieving a specific viscosity (η V = 1 × 10 8 Pa·s), as well as for a commonly reported higher viscosity extrapolation (η V = 1 × 10 12 Pa·s). The evidence suggests that vitrimers may follow universal curves for E A vs T V , as a function of ρ. Finally, this study achieves the first of these comparisons of molecular simulations to experiments and reveals critical insights toward creating vitrimers by design, while providing a route for the prediction of T V from kinetic Monte Carlo molecular dynamics simulations.

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

Scaling kinetic Monte-Carlo simulations of grain growth with combined convolutional and graph neural networks

Graph neural networks (GNN) have emerged as a promising machine learning method for microstructure simulations such as grain growth. However, accurate modeling of realistic grain boundary networks requires large simulation cells, which GNN has difficulty scaling up to. To alleviate the computational costs and memory footprint of GNN, we suggest a hybrid architecture combining a convolutional neural network (CNN) based bijective autoencoder to compress the spatial dimensions, and a GNN that evolves the microstructure in the latent space of reduced spatial sizes. Our results demonstrate that the new design significantly reduces computational costs with using fewer message passing layer (from 12 down to 3) compared with GNN alone. The reduction in computational cost becomes more pronounced as the spatial size increases, indicating strong computational scalability. For the largest mesh evaluated (160 3 ), our method reduces memory usage and runtime in inference by 117× and 115×, respectively, compared with GNN-only baseline. More importantly, it shows higher accuracy and stronger spatiotemporal capability than the GNN-only baseline, especially in long-term testing. Such combination of scalability and accuracy is essential for simulating realistic material microstructures over extended time scales. The improvements can be attributed to the bijective autoencoder’s ability to compress information losslessly from spatial domain into a high dimensional feature space, thereby producing more expressive latent features for the GNN to learn from, while also contributing its own spatiotemporal modeling capability. Training data are generated from stochastic grain growth simulations, providing realistic variability for learning robust microstructure evolution. Comprehensive system validation confirms that the model is accurate, robust, and scalable.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Catalytic resonance theory for the kinetic signatures of multiple wavelength photocatalysis

Multi-step surface chemistry promoted with one or two wavelengths of light exhibited distinct kinetic signatures indicative of the light-adsorbate interactions and the number of photon-sensitive elementary steps. In this work, kinetic Monte Carlo simulations identified the kinetic response of general surface mechanisms to variation in the per-site photon flux of one or two wavelengths of incident light. Photocatalytic rates were described via a non-dimensional photon flux, identifying multiple kinetic regimes with unique degrees of rate control. Under photon-controlled kinetic conditions, maximum quantum yield and turnover frequency were obtained under constant illumination, while pulsing of one or more light sources was shown to exhibit slower rates and less efficient photocatalytic promotion due to the inherent dynamic nature of light, which comprises a stream of photons. Simulations provided distinct kinetic signatures in catalytic rates and Arrhenius kinetics for specific light-surface-adsorbate interactions corresponding to photocatalytic promotion of specific steps in surface chemistry.

catalytic resonance↗
DOE OSTI2023

Understanding grain boundary segregation under irradiation: A mesoscale modeling study

Grain boundaries (GBs) are an important component in the design of irradiation resistant materials as they play a crucial role as point defect sinks for non-equilibrium vacancies and self-interstitials. However, in the presence of sustained point defect fluxes, the coupled transport of atomic components with inherently different diffusivities often leads to radiation-induced segregation (RIS), with undesired effects on properties such as intergranular corrosion and stress-assisted cracking. In this talk, we investigate the fundamental mechanisms leading to GB segregation in BCC FeCr-based alloys using multiscale modeling and experiments. A phase-field model is employed to account for both thermodynamic and kinetic segregation mechanisms. The full Onsager transport matrix is derived using kinetic Monte Carlo, molecular dynamics, and kinetic cluster expansion (KineCluE) calculations. The thermodynamics of GB segregation is described via an atomic density-based modification to the CALPHAD free energy. The modeling predictions are then compared against segregation data from atom probe characterization for various temperatures, compositions, and neutron irradiation conditions. We will conclude by discussing the utility of the multiscale modeling approach to shed light on the fundamental mechanisms leading to solute-GB interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Thermodynamics and kinetics of interconversion between platinum nanoparticles and cations in zeolites

In metal-containing zeolites, sintering and redispersion processes exercise control over the identity of metal site structures, but the thermodynamic and kinetic factors that influence these molecular-level processes are not completely understood. Furthermore, we assess the ability of first principles informed free energy models (for supported and unsupported nanoparticles) and kinetic models integrating Ostwald ripening with atom trapping to describe the interconversion between Pt cations and nanoparticles encapsulated in chabazite (CHA) zeolites. Density functional theory-derived thermodynamic phase diagrams show that the interconversion between cations, favored in oxidizing environments, and particles, favored in reducing environments, is fully reversible within a wide range of their respective conditions (temperatures and pressures) for CHA and several other zeolite topologies. Kinetic Monte Carlo simulations of Pt redispersion are consistent with experimentally observed redispersion kinetics of encapsulated Pt nanoparticles in CHA zeolites, and model results suggest the zeolite host imparts additional stability for Pt nanoparticles. We envision our thermodynamic and kinetic models for Pt-CHA are also capable of describing nanoparticle and cation interconversions for other zeolite frameworks under similar conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Computational Investigation of the Role of Active Site Heterogeneity for aSupported Organovanadium(III) Hydrogenation Catalyst

A crucial consideration for supported heterogeneous catalysts is the nonuniformity of the active sites, particularly for supported organometallic catalysts. Standard spectroscopic techniques, such as X-ray absorption spectroscopy, reflect the nature of the most populated sites, which are often intrinsically structurally distinct from the most active catalytic sites. Additionally, with computational models, often, only a few representative structures are used to depict catalytic active sites on a surface, even though there are numerous observable factors of surface heterogeneity that contribute to the kinetically favorable active species. A previously reported study on the mechanism of a surface organovanadium-(III) catalyst [(SiO 2 )V III (Mes)(THF)] for styrene hydrogenation yielded two possible mechanisms: heterolytic cleavage and redox cycling. These two mechanistic scenarios are challenging to differentiate experimentally since the kinetic readouts of the catalyst are identical. To showcase the importance of modeling surface heterogeneity and its effect on catalytic activity, density functional theory (DFT) computational models of a series of potential active sites of [(SiO 2 )V III (Mes)(THF)] for the reaction pathways are applied in combination with kinetic Monte Carlo (kMC) simulations. Computed results were then compared to the previously reported experimental kinetic study: (1) DFT free-energy reaction pathways indicated the likely active site and pathway for styrene hydrogenation, a heterolytic cleavage pathway requiring a bare tripodal vanadium site. (2) From the kMC simulations, a mixture of different bond lengths from the support oxygen to the metal center was required to qualitatively describe the experimentally observed kinetic aspects of a supported organovanadium(III) catalyst for olefin hydrogenation. This work underscores the importance of modeling surface heterogeneity in computational catalysis.

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

Catalytic resonance theory for the kinetics of photon-promoted catalysis

The illumination of catalytic surfaces with a continuous or pulsed stream of photons dynamically modulates surface chemistry for faster rates, non-equilibrium conversion, or product selectivity control. To establish fundamental principles of dynamic photon-modulated catalysis, the photocatalytic conversion of a generic surface reaction was simulated using the kinetic Monte Carlo method to understand the kinetic implications of an independent stream of photons that promotes surface product desorption. The time-averaged photocatalytic rate at differential conditions for varying photon flux and temperatures indicated three kinetic regimes described by product thermal desorption control, surface reaction control, and an intermediate kinetic regime with a zero slope Arrhenius plot, consistent with a degree of rate control dominated by the photon arrival frequency (i.e., per-site photon flux). Here, the maximum photocatalytic rate occurred orders of magnitude above the Sabatier limit at the resonance frequency, identified as the photon arrival frequency matching the surface reaction rate constant.

Canavan, Jesse R. [University of Minnesota, Minnea↗
DOE OSTI2024

Dynamic data-driven multiscale modeling for predicting the degradation of a 316L stainless steel nuclear cladding material

Here, we have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamic Data-Driven Multiscale Modeling for Predicting Microstructure and Mechanical Properties of 316L Stainless Steel

We have developed a long short-term memory stacked ensemble (LSTM-SE) surrogate modeling approach that can provide rapid predictions of microstructural evolution and the resultant mechanical properties of American Iron and Steel Institute (AISI) 316L series stainless steel (316LSS) fuel cladding under conditions of varying temperature and radiation dose rate. To acquire training data, we developed and implemented a kinetic Monte Carlo (KMC) model to simulate precipitation kinetics of M 23 C 6 , γ', and G phases within SS316L cladding. Experimentally reported precipitation kinetics of SS316L in literature were linked to the kinetic parameters of the simulated precipitation in our KMC model. The model was then used to simulate microstructure evolution under synthetically generated treatments of varying temperature and radiation dose rate, for periods of up to 3000 hours. Changes in volume fraction, number density, and particle size of precipitates were recorded, and particle area fractions were correlated using statistical methods to develop the surrogate model. Simultaneously, the mechanical properties of the simulated microstructures were evaluated using microstructure-based finite element method (FEM) analysis to determine the elastic modulus, yield stress, ultimate tensile strength, and elongation to failure of the aged microstructures. Using this approach, our surrogate model can predict precipitation behavior within 0.25% volume fraction and mechanical properties within 6% relative error from the values predicted by the KMC and FEM models using 50 training simulations as input. The trained recurrent neural network-based model can return estimations of precipitation kinetics and mechanical properties ~1000 times faster than the physics-based codes. This work demonstrates, as a proof of concept, that reactor material service lifetimes under variable service conditions can be predicted for a statistics-based model from a practicably obtainable dataset.

36 MATERIALS SCIENCE↗
DOE OSTI2024

New insights on the effects of chemistry and temperature on α’ precipitation during aging of FeCrAl alloys

Here, the effect of Al on α' precipitation in FeCrAl alloys was studied through thermal aging of several model binary FeCr and FeCrAl alloys with Cr content of 13 wt.%, 17 wt.%, and 25 wt.% with and without 5.5 wt.% Al. Aging was performed up to 1,000h at varying temperatures from 400-500°C. At 400°C no age hardening was observed due to slow kinetics at this temperature. For both the 17Cr and 25Cr alloys, the addition of Al shows a lowering of the miscibility gap, consistent with other reports in literature. Interestingly, however, for the 25Cr alloys the addition of Al in the FeCrAl ternary alloy accelerated α' precipitation below the miscibility gap. Such enhanced precipitation has also been predicted by our atomistic kinetic Monte Carlo (AKMC) simulations. While previous literature has often focused on Al suppressing α', here we show that while Al can lower the miscibility gap in the Fe-Cr ferritic system thermodynamically, it may also enhance the kinetics of precipitation.

36 MATERIALS SCIENCE↗