Molecular dynamics simulation of metallic impurity diffusion in liquid lead-bismuth eutectic (LBE)

Corrosion of stainless steels by lead-bismuth eutectic (LBE) is an important problem which depends, amongst other things, on the diffusion of the steel components inside this liquid alloy. Here we present the results of classical molecular dynamics simulations of the diffusion of Fe and Ni within LBE. The simulations complement experimental studies of impurity diffusion by our group and provide an atomic-level understanding of the relevant diffusion phenomena. They are based on the embedded atom method (EAM) to represent many-body interactions among atoms. The EAM potentials employed in our simulations have been validated against ab initio density functional calculations. We show that the experimental and simulation results for the temperature-dependent viscosity of LBE and the impurity diffusion coefficients can be reconciled by assuming that the Ni and Fe diffuse mainly as nanoscopic clusters below 1300 K. The average Fe and Ni cluster sizes decrease with increasing the temperature and there is essentially single-atom diffusion at higher temperatures.