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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Bulk and surface diffusion of neodymium in alpha-uranium: Ab initio calculations and kinetic Monte Carlo simulations

A fundamental understanding of lanthanide transport in metallic fuels is critical for high fidelity modeling of the fuel-cladding chemical interaction (FCCI) phenomenon, which can lead to the formation of brittle intermetallic compounds and premature failure of the cladding. Here we report a combined ab initio density functional theory (DFT) and kinetic Monte Carlo (KMC) study of the bulk diffusivity of Nd in α-U, fully taking into account the effect of radiation enhanced diffusion. In this work, the vacancy mechanism is considered to be the dominant mechanism for the bulk diffusion of Nd since a Nd interstitial is found to be intrinsically unstable in α-U. The surface diffusivity of a Nd adatom on α-U (001) surface has been further predicted using KMC simulations parameterized by DFT calculations. The present study suggests that Nd transport via the surface diffusion mechanism can be many orders of magnitude faster than bulk diffusion. Furthermore, the results from the present lower length scale study can be used to inform mesoscale phase-field simulations to determine the effective diffusion coefficient of Nd through α-U with a porous microstructure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Exploring driving forces for length growth in graphene nanoribbons during chemical vapor deposition of hydrocarbons on Ge(0 0 1) via kinetic Monte Carlo simulations

Graphene grown slowly on Ge(0 0 1) using chemical vapor deposition of hydrocarbons leads to high-aspect ratio graphene nanoribbons with smooth edges and a technologically relevant band gap at room temperature; however, the driving forces leading to one-dimensional growth of such graphene crystals are not well understood. Here, we combine a lattice kinetic Monte Carlo approach based on steps in graphene growth and experimental measurements to study the growth of graphene nanoribbons via chemical vapor deposition on Ge(0 0 1). To identify potential reasons for growth of graphene as anisotropic ribbons, we study the impact of anisotropy in various growth parameters on the resulting graphene crystals. Comparing our model with experimental measurements indicates that anisotropy in the stabilization of a graphene precursor species bound to the graphene edge is the most likely reason why high aspect ratio graphene ribbons with smooth edges grow on Ge(0 0 1). Using the growth model developed here, we reproduce experimental trends in the synthesis of graphene nanoribbons on Ge(001) and arrive at an intuitive picture for their growth. These insights shed light on the driving forces governing this highly anisotropic regime of crystal growth.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of Accelerated Kinetic Monte Carlo Code for Simulation of Helium Bubble Evolution

A mesoscale model to predict helium bubble evolution is needed for tritium applications. Such a model requires that the conventional kinetic Monte Carlo (kMC) simulations be significantly accelerated. The objective of this report is to (a) highlight the concepts and mathematical expressions of the accelerated method for defect implementation that have not been published, (b) show an example input file to run the kMC code, and (c) provide suggestions on future improvement following my retirement.

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

Self-Learning Kinetic Monte Carlo Simulations of Radiation Damage in Nuclear Fuels

Understanding how irradiation affects the thermo-physical and mechanical properties of nuclear materials, such as thermal conductivity degradation in fuels and embrittlement of structural components, is critical to the safety and efficiency of nuclear reactors. These effects are largely governed by the formation and evolution of atomic-scale point defects and defect clusters. Due to their small sizes, however, these defects are invisible under high-resolution scanning transmission electron microscopy. This project aims to fill this experimental knowledge gap by integrating density functional theory (DFT), machine learning interatomic potential (MLIP), and kinetic Monte Carlo (KMC) techniques to predict longtime evolution of irradiation-induced defects in nuclear fuels.

36 - MATERIALS SCIENCE↗
DOE OSTI2023

Electrochemical Degradation of Pt 3 Co Nanoparticles Investigated by Off-Lattice Kinetic Monte Carlo Simulations with Machine-Learned Potentials

