SCFMO calculations of heteroatomic systems with the variableβ approximation