Search NASA⌕ Search

SEARCH · Search NASA

Results for “kMC”

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 OSTI2023

Nucleation-mediated reshaping of facetted metallic nanocrystals: Breakdown of the classical free energy picture

Shape stability is key to avoiding degradation of performance for metallic nanocrystals synthesized with facetted non-equilibrium shapes to optimize properties for catalysis, plasmonics, and so on. Reshaping of facetted nanocrystals is controlled by the surface diffusion-mediated nucleation and growth of new outer layers of atoms. Here, Kinetic Monte Carlo (KMC) simulation of a realistic stochastic atomistic-level model is applied to precisely track the reshaping of Pd octahedra and nanocubes. Unexpectedly, separate constrained equilibrium Monte Carlo analysis of the free energy profile during reshaping reveals a fundamental failure of the classical nucleation theory (CNT) prediction for the reshaping barrier and rate. Why? Nucleation barriers can be relatively low for these processes, so the system is not locally equilibrated before crossing the barrier, as assumed in CNT. This claim is supported by an analysis of a first-passage problem for reshaping within a master equation framework for the model that reasonably captures the behavior in KMC simulations.

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

Probing exciton diffusion dynamics in photosynthetic supercomplexes via exciton–exciton annihilation

Photosynthesis converts solar energy into chemical energy through coordinated energy transfer between light-harvesting complexes and reaction centers (RCs). Understanding exciton motion, particularly the exciton diffusion length, is essential for optimizing energy efficiency in photosystems. In this work, we combine intensity-cycling transient absorption spectroscopy with kinetic Monte Carlo (kMC) simulation to investigate exciton motion in the C2S2 photosystem II supercomplex of spinach. Using exciton–exciton annihilation, revealed in the fifth-order response, we experimentally estimate an exciton diffusion length of 10.9 nm based on a 3D normal diffusion model, suggesting the ability of excitons to traverse the supercomplex. However, kMC simulations reveal that exciton motion is sub-diffusive because of spatial constraints and the strong RC traps. An anomalous diffusion model analysis of the experimental data yields a diffusion length of 9.7 nm, while the simulated diffusion length is 7.4 nm. The variable exciton residence time across subunits, partly influenced by their connectivity to the trap, indicates inhomogeneous annihilation probability and suggests how plants balance efficient light harvesting with photoprotection. We also explore the influence of specific assumptions in the annihilation simulation, which are challenging to access in more complex environments, such as the thylakoid membrane. Our study provides a framework for studying exciton dynamics using exciton–exciton annihilation, which can be extended to understand the light-harvesting efficiencies of larger, more complex photosynthetic assemblies.

Zhang, Kunyan (ORCID:000000026830409X)↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗
DOE OSTI2024

Structural and mechanical properties of monolayer amorphous carbon and boron nitride

Amorphous materials exhibit various characteristics that are not featured by crystals and can sometimes be tuned by their degree of disorder (DOD). Here, we report results on the mechanical properties of monolayer amorphous carbon (MAC) and monolayer amorphous boron nitride (maBN) with different DOD. The pertinent structures are obtained by kinetic-Monte-Carlo (kMC) simulations using machine-learning potentials (MLP) with density-functional-theory (DFT)-level accuracy. An intuitive order parameter, namely the areal fraction F x occupied by crystallites within the continuous random network, is proposed to describe the DOD. We find that F x captures the essence of the DOD: Samples with the same F x but different sizes and arrangements of crystallites, obtained using two distinct kMC procedures, have virtually identical radial distributions functions as well as bond-length and bond-angle distributions. Furthermore, by simulating the fracture process with molecular dynamics, we found that the mechanical responses of MAC and maBN before fracture are mainly determined by F x and are insensitive to the sizes and specific arrangements and to some extent the numbers and area distributions of the crystallites. The behavior of cracks in the two materials is analyzed and found to mainly propagate in meandering paths in the CRN region and to be influenced by crystallites in distinct ways that toughen the material. Furthermore, the present results reveal the relation between structure and mechanical properties in amorphous monolayers and may provide a universal toughening strategy for 2D materials.

2-dimensional systems↗
DOE OSTI2023

