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 longrange 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 onlattice diffusion in hcp metals, or systems with similar trapping basins.
 Authors:

 Idaho National Lab. (INL), Idaho Falls, ID (United States)
 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/JOU1640134
Journal ID: ISSN 20452322; srep41033
 Grant/Contract Number:
 AC0705ID14517
 Resource Type:
 Accepted Manuscript
 Journal Name:
 Scientific Reports
 Additional Journal Information:
 Journal Volume: 7; Journal ID: ISSN 20452322
 Publisher:
 Nature Publishing Group
 Country of Publication:
 United States
 Language:
 English
 Subject:
 36 MATERIALS SCIENCE; accelerated kinetic Monte Carlo; hydrogen diffusion; atomistic models; condensedmatter 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 longrange 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 onlattice diffusion in hcp metals, or systems with similar trapping basins.},
doi = {10.1038/srep41033},
journal = {Scientific Reports},
number = ,
volume = 7,
place = {United States},
year = {2017},
month = {1}
}
Web of Science
Works referenced in this record:
Metal hydride materials for solid hydrogen storage: A review☆
journal, June 2007
 Sakintuna, B.; Lamaridarkrim, F.; Hirscher, M.
 International Journal of Hydrogen Energy, Vol. 32, Issue 9
Generalized Gradient Approximation Made Simple
journal, October 1996
 Perdew, John P.; Burke, Kieron; Ernzerhof, Matthias
 Physical Review Letters, Vol. 77, Issue 18, p. 38653868
Lattice Parameters, Thermal Expansions, and Grüneisen Coefficients of Zirconium, 4.2 to 1130°K
journal, April 1966
 Goldak, J.; Lloyd, L. T.; Barrett, C. S.
 Physical Review, Vol. 144, Issue 2
Determination of fracture strength of δzirconium hydrides embedded in zirconium matrix at high temperatures
journal, April 2013
 Kubo, T.; Kobayashi, Y.; Uchikoshi, H.
 Journal of Nuclear Materials, Vol. 435, Issue 13
Diffusion coefficient of hydrogen in alpha zirconium, Zircaloy2 and Zircaloy4
journal, June 1972
 Kearns, J. J.
 Journal of Nuclear Materials, Vol. 43, Issue 3
A review on hydride precipitation in zirconium alloys
journal, November 2015
 Bair, Jacob; Asle Zaeem, Mohsen; Tonks, Michael
 Journal of Nuclear Materials, Vol. 466
The HZr (hydrogenzirconium) system
journal, August 1990
 Zuzek, E.; Abriata, J. P.; SanMartin, A.
 Bulletin of Alloy Phase Diagrams, Vol. 11, Issue 4
Efficient iterative schemes for ab initio totalenergy calculations using a planewave basis set
journal, October 1996
 Kresse, G.; Furthmüller, J.
 Physical Review B, Vol. 54, Issue 16, p. 1116911186
Using density functional theory to study hydrogen diffusion in metals: A brief overview
journal, October 2007
 Sholl, David S.
 Journal of Alloys and Compounds, Vol. 446447
Diffusion of ${\mathrm{H}}_{2}$ , ${\mathrm{D}}_{2}$ , and ${\mathrm{T}}_{2}$ in SingleCrystal Ni and Cu
journal, July 1971
 Katz, L.; Guinan, M.; Borg, R. J.
 Physical Review B, Vol. 4, Issue 2
Adaptive kinetic Monte Carlo for firstprinciples accelerated dynamics
journal, September 2008
 Xu, Lijun; Henkelman, Graeme
 The Journal of Chemical Physics, Vol. 129, Issue 11
Hydrogen accommodation in Zr second phase particles: Implications for H pickup and hydriding of Zircaloy2 and Zircaloy4
journal, April 2013
 Burr, P. A.; Murphy, S. T.; Lumley, S. C.
 Corrosion Science, Vol. 69
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. 99019904
Homogeneous hydride formation path in αZr: Molecular dynamics simulations with the chargeoptimized manybody potential
journal, June 2016
 Zhang, Yongfeng; Bai, XianMing; Yu, Jianguo
 Acta Materialia, Vol. 111
An energy basin finding algorithm for kinetic Monte Carlo acceleration
journal, April 2010
 Puchala, Brian; Falk, Michael L.; Garikipati, Krishna
 The Journal of Chemical Physics, Vol. 132, Issue 13
Hydride Formation in Zirconium Alloys
journal, October 2012
 Motta, Arthur T.; Chen, LongQing
 JOM, Vol. 64, Issue 12
The diffusion and trapping of hydrogen in steel
journal, January 1970
 Oriani, R. A.
 Acta Metallurgica, Vol. 18, Issue 1
Atomicscale Abinitio study of the ZrH system: I. Bulk properties
journal, August 2002
 Domain, C.; Besson, R.; Legris, A.
 Acta Materialia, Vol. 50, Issue 13
Modeling proton mobility in acidic zeolite clusters: II. Room temperature tunneling effects from semiclassical rate theory
journal, April 2000
 Fermann, Justin T.; Auerbach, Scott
 The Journal of Chemical Physics, Vol. 112, Issue 15
A firstprinciples study of surface and subsurface H on and in Ni(111): diffusional properties and coveragedependent behavior
journal, August 2003
 Greeley, Jeff; Mavrikakis, Manos
 Surface Science, Vol. 540, Issue 23