Molecular Dynamics Study of Hydrogen in α-Zirconium

Abstract

Molecular dynamics approach is used to simulate hydrogen (H) diffusion in zirconium. Zirconium alloys are used in fuel channels of many nuclear reactors. Previously developed embedded atom method (EAM) and modified embedded atom method (MEAM) are tested and a good agreement with experimental data for lattice parameters, cohesive energy, and mechanical properties is obtained. Both EAM and MEAM are used to calculate hydrogen diffusion in zirconium. At higher temperatures and in the presence of hydrogen, MEAM calculation predicts an unstable zirconium structure and low diffusion coefficients. Mean square displacement (MSD) of hydrogen in bulk zirconium is calculated at a temperature range of 500–1200 K with diffusion coefficient at 500 K equals 1.92 * 10−7 cm2/sec and at 1200 K has a value 1.47 * 10−4 cm2/sec. Activation energy of hydrogen diffusion calculated using Arrhenius plot was found to be 11.3 kcal/mol which is in agreement with published experimental results. Hydrogen diffusion is the highest along basal planes of hexagonal close packed zirconium.