Aromatic borozene.

Based on our comprehensive theoretical investigation and known experimental results for small boron clusters, we predict the existence of a novel aromatic inorganic molecule, B12H6. This molecule, which we refer to as borozene, has remarkably similar properties to the well-known benzene. Borozene is planar, possesses a large first excitation energy, D3hsymmetry, and more importantly is aromatic. Furthermore, the calculated anisotropy of the magnetic susceptibility of borozene is three times larger in absolute value than for benzene. Finally, we show that borozene molecules may be fused together to give larger aromatic compounds with even larger anisotropic susceptibilities.

Introduction

Why certain molecules are more stable than others is not always easy to understand. Nature’s diversity does not always permit a simple answer for the structure of all compounds. However, a very useful concept in structural stability is aromaticity [1,2], which was first developed to account for the properties of organic compounds involving ring structures such as benzene (C6H6) and more recently extended to inorganic systems [3]. However, the question arises whether aromatic hydrocarbons are the only structures where an “aromatic ring” acts as a building block and plays a key role in their stability. In this study, we predict the existence of a novel inorganic molecule, B12H6, that has remarkably similar properties to benzene. This molecule, which we call borozene for brevity, is planar, possesses a large first excitation energy, exhibits a highly aromatic character, and similar to benzene is a building block of possibly much larger aromatic compounds.

Small all-boron clusters, B (n < 20) have been widely investigated both experimentally and theoretically [4-9]. All these studies indicate that small boron clusters assume in most cases quasi-planar structures and in some special cases even perfectly planar structures. In contrast, neutral and anionic boron hydrides, BH, are all known to have three-dimensional deltahedral structures [10]. There is yet little known about the structure of small boron hydrides where the number of hydrogen atoms is smaller than the number of boron atoms (see ref. [11] and references therein). One such example is the recently studied σ-aromatic and π-antiaromatic cluster, which is fully planar [11].

It was experimentally established that one of the most stable all-boron clusters is made up of twelve boron atoms, is quasi-planar in shape, and possesses a large first excitation energy of 2.0 eV [4]. The B12 structure consists of 13 B3 triangles with 12 outer triangles surrounding a central one; the atoms forming the central triangle are situated above a nine-member boron ring making B12 a convex structure of C3v symmetry [4-6,9]. Our calculations revealed that the B12 cluster has three outer boron pairs that are 5% shorter than the average B–B bond lengths between the rest of the boron atoms. This suggests the presence of strong covalent bonds between those atoms. However, it is possible to increase the B–B bond lengths of the three outer pairs by 17% by attaching hydrogen atoms to the outer boron atoms (see Fig. 1a). The unexpected consequence of this is that the molecule becomes perfectly planar. This finding motivated us to investigate the interaction between B12 cluster and up to four hydrogen molecules.

Figure 1

aPlot of the structures and total electronic densities of B12H6. Note that the density of electrons is weaker at the center of the “boron ring.”bQuasi-planarCsstructure of B12H8, which is the energetically preferred configuration for B12with 4 H2molecules, attached to it.cA contour plot of NICS(x,y) for B12H6in plane (left) and at 1 Å above the planar molecule (right)

Computational Details

The structure and electronic properties of all clusters were obtained at the X3LYP/6-311G(d,p) level of theory using tight convergence criteria as implemented in FreeON (formerly known as MondoSCF), a suite of programs using Gaussian basis sets and all-electron Hartree–Fock, density functional theory or hybrid approach for self-consistent electronic structure calculations [12-17]. The initial search for the most stable structures of the boron hydride B12H have been done at the X3LYP/6-31G(d,p) level of theory starting from the energy-minimum structure of B12H and the low-lying isomers in each case have been reoptimized using the 6-311G(d,p) basis set. The obtained local energy-minimum structures are well separated in energy from its higher isomers by at least 22 kcal/mol in the case of B12H, where n ≤ 6, and 17 kcal/mol in the case of the B12H8 cluster. The B22H8 and B60H12 clusters, which result from the fusion of two and six B12H6 molecules, respectively, were fully optimized using symmetry-unrestricted calculations.

To ensure that the structures for the neutral and negatively charged B12H6molecules correspond to a local minima of energy, the nature of the stationary points where checked by vibrational frequency calculations. In addition, several tests have been done for neutral B12H6optimizing this molecule starting from structures with randomly displaced atoms. For all tested cases, the substructure of boron atoms became planar after optimization to within 0.012 and 0.014 Å at the X3LYP/6-311++G(d,p) and UHF-MP2/6-311++G(d,p) levels of theory, respectively. Precise calculations at the X3LYP/6-311G(d,p) level of theory and using very tight convergence criteria revel, however, that the fully planar B12H6molecule is the most energetically favorable structure. Similar tests were done for B22H8and B60H12, and after optimization at the X3LYP/6-311G(d,p) level of theory both substructures of boron atoms became planar to within 0.016 Å.

