Molecular dynamics calculation of thermal expansion coefficient of a series of rare-earth zirconates

To develop more reliable and stable thermal barrier coatings, low thermal conductivity materials with higher thermal expansion coefficients (TECs), like rare-earth zirconates (A 2Zr 2O 7, A = La, Nd, Sm, Eu, Gd, Er, Yb and Lu), have been the focus of a great deal of attention in recent years. But the mechanisms underlying the values of TECs are still poorly understood. In this paper, the TECs of a series of rare-earth zirconates are calculated using molecular dynamics. Two typical crystal structures, pyrochlore and fluorite, are considered. By investigating the potential functions and the equilibrium-location deviations between the atoms, it is found that the Zr–O bond is the most important factor that determines the overall TECs, and the A–O bond plays a secondary role. However, O–O has little effect on the TECs. In addition, the fluorite structure has a relatively higher TEC due to a weaker Zr–O bond. Calculated values are consistent with the experimental observations. Using the method presented in this study, two doped zirconates (Gd 0.4Sm 0.5Yb 0.1) 2Zr 2O 7 and Sm 2(Ce 0.3Zr 0.7) 2O 7 are designed and compared. As expected, both of them have higher TECs than undoped Sm 2Zr 2O 7, and doping at the A-site is more efficient than doping at the Zr-site.