Computational identification of Ga-vacancy related electron paramagnetic resonance centers in β-Ga2O3