The first singlet excitation energy, nuclear magnetic resonance (NMR) shielding tensors, and magnetic susceptibility tensors were calculated using the Gaussian03 package [18]. This package was also used to further optimize the B12H6 molecule using the 6-311++G(d,p) basis set. To obtain the nucleus independent chemical shift (NICS) values (from the NMR shielding tensors), we have used the GIAO (gauge-independent atomic orbital) method and the magnetic susceptibility tensors were calculated using the CSGT (continuous set of gauge transformations) method. All computations have been performed at the X3LYP/6-311++G(d,p) level of theory except for the B60H12 cluster, for which we have used the RHF/6-31G(d,p) level of theory. The anisotropy of magnetic susceptibility is defined as the difference between out-of-plane and the average in-plane components of the susceptibility tensor.

The MOs of B12H6 and C6H6 were calculated at the RHF/6-311++G(d,p) level of theory using the GAMESS-US package [19]. The same package was used to calculate the π–π interaction between molecules in borozene and benzene dimmers at the RHF-MP2/6-311G(d,p) and RHF-MP2/6-311++G(d,p) levels of theory, respectively. The counterpoise correction was applied to account for the basis set superposition error.

Results and Discussion

The search for the stable structures of B12H, with n ≤ 6 is simplified by the fact that the most energetically favorable configurations, as far as we can determine, are those where the hydrogen atoms are directly attached to the outer boron atoms of the molecule. We have established that the most likely stable configuration for B12H2 is when the hydrogen atoms are attached to one of the outer short-bonded boron pairs of the B12 cluster. The energetically preferred configuration for B12H4 is when the hydrogen atoms are attached to one of the two remaining outer short-bonded boron pairs in the B12H2. Finally, the B12H6 cluster has all short B–B pairs, from B12, with hydrogen atoms attached to them. Only B12H6 is a fully planar molecule, whereas B12H2 and B12H4 are quasi-planar with Cs symmetry. In Fig. 1a, we have shown the structure of B12H6. The hydrogenation energy, defined as ΔE = E(B12H) − E(B12H) − E(H2) where E is the total energy, is −44 kcal/mol for n = 2, −45 kcal/mol for n = 4, and −51 kcal/mol for n = 6. We have found, however, that if a fourth H2 molecule is attached to B12H6 the hydrogenation energy increases to −2 kcal/mol (i.e., the H2 molecule is weakly bound to B12H6). It is also important to mention that our ΔE values are about two times larger than the predicted energy of hydrogenation of the cluster [11], which is an additional indication of unusual stability of the B12 structure. The B12H8 molecule is shown in Fig. 1b and can be described as a distorted B12H6 cluster with two extra (one terminal and one bridging) hydrogen atoms attached to it. The B–H bond lengths are 1.36 and 1.21 Å for the bridging and terminal hydrogen atoms, respectively, whereas the remaining B–H distances in B12H8 and in all other described earlier B12H (n ≤ 6) clusters are the same and equal to 1.18 Å. The last value is very close to the calculated bond lengths B–H = 1.19 Å in borane, BH3.

Although no single measure of aromaticity is without limitations, the anisotropy of the magnetic susceptibility (AMS) is an important indicator of diatropicity [20]. The B12H6 molecule has very important properties: it is planar with D3h symmetry; it possesses a large first excitation energy of 2.6 eV and a large AMS of −208.2 cgs-ppm. Also, the B12H6 molecule can be a building block of larger planar molecules with similar structural and physical characteristics. In Fig. 2a, b are shown what we call boron analogs of naphthalene (B22H8) and coronene (B60H12), which are fusions of two and six B12H6, clusters, respectively (It should be pointed out that the idea of boron analogs of both benzene and naphthalene are not new, since the all-boron B12 molecule was previously indicated [4] to be an analog of benzene and very recently also an all-boron analog of naphthalene was proposed [7].). Also of interest, the HOMO–LUMO (HOMO, highest occupied molecular orbital; LUMO, lowest unoccupied molecular orbital) gap decreases with cluster size. The gap values are 3.6, 2.4, 1.3 eV for B12H6, B22H8, B60H12, respectively. Furthermore, the absolute AMS value increases with cluster size; the AMS values are −294.3, −454.7 for B22H8, B60H12, respectively. This and other values for the studied boron hydrides are summarized in Table 1. In this table, we have also included the calculated values for B12 with enforced planarity (D3h) and the fully relaxed structure (C3v). We have in addition included the values for three hydrocarbons for comparison. It should be noted that the absolute value of AMS for B12H6 is three times larger than our value for benzene (−67.5 cgs-ppm) and 7% larger than the value for the C3v B12 cluster (−192.9 cgs-ppm).

Figure 2