In fuel cell applications, the durability of catalysts is critical for large-scale industrial implementation. However, limited synthesis controllability and spectroscopic resolution impede a comprehensive understanding of degradation mechanisms at the atomic level. In this study, we develop a machine-learned potential (MLP) to simulate the degradation processes for Pt3Co nanoparticles. The precision of MLP is determined to be comparable to that of density functional theory calculations. Using off-lattice kinetic Monte Carlo simulations with MLP, we successfully replicate established experimental trends and offer a logical resolution to ongoing debates regarding atomic orderings. Based on the simulation results, we suggest design principles for Pt3Co nanoparticles that combine high activity and durability. Finally, we validate the wide applicability of our method by successfully applying it to Pt3Ni and Pt3Co0.5Ni0.5 nanoparticles. The research serves as a guideline for developing MLPs for alloy electrochemical catalysts and lays the foundation for designing more durable and active fuel-cell catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Surface Dynamics of Selective Hydration of Rutile TiO 2 : A Kinetic Monte Carlo Approach

The utility of atomic layer deposition (ALD) continues to expand beyond conformal thin-film fabrication to include area- or site-selective ALD. We previously identified a strategy for site-selective ALD (SS-ALD) via the evaluation of thermodynamically preferred hydration of rutile TiO 2 surfaces, as deduced from electronic structure calculations. Here, we report a novel kinetic Monte Carlo (KMC) model that allows for the investigation of surface dynamics and kinetics that improves our understanding of and intuition for the selective hydration strategy. Here we demonstrate the validity of the strategy with respect to step-edge defects for the lowest energy (110) facet as well as report results for the other common facets which agree with experimental STM observations. The results here indicate that the selective hydration strategy is feasible both thermodynamically (evaluated in our previous publication) and kinetically (from the KMC model). Because diffusion has a slower rate than others, we find that any proximity effects between terrace and defect sites are unlikely to affect the selective hydration strategy for rutile TiO 2 . The KMC model further provides relevant timescales for achieving selectivity experimentally and establishes the kinetic viability of the selective hydration approach to SS-ALD.

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

SEI Formation and Lithium-Ion Electrodeposition Dynamics in Lithium Metal Batteries via First-Principles Kinetic Monte Carlo Modeling

The stabilization and enhanced performance of lithium metal batteries (LMBs) depend on the formation and evolution of the Solid Electrolyte Interphase (SEI) layer as a critical component for regulating the Li metal electrodeposition processes. This study employs a first-principles kinetic Monte Carlo (kMC) model to simulate the SEI formation and Li + electrodeposition processes on a lithium metal anode, integrating both the electrochemical electrolyte reduction reactions and the diffusion events giving place to the SEI aggregation processes during battery charge and discharge processes. The model replicates the competitive interactions between organic and inorganic SEI components, emphasizing the influence of the cycling regime. Results indicate that grain boundaries within the SEI facilitate faster lithium-ion transport compared to crystalline regions, crucial for improving the performance and stability of LMBs. The findings underscore the importance of dynamic SEI modeling for further development of next-generation high-energy-density batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Accelerated kinetic Monte Carlo: A case study; vacancy and dumbbell interstitial diffusion traps in concentrated solid solution alloys

Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A DFT-based kinetic Monte Carlo simulation of multiphase oxide-metal thin film growth

Functional thin films of nanoscale metal pillars in oxide or nitride matrices known as vertically aligned nanocomposite (VAN) have gained much interest owing to their unique strain-coupled and highly anisotropic properties. So far, the deposition of these films has been explored mostly experimentally. In this work, a density functional theory (DFT)-based kinetic Monte Carlo simulation model using Bortz–Kalos–Lebowitz algorithm was developed to understand the growth of VAN films deposited by pulsed laser technique on mismatching substrates. The model has been parameterized and applied to understand the kinetics of growth thin films consisting of Au pillars in CeO2 matrix deposited on SrTiO3 substrates. The effects of pulsed laser deposition (PLD) conditions including the pulse frequency, deposition flux, and substrate temperature were explored. The simulations indicate that the Au pillar size and shape exhibit significant dependence on the PLD conditions. Namely, increasing the temperature increases the average pillar size and lowers the pillar density, and vice versa. In addition, the simulations revealed that increasing the deposition rate results in lowering the average pillar size and increasing the density. Particularly, the DFT results suggest that Au pillar size can be tuned during the initial growth of the first monolayer due to the significantly low activation barrier. Our analysis showed that the relationship between the average pillar size and pillar density is influenced by the kinetics. Furthermore, autocorrelation analysis showed that pillars self-organize in quasi-ordered patterns at certain windows of the deposition conditions, which is attributed to the complex nature of the chemical interactions in the system, the kinetics, and the deposition parameters.

