I am investigating the formation energy of defects in hybrid perovskite (MAPbI3), so I need to calculate the chemical potentials of Methylammonium, lead, and iodine.
I am using FHI-aims package in my project.
How to calculate these potentials?
(Note: I am still new in the computational field)