Equations of state for crystalline zirconium iodide: The role of dispersion

We present the first-principle equations of state of several zirconium iodides, ZrI2, ZrI3, and ZrI4, computed using density functional theory methods that apply various methods for introducing the dispersion correction. Iodides formed due to reaction of molecular or atomic iodine with zirconium and zircaloys are of particular interest due to their application to the cladding material used in the fabrication of nuclear fuel rods. Stress corrosion cracking (SCC), associated with fission product chemistry with the clad material, is a major concern in the life cycle of nuclear fuels, as many of the observed rod failures have occurred due to pellet–cladding chemical interactions (PCCI) [A. Atrens, G. Dannhäuser, G. Bäro, Stress-corrosion-cracking of zircaloy-4 cladding tubes, Journal of Nuclear Materials 126 (1984) 91–102; P. Rudling, R. Adamson, B. Cox, F. Garzarolli, A. Strasser, High burn-up fuel issues, Nuclear Engineering and Technology 40 (2008) 1–8]. A proper understanding of the physical properties of the corrosion products is, therefore, required for the development of a comprehensive SCC model. In this particular work, we emphasize that, while existing modeling techniques include methods to compute crystal structures and associated properties, it is important to capture intermolecular forces not traditionally included, such as van der Waals (dispersion) correction. Furthermore, crystal structures with stoichiometries favoring a high I:Zr ratio are found to be particularly sensitive, such that traditional density functional theory approaches that do not incorporate dispersion incorrectly predict significantly larger volumes of the lattice. This latter point is related to the diffuse nature of the iodide electron cloud.