Analysis of the relative stabilities of ortho, meta, and para MClY(XC 4H4)(PH3)2 heterometallabenzenes (M = Rh, Ir; X = N, P; Y = Cl and M = Ru, Os; X = N, P; Y = CO)