DOI:
10.1039/B512700A
(Paper)
Dalton Trans., 2006, 686-692
The relative stabilities of 11-vertex nido- and 12-vertex closo-heteroboranes and -borates: facile estimation by structural or connection increments†
Received
7th September 2005
, Accepted 27th October 2005
First published on 22nd November 2005
Abstract
The relative thermodynamic stabilities of ortho-, meta- and para-isomers of 12-vertex closo-heteroboranes and -borates with different p-block heteroatoms were determined using density functional theory. More electronegative (smaller) heteroatoms tend to occupy non-adjacent, whereas less electronegative (larger) heteroatoms tend to occupy adjacent vertices in the thermodynamically most stable closo-diheterododecaborane isomers. The computed relative stabilities agree perfectly with experimental observations. The energy differences of para- and meta- relative to ortho-isomers of 12-vertex closo-heteroboranes generally depend on the extent of electron localization by a given heteroatom and show highly periodic trends, i.e. increase along the period and decrease down the group. The energy penalties for the HetHet structural feature (two heteroatoms adjacent to each other) for the 12-vertex closo-cluster are apparently significantly different from those for the 11-vertex nido-cluster. Reformulating two 11-vertex nido-structural features, i.e. Het5k(2) and HetHet, in terms of connection increments along with the additional structural feature HetHetm give the relative stabilities of various isomeric 11-vertex nido- as well as 12-vertex closo-heteroboranes and -borates, using one unique set of increments.
1 Introduction
Beside their potential use in medical applications,1,2 heteroboranes are of interest due to their unusual non-classical structure and their key role in stimulating new concepts and a general picture of chemical bonding.3 Simple qualitative rules presented by Williams4–6 and Wade7,8 are helpful for both understanding the building principles and to identify possible synthetic targets. More sophisticated quantitative rules also called structural increments arising due to disfavoring structural features in heteroboranes3,9–12 have also been presented. The additive nature of these structural increments, i.e. their ability to accurately reproduce the DFT computed relative stabilities is now proven to be generally true for various nido-heteroboranes and -borates with diverse numbers, connectivities and types of heteroatoms. Sets of structural features and corresponding energy penalties were reported for 6-,9 10-10 and 11-vertex nido-heteroboranes and -borates.3,11,12 They allow to easily derive the relative stabilities of various isomers with good accuracy, once corresponding energy increments “penalties” are assigned to relevant structural features based on DFT results computed for a selected set of structures. The energy penalties for the structural features HetHet (two adjacent heteroatoms) and HetC (a heteroatom adjacent to a carbon atom) possess periodic trends (increase along the period and decrease down the group) for various p-block heteroatoms in 11-vertex nido-heteroboranes and -borates.3 Here, we extend the study of HetHet and HetC structural features to the 12-vertex closo-cluster for comparison with those from the 11-vertex nido-cluster. Both clusters have identical numbers of skeletal electrons (26) and hence have allied cluster shapes: the 11-vertex nido-cluster is derived from the 12-vertex closo-icosahedron by elimination of one vertex. Therefore, the influence of heteroatoms can be expected to be similar in both cases, particularly the positions of various heteroatoms relative to each other. For closo-compounds the number of possible isomers is very limited due to the high symmetry of the spherical cluster shapes. In 12-vertex closo-diheteroboranes, each structural feature is present only in one isomer (e.g. 1,2-C2B10H12 has ortho carbon atoms, its 1,7- and 1,12-isomers have carbon atoms in meta and para positions, respectively). On the basis of these, it is therefore impossible to decide if the destabilization due to certain structural features is special for one isomer or has a general meaning. For the 11-vertex nido-clusters, the number of possible isomeric structures is vast due to the less symmetric cluster shape and due to the presence of additional hydrogen atoms. It has already been shown that the relative energies can, nevertheless, be estimated quite accurately by a small number of energy penalties related to certain structural features. Here we try to find quantitative rules that apply both for 12-vertex closo- and 11-vertex nido-compounds.
2 Computational details
All structures were optimized at the RB3LYP13,14/LANL2DZ level of density functional theory using the Gaussian 03 program15 with d-type polarization functions16 for B (ζ = 0.388), Al (ζ = 0.198), Ga (ζ = 0.207), In (ζ = 0.160), Tl (ζ = 0.146), C (ζ = 0.600), Si (ζ = 0.262), Ge (ζ = 0.207), Sn (ζ = 0.183), Pb (ζ = 0.164), N (ζ = 0.864), P (ζ = 0.340), As (ζ = 0.293), Sb (ζ = 0.211) and Bi (ζ = 0.185) followed by frequency calculations and zero-point energy calculations at the same level. Single-point energies of the optimized geometries were computed at the RB3LYP/SDD level with additional d-type polarization functions (see above), p-type polarization function for hydrogen atoms (ζ = 1.000) and an sp set of diffuse functions for B, C and N (ζ = 0.0315, 0.0438 and 0.0639, respectively),17 for Al, Si and P (0.0318, 0.0331 and 0.0348, respectively)18 and for Ga, Ge, As, In, Sn, Sb, Tl, Pb and Bi (0.0205, 0.0222, 0.0287, 0.0223, 0.0231, 0.0259, 0.0170, 0.0171 and 0.0215, respectively).19 The relative energies reported for isomeric structures correspond to the RB3LYP/SDD//RB3LYP/LANL2DZ + ZPE level.
3 Results and discussion
3.1 Periodic trends of heteroatom energy penalties for the 12-vertex closo-cluster
Thermodynamic stabilities of Het2B10H10 isomers have been determined where Het = two-electron donating group 13 heterogroup, i.e. H–Al, H–Ga, H–In, H–Tl, three-electron donating exo-substituted group 14 heteroatoms, i.e. H–C, H–Si, H–Ge, H–Sn and H–Pb, and three-electron donating group 15 heteroatoms, i.e. N, P, As, Sb and Bi. The thermodynamic stabilities of the diheterododecaboranes and -borates depend on the position of the heteroatoms relative to each other and can be expressed in the following terms: HetHeto is the energy difference of the ortho
(1,2)- with respect to the para
(1,12)-isomer (Fig. 1) while HetHetm is the energy difference of the meta
(1,7)- with respect to the para (1,12)-isomer. HetHetm
→
o, the energy difference of the ortho
(1,2)- with respect to the meta
(1,7)-isomer results as the difference of HetHetm and HetHeto (Chart 1). Generally, HetHeto and HetHetm
→
o have larger values than HetHetm for diheteroboranes and -borates. The thermodynamic stabilities of the diheterododecaboranes and -borates also depend upon the position of heteroatoms in the Periodic Table. Energy penalties for HetHeto (two heteroatoms adjacent to each other) decrease down the group. For example, among group 14 heteroatoms (Fig. 2, Chart 1), preference of paravs.ortho of C2B10H12 (+18 kcal mol−1) is diminished in Si2B10H12 (+4 kcal mol−1) and reversed for the Ge homolog (−2.4 kcal mol−1). The ortho preference is even more pronounced for Sn (−13 kcal mol−1) and Pb (−26.2 kcal mol−1). The same trend is observed for N (+25.0 kcal mol−1) → Bi (−12.6 kcal mol−1). Among exo-substituted group 13 heteroatoms, aluminium has the largest HetHeto energy penalties while thalium has the smallest.
 |
