Thermodynamic modeling of native defects in ZnO