I like to use PyDII for defect density calculation using VASP, In my case we like to study the defect formation energy of Zn4Sb3 and Ce doped Zn4Sb3, since we believe doping induced defects enhances the thermopower. Can you please comment on the possibility of studying the defect formation energy ??