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: 3 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Generalized Gradient Approximation Made Simple
journal, October 1996

  • Perdew, John P.; Burke, Kieron; Ernzerhof, Matthias
  • Physical Review Letters, Vol. 77, Issue 18, p. 3865-3868
  • DOI: 10.1103/PhysRevLett.77.3865

Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set
journal, October 1996


A climbing image nudged elastic band method for finding saddle points and minimum energy paths
journal, December 2000

  • Henkelman, Graeme; Uberuaga, Blas P.; Jónsson, Hannes
  • The Journal of Chemical Physics, Vol. 113, Issue 22, p. 9901-9904
  • DOI: 10.1063/1.1329672