Diffusion of interstitials in metallic systems, illustration of a complex study case: aluminum

While diffusion mechanisms of interstitial elements in fcc systems are generally well-known, especially in the case of H atoms, we show in this work that even in the case of a simple metallic system (aluminum), the diffusion of interstitials exhibits a wide variety of paths and mechanisms that depend on the specie. We used an approach based on first-principles calculations associated with kinetic Monte-Carlo simulations and a multi-state diffusion formalism to compute the diffusion coefficients of five interstitial elements: hydrogen, boron, carbon, nitrogen and oxygen. For instance, at the atomic scale, whilst we find that C atoms prefer to be located in octahedral sites (labeled o) rather than in tetrahedral positions (labeled t), we find one additional stable position in the lattice (M). The diffusion through these three stable positions are thus studied in detail. In the case of B atoms, for which the tetrahedral site is found unstable, the diffusion path is between o-o sites. Similarly, in the case of oxygen, t positions are found to be the only stable positions (o are unstable) and the path of migration, along t-t direction, is found through a twice degenerated asymmetric transition state. In the case of H and N atoms for which t and o sites are stable, we explain why the only path is along the t-o direction. Finally, we discuss explicit formulas to compute coefficients of diffusion of interstitials in fcc structures.