Molecular dynamics simulation of interaction of H with vacancy in W

Molecular dynamics simulations were performed to investigate the interaction between H and vacancy in W using an analytical bond-order potential to describe the interactions between W-W, W-H and H-H. The most stable configuration for H in W is the tetrahedron interstitial site. We calculated the binding energies of an H and a vacancy to an H-vacancy cluster (HnVm) in W, respectively, where n and m ranged from 0 to 10. The binding energy was almost unchanged. The binding energy of a vacancy to H-vacancy cluster is about 0.4 eV, which is higher than the binding energy of an H to H-vacancy cluster. Vacancy is much easier to bond with H-vacancy cluster than H. And H is easier to stay in the tetrahedron interstitial site or octahedron interstitial site in bcc W. (C) 2009 Elsevier B.V. All rights reserved.