skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Anisotropic hydrogen diffusion in α-Zr and Zircaloy predicted by accelerated kinetic Monte Carlo simulations

Abstract

Here, this report presents an accelerated kinetic Monte Carlo (KMC) method to compute the diffusivity of hydrogen in hcp metals and alloys, considering both thermally activated hopping and quantum tunneling. The acceleration is achieved by replacing regular KMC jumps in trapping energy basins formed by neighboring tetrahedral interstitial sites, with analytical solutions for basin exiting time and probability. Parameterized by density functional theory (DFT) calculations, the accelerated KMC method is shown to be capable of efficiently calculating hydrogen diffusivity in α-Zr and Zircaloy, without altering the kinetics of long-range diffusion. Above room temperature, hydrogen diffusion in α-Zr and Zircaloy is dominated by thermal hopping, with negligible contribution from quantum tunneling. The diffusivity predicted by this DFT + KMC approach agrees well with that from previous independent experiments and theories, without using any data fitting. The diffusivity along < c > is found to be slightly higher than that along < a >, with the anisotropy saturated at about 1.20 at high temperatures, resolving contradictory results in previous experiments. Demonstrated using hydrogen diffusion in α-Zr, the same method can be extended for on-lattice diffusion in hcp metals, or systems with similar trapping basins.

Authors:
 [1];  [1];  [2]
  1. Idaho National Lab. (INL), Idaho Falls, ID (United States)
  2. Virginia Polytechnic Inst. and State Univ. (Virginia Tech), Blacksburg, VA (United States)
Publication Date:
Research Org.:
Idaho National Lab. (INL), Idaho Falls, ID (United States)
Sponsoring Org.:
USDOE Office of Nuclear Energy (NE)
OSTI Identifier:
1357778
Report Number(s):
INL/JOU-16-40134
Journal ID: ISSN 2045-2322; srep41033
Grant/Contract Number:
AC07-05ID14517
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Scientific Reports
Additional Journal Information:
Journal Volume: 7; Journal ID: ISSN 2045-2322
Publisher:
Nature Publishing Group
Country of Publication:
United States
Language:
English
Subject:
36 MATERIALS SCIENCE; accelerated kinetic Monte Carlo; hydrogen diffusion; atomistic models; condensed-matter physics; metals and alloys

Citation Formats

Zhang, Yongfeng, Jiang, Chao, and Bai, Xianming. Anisotropic hydrogen diffusion in α-Zr and Zircaloy predicted by accelerated kinetic Monte Carlo simulations. United States: N. p., 2017. Web. doi:10.1038/srep41033.
Zhang, Yongfeng, Jiang, Chao, & Bai, Xianming. Anisotropic hydrogen diffusion in α-Zr and Zircaloy predicted by accelerated kinetic Monte Carlo simulations. United States. doi:10.1038/srep41033.
Zhang, Yongfeng, Jiang, Chao, and Bai, Xianming. Fri . "Anisotropic hydrogen diffusion in α-Zr and Zircaloy predicted by accelerated kinetic Monte Carlo simulations". United States. doi:10.1038/srep41033. https://www.osti.gov/servlets/purl/1357778.
@article{osti_1357778,
title = {Anisotropic hydrogen diffusion in α-Zr and Zircaloy predicted by accelerated kinetic Monte Carlo simulations},
author = {Zhang, Yongfeng and Jiang, Chao and Bai, Xianming},
abstractNote = {Here, this report presents an accelerated kinetic Monte Carlo (KMC) method to compute the diffusivity of hydrogen in hcp metals and alloys, considering both thermally activated hopping and quantum tunneling. The acceleration is achieved by replacing regular KMC jumps in trapping energy basins formed by neighboring tetrahedral interstitial sites, with analytical solutions for basin exiting time and probability. Parameterized by density functional theory (DFT) calculations, the accelerated KMC method is shown to be capable of efficiently calculating hydrogen diffusivity in α-Zr and Zircaloy, without altering the kinetics of long-range diffusion. Above room temperature, hydrogen diffusion in α-Zr and Zircaloy is dominated by thermal hopping, with negligible contribution from quantum tunneling. The diffusivity predicted by this DFT + KMC approach agrees well with that from previous independent experiments and theories, without using any data fitting. The diffusivity along < c > is found to be slightly higher than that along < a >, with the anisotropy saturated at about 1.20 at high temperatures, resolving contradictory results in previous experiments. Demonstrated using hydrogen diffusion in α-Zr, the same method can be extended for on-lattice diffusion in hcp metals, or systems with similar trapping basins.},
doi = {10.1038/srep41033},
journal = {Scientific Reports},
number = ,
volume = 7,
place = {United States},
year = {Fri Jan 20 00:00:00 EST 2017},
month = {Fri Jan 20 00:00:00 EST 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 1work
Citation information provided by
Web of Science

Save / Share:
  • 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
  • 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
  • The diffusion of dilute hydrogen in fcc Ni–Al and Ni–Fe binary alloys was examined using kinetic Monte Carlo method with input kinetic parameters obtained from first-principles density functional theory. The simulation involves the implementation of computationally efficient energy barrier model that describes the configuration dependence of the hydrogen hopping. The predicted hydrogen diffusion coefficients in Ni and Ni 89.4Fe 10.6 are compared well with the available experimental data. In Ni–Al, the model predicts lower hydrogen diffusivity compared to that in Ni. Overall, diffusion prefactors and the effective activation energies of H in Ni–Fe and Ni–Al are concentration dependent of themore » alloying element. Furthermore, the changes in their values are the results of the short-range order (nearest-neighbor) effect on the interstitial diffusion of hydrogen in fcc Ni-based alloys.« less