J . Nmn. Sci. Coun. Sri Lanka 1987 15 (1): 47 -59 SELF-CONSISTENT CHARGE AND CONFIGURATION (SCCC) CALCU- LATIONS ON 1,6-DICARBA-CLOSO-HEXABORANE (6), 2,IDICARBA- CLOSO-HEPTABORANE (7) AND THEIR METALLO-DERIVATIVES RAJESWARY MAGESWARAN Department o f Chemistry, University o f Jaffna, Jaffna. and N. J . FITZPATRICK Department of Chemistry, University College, Dublin, Dublin 4, Ireland (Date o f receipt : 14 May 1985)- (Date or acceptauce : 2 March 1987) Abstract : The bonding in 1, 6-dicarba-closo-hexaborane [61 [I , [C B H {P~(co) 11 [I11 .2,4-dicaiba-closo-heptaborane[7] [IIII . [C2B4H6{Pe(Cb):~ IN], ana [C B H { Fe(C0) ) ] [V] was investigated by c a y i n s out self consistent andZc~fiflguration [ ~ C C ] calculations. These investigations show that an apical boroh atom pses 2px, 2p atomic orbitals and a [2s, 2pZl hybrid orbital for bonding with rest of thy cluster. In cluster compounds containing Fe(C0) units, the iron atom uses predominantly [3dz. 4p I and [3d 4p I 3 combinations for bonding with the rest of the cluster and t h s e is no eevigce &r the involvement of a [4s, 4pZ, 3dZ21 hybrid orbital similar t o the [2s, 2pZl hybrid observed in the boron analogue. Analysis of the total overlap population associated with metallo.group and the rest of the cluster shows that the Fe(C0)3 group is weakly bonded ifito the cluster when compared with BH group of compounds in which Fe(C0) is replaced by BH. This analysis also indicates that Fe(C0) .group uses less than 2wo electrons for cluster bonding. 3 1. Introduction . . ' . Understanding of the electronic structure and bonding of metal carbonyl clusters is of present interest. The number of orbitals invoIved in bonding [ > 100 atomic orbitals] makes it difficult to perform detailed molecular orbital calculations on metal carbonyl clusters. But it is believed that certain metallo groups such as Fe(C0)3 and CO(CO)~ are isolobal and isoelectronic h t h the BH and CH groups respectively. In terms of Wade's approacd9 a 3H or a Fe(C0) group formally supplies two electrons for cluster bonding whereas CH and Co(C0) groups supply three each. In addition, many boranes, carboranes, metadoboranes, rnetallocarboranes, and metal carbonyl clusters have similar geometries? 3919 Thus the study of the electronic structure and bgnding of boranes, carboranes, metallo-boranes and metallo- carboranes cauld be used to compare the bonding capabilities of transition metal carbonyl fragments, M(C0)3 with those of the BH or CH Units in cluster compounds. In our study, self-consistent charge and configuration (SCCC) calculations were used to investigate the bonding in 1,6-dicarba- 48 ~ a j e s w a s j Mageswaran and N. J. Fitzpatrick closo-hexaborane [6] [I] , [C2 BI HI {F~(co) ] [11] , 2,4-dicarba-close heptaborane [7] [111] , [CIB4HI { F ~ ( c o ) ~ ~ ] [IV] and [C,B3H, {F~(CO) 3) 1 [Vl 2. Computational Method The method used was the FORT ICON^ of Hoffman and co- worker^?^ This calculation is an all valence electron, extended Huckel calculation [3d, 4s and 4p atomic orbitals ontransition metal atoms, 2s and 2p on second row atoms and on hydrogen atom are used] and was used in its charge iteration mode according to Hii = VSIE [QJ where Hii is the diagonal Hamiltonian matrix element and VSIE[OJ is the valence state ionisation energy of drbital -i when the atom has total charge Q. The off diagonal Hamiltonian matrix elements are calculated in the normal manner9 using the expression, where ki is a constant. The VSIE [Q] functions are assumed t o be of the form, VSIE [Q] = A Q ~ + BQ + C where A,B and C are parameters~6~z0 which depend on the atom and the orbital. Iterations are continued hntd successive ones producedd 1 0 - ~ e change in the charge on any atom [e = electro'hic charge]. The well known problem of obtaining convergence with such self-cornistent charge calcula- tions was eased by an improved damping scheme and by the use of a Made: lung correction. The orbital occupations are summed separately over the p and d orbitals on each atom. The resulting p and d occupations, along with s occupations, are damped according to the equation input i; output k ' inputk = ! + - P ) r where p is the summed orbital occupation of a given type, s,p^ar d, indexed by r on the kth cycle andxis a damping parameter. The geometries of carboranes used were taken from reported,17f18 electron diffraction data. The same geometries were used for metallocar- boranes. The geometry of the FelCO) unit was taken from that of [Fe(CO) $. -C4H I ' O -Fe - C and C-d bond lengths of. 0.182 nm and 0.1 15 nm respectlvefy were used. Changing the position of the metallo-units by 2 10% with respect to the carborane base produced no significant change in the results. - - T ab le 1 . G ro ss a to m ic c ha rg es in m ol ec ul es [ I] - [ V l 6 M ol ec ul es A to m s it es a s sh ow n in d ia gr am s [I ] to [ V ] S 2. 1 2 3 4 5 6 7 C ar bo ny l g ro up 's 4 3 C ar bo n O xy ge n K o d e d e - 3 n [I ] -0 .0 14 0. 09 4 [I1 1 -0 .0 55 0. 05 7 [I1 11 0. 12 4 -0 .0 22 II va l 0. 82 4 -0 .0 62 [I V b] ' 0. 07 8 -0 .0 37 [V al 0. 72 3 -0 .0 74 [V bl -0 .8 48 0. 01 0 d - C ha rg e on t he c ar bo nl ox yg en a to m o f ca rb on yl g ro up w he n th e F e( C 0) un it is a pi ca l. 3 e - C ha rg e on th e ca rb on lo xy ge n at om o f ca rb on yl g ro up w he n th e F e( C 0) un it is b as al . N ot e: 1 ,2 ,3 e tc . c or re sp on ds to th e at om p os it io ns la be lle d 1, 2, 3 et c. in t he ' )d ia gr am s [I ] to [ V ] . R zjes wary Mageswaran and N. J. Fitzpatrick Table 2. Charges localised on BH2 Ch and Fe(C0)3 units in molecules [ I ] - [V] Atom sites as shown in diagrams [I] - [IV] Molecule 1 2 3 4 5 6 7 [Vb] -1.204 -0.033 0.007 -0.118 1.447 -0.060 -0.038 Table 3. Overlap population in molecules [I] and [I11 Bond I I1 3. Results and Discussion The details of the cluster bonding in carbsranes and metallocarboranes can be considered from three pieces of information provided by the SCCC calculation. They are [i] the gross atomic charges [ii] the overlap popula- . tions and [iii] the molecular orbital energy levels and their LCAO expan- sions. 3.1 1,6-C2B,H6 [I] and 2,4-C2B, H, [111] The gross atomic charges (Table 1) and the overlap populations (Tiable 3) indicate that the four BI-I units in 1,6-C,B4H6 are equivalent and so are the two CH units. This is in agreement with the predicted1 D4,, symmetfy of the molecule. The presence of three distinct type of BH units in 2 , 4 Self-consistent Charge and Configu~ation C2 B5 33, is shown by the localised charges [Tables 1 and 21, the boron atom in apical position being more positively charged. The overlap populations [Table 41 bonding these BH units into the cluster [in 2,4C2B, H,] are slightly different and reflect the positive charges on boron atoms! Thus &e two apical BH units ate less strongly bonded into the cluster when compared with the three basal BI-I units. The gross atomic charges on carbons and borons agree with the more electronegative character of carbon compared with boron. Our values for located charges on C--H and B-H moieties [-0.045 and 0.022 respectively] of 1,6-C2B4H, arek comparable with the ' values [-0.05 and 0.03 respectively] reported1 eadier by a different SCF Table 4. Overlap populations in Molecules [III] - [V] Bond [I111 [IVal DVbI [Val ' [Vbl B1[71 - H 0.836 0.828 0.862 - 7 0.856 C' - H 3 0.827 0.861 0.830 0.881 0.852 B -H 4 8.836 0.858 0.838 0.871 0.848 C - H 0.827 0.861 0.856 0.881 0.891 B6[5] - H 0.839 0.860 9.860 0.871 0.881 , B7-H 0.836 0.828 0.860 - 0.856 B 7 ~ 1 ~ - c 2 0.469 0.536 0.518 - 3 2 0.558 B - C 0.776 0.865 0.823 0.985 0.879 B3 - B7[l1 0.355 0.414 0.3 94 - 0.427 B6[51 - *7[11 0.447 0.538 0.5 84 - 6 2 0.657 B - C 0.726 0.878 0.785 0.994 0.908 {B5 - c41 6 5 B - B- 0.693 0.75 1 - 0.842 B7 - c4 0.469 0.536 0.65 1 - 0.712 ~ e ' - c2 - 0.359 - 0.245 0.183 1 3 Fe - B - 0,294 - 0.235 0.147 1 4 Fe - C - 0.359 - 0.245 0.141 Fel - B6t51 - 0.290 - 0.212 0.115 5 7111 - - 0.3 15 - 0.293 Fe - B I7e5 - c2 - - 0.340 - 0.299 ~ e ' - B6 - - 0.285 - 0.217 1 5 Fe - Fe - - - - 0.777 calculation. Three cluster bonding molecular orbitals are found in %@.,B,H6. The HOMO[2b2g] is a B4 combination [E = -11.75eVl followed by a degenerate pa r -[3eu] [E = -12.38eVl which has CZB4 combination. In 2,4C2B5H,, two cluster bonding molecular orbitals ' [E = - 25.95 and -12.81eVI of C B5 combination have been found in addition to the HOMO [E = -11.2e~f which is a combination of mainly the planar carbon and boron atoms. Our results [Table 41 also show that in 2,4C,B, H,, the site tb site fBH] b,d - [BH] overlap population is SCHEME - [I] Raieswary Magestuaran and N. J. Fitzpatrick ( I l l ) X : Y -- Z : BH ( I h ) X : Fe(CO),, Y : Z : 8 H (IVb) X = Y : B H , Z : H(COl, ( M I X = Y = k ( C O \ , 2 : BH (Vb) X - Z F e ( C O 1 , Y : BH Self-consistent Charge and Conjipration 53 - - -- - - - greater than that of [BH] basal - [BH], ic J population and this is in agreement with the difference in observe818 bond lengths [BH] basal - IBH] basal = 0.1659nm, [B'H] - [B'H] = 0.1772nm and [B'H] - [B'H] = 0.1852mml. Let us commence the discussion of metallocarboranes by comparing the orbitals available on pH an$ Fe(C0) units for cluster bonding. In terms of Wade's approacH9 both the BH and the Fe(CO), units use two electrons and three orbitals for cluster bonding. Iron in the Fe(CO), unit has nine valence shell orbitals [c.f. boron which has four valence shell orbitals], three of these are used to form M-C 6 -bonds, three more as well as three electron pairs. are used to form n -b onds to the carbonyl groups. This leaves the Fe atom with three valence shell orbitals and two electrons. SCCC calculations show that in BH the HOMO is a 6 4evel composed of 2s and 2p while LUMOs[e] [2p4 and 2py] are the n-levels. In the case of Fe(CO)', the HOMO is a hybnd [a l ] orbital having 4s, 4pZ and 3dZ2 components and the LUMOs[e] orbitals are formed from [4px, 3dxZ] and [4p , 3dy z ] combination as has been reported!' This is in agreement with waJee's skeletal electron counting technique. Thus the orbitals* available for cluster bonding on BH and Fe(C0) units, as shown in Figure 1, are similar ~JI symmetry as reported1' earlier. [1'~he 2D projections of the HOMO and LUMO orbitals of BH and Fe(C0)3 units obtained by z plot progra- mme1-4,9,16,20 al so had similar shapes] . 3.2 [C,B, H, (F~(co),}I I Replacement of the BH unit is 1,6-C2B4H6 by the Fe(CO), unit, as in [11], v a s e s [Table 11 the net negative charge on carbons by 0.041 and decreases the positive charge on boron by values in the range 0.020-0.037. Table 3 shows that the overlap population bonding the metallo-unit to the C2B H base is about 5 7% of that bonding the BH in 1,6-C B4H6. In other words t i e replacement of a BH unit by a Fe(CO)3 group maftes the resulting compound even more electron-deficient than the parent carborane. A similar result has been observed in borane c1uste1-s.~ Table [I] shows that the iron atom of ferracarborane [II] cairies a positive charge 10.7331 and this may be the reason for the reduction in overlap population betwen the metallo- unit and carborane base. The positive charge on iron atom may be due to the valence electrons in iron atom being largely delocalised over the carbonyl ligands. Replacement of a BH unit in 1,6-C2B4H6. by a Fe(C0)' unit did not change the energy of the cluster orbitals appreciably. The comparison of the cluster bonding molecular orbitals shows that the roles of the boron 2p and 2py. atomic orbitals in carboranes are taken up by the iron 3dxz an3 - 3dyz mth some admixture of 4px and 4p . According to the isolobd principle the 3dz2, 4s and 4pz metal atomic orsitals should combine to form a hybrid orbital which takes on the role of the boron 2s and Zp, orbitals. P In Fe(CO),, the 3dZ2 atomic orbital is involved in bonding but plays no ( a ) BH cluster orbitals ( b , Fa (CO)., c lus te r orbitals Rajeswary Mageswaran and N. J. Fitzpatrick (3d, ,4p,) Combination Figure 1. The orbitals available for cluster bonding on BH and Fe(C0)3 units. Self-consistent Charge and Configuration 5 5 significant part in cluster bonding, whereas, the al cluster molecular orbital of BH unit (which is equivalent to 2s, 2pl hybnd of boron atom] plays a significant part in cluster bonding. A simllar result has been reported8 in metdoboranes. LCAO expansion of HOMO of [11] shows that it has a small contribution from the 3d 2 2 of Fe and carbonyl carbon and oxygen x - Y atomic orbitals. Table 5. Totd cnergy of Valence Electrons . Molecule EnergyIeV [I] -403.8 [I11 -940.0 [11fI -462.3 [ IVal -992.4 [IVbI -987.9 [Val -1538.1 [Vbl -1521.5 3.3 [C2B4H, (F~(co)$] [IV] The replacement of an apical BH unit in 2,4-C2B5H7 as in [IVa] , [il increases the net negative charge on carbon atoms by 0.040 [ii] decreases the positive charges on boron atoms [Table 11 and [iiil decreases the overlap population bonding the replaced site Fe(CO), , c.f.BH, to the C2B H5 base by about 27%. The replacement of a basal BH unit by a ~ e ( d o ) , unit in 2,4-C2 B5 11, as in [IVb] , also increases the net negative charges on carbon atoms but the increase is low [0.015 on c2 and 0.022 on C4] and decrease the positive charges on boron atoms [again the decrease is low, Table 11. The total overlap popdation associated with the metallo- group in [IVb] is markedly less than that in [IVa] and is much smaller than that of the BII replaced [Ca 54%]. There is a corresponding increase in overlap population between adjacent boron and carbon atoms in [IVa] and [IVb] [Table 41 compared with those in 2,4 C2B5H7. The total overlap population associated with the metallo-group and the total energy of valence electrons [Table 51 indicate that the l-isomer [IVa] is more stable than the 5-isomer [IVb] . As mentioned earlier, the 2,4-C B5H7 has three cluster bonding molecular orbitals. two of C2B5 cornginations and the third VOMO] of basal carbon and boron atoms combination. The replacement of either the apicaI BIH or basal B5H by a Fe(C0)3 unit did not change the energy of any of these orbitals appreciably. The comparison of the LCAO expansion of these orbitals shows that there is no significant contribution from the metal 3dZ2, 4s and 4pZ atomic orbitals to these cluster bonding, molecular orbitals and the roles of 2px and 2py atomic orbitals of boron are taken up by iron 3dxz and 3dyz atomic orb~tals with small contributions 5 6 Rajeswary Mageswaran and N.J. Fitzpatrick from 4px and 4p . LCAO expansion of HOMOS of both the ferracarboranes PVa] and [ ~ V b r skiow that they have a small contribution from the 3dx:- y~ atomic orbital of iron and carbonyl carbon and oxygen atomic orbitals. . . The replacement of the second apical BH unit by a Fe(C0) unit as in [Val, increases the negative charge on carbon atoms further but the increase is low [0.012 compared with 0.040 in IVa] and decreases the positive charge on basal boron atoms by similar ratios [Table I]. The total overlap population associated with bonding of Fe(C0)3 unit with C2B3H, base had also been decreased further by about the same ratio [Ca 27%]. The replacement of the basal B" H unit by the Fe(C0)3 unit, as in [Vb] , changes the charge on c2 from -0.062 to +0.010 [but the charge on c4 has changed only slightly from -01062 t o 0.0651 and increases the positive charges on boron atoms [Table 11. Also, it should be noted in [Vb], Fel atom has a negative net charge whereas I?e5 atom has a high positive net charge [Table 11. The total overlap populations associated with the bonding of metallo-groups with the rest of the cluster units in [Vb] are markedly less than those.in [Val and shows that the basal Fe5 (CO) unit is weakly bonded into the cluster unit when compared with t i e apical Fe(C0) unit. The total overlap population associated with the bonding of metallo-groups with the rest of the cluster unit and the total energy of valence electrons [Table 51 suggest that the compound [Val is more stable than compound [Vb] . SCHEME [ Z ] Self-consistent Charge and Configuration It is interesting to note that the compound 1 ,7 ,2 ,4- [ (~5-C,H~)~ Co C2 D H, ] [VI] which is analogous to [Val had been synthesised6 and is stahe. dur calculations also show that the replacement of either B'H or B~ H or B1 H and B~ H by Fe(C0) units did not change the energy of cluster bonding molecular orbitals appreciably. Analysis of LCAO expansion of molecular orbitals for [Val and [Vb] is too complex because of the two Fe(C0)3 units present in these molecules; But LCAO expansion of HOMOS of [Val and [Vb] show that there are contributions from 3d 2 2 atomic' orbitals of iron atoms and atomic orbitals of the carbonyf LLbon and oxygen. 3.5 Coinparison of bonding by BH and Fe(C0)3 units As mentioned earlier in terms of Wade's approach the BH and the Fe(C0)3 units use two electrons each and the CH unit uses 3 electrons for cluster bonding. The overlap population and the atomic orbital contributions in cluster molecular orbitals in carboranes and metallocarboranes studied show that the BH and the CH units use approximately two electrons each and that the Fe(C0)3 unit uses less than two electrons [varies with the position and the number of BH units replaced by Fe(C0)3 units] for cluster bonding. Our results [Tables 3 and 41 also show that the CH unit uses only one electron from carbon for the C-H bond. Hence the fourth electron froni carbon must have become delocalised into the cluster bonding. The total overlap popula- tions associated with the bonding of Fe(C0)3 groups into the cluster in all ferracarboranes studied are mudh smaller than those of the BH replaced. Table 1 shows that in ferracarboranes [except the Fel of (Vb)] the iron atoms cany net gross positive atomic charges and these positive charges on iron atoms may be the reason for the reduction in total overlap population between the Fe(CO)3 unit and the rest of the cluster unit. The analysis of localised charges on the Fe(C0)3 units [Table 21 and the gross atomic charges on the carbonyl carbon and the oxygen atoms [Table I ] suggest that the high positive charges on iron atoms may be due to the valence electrons in iron atom being largely delocalised over the carbonyl ligands and this may be the reason for weak bonding between the Fe(C0)3 unit and the rest of the cluster unit. As discussed earlier, the analysis of LCAO of the cluster bonding molecular orbitals shows that only 3dx,, 3dyZ, 4px and 4p atomic orbitals Y of iron atom are used in cluster bonding but aU Zs, 2p,, 2p and 2pZ orbitals Y of both boron and carbon atoms are used. 4. Conclusion SCCC calculations on carbones [I] and [111] and ferracarboranes [11] [IVa] , [IVb] , [V] and [Vb] show that even though in terms of Wade's approach, the BH and the Fe(C0) units are considered to be isolobal and isoelectronic, 5 8 Rajeswary Mageswaran and N.J. Fitzpatrick the Fe(C0)3 unit is weakIy bonded into the cluster when compared with the BPI unit replaced and this may be due to the valence electrons in iron atom being largely delocalised over the carbonyl ligands in ferracarboranes. Our calculations also show that even though the CH unit appears to use only two electrons for cluster bonding the third electron may be delocalised into the cluster. The comparison of the LCAO of the cluster bonding molecular orbitals shows that even though the roles of 2px and 2p atomic orbitals of Y boron are taken up by the iron 3dxZ and 3dyZ orbitals w t h some admixture of 4px and 4p atomic orbitds there is no evidence for the involvement of [4s, 4pz, 3dz2f hybrid orbital which is analogous to boron [2s, 2pz] hybrid. References 1. BALLHAUSEN, C. J. & GRAY, H. B. (1964). Molecular Orbital Theory, - . W. A. Benjamin Inc. New York, 1.20. 2. BASCH, H. & GRAY, H. B. (1967). Tbeoret. Chem. Acta, 4 : 367 3. BASCH, H., VISTE, A. & GRAY, H. B. (1960) J. Chem. Phys, 44 : 10. 4. BASCH, II., VISTE, A. & GRAY, H. B. (1965). Theoret. Chem. Acta, 3 : 33. 5. BEER, D. c., GRIMES, R. N., SNEDDON, L. G., MILLER, V. R. & WEISS, R. (1974). Inorg. Cbem., 1 3 : 1138. 6 . BEER, D. C. MILLER, V. R., SNEDDON, L. G., GRIMES, R. N., MATHEW, M. & BALENIK, G. J. (1973). J.~mer.Chem. Soc. 95 : 3046. 7. BRINT, P. & SPALDING, T. R. (1979). Inorg. Nuclear Chem. Letters, 1-5 : 355. 8. BRINT, P. & SPALDING, T. R. (1980). J. Chem. Soc. Dalton Trans., 1236. 9. BROWN, D. A., CHAMBERS, W. J. & FITZPATRICK, N. J. (1972). 'inorganic Chemica Acta, 7 and references cited there in. 10. DAVIS, N. T. & SPEED, C. S. (1970). J. Organometal, Chem. 21 : 401'. \11. ELIAN, M., @HEN, M. M. L., MINGOS, D. M. P. &HOFFMAN, R. (1976), Inorg. Chem., 15 : 1148. 12. EPSTEIN, I. R., KOETZE, T. F., STEVENS, R. M. & LIPSCOMB, W. N. (1970). J. Amer. Chem. Soc., 92 : 7019. 13. GRIMES, R. N. (1978). Acc. of Chem. Res., 11 : 420. 14. GRIMES, R. N. & WEISS, R. (1976). J. Organometal. Chem. 27 : 113. Selfconsistent Charge and Confipration . 59 15. HOWELL, J., ROSSI, A., WALLACE, D., HARAKI, K. & HOFFMAN, R. (1977). FORTICON8: Quantum Chemistry Programme Exchange No. 334: 16. McGLYNN, S. P., VANQUIKENBORNE, L. G., KINOSHITA, M. & CARROLL, D. G . Introduction t o Applied Quantum Chemistry Edited by Holt, Rinehart & Winston, Inc., New York, 423-431. 17. McNEILL, E. A., GALLABER, K. L., SCHOLER, F. R. & BAUER, S. H. (1973). Inorg. Chem., 12 : 2108. 18. McNEILL, E. A. & SCHOLER, F. R. (1975). J. Mol. Stmct . , 27 : 151. 19. WADE, K. (1976). Adv. Inorg. Chem. Radiochem., 18 : 1. 20. WHITHEAD, M. A. (1970). Semi-empirical All-ualence-electron SCF-CNDO Theory, edited by Sinanoglu, 0. & Wiberg, K. Yale Press, 50. JNSF 15_1_47.pdf JNSF 15_1_47 (2).pdf JNSF 15_1_47 (3).pdf JNSF 15_1_47 (4).pdf JNSF 15_1_47 (5).pdf JNSF 15_1_47 (6).pdf JNSF 15_1_47 (7).pdf JNSF 15_1_47 (8).pdf JNSF 15_1_47 (9).pdf JNSF 15_1_47 (10).pdf JNSF 15_1_47 (11).pdf JNSF 15_1_47 (12).pdf JNSF 15_1_47 (13).pdf