| Fig. 1 Structural features HetHet′o, HetHet′m and HetHet′m
→
o for 12-vertex closo-HetHet′B10H10. | |
 |
| Fig. 2 HetHeto energy penalties for 12-vertex closo-clusters generally decrease down the group and increase along the period. The energy penalty for NNo (included in parenthesis) was obtained by fixing the N–N bond distance to be 1.775 Å in the ortho-N2B10H10 isomer because an icosahedral structure without fixed bond distances did not survive geometry optimization. | |
 |
| Chart 1 Periodic trends of HetHet and HetC energy penalties in 12-vertex closo-hetero(carba)boranes and -borates. | |
HetHeto energy penalties generally increase (become more positive) along one period; e.g. along the second period, energy penalties increase from zero (for two adjacent boron atoms) to 15.9 kcal mol−1 for two adjacent carbon atoms. The energy penalty for NN (two adjacent nitrogen atoms) was obtained by fixing N–N bond distance in the ortho isomer to be 1.775 Å, as otherwise the cluster distorted upon full optimization to have a four-membered open face. HetHeto energy penalties increase along Al to P in the 3rd period, Ga to As in the 4th, In to Sb in the 5th and Tl to Bi in the 6th period (Fig. 2, Chart 1).
Thermodynamic stabilities of HetCB10H11 isomers (where Het = heteroatom) were also computed and corresponding values for the relative placement of a heteroatom adjacent to a carbon atom (HetCo) are also included in Chart 1 (also see Appendix I, ESI†). The energy range of HetCo energy penalties is much smaller and the values generally diminish down the 13th, 14th and 15th group, although the trends are less strictly followed. For example, Ga has a higher energy penalty than Al in group 13, and Ge has a higher HetCo energy penalty than Si in group 14. HetCo values, however, always become more positive along the 3rd, 4th, 5th and 6th period, without any exception.
HetHetm
→
o and HetCm
→
o energy penalties for 12-vertex closo-Het2B10H10 and HetCB10H11 clusters also decrease down the group and increase along the period (Chart 1).
Periodic trends of energy penalties can be explained on the basis of the extent of electron localization due to a heteroatom. Within the same group, all heteroatoms formally donate the same number of electrons to the cluster. However, the larger electronegativities of smaller heteroatoms, generally result in enhanced electron localization and hence larger energy penalties. There is a steady increase in HetHet and HetC energy penalties of group 14 members as compared to group 13 members but a less pronounced increase is observed for group 15 members (see the HetHeto and HetCo curves in Fig. 2 and Appendix I of ESI†). This is because group 13 heteroatoms formally localize two electrons at one vertex but group 14 and 15 heteroatoms localize three electrons. The increase in energy penalties for group 15 heteroatoms as compared to group 14 heteroatoms is, however, due to higher electronegativity, that further enhances the extent of electron localization.
3.2 Redefining Het5k(2) and HetHet in terms of connection increments
The HetHetm
→
o energy penalties for group 14 heteroatoms in 11-vertex nido- and equivalent HetHetm
→
o in 12-vertex closo-clusters are listed next to each other in Table 1 (also see Fig. 3). The CCm
→
o energy penalty for two neighboring carbon atoms is very similar for both closo- and nido-clusters. For all other heteroatoms, HetHetm
→
o values are surprisingly smaller for the 12-vertex closo-cluster as compared to HetHetm
→
o for the 11-vertex nido-cluster. HetHetm
→
o values reported for the nido-cluster are all positive (meaning heteroatoms apart are generally more favorable) but a number of heteroatoms (especially large ones) have negative HetHetm
→
o (and also HetHeto) values for the 12-vertex closo-cluster (meaning ortho isomers are more stable). Furthermore, the difference of HetHetm
→
o energy penalties between the 11-vertex nido- and 12-vertex closo-cluster (Δ) increases with the size of the heteroatom (Sn > Ge > Si > C, Fig. 3). HetHet energy penalties are also significantly different for the 11-vertex nido- and the 12-vertex closo-cluster for group 15 heteroatoms (Table 1, Fig. 3). The difference of energy penalties of a given heteroatom for the two clusters increases from P to Sb (Fig. 3). This apparently means different energy penalties for the same structural features in 11-vertex nido- and 12-vertex closo-cluster, in spite of very similar cluster structure and the same electronic requirement of 13 skeletal electron pairs for each cluster. This prompted us to consider increments for individual “bonds” (or better bonding connections) rather than for structural features for a universal treatment of both cluster types.
Table 1 Direct comparison of HetHet energy penalties (kcal mol−1) for exo-substituted group 14 heteroatoms in 12-vertex closo- and 11-vertex nido-clusters
Heteroatom |
12-Vertex closo- |
11-Vertex nido-a |
Δ
|
Group 13 heteroatoms |
H–C |
16 |
16 |
0 |
H–Si |
3 |
9 |
6 |
H–Ge |
−3 |
4 |
7 |
H–Sn |
−7 |
3 |
10 |
Group 14 heteroatoms |
Values taken from ref. 12.
The NNm
→
o energy penalty was obtained by fixing the N–N bond distance in ortho-N2B10H10 to 1.775 Å. The full optimization of ortho-N2B10H10 led to a distorted structure with a four-membered open face where the nitrogen atoms occupy opposite sites (reminiscent of a diamond–square–diamond intermediate).
|
Nb |
25 |
41 |
32 |
P |
4 |
11 |
7 |
As |
−1 |
7 |
8 |
Sb |
−7 |
4 |
11 |
 |
