First-principles Calculations of Atomic Diffusion in Crystalline Solids
This approach (nudged elastic band method) requires knowledge of diffusion mechanisms. A complementary approach using molecular dynamics can help explore diffusion pathways.
Written by Cheng-Wei Lee (clee2 [at] mines [dot] edu)
Atomic diffusion mechanisms
There are multiple textbook diffusion mechanisms and the vacancy-mediated diffusion is one of the most common diffusion mechanisms in ionic compounds. The associated diffusion coefficient can be written as:
D=61​fZf​Zm​l2CD​Γ where is for 3D system ( and ) and , , and are geometry factors, which depend on the crystal structure. is the distance the vacancy travel for each jump, is the concentration of participating defects, e.g. vacancies. Lastly, is the successful jump frequency,
where is the attempted jump frequency, is the migration barrier of a diffusing atom (see figure below)
Therefore, the diffusion coefficient for vacancy-mediated diffusion in crystalline solids can be calculated from first principles:
The attempt frequency can be approximated by the local dynamic matrix of the moving atom and the force due to smaller change in moving atom can be calculated using DFT.
The defect concentration can be calculated using the supercell approach (see defect calculations)
The migration barrier can be calculated using