Generic two-phase coexistence in a type-2 Schloegl model for autocatalysis on a square lattice: Analysis via heterogeneous master equations

Schloegl's second model (also known as the quadratic contact process) on a square lattice involves spontaneous annihilation of particles at lattice sites at rate p, and their autocatalytic creation at unoccupied sites with n ≥ 2 occupied neighbors at rate k n . Here, kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition with generic two-phase coexistence: the p value for equistability of coexisting populated and vacuum states, p eq (S), depends on the orientation or slope, S, of a planar interface separating those phases. The vacuum state displaces the populated state for p > eq (S), and the opposite applies for p < p eq (S) for 0 < S < ∞. The special “combinatorial” rate choice k n = n(n–1)/12 facilitates an appealing simplification of the exact master equations for the evolution of spatially heterogeneous states in the model, which aids analytic investigation of these equations via hierarchical truncation approximations. Truncation produces coupled sets of lattice differential equations which can describe orientation-dependent interface propagation and equistability. The pair approximation predicts that p eq (max) = p eq (S = 1) = 0.09645 and p eq (min) = p eq (S → ∞) = 0.08827, values deviating less than 15% from KMC predictions. In the pair approximation, a perfect vertical interface is stationary for all p < p eq (S = ∞) = 0.08907, a value exceeding p eq (S → ∞). One can regard an interface for large S → ∞ as a vertical interface decorated with isolated kinks. For p < p eq (S = ∞), the kink can move in either direction along this otherwise stationary interface depending upon p, but for p = p eq (min) the kink is also stationary.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Trends in spatial correlations, two-phase coexistence, and criticality in a class of type-2 Schloegl models for autocatalysis

A class of type-2 Schloegl models is considered for particles on a square lattice with variable-range cooperativity. These models involve: (i) spontaneous particle annihilation at rate p; (ii) autocatalytic particle creation at unoccupied sites (i, j) with n ⩾ 2 particles within a specified neighborhood, Ω 𝑁 (i, j), of sites at rate $k_n$ = $\frac{^{(^n_2)}}{_{(^N_2)}}$ = $\frac{n{(n-1)}}{_{N(N-1)}}$; and (iii) possible spontaneous particle creation at unoccupied sites at “small” rate ɛ ⩾ 0. In some cases, Ω 𝑁 just includes all symmetry-equivalent sites at a single specific distance 𝑑 (in units of lattice constants) from the unoccupied site, e.g., 𝑑 = 1 (nearest-neighbor sites) where 𝑁 = 4, or 𝑑 = √5 (or √13 or…) where 𝑁 = 8. In other cases, Ω 𝑁 includes sites multiple distances from the unoccupied site, e.g., 𝑑 = {1,√2}, where 𝑁 = 8. Kinetic Monte Carlo (KMC) simulation reveals that these models exhibit a nonequilibrium discontinuous phase transition between high- and low-density states below a critical point, ɛ < ɛ c , with generic two-phase coexistence (2PC) at least for smaller 𝑁. With some exceptions, there is an approach toward mean-field behavior with increasing 𝑁 (so the regime of generic 2PC shrinks, and ɛ c approaches the mean-field value of 1/27). Additional insight into trends is provided by analysis of the exact master equations for the models via hierarchical truncation. These truncations utilize suitably tailored pair approximations which reflect the dominant nonequilibrium spatial correlations. These correlations in turn are shown to reflect the details of the autocatalytic particle creation process. For spatially heterogeneous states, the truncations produce coupled sets of lattice differential equations (LDE) which can describe orientation-dependent propagation of an interface between high- and low-density steady states for ɛ < ɛ c . Pair approximation values of 𝑝 = 𝑝 eq where the interface is stationary, and its orientation-dependence, are in semiquantitative agreement with KMC results. In conclusion, this comparison accounts for propagation failure in the LDE which complicates interpretation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The Critical Stack Pressure to Alter Void Generation at Li/Solid-Electrolyte Interfaces during Stripping

