I am studying the ionic conductivity of LiGe2(PO4)2 computationally by doping using different dopants. To do this, I need to calculate the defect energy after doping LiGe2(PO4)2. The defect energy is calculated as:
Edefect = E pure – Edoped + niμi
Where Edefect and Epure are the total formation energy of the LiGe2(PO4)3 with and without the dopants, respectively; ni is the number of atoms of element i added to (or removed from) LiGe2(PO4)3 to create and charge balance the supercell; and µi is the chemical potential of element i.
So could you help me how to calculate µi using VASP? Thank you.