| Fig. 3 HetHetm
→
o energy penalties for 11-vertex nido- and 12-vertex closo-clusters of group 14 and 15 heteroatoms. The energy penalties for NNm
→
o (included in parenthesis) was obtained by fixing the N–N bond distance to be 1.775 Å in ortho-N2B10H10 isomer as an icosahedral structure without fixed bond distances did not survive geometry optimization. | |
Fig. 4(a) shows two isomeric [CB10H11]3− structures, i.e. 7-[CB10H11]3− (A) and 2-[CB10H11]3− (B). The former has a carbon atom at the open face while the latter has the carbon atom at an unfavorable five-coordinate vertex (vertex number 2). Hence, the difference may be described by the disfavoring structural feature C5k(2) with an increment of 28 kcal mol−1 in the latter (B), while the former (A) is a structure without any disfavoring structural feature (carbon at the lowest possible coordinated vertex).11 The first structure, 7-[CB10H11]3−, described in terms of heteroatom cluster connections has four C–B bonding connections while the latter has five C–B bonding connections. The presence of one additional C–B connection in B as compared to A has to reflect the stability difference of 28 kcal mol−1. Hence, the C–B connection increment is attributed a (destabilizing) value of 28 kcal mol−1.20
![Comparison of (a) 7-[CB10H11]3− and 2-[CB10H11]3− as well as (b) 7,9-[C2B9H11]2− and 7,8-[C2B9H11]2− using (I) the structural increment approach and (II) the connection increment approach (all energy values in kcal mol−1). The values reported in Table 2 are slightly different due to statistical fitting procedure.](/image/article/2006/DT/b512700a/b512700a-f4.gif) |
| Fig. 4 Comparison of (a) 7-[CB10H11]3− and 2-[CB10H11]3− as well as (b) 7,9-[C2B9H11]2− and 7,8-[C2B9H11]2− using (I) the structural increment approach and (II) the connection increment approach (all energy values in kcal mol−1). The values reported in Table 2 are slightly different due to statistical fitting procedure. | |
The 16 kcal mol−1 energy difference of 7,9-[C2B9H11]2− (C) and 7,8-[C2B9H11]2− (D) (Fig. 4(b)) is due to the presence of adjacent carbon atoms (structural feature CC) in the latter with an energy penalty of 16 kcal mol−1.11 These two structures can also be distinguished on the basis of connection increments: 7,9-[C2B9H11]2− has eight C–B bonding connections while 7,8-[C2B9H11]2− has six C–B bonding connections in addition to one C–C bonding connection. 7,9-[C2B9H11]2− accumulates 224 kcal mol−1 due to eight C–B (8 × 28) and is 16 kcal mol−1 less stable than the 7,8-isomer. In order to reproduce the relative energy, the latter must have a total of 240 kcal mol−1 of connection increment energy. Six C–B contribute 168 kcal mol−1, the remaining 72 kcal mol−1 are deduced as connection increment for C–C.
Two structures differing in para and meta positions of two carbon atoms may be differentiated on the basis of CCm (Fig. 5). It cannot be treated in terms of connection increments, due to the same number of C–B connections in the para- and meta-isomers (Fig. 5). However, this structural feature has a very small energy penalty (2 kcal mol−1) for both 11-vertex nido- and 12-vertex closo-clusters.21 Both 2,9-[C2B9H11]2− (E) and 2,8-[C2B9H11]2− (F) have nine C–B connections. Four C–B connections are due to a carbon atom at the open face and the remaining five C–B connections involve a carbon atom in the middle belt. Likewise closo-C2B10H12 has ten C–B connections for both the para (G) and the meta-isomer (H).
![Structural feature CCm differentiates between 2,9- (E) and 2,8-[C2B9H11]2− (F) and 1,12- (G) and 1,7-C2B10H12 (H) (all values in kcal mol−1).](/image/article/2006/DT/b512700a/b512700a-f5.gif) |
| Fig. 5 Structural feature CCm differentiates between 2,9- (E) and 2,8-[C2B9H11]2− (F) and 1,12- (G) and 1,7-C2B10H12 (H) (all values in kcal mol−1). | |
Upon statistical fitting to Ecalc of the isomers reported in Table 2, the energy penalties for C–B, C–C and CCm change slightly to 27, 71 and 1 kcal mol−1, respectively. These energy penalties suffice to reproduce the relative energies of various 11-vertex nido-carborates and 12-vertex closo-carboranes and -borates (Table 2). For instance, 7,8,10-[C3B8H11]− has ten C–B connections, one C–C connection and twice the structural feature CCm. ∑Einc (the sum of increments) is 343 kcal mol−1. 7,8,9-[C3B8H11]− has eight C–B connections, two C–C connections and one CCm structural feature. The sum of increments (∑Einc) is 359 kcal mol−1 with Eincrel (the relative energy of the 7,8,9-isomer as compared to the 7,8,10-isomer) is 16 kcal mol−1 as predicted by the increment scheme. Ecalc (the DFT computed relative energy) comes out to be 19 kcal mol−1 for the 7,8,9-isomer as compared to the 7,8,10-isomer. The difference in relative energy between the two methods is 3 kcal mol−1. The same set of increments has been applied to 11-vertex nido-carborates with one, two and three carbon atoms and also to closo-carboranes with two and three carbon atoms (Table 2).
Table 2 The relative stabilities (kcal mol−1) of 11-vertex nido-a and 12-vertex closo-carborane isomers based on a common set of increments (values in kcal mol−1 are specified in the column headers)
|
|
C–B |
C–C |
CCm |
|
|
|
|
|
Formula |
27 |
71 |
1 |
∑Einc |
E
inc
rel
|
E
calc
b
|
ΔEc |
If endo-hydrogen atoms are present, additional increments as reported in ref. 11 may be needed. These can be used together with the bonding connection increments reported here.
Computed relative energies for 11-vertex nido-carboranes and -borates are taken from ref. 11.
ΔE is the energy difference of Eincrel and Ecalc.
|
nido
|
7-[CB10H11]3− |
4 |
|
|
108 |
0 |
0 |
0 |
|
2-[CB10H11]3− |
5 |
|
|
135 |
28 |
26 |
2 |
|
7,9-[C2B9H11]2− |
8 |
|
1 |
136 |
0 |
0 |
0 |
|
7,8-[C2B9H11]2− |
6 |
1 |
|
216 |
14 |
16 |
−2 |
|
2,9-[C2B9H11]2− |
9 |
|
|
233 |
26 |
27 |
−1 |
|
2,8-[C2B9H11]2− |
9 |
|
1 |
260 |
28 |
29 |
−1 |
|
2,7-[C2B9H11]2− |
7 |
1 |
|
244 |
42 |
43 |
−1 |
|
7,8,10-[C3B8H11]− |
10 |
1 |
2 |
343 |
0 |
0 |
0 |
|
7,8,9-[C3B8H11]− |
8 |
2 |
1 |
359 |
16 |
19 |
−3 |
|
closo
|
1,12-C2B10H12 |
10 |
|
|
270 |
0 |
0 |
0 |
|
1,7-C2B10H12 |
10 |
|
1 |
270 |
1 |
2 |
−1 |
|
1,2-C2B10H12 |
8 |
1 |
|
287 |
17 |
18 |
−1 |
|
1,7,9-[C3B9H12]+ |
15 |
|
3 |
408 |
0 |
0 |
0 |
|
1,2,12-[C3B9H12]+ |
13 |
1 |
1 |
423 |
15 |
14 |
1 |
|
1,2,8-[C3B9H12]+ |
13 |
1 |
2 |
424 |
16 |
17 |
−1 |
|
1,2,4-[C3B9H12]+ |
11 |
2 |
1 |
440 |
32 |
33 |
−1 |
|
1,2,3-[C3B9H12]+ |
9 |
3 |
|
456 |
48 |
49 |
−1 |
Following the same approach, the relative stabilities of various other 11-vertex nido- and 12-vertex closo-heteroboranes and -borates (Appendices II through V, ESI†) can successfully be estimated by a single set of increments.
This connection increment approach can be further expanded to 11-vertex nido- and 12-vertex closo-heterocarbaboranes and -borates with two different heteroatoms (i.e. Het and C). Table 3 indicates how the relative stabilities of germaboranes and -borates, germacarbaboranes and -borates can be reproduced using the additional connection increments, Ge–C, Ge–B, Ge–Ge, and structural increments CGem and GeGem. The method can be used to get the relative stabilities of other heterocarbaboranes and -borates.
Table 3 The relative stabilities (kcal mol−1) of 11-vertex nido- and 12-vertex closo-germaborane and germacarbaborane isomers based on a common set of increments (values in kcal mol−1 are specified in the column headers)
|
|
C–B |
Ge–B |
Ge–Ge |
GeGem |
Ge–C |
GeCm |
|
|
|
|
|
Formula |
28 |
44 |
90 |
2 |
82 |
1 |
∑Einc |
E
inc
rel
|
E
calc
a
|
ΔE |
E
calc values for 11-vertex nido-germaboranes and -borates and germacarbaboranes and -borates are taken from ref. 12.
|
nido
|
7-[GeB10H11]3− |
|
4 |
|
|
|
|
176 |
0 |
0 |
0 |
|
2-[GeB10H11]3− |
|
5 |
|
|
|
|
220 |
44 |
44 |
0 |
|
7,9-[Ge2B9H11]2− |
|
8 |
|
1 |
|
|
354 |
0 |
0 |
0 |
|
7,8-[Ge2B9H11]2− |
|
6 |
1 |
|
|
|
354 |
1 |
4 |
−4 |
|
7,9-[GeCB9H11]2− |
|
4 |
|
|
|
1 |
289 |
0 |
0 |
0 |
|
7,8-[GeCB9H11]2− |
|
3 |
|
|
1 |
|
298 |
9 |
8 |
1 |
|
closo
|
1,12-GeCB10H12 |
5 |
5 |
|
|
|
|
360 |
0 |
0 |
0 |
|
1,7-GeCB10H12 |
5 |
5 |
|
|
|
1 |
361 |
1 |
4 |
−3 |
|
1,2-GeCB10H12 |
4 |
4 |
|
|
1 |
|
370 |
10 |
15 |
−5 |
|
1,12-Ge2B10H12 |
|
10 |
|
|
|
|
440 |
0 |
0 |
0 |
|
1,7-Ge2B10H12 |
|
10 |
|
1 |
|
|
442 |
2 |
0 |
1 |
|
1,2-Ge2B10H12 |
|
8 |
1 |
|
|
|
442 |
2 |
−2 |
4 |
Some 11-vertex nido-heteroboranes and -borates like 2-[NB10H11]2−,3 2-[SnB10H11]3− and 2-[AsB10H11]2−,12 did not optimize to a regular 11-vertex nido-cluster due to cluster distortion.3 Likewise differences of the estimated against the DFT computed relative energy larger than 5 kcal mol−1 result for 2,7-, 2,8- and 2,9-[Sn2B9H11]2− structures. They all have a tin atom at vertex number 2 and show significant cluster distortion (Appendix VII, ESI†).
The strength of bonding interaction of one heteroatom with its cluster neighbors should depend on its overall connectivity. Bond increments therefore should also depend on the connectivity. However, as we see here, at least for fourfold vs. fivefold cluster connections, the difference seems to be small enough for the approach to work well.
3.3 Heavy heteroatoms at adjacent positions in the thermodynamically most stable 12-vertex closo-isomer; a phenomenon supported by experimental results
The qualitative heteroatom placement rules of Williams4–6 have long remained the only theoretical tool to quickly determine the positions of heteroatoms in the thermodynamically most stable heteroboranes and -borates. These rules suggest least connected, non-adjacent vertices for heteroatoms in the thermodynamically most stable heteroborane isomers. DFT computations coupled with our quantitative structural increment approach, however, helped to rationalize the presence of heteroatoms at highly connected vertices in the thermodynamically most stable isomers for a number of structures in 10-10 and 11-vertex nido-heteroboranes and -borates.11,3 Here, our computations demonstrate that if equivalently connected vertices are available (as in the 12-vertex closo-cluster), heavy heteroatoms tend to occupy adjacent vertices in the thermodynamically most stable isomer. Chart 1 shows larger (more positive) HetHet energy penalties for smaller heteroatoms. As a consequence, they rearrange to the more stable meta and para isomers upon heating (Fig. 6) as is known from experiments.22 For example, 2,1-PCB10H11 undergoes thermal rearrangement at 485 °C to form the 1,7-isomer,23 while at higher temperature (650 °C) significant amounts of the 1,12-isomer are formed.24 However, numerous experimental structures were reported with two heavy heteroatoms in the ortho position such as 1,2-As2B10H10,25 1,2-Sb2B10H10,25,26 1,2-SbAsB10H1026 or 1,2-HetBiB10H10 (Het = P, As, Sb, Bi)27 but no meta or para-rearrangements were reported. As seen from smaller (more negative) HetHet energy values, the ortho-isomers are favored as compared to meta- and para-isomers for larger heteroatoms (Fig. 6) and no isomerization takes place upon heating the ortho compounds.
 |