The lithium stripping process generates vacancies, which may accumulate as voids and lead to uneven current distribution and dendrite growth in the following plating cycles. A stack pressure is typically required during stripping, but how to optimize the stack pressure is not clear. In this work, extremely lithiophilic Li/Li 2 O and lithiophobic Li/LiF interfaces were used to reveal the combining effect of interface interaction and stack pressure induced lithium creep on the stripping critical current density (CCD). A multiscale simulation scheme with Density Functional Theory (DFT), kinetic Monte Carlo (KMC) simulations, and an analytical model was developed. The analytical model predicted lithiophobic interfaces require a higher stack pressure than lithiophilic interfaces to reach the same CCD. The KMC simulations also showed higher stack pressure is needed at lithiophobic interfaces to accelerate Li vacancy diffusion into the bulk and maintain a flat surface. This stack pressure needs to be high enough to alter the Li forward-and-backward hopping barriers at the interface. This multiscale simulation scheme illustrates the importance to include the chemical-mechanical effects during Li stripping morphology evolution. It can be used to design ideal interlayer coating materials to maintain a flat Li surface during cycling.

25 ENERGY STORAGE↗
DOE OSTI2023

KinCat v.1.0

SAND2024-02099O The software is designed to allow researchers to perform kinetic Monte Carlo (KMC) simulations of catalytic reactions on a 2D lattice. The code is written efficiently to run on a variety of shared memory computing architectures (e.g. GPU, multi-core) and to natively express the full complexity of lateral interactions on reaction rates. The software allows researchers to perform KMC simulations of catalytic reactions on a 2D lattice. It uses parallel shared-memory computing architectures to reduce run-times and allows for simultaneous simulation of multiple independent runs. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Najm, Habib↗
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

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

An Integrated Simulation of Multiple-Pass U-10Mo Alloy Hot Rolling and Static Recrystallization

To achieve a desired microstructure and minimize the thickness variation in rolled foils, researchers must understand the effects of foil fabrication process variables on microstructure evolution. We developed an integrated simulation of deformation and recrystallization that employs the finite element method (FEM) and the kinetic Monte Carlo (KMC) Potts model, respectively, to investigate microstructure evolution during multiple-pass hot rolling and heat treatment in polycrystalline U-10Mo fuel. Scanning electron microscopy and electron backscatter diffraction images of microstructures were directly used as input in FEM calculation of deformation, and the calculated strains were used to determine the driving force of nucleation and growth of recrystallized grains in the Potts model. Grain structures predicted by the Potts model were used to update the grain structure and material properties for FEM. Simulation alternated between FEM and the Potts model to simulate grain structure evolution during multiple rolling and heat treatments. The initial model parameters were determined by benchmarking the recrystallization kinetics against experimental data. Then, the model was applied to predict the grain structure evolution. Results showed that our model can capture the coupling between deformation and recrystallization and can quantitatively reproduce the observed U-10Mo recrystallization and grain growth kinetics. The simulation results demonstrated that the developed model can predict U-10Mo grain structures as a function of initial microstructure and foil fabrication parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Anomaly Detection in Materials Digital Twins with Multiscale ICME for Additive Manufacturing

Detecting anomaly in fatigue and fracture experimental materials science is an interesting yet challenging topic. The reasons are threefold. First, the anomalous microstructure feature that gives rise to structural failure is small, sometimes in the order of 10 -7 of the interrogated volume. This, in turn, results in a highly imbalanced classification problem in machine learning (ML). Second, the consequence is high, in the sense that the test specimen is destructed in such case. Third, the convolution between microstructure stochasticity and the small probability of void nucleation, growth, and coalescence makes failure and fracture a hard-to-predict and challenging problem in materials science due to its irreproducibility, even experimentally. In this paper, we developed a materials digital twin and applied anomaly detection methods to detect voids and anomaly in additive manufacturing (AM). The materials digital twin is driven by two integrated computational materials engineering (ICME) models, which are kinetic Monte Carlo (kMC) and crystal plasticity finite element method (CPFEM). In conclusion, we demonstrated that by using anomaly detection, it is possible to detect voids and other defects in materials digital twin, which paves way for future research in integrating materials digital twin with its physical counterpart.

ICME↗
DOE OSTI2020

An active learning high-throughput microstructure calibration framework for solving inverse structure–process problems in materials informatics

