First-principles calculations of impurity diffusion coefficients in dilute Mg alloys using the 8-frequency model

Diffusion in dilute Mg–X alloys, where X denotes Al, Zn, Sn and Ca impurities, was investigated with first-principles density functional theory in the local density approximation. Impurity diffusion coefficients were computed as a function of temperature using the 8-frequency model which provided the relevant impurity and solvent (Mg) jump frequencies and correlation factors. Minimum energy pathways for impurity diffusion and associated saddle point structures were computed with the climbing image nudged elastic band method. Vibrational properties were obtained with the supercell (direct) method for lattice dynamics. Calculated diffusion coefficients were compared with available experimental data. For diffusion between basal planes, we find DMg–Ca > DMg–Zn > DMg–Sn > DMg–Al, where D is the diffusion coefficient. For diffusion within a basal plane, the same trend holds except that DMg–Zn overlaps with DMg–Al at high temperatures and DMg–Sn at low temperatures. These trends were explored with charge density contours in selected planes of each Mg–X alloy, the variation of the activation energy for diffusion with the atomic radius of each impurity and the electronic density of states. The theoretical methodology developed herein can be applied to impurity diffusion in other hexagonal materials.

Research highlights
► Implemented the eight frequency model for impurity diffusion in hexagonal metals.
► Model inputs were energetics/vibrational properties from first princples.
► Predicted diffusion coefficients for Al, Ca, Zn and Sn impurity diffusion in Mg.
► Successful prediction of partial correlation factors and jump frequencies.
► Good agreement between calculated and experimental results.