| Fig. 6 (a) ortho-2,1-PCB10H11 undergoes thermal rearrangement to give meta and para isomers due to a PCo energy penalty of 15.6 kcal mol−1. (b) ortho-Bi2B10H10 does not undergo any rearrangement: it is the most stable isomer as reflected by a negative BiBio energy penalty (−12.6 kcal mol−1). | |
The stability order of para > meta > ortho is not generally valid! Chart 1 shows that this is only the case for Het = C, Si, N, P.
Note that in the third column in Chart 1 (group 14 members) there is no parallel trend of para
→
ortho and meta
→
ortho quantitatively, i.e., both become more negative when going down, but para
→
ortho much more. As a consequence para
→
meta becomes really significant. In all other cases, it remains small.
4 Conclusion
A single connection increment scheme is presented that can give the relative stabilities of 11-vertex nido- as well as 12-vertex closo-clusters with different heteroatom substitution patterns. The Het5k and HetHet structural increments proposed for the 11-vertex nido-cluster3,11 may be transformed to bonding connection increments that can be used for 12-vertex closo-clusters. Connection increments are still applicable to 11-vertex nido-heteroboranes and -borates along with other structural features. Hence, for the 11-vertex nido-cluster, the resulting estimated relative energies have exactly the previously reported values. The energy differences between para- to ortho-, and meta- to ortho-12-vertex closo-diheteroborane clusters follow strictly periodic trends. Smaller heteroatoms tend to occupy non-adjacent, while larger heteroatoms tend to occupy adjacent vertices in thermodynamically most stable diheteroborane isomers.
Acknowledgements
Financial support by DFG is gratefully acknowledged.
References and notes
- M. F. Hawthorne and A. Maderna, Chem. Rev., 1999, 99, 3421–3434 CrossRef CAS.
- E. A. Mizusawa, H. L. Dahlman, S. J. Bennet and M. F. Hawthorne, Proc. Natl. Acad. Sci. USA, 1982, 79, 3011–3014 CAS.
- F. A. Kiani and M. Hofmann, Inorg. Chem., 2005, 44, 3746-3754.
- R. E. Williams, J. Am. Chem. Soc., 1965, 87, 3513–3515 CrossRef CAS.
-
R. E. Williams, in Progress in Boron Chemistry, ed. R. J. Brotherton and H. Steinberg, Pergamon Press, England, 1970, vol. 2, ch. 2, p. 57 Search PubMed.
- R. E. Williams, Chem. Rev., 1992, 92, 177–207 CrossRef CAS and references therein.
- K. Wade, Adv. Inorg. Chem. Radiochem., 1976, 18, 1–66 CAS.
-
K. Wade, in Metal Interactions with Boron Clusters, ed. R. N. Grimes, Plenum Press, New York, 1982, ch. 1, pp. 1–41 Search PubMed.
- M. Hofmann, M. A. Fox, R. Greatrex, P. v. R. Schleyer and R. E. Williams, Inorg. Chem., 2001, 40, 1790–1801 CrossRef CAS.
- F. A. Kiani and M. Hofmann, Eur. J. Inorg. Chem., 2005, 12, 2545–2553 CrossRef.
- F. A. Kiani and M. Hofmann, Inorg. Chem., 2004, 43, 8561–8571 CrossRef CAS.
- F. A. Kiani and M. Hofmann, J. Mol. Model., 2005 DOI:10.1007/s00894-005-0037-3.
- A. D. Becke, J. Chem. Phys., 1993, 98, 5648–5652 CrossRef CAS.
- C. Lee, W. Yang and R. G. Parr, Phys. Rev. B, 1988, 37, 785–789 CrossRef CAS.
-
M. J. Frisch, G. W. Trucks, H. B. Schlegel, G. E. Scuseria, M. A. Robb, J. R. Cheeseman, J. A. Montgomery, Jr., T. Vreven, K. N. Kudin, J. C. Burant, J. M. Millam, S. S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G. A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J. E. Knox, H. P. Hratchian, J. B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R. E. Stratmann, O. Yazyev, A. J. Austin, R. Cammi, C. Pomelli, J. Ochterski, P. Y. Ayala, K. Morokuma, G. A. Voth, P. Salvador, J. J. Dannenberg, V. G. Zakrzewski, S. Dapprich, A. D. Daniels, M. C. Strain, O. Farkas, D. K. Malick, A. D. Rabuck, K. Raghavachari, J. B. Foresman, J. V. Ortiz, Q. Cui, A. G. Baboul, S. Clifford, J. Cioslowski, B. B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R. L. Martin, D. J. Fox, T. Keith, M. A. Al-Laham, C. Y. Peng, A. Nanayakkara, M. Challacombe, P. M. W. Gill, B. G. Johnson, W. Chen, M. W. Wong, C. Gonzalez and J. A. Pople, GAUSSIAN 03 (Revision C.02), Gaussian,
Inc., Wallingford, CT, 2004 Search PubMed.
-
S. Huzinaga and J. Andzelm, Gaussian Basis Sets for Molecular Calculations, Elsevier, Amsterdam, 1984, pp. 23–25 Search PubMed.
- T. Clark, J. Chandrasekhar, G. W. Spitznagel and P. v. R. Schleyer, J. Comput. Chem., 1983, 4, 294–301 CrossRef CAS.
-
G. W. Spitznagel, Diplomarbeit, Erlangen, 1982.
- These values were optimized for the atomic ground state anion, using ROHF with a flexible ECP basis set, by Ted Packwood at NDSU. Diffuse functions for these and other heteroatoms may be found at the website: http://phoenix.liu.edu/%E2%88%BCnmatsuna/gamess/refs/basis.refs.html.
- As there are no classical 2c–2e bonds in the clusters considered we use bonding “connection increments” rather than bond increments to refer to bonding connections in clusters.
- As the preference of para over meta positions is only minor and in order to make the overall structural increment scheme simpler, a CCm increment was not included for the 11-vertex nido-cluster.
- B. Stibr, Collect. Czech. Chem. Commun., 2002, 67, 843–868 CrossRef CAS and references therein.
- J. L. Little, J. T. Moran and L. J. Todd, J. Am. Chem. Soc., 1967, 89, 5495–5496 CrossRef CAS.
- L. J. Todd, J. L. Little and H. T. Silverstein, Inorg. Chem., 1969, 8, 1698–1703 CrossRef CAS.
- X. L. R. Fontaine, J. D. Kennedy, M. McGrath and T. R. Spalding, Magn. Reson. Chem., 1991, 29, 711–720 CAS.
- J. L. Little, Inorg. Chem., 1979, 18, 1598–1600 CrossRef CAS.
- J. L. Little, M. A. Whitesell, J. G. Kester, K. Folting and L. J. Todd, Inorg. Chem., 1990, 29, 804–808 CrossRef CAS.
Footnote |
† Electronic supplementary information (ESI) available: A plot of HetCo energy penalties for 12-vertex closo-clusters (Appendix I), a comparison of relative thermodynamic stabilities of various stanna-, phospha-, arsa- and stibaboranes and -borates (Appendices II–V), Cartesian coordinates from RB3LYP/SDD//LANL2DZ optimizations of various Het2B10H10 and HetCB10H11 structures (Appendix VI) and optimized geometries of three 11-vertex nido distannaboranes, i.e., 2,7-, 2,8- and 2,9-isomers of [Sn2B9H11]2− (which distorted upon optimization, Appendix VII). See DOI: 10.1039/b512700a |
|
This journal is © The Royal Society of Chemistry 2006 |
Click here to see how this site uses Cookies. View our privacy policy here.