A reaction model explaining (i) the hydrogen-mediated transformation of the vacancy-oxygen (V O) center into a vacancy-oxygen-hydrogen center (V OH * ), with an energy level at 0.37 eV below the conduction-band edge (E c ), and (ii) the passivation of the divacancy center is presented. V OH * dissociates with a rate of 2 × 10 −5 s −1 at 195 • C, causing V O to recover after long duration (>10 4 min), while a similar evolution occurs at 300 • C on a time scale of the order of 10 min. The diffusivity of the monatomic hydrogen used in the model agrees closely with the established values for the diffusivity of protons. After the recovery of V O, further annealing at higher temperatures and/or longer durations transforms V O into the "ordinary" vacancy-oxygen-hydrogen center with an energy level at E c − 0.32 eV (V OH). V OH is subsequently transformed into V OH 2 . For temperatures above 250 • C, two additional hydrogen-related levels occur (∼ 0.17 and ∼0.58 eV below E c ) with a one-to-one ratio and a possible association with different charge states of a V 2 OH center is discussed.