First principles based phase diagram calculations for the octahedral-interstitial solid solution system ZrOx