skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information
  1. Self-learning kinetic Monte Carlo simulations of diffusion in ferromagnetic α -Fe–Si alloys

    Diffusion in α-Fe-Si alloys is studied using AKSOME, an on-lattice self-learning KMC code, in the ferromagnetic state. Si diffusivity in the α-Fe matrix were obtained with and without the magnetic disorder in various temperature ranges. In addition we studied vacancy diffusivity in ferromagnetic α-Fe at various Si concentrations up to 12.5at.% in the temperature range of 350–550 K. The results were compared with available experimental and theoretical values in the literature. Local Si-atom dependent activation energies for vacancy hops were calculated using a broken-model and were stored in a database. The migration barrier and prefactors for Si-diffusivity were found tomore » be in reasonable agreement with available modeling results in the literature. Magnetic disorder has a larger effect on the prefactor than on the migration barrier. Prefactor was approximately an order of magnitude and the migration barrier a tenth of an electron-volt higher with magnetic disorder when compared to a fully ferromagnetic ordered state. In addition, the correlation between various have a larger effect on the Si-diffusivity extracted in various temperature range than the magnetic disorder. In the case of vacancy diffusivity, the migration barrier more or less remained constant while the prefactor decreased with increasing Si concentration in the disordered or A2-phase of Fe-Si alloy. Important vacancy-Si/Fe atom exchange processes and their activation barriers were also identified and discuss the effect of energetics on the formation of ordered phases in Fe-Si alloys.« less
  2. DEVELOPMENT OF INTERATOMIC POTENTIALS IN TUNGSTEN-RHENIUM SYSTEMS

    Reference data are generated using the ab initio method to fit interatomic potentials for the W-Re system. The reference data include single phases of W and Re, strained structures, slabs, systems containing several concentrations of vacancies, systems containing various types of interstitial defects, melt structures, structures in the σ and χ phases, and structures containing several concentrations of solid solutions of Re in bcc W and W in hcp Re. Future work will start the fitting iterations.
  3. INTERACTION OF INTERSTITIAL CLUSTERS WITH RHENIUM, OSMIUM, AND TANTALUM IN TUNGSTEN

    In the previous semi annual report, we explored the stability of interstitial clusters in W up to size seven. In this report, we study the binding of those clusters to Re, Os, and Ta atoms. For each cluster size, the three most stable configurations are considered to average the binding property. The average binding energy to a Re decreases from 0.79 eV for a size-1 cluster (a [111] dumbbell) to 0.65 eV for a size-7 cluster. For Os, the binding decreases from 1.61 eV for a [111] dumbbell to 1.34 eV for a size-7 cluster. Tantalum is repulsive to interstitialmore » clusters with binding energy ranges from -0.61 eV for a [111] dumbbell to -0.5 eV for a size-7 cluster.« less
  4. Ab initio study of interstitial cluster interaction with Re, Os, and Ta in W

    The stability of tungsten self-interstitial atom (SIA) clusters is studied using first-principles methods. Clusters from one to seven SIAs are systematically explored from 1264 unique configurations. Finite-size effect of the simulation cell is corrected based on the scaling of formation energy versus inverse volume cell. Furthermore, the accuracy of the calculations is improved by treating the 5p semicore states as valence states. Configurations of the three most stable clusters in each cluster size n are presented, which consist of parallel [111] dumbbells. The evolution of these clusters leading to small dislocation loops is discussed. The binding energy of size-n clustersmore » is analyzed relative to an n → (n-1) + 1 dissociation and is shown to increase with size. Extrapolation for n > 7 is presented using a dislocation loop model. In addition, the interaction of these clusters with a substitutional Re, Os, or Ta solute is explored by replacing one of the dumbbells with the solute. Re and Os strongly attract these clusters, but Ta strongly repels. The strongest interaction is found when the solute is located on the periphery of the cluster rather than in the middle. The magnitude of this interaction decreases with cluster size. Empirical fits to describe the trend of the solute binding energy are presented.« less
  5. STABILITY OF SMALL SELF-INTERSTITIAL CLUSTERS IN TUNGSTEN

    Density functional theory was employed to explore the stability of interstitial clusters in W up to size seven. For each cluster size, the most stable configuration consists of parallel dumbbells. For clusters larger than size three, parallel dumbbells prefer to form in a multilayer fashion, instead of a planar structure. For size-7 clusters, the most stable configuration is a complete octahedron. The binding energy of a [111] dumbbell to the most stable cluster increases with cluster size, namely 2.49, 3.68, 4.76, 4.82, 5.47, and 6.85 eV for clusters of size 1, 2, 3, 4, 5, and 6, respectively. For amore » size-2 cluster, collinear dumbbells are still repulsive at the maximum allowable distance of 13.8 Å (the fifth neighbor along [111]). On the other hand, parallel dumbbells are strongly bound together. Two parallel dumbbells in which the axis-to-axis distance is within a cylindrical radius of 5.2 Å still exhibit a considerable binding of 0.28 eV. The most stable cluster in each size will be used to explore interactions with transmutation products.« less
  6. Ab initio study of interstitial cluster interaction with Re, Os, and Ta in W

  7. Self-learning kinetic Monte Carlo simulations of Al diffusion in Mg

    Atomistic on-lattice self-learning kinetic Monte Carlo (SLKMC) method was used to examine the vacancy-mediated diffusion of an Al atom in pure hcp Mg. Local atomic environment dependent activation barriers for vacancy-atom exchange processes were calculated on-the-fly using climbing image nudged-elastic band method (CI-NEB) and using a Mg-Al binary modified embedded-atom method (MEAM) interatomic potential. Diffusivities of vacancy and Al atom in pure Mg were obtained from SLKMC simulations and are compared with values available in the literature that are obtained from experiments and first-principle calculations. Al Diffusivities obtained from SLKMC simulations are lower, due to larger activation barriers and lowermore » diffusivity prefactors, than those available in the literature but have same order of magnitude. We present all vacancy-Mg and vacancy-Al atom exchange processes and their activation barriers that were identified in SLKMC simulations. We will describe a simple mapping scheme to map a hcp lattice on to a simple cubic lattice that would enable hcp lattices to be simulated in an on-lattice KMC framework. We also present the pattern recognition scheme used in SLKMC simulations.« less
  8. OBJECT KINETIC MONTE CARLO SIMULATIONS OF RADIATION DAMAGE ACCUMULATION IN TUNGSTEN

    The objective of this work is to understand the accumulation of radiation damage created by primary knock-on atoms (PKAs) of various energies, at 300 K and for a dose rate of 10-4 dpa/s in bulk tungsten using the object kinetic Monte Carlo (OKMC) method.
  9. OBJECT KINETIC MONTE CARLO SIMULATIONS OF RADIATION DAMAGE IN TUNGSTEN SUBJECTED TO NEUTRON FLUX WITH PKA SPECTRUM CORRESPONDING TO THE HFIR

    The objective of this work is to study the damage accumulation in pure tungsten (W) subjected to neutron bombardment with a primary knock-on atom (PKA) spectrum corresponding to the High Flux Isotope Reactor (HFIR), using the object kinetic Monte Carlo (OKMC) method.
  10. OBJECT KINETIC MONTE CARLO SIMULATIONS OF RADIATION DAMAGE IN BULK TUNGSTEN

    We used our recently developed lattice based OKMC code; KSOME [1] to carryout simulations of radiation damage in bulk W. We study the effect of dimensionality of self interstitial atom (SIA) diffusion i.e. 1D versus 3D on the defect accumulation during irradiation with a primary knock-on atom (PKA) energy of 100 keV at 300 K for the dose rates of 10-5 and 10-6 dpa/s. As expected 3D SIA diffusion significantly reduces damage accumulation due to increased probability of recombination events. In addition, dose rate, over the limited range examined here, appears to have no effect in both cases of SIAmore » diffusion.« less
...

Search for:
All Records
Creator / Author
"Nandipati, Giridhar"

Refine by:
Resource Type
Availability
Publication Date
Creator / Author
Research Organization