where the net charge of a system takes into account all defect species (X) with charge q, free electrons (n e ) and free holes (n h ). The free carrier concentrations are determined according to the Fermi-Dirac distribution function. The concentration c of a defect (in thermodynamic equilibrium) at temperature T is related to its formation energy (E f ):