Cu diffusion in α-Fe: determination of solute diffusivities using atomic-scale simulations

The nucleation and growth of radiation-induced Cu-rich clusters have proved an important topic in the study of microstructural and chemical evolution responsible for property changes in pressure vessel steels subjected to neutron irradiation. Small Cu clusters can act as precipitate precursors and contribute to the hardening of the material by impeding dislocation glide. Reliable Cu diffusion coefficients provide critical input to Monte Carlo and rate theory approaches to study the kinetics of diffusion and clustering that lead to the formation of such Cu-rich features. In this paper we report the results of a molecular dynamics study of Cu solute-atom diffusion in a Fe–0.9at%Cu alloy, using high-temperatures and high-vacancy concentrations to hasten the convergence of the calculation. We find a value of 0.61 eV for the migration energy of Cu in the b.c.c Fe matrix, in good agreement with available experimental results. The main diffusion mechanisms under the conditions simulated are identified and discussed, and the effect of vacancy and solute clustering on Cu migration is assessed. In addition, a rigorous computation of the solute enhancement factor, correlation factor and solute diffusivity is performed within the five-frequency theoretical framework and the harmonic approximation. The results obtained are compared with those yielded by the direct MD simulations and the observed differences discussed.