Determining a process–structure–property relationship is the holy grail of materials science, where both computational prediction in the forward direction and materials design in the inverse direction are essential. Problems in materials design are often considered in the context of process–property linkage by bypassing the materials structure, or in the context of structure–property linkage as in microstructure-sensitive design problems. However, there is a lack of research effort in studying materials design problems in the context of process–structure linkage, which has a great implication in reverse engineering. In this paper, given a target microstructure, we propose an active learning high-throughput microstructure calibration framework to derive a set of processing parameters, which can produce an optimal microstructure that is statistically equivalent to the target microstructure. The proposed framework is formulated as a noisy multi-objective optimization problem, where each objective function measures a deterministic or statistical difference of the same microstructure descriptor between a candidate microstructure and a target microstructure. Furthermore, to significantly reduce the physical waiting wall-time, we enable the high-throughput feature of the microstructure calibration framework by adopting an asynchronously parallel Bayesian optimization to exploit high-performance computing resources. Case studies in additive manufacturing and grain growth are used to demonstrate the applicability of the proposed framework, where kinetic Monte Carlo (kMC) simulation is used as a forward predictive model, such that for a given target microstructure, the target processing parameters that produced this microstructure are successfully recovered.

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 OSTI2023

Understanding structure-processing relationships in metal additive manufacturing via featurization of microstructural images

Understanding and predicting accurate property-structure-processing relationships for additively manufactured components is important for both forward and inverse design of robust, reliable parts and assemblies. While direct mapping of process parameters to properties is sometimes plausible, it is often rendered difficult due to poor microstructural control. Exploring the direct relationship between processing conditions and microstructural features can thus provide significant physical insights and aid the overall design process. Here, in this study, we develop an automated high-throughput framework to simulate an uncertainty-aware additive manufacturing (AM) process, characterize microstructural images, and extract meaningful features/descriptors. A kinetic Monte Carlo (KMC) based model of the AM process is used to simulate microstructural evolution for a diverse set of experimentally relevant processing conditions. We perform a parametric study to explore the relationship between microstructural features and processing conditions. Our results indicate that a many-to-one mapping can exist between processing conditions and typical descriptors; therefore, multiple descriptors are thus necessary to unambiguously represent microstructural images. Our work provides crucial quantitative and qualitative in-formation that would aid in the selection of features for microstructural images. Featurized microstructures could then be utilized to build data-driven models for predictive control of microstructures and thereby properties of additively manufactured components.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗
DOE OSTI2022

Engineering defect energy landscape of CoCrFeNi high-entropy alloys by the introduction of additional dopants

The concept of high-entropy alloys (HEAs) focusing on tuning the overall chemical complexity represents a novel alloy design strategy. In contrast, alloying of a metallic matrix with minor doping elements with limited and localized tunability has been a common practice to improve material performance. Combining the idea of globally engineering defect energy landscape in HEAs and the localized doping strategy in dilute alloys, in this work, we explore doping effects of minor elements in a HEA matrix to further enhance the overall and localized chemical tunability, aiming to improve its irradiation resistance. Specifically, we study the influence of minor Al, Cu, Ti, and Pd substitutional doping elements on defect energetics in a CoCrFeNi model HEA based on density-functional theory (DFT) calculations. The DFT results indicate that the formation and migration energies of vacancies can be strongly influenced when a dopant is introduced at the first nearest neighbor shells around a vacancy. On the other hand, interstitial energetics are only slightly affected. Among the four elements, Ti and Pd generally decrease vacancy formation energies and increase vacancy migration energies more significantly than Al and Cu. The doping effects become more pronounced when the concentration of the substitutional dopants increases. Based on the energy distributions obtained from DFT, we build a kinetic Monte Carlo (kMC) model to assess the impact of dopants on vacancy-mediated diffusivity in the doped HEAs. Our results suggest that Ti and Pd can lower the tracer diffusivity in the considered HEAs and act as trapping sites, whereas Cu may enhance the atomic transport. This work indicates that substitutional doping in HEAs is an effective strategy in metallurgy to further tune the defect and transport properties of complex alloys.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