Plot of the structures and total electronic densities ofaB22H8, andbB60H12. In both molecules, all B–H distances are the same and equal to 1.18 Å

Table 1

Molecular symmetry, HOMO–LUMO energy gaps, and the isotropic and anisotropic values of magnetic susceptibility for the studied planar boranes and hydrocarbons

Structure	Symmetry	HOMO–LUMO gap (eV)	Magnetic susceptibility (cgs-ppm)
			Isotropy	Anisotropy
B12	C3v	3.73	−105.4	−192.9
B12	D3h	3.58	−107.6	−213.6
B12H6	D3h	3.67	−92.0	−208.2
B22H8	D2h	2.38	−147.9	−294.3
B60H12	D6h	1.30	−286.5	−454.7
C6H6	D6h	6.86	−53.0	−67.5
C10H8	D2h	4.93	−90.4	−128.5
C24H12	D6h	4.13	−251.1	−474.5

For comparison, we have also included our results for B12with enforced planarity (D3h) and the fully relaxed structure (C3v)

To gain information about the individual contributions of the B3triangles to the overall aromaticity of the B12H6molecule, we have studied its two-dimensional NICS map. In Fig. 1c, we have shown the contour plot of NICS(x, y) in plane (left) and at 1 Å above the B12H6molecule (right). It is clearly seen from the left part of the figure that the NICS values are negative inside the twelve outer B3triangles of the molecule, suggesting a flow of a global diatropic current around the central triangle. The central part of the molecule has a paratropic current flowing inside the inner B3triangle which is not overwhelmed by a diatropic current due to an electron charge transfer from the center of the structure towards the outer boron triangles after hydrogenation. Also, this antiaromatic region is spatially localized since the NICS values are negative at 1 Å (see the right part of Fig. 1c) and 2 Å above and below the center of the B12H6molecule (NICS(1) = −3.9 ppm; NICS(2) = −5.2 ppm).

The benzene dimer is the simplest prototype of the aromatic π–π interactions which is an important weak interaction present in aromatic supramolecular systems [12]. Since the number of π electrons in borozene and benzene is the same and the molecular orbital (MO) picture for these electors is very similar (see Fig. 3), we may expect that the strength of the aromatic–aromatic interaction in a borozene dimer is comparable to that of the benzene dimer. To investigate this, we have considered the simplest case where the molecules in the dimer have the parallel “sandwich” configuration. In Fig. 4, we plotted the association energy versus the distance between the molecules in the B12H6 and C6H6 dimmers. From this figure, we can see that the association energy for the borozene dimmer, in its equilibrium position, is about five times larger than the corresponding energy for the benzene dimer. This result suggests a stronger polarization contribution from borozene’s π-MOs, which we theorize is a consequence of more delocalized π-electrons in the borozene dimer with respect to the benzene dimer.

Figure 3

Comparison of the π molecular orbitals of benzene with the corresponding π molecular orbitals of B12H6

Figure 4

Potential energy curves for B12H6and benzene dimers versus the center-to-center distance between the monomers. The association energies are −1.99 and −9.81 kcal/mol and the equilibrium distances are 3.8 and 3.9 Å for C6H6and B12H6, respectively

Finally, we have also calculated the electron affinity (EA) and the ionization energy (IE) for B12H6. The values are summarized in Table 2calculated at two levels of theory. The agreement between the values obtained using the hybrid functional and MP2 is quite good. The EA values are positive what means that the molecule is stable with respect to B12H6and a free electron. However, an addition of an electron to costs energy since the EA is negative and large (see Table 2). A removal of one electron from B12H6costs a lot of energy since the vertical IE values are very large: 8.52 and 8.72 eV for X3LYP and MP2, respectively. The structure of is shown in Fig. 5. This molecule hasC2symmetry and the axis of rotation is shown in the figure. TheD3hsymmetry of B12H6is reduced here toC2since the charged molecule is “twisted” along the shown axis.

Table 2

Vertical electron affinities and ionization energies for B12H6 and calculated using the 6-311++G(d,p) basis set

	Electron affinity (eV)		Ionization energy (eV)
	B12H6		B12H6
X3LYP	1.69 (1.80)	−2.57	8.52	1.97 (1.80)
UHF-MP2	1.52 (1.71)	−2.18	8.72	1.89 (1.71)

The numbers in brackets correspond to the adiabatic values for EA and IE. All values have been obtained from differences in total energies

Figure 5

Front (left) and side (right) views of the negatively charged B12H6molecule. TheC2rotation axis is shown

Although the specific route for the synthesis of the B12H6(or ) structure is not yet known, it is clear from our investigation that to some extent the chemistries of B12H6and benzene may be very similar, suggesting that similar methods could be employed to synthesize this and related compounds. Given the technological importance of benzene and its derivatives, we believe that the B12H6molecule will have a significant technological impact and deserves further extensive study.