Physics↗
NASA NTRS2015

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Ytterbium disilicate is of interest as a potential environmental barrier coating for aerospace applications, notably for use in next generation jet turbine engines. In such applications, the transport of oxygen and water vapor through these coatings to the ceramic substrate is undesirable if high temperature oxidation is to be avoided. In an effort to understand the diffusion process in these materials, we have performed kinetic Monte Carlo simulations of vacancy-mediated and interstitial oxygen diffusion in Ytterbium disilicate. Oxygen vacancy and interstitial site energies, vacancy and interstitial formation energies, and migration barrier energies were computed using Density Functional Theory. We have found that, in the case of vacancy-mediated diffusion, many potential diffusion paths involve large barrier energies, but some paths have barrier energies smaller than one electron volt. However, computed vacancy formation energies suggest that the intrinsic vacancy concentration is small. In the case of interstitial diffusion, migration barrier energies are typically around one electron volt, but the interstitial defect formation energies are positive, with the result that the disilicate is unlikely to exhibit experience significant oxygen permeability except at very high temperature.

Coatings↗
DOE OSTI2022

Unraveling small-scale defects in irradiated ThO 2 using kinetic Monte Carlo simulations

Point defects and their clusters generated through irradiation can have significant impact on the physical and mechanical properties of materials. However, direct experimental visualization of these small-scale defects using high-resolution scanning transmission electron microscopy remains a challenging task. In this study, using thorium dioxide (ThO 2 ) with the fluorite structure as a model system, we demonstrate the use of ab initio basin-hopping simulations in synergy with object kinetic Monte Carlo simulations as a powerful tool for identifying small defect complexes in irradiated materials. In addition to providing quantitative insights into defect evolution in ThO 2 under irradiation, our study reveals an unexpected role of bound anti-Schottky defect clusters in mediating defect transport. Remarkably, despite their poor thermal stability against dissociation at high temperatures, the transient formation of bound anti-Schottky defects under irradiation and their subsequent migration provide the dominant mechanism for the growth of large interstitial loops that have been experimentally observed in ThO 2 .

36 MATERIALS SCIENCE↗
NASA NTRS2017

Kinetic Monte Carlo Simulations of Oxygen Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Composite (CMC) materials are of interest for use in next-generation turbine engine components, offering a number of significant advantages, including reduced weight and high operating temperatures. However, in the hot environment in which such components operate, the presence of water vapor can lead to corrosion and recession, limiting the useful life of the components. Such degradation can be reduced through the use of Environmental Barrier Coatings (EBCs) that limit the amount of oxygen and water vapor reaching the component. Candidate EBC materials include Yttrium and Ytterbium silicates. In this work we present results of kinetic Monte Carlo (kMC) simulations of oxygen diffusion, via the vacancy mechanism, in Yttrium and Ytterbium disilicates, along with a brief discussion of interstitial diffusion.

Good, Brian S.↗
NASA NTRS2017

Kinetic Monte Carlo Simulations of Oxygen Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Composite (CMC) materials are of interest for use in next-generation turbine engine components, offering a number of significant advantages, including reduced weight and high operating temperatures. However, in the hot environment in which such components operate, the presence of water vapor can lead to corrosion and recession, limiting the useful life of the components. Such degradation can be reduced through the use of Environmental Barrier Coatings (EBCs) that limit the amount of oxygen and water vapor reaching the component. Candidate EBC materials include Yttrium and Ytterbium silicates. In this work we present results of kinetic Monte Carlo (kMC) simulations of oxygen diffusion, via the vacancy mechanism, in Yttrium and Ytterbium disilicates, along with a brief discussion of interstitial diffusion. An EBC system typically includes a bond coat located between the EBC and the component surface. Bond coat materials are generally chosen for properties other than low oxygen diffusivity, but low oxygen diffusivity is nevertheless a desirable characteristic, as the bond coat could provide some additional component protection, particularly in the case where cracks in the coating system provide a direct path from the environment to the bond coat interface. We have therefore performed similar kMC simulations of oxygen diffusion in this material.

