Thermodynamics of impurities in semiconductors

First-principles electronic-structure calculations of impurities in semiconductors are extended to finite temperatures. We first calculate the vibrational free energy in cubic C, Si, and Ge and hexagonal GaN and show that reliable phonon densities of state, specific heats, and other thermodynamic quantities can be obtained in the same (64-atom) supercells commonly used to study impurities. Then, we use examples in Si to quantify various free energy contributions. Only the vibrational free energy plays a role in the temperature dependence of the relative energy of the two CH₂* complexes, and the configurational entropy dominates when calculating the binding energy of copper pairs, while the rotational free energy is critical to interstitial H₂, HD, and D₂. The contribution to the free energy of electrons (holes) originating from occupied (empty) localized impurity levels in the gap is estimated and found to be very small.