This work systematically explores 19 unique configurations of the close-associate Ga-O divacancies (V Ga V O ) in β-Ga 2 O 3 , including their complexes with H impurities, using hybrid functional calculations. Interestingly, most configurations are found to retain the negative-U behavior of V O , as they exhibit a thermodynamic (−/3−) charge-state transition level energetically located in the upper part of the band gap, where the 3− charge state is associated with the formation of a Ga-Ga dimer. The energy positions of the thermodynamic (−/3−) chargestate transition levels divide the divacancy configurations into three different groups, which can be understood from the three possible Ga-Ga dimerizations resulting from the tetrahedral and octahedral Ga sites. The relative formation energies of the different divacancy configurations, and hence the electrical activity of the divacancies, is found to depend on the Fermi-level position, and the energy barriers for transformation between different divacancy configurations are explored from nudged elastic band calculations. Hydrogenation of the divacancies is found to either passivate their negative-U charge-state transition levels or shift them down in Fermi level position, depending on whether the H resides at V O or forms an O-H bond at V Ga , respectively. Finally, the divacancy is discussed as a potential origin of the so-called E * 2 center previously observed by deep-level transient spectroscopy.