Determining the site preference of trivalent dopants in bixbyite sesquioxides by atomic-scale simulations