Good, Brian S.↗
NASA NTRS2017

Kinetic Monte Carlo Simulations of Diffusion in Environmental Barrier Coating Materials

Ceramic Matrix Components (CMC) components for use in turbine engines offer a number of advantages compared with current practice. However, such components are subject to degradation through a variety of mechanisms. In particular, in the hot environment inside a turbine in operation a considerable amount of water vapor is present, and this can lead to corrosion and recession. Environmental Barrier Coating (EBC) systems that limit the amount of oxygen and water reaching the component are required to reduce this degradation and extend component life. A number of silicate-based materials are under consideration for use in such coating systems, including Yttterbium and Yttrium di- and monosilicates. In this work, we present results of kinetic Monte Carlo computer simulations of oxygen diffusion in Yttrium disilicate, and compare with previous work on Yttterbium disilicate. Coatings may also exhibit cracking, and the cracks can provide a direct path for oxygen to reach the component. There is typically a bond coat between the coating and component surface, but the bond coat material is generally chosen for properties other than low oxygen diffusivity. Nevertheless, the degree to which the bond coat can inhibit oxygen diffusion is of interest, as it may form the final defense against oxygen impingement on the component. We have therefore performed similar simulations of oxygen diffusion through HfSiO4, a proposed bond coat material.

Computer Simulation↗
NASA NTRS2015

Kinetic Monte Carlo Simulation of Oxygen Diffusion in Ytterbium Disilicate

Silicon-based ceramic components for next-generation jet turbine engines offer potential weight savings, as well as higher operating temperatures, both of which lead to increased efficiency and lower fuel costs. Silicon carbide (SiC), in particular, offers low density, good strength at high temperatures, and good oxidation resistance in dry air. However, reaction of SiC with high-temperature water vapor, as found in the hot section of jet turbine engines in operation, can cause rapid surface recession, which limits the lifetime of such components. Environmental Barrier Coatings (EBCs) are therefore needed if long component lifetime is to be achieved. Rare earth silicates such as Yb2Si2O7 and Yb2SiO5 have been proposed for such applications; in an effort to better understand diffusion in such materials, we have performed kinetic Monte Carlo (kMC) simulations of oxygen diffusion in Ytterbium disilicate, Yb2- Si2O7. The diffusive process is assumed to take place via the thermally activated hopping of oxygen atoms among oxygen vacancy sites or among interstitial sites. Migration barrier energies are computed using density functional theory (DFT).

Computer Simulation↗
NASA NTRS2011

Kinetic Monte Carlo Investigation of the Effects of Vacancy Pairing on Oxygen Diffusivity in Yttria-Stabilized Zirconia

Yttria-stabilized zirconia s high oxygen diffusivity and corresponding high ionic conductivity, and its structural stability over a broad range of temperatures, have made the material of interest for use in a number of applications, for example, as solid electrolytes in fuel cells. At low concentrations, the stabilizing yttria also serves to increase the oxygen diffusivity through the presence of corresponding oxygen vacancies, needed to maintain charge neutrality. At higher yttria concentration, however, diffusivity is impeded by the larger number of relatively high energy migration barriers associated with yttrium cations. In addition, there is evidence that oxygen vacancies preferentially occupy nearest-neighbor sites around either dopant or Zr cations, further affecting vacancy diffusion. We present the results of ab initio calculations that indicate that it is energetically favorable for oxygen vacancies to occupy nearest-neighbor sites adjacent to Y ions, and that the presence of vacancies near either species of cation lowers the migration barriers. Kinetic Monte Carlo results from simulations incorporating this effect are presented and compared with results from simulations in which the effect is not present.

Good, Brian S.↗
DOE OSTI2021

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

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

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