Abstract
Using density functional theory and quantum Monte Carlo calculations, we show that B_{12}H_{n }and B_{12}F_{n }(n = 0 to 4) quasiplanar structures are energetically more favorable than the corresponding icosahedral clusters. Moreover, we show that the fully planar B_{12}F_{6 }cluster is more stable than the threedimensional counterpart. These results open up the possibility of designing larger boronbased nanostructures starting from quasiplanar or fully planar building blocks.
Background
The icosahedral B_{12}H_{12}^{2 }cluster is the most stable molecule among the number of polyhedral boranes synthesized so far [1]. A largescale and efficient synthesis of fully fluorinated boron hydrides, e.g., icosahedral B_{12}F_{12}^{2}, has been also reported [2]. On the other hand, the allboron C_{3v}B_{12 }cluster is quasiplanar, and it was reported to be one of the most stable allboron clusters. It was also established by extensive computations that the quasiplanar B_{12 }cluster is much lower in energy than the allboron icosahedral B_{12 }cluster. This was reported not only for the neutral clusters [3], but also for the charged ones [4]. It is then interesting to investigate what happens with the relative stability of the two (quasiplanar and threedimensional (3D)) allboron structures upon addition of hydrogen or fluorine atoms. This is the purpose of this study.
Quasiplanar and 3D boron clusters with the number of hydrogen atoms smaller than the number of boron atoms have been studied both theoretically [511] and experimentally [1214]. Ohishi et al. [12] reported the formation of B_{12}H_{n}^{+ }(n = 0 to 12) cationic clusters through ionmolecule reactions of the decaborane ions (B_{10}H_{n}^{+}, n = 6 to 14) with diborane molecules (B_{2}H_{6}) in an external quadrupole static attraction ion trap. The mass spectrum analysis revealed that among the B_{12}H_{n}^{+ }clusters with different hydrogen content n, the B_{12}H_{8}^{+ }molecule was the main product. In the same study, using first principle calculations with the Becke 3parameter LeeYangParr (B3LYP) hybrid functional and the 631G(d) basis set, the authors compared the relative energies of quasiplanar and 3D B_{12}H_{n}^{+ }clusters with n varying from 0 to 12. According to that study, twodimensional (2D) clusters with n = 0 to 5 are energetically preferred over the 3D structures, whereas 3D clusters are energetically favored for n ≥ 6. In a more recent combined experimental/theoretical study, Ohishi et al. [14] suggested that quasiplanar B_{12}H_{n}^{+ }with n = 0 to 3 clusters can be obtained by further removal of H atoms from the decaborane ions. This opens up the possibility of changing the structure of the B_{12}H_{n}^{+ }cluster by controlling the number of hydrogen atoms in the cluster.
To our knowledge, there are no previous studies on the structure and properties of quasiplanar B_{12}F_{m }clusters. However, the structures of two polyboron fluorides, B_{8}F_{12 }and B_{10}F_{12}, revealing unusual open structures were recently determined [15].
Methods
The initial search for the most stable structures of the boron hydrides B_{12}H_{n }and boron fluorides B_{12}F_{n }was done at the B3LYP/631G(d) level of theory using the FreeON code [16] with no symmetry restrictions. For clusters with an even number of hydrogen or fluorine atoms (even number of electrons), the computations were performed for the singlet multiplicity only, whereas doublet and quartet multiplicities were considered for clusters with an odd number of hydrogen or fluorine atoms (odd number of electrons). In the later case, structures with lower multiplicity were energetically more favorable. For the charged structures, a similar analysis has been done, and the ground state was found to have the lowest multiplicity. Next, the lowlying isomers of B_{12}H_{n }and B_{12}F_{n }have been reoptimized using the GAMESSUS code [17] at the B3LYP/6311++G(d, p) level of theory, and for the resulting structures, the vibrational analysis has been done to identify true local minima. It is important to mention that the 'bare' B_{12 }icosahedron undergoes distortions after structural optimization and that its symmetry is S_{2 }[3], not I_{h}. However, we will refer to that structure and its derivatives as icosahedral or 3D. The quasiplanar or fully planar clusters, for a change, will be often labeled as 2D structures.
The nucleus independent chemical shift (NICS) values and magnetic susceptibility tensors were calculated using the Gaussian 03 package [18] at the B3LYP/6311++G(d, p) level of theory. To obtain the NICS values, we have used the gaugeindependent atomic orbital method, and the magnetic susceptibility tensors were calculated using the continuous set of gauge transformations method.
The quantum Monte Carlo calculations have been done using the QWalk [19] package in two steps. The first step involved optimizing the trial manybody wave function by doing variational Monte Carlo calculations. The trial wave function was of the SlaterJastrow form. The Slater determinants were constructed using B3LYP orbitals, generated using the GAMESSUS code with the previously optimized geometries within the B3LYP/6311++G(d, p) level of theory. For the calculations, we have used Gaussian basis sets with effective core potentials [20]. In the second step, we have done fixednode diffusion Monte Carlo (DMC) calculations with the previously optimized trial wave functions. In the computations, we have used a time step of 0.005 a.u. The DMC error bars are about 0.1 eV.
Results and discussion
Planar vs icosahedral structures
The procedure for determining the most stable isomers of B_{12}H_{n }was very similar to that reported in [12], namely, we started with optimized icosahedral and quasiplanar B_{12 }clusters, and for a given n, we have calculated the total energies of all possible clusters that resulted from adding hydrogen atoms to the vertices of the distorted icosahedron or to the outer boron atoms of the quasiplanar structure; 2D clusters with an even n have been considered in our previous work [6], and here, we have extended the investigation to an odd n. The energetically most favorable 2D and 3D B_{12}H_{n }structures are shown in Figure 1. The minimumenergy cluster structures of B_{12}F_{n}, shown in Figure 2, have been found by replacing the hydrogen atoms of the lowlying B_{12}H_{n }isomers by fluorine atoms. Interestingly enough, the resulting structures are similar to those found for B_{12}H_{n}. One of the small differences is that the BF bonds are on average 13% longer than the BH bonds.
Figure 1. The most stable structures of 3D and 2D B_{12}H_{n }(n = 0 to 8) clusters. The symmetry of each cluster is given.
Figure 2. The most stable structures of 3D and 2D B_{12}F_{n }(n = 1 to 8) clusters. The symmetry of each cluster is given.
In Figure 3, we have plotted the total energy difference between quasiplanar (or fully planar) and icosahedral B_{12}X_{n }(X = H, F) clusters as a function of n, the number of H or F atoms in the cluster. As can be seen from the figure, the quasiplanar clusters with up to four hydrogen atoms are more stable than the corresponding icosahedral structures (a similar result has been recently reported [11] for B_{12}H_{n}^{0/ }clusters). The same is true for the fully planar B_{12}F_{6 }molecule, which is 0.63 eV lower in energy than the 3D cluster. The 2D and 3D B_{12}F_{5 }isomers are almost degenerated in energy. From Figure 3, we can also see that the energy difference, E_{2D } E_{3D}, increases monotonically with n, with the exception of the two 'minima' for B_{12}F_{4 }and B_{12}F_{6}. These two minima may suggest an additional stabilization of the 2D structures over the 3D counterparts due to the presence of aromatic stabilization energy.
Figure 3. Energy difference between quasiplanar and icosahedral B_{12}X_{n }(X = H, F) clusters as a function of number atoms.
Similar results to those presented in Figure 3 were reported for the icosahedral and quasiplanar B_{12}H_{n}^{+ }structures [12]. However, in their recent work, Ohishi et al. [14] have used the PBE0 functional instead of the B3LYP functional to determine the energies of the B_{12}H_{n}^{+ }clusters. The authors' choice was motivated by the fact that the B3LYP functional may overestimate the energy difference between 2D and 3D structures. To address this problem, we calculated the energy difference between the 2D and 3D structures of B_{12}, B_{12}H_{6}, and B_{12}F_{6 }using the very accurate DMC approach. The DMC E_{2D } E_{3D }values are 5.13, 0.79, and 0.47 eV for B_{12}, B_{12}H_{6}, and B_{12}F_{6}, respectively, whereas the corresponding B3LYP values are 5.34, 0.41, and 0.63 eV, respectively (see Figure 3). This means that the DMC values are shifted up by a value not larger than about 0.4 eV with respect to the B3LYP values. This, however, does not affect the conclusions that are drawn from Figure 3, since even if we shift up the curves by 0.4 eV, the quasiplanar B_{12}H_{n }and B_{12}F_{n }(n = 0 to 4), and the fully planar B_{12}F_{6 }clusters still remain energetically favorable.
Fully planar clusters: B_{12}H_{6 }vs B_{12}F_{6}
As calculated here and also reported in [11], the fully planar B_{12}H_{6 }cluster corresponds to a local minimum of energy, whereas the D_{3h}B_{12}F_{6 }structure wins the competition with other 2D and 3D isomers and corresponds to a global minimum of energy. Many properties of the B_{12}H_{6 }cluster have been previously described in [6], but for consistency purposes, we have repeated some of those calculations at the B3LYP/6311++G(d, p) level of theory. The highest occupied molecular orbitallowest unoccupied molecular orbital (HOMOLUMO) gaps of the planar B_{12}H_{6 }and B_{12}F_{6 }structures are 3.54 and 4.39 eV, respectively, whereas the HOMOLUMO gaps of the corresponding 3D clusters are the same and equal to 2.73 eV. The BH and BF interatomic distances are 1.179 and 1.326 Å in B_{12}H_{6 }and B_{12}F_{6}, respectively. For comparison, the computed BH and BF bond lengths in borane (BH_{3}) and boron trifluoride (BF_{3}) are 1.190 and 1.318 Å, respectively.
While both 2D structures, B_{12}H_{6 }and B_{12}F_{6}, have similar shape and size, they exhibit quite different magnetic properties that are directly related to aromaticity. First, we have computed the anisotropy of magnetic susceptibility. The values for B_{12}H_{6 }and B_{12}F_{6 }are 208.1 and 125.8 cgs ppm, respectively. The isotropic values of the magnetic susceptibility are 91.9 and 118.2 cgs ppm for B_{12}H_{6 }and B_{12}F_{6}, respectively. These results suggest that the induced ring current is stronger for B_{12}H_{6 }than for B_{12}F_{6}. Similarly, as reported in [6] for B_{12}H_{6}, the central part of the B_{12}F_{6 }molecule has a paratropic current flowing inside the inner B_{3 }triangle. The antiaromaticity of the inner triangle is, however, smaller for B_{12}F_{6 }than for B_{12}H_{6 }since the NICS(0) values are 3.9 and 13.3 ppm, respectively. A global aromatic current is dominant above the B_{12}F_{6 }(B_{12}H_{6}) molecule since the NICS values are negative, NICS(1) = 5.5 ppm (3.6 ppm) and NICS(2) = 4.8 ppm (5.0 ppm), 1 and 2 Å above and below the center of the cluster.
To examine the influence of charge on the structure and stability of the fully planar clusters, we have also studied charged 2D and 3D B_{12}X_{6 }(X = H, F) structures. The lowest energy 2D structures identified for B_{12}H_{6 }and B_{12}F_{6 }were used as initial structures for the structural optimization at a given charge state. For the 3D structures, we have made a search over all possible configurations of the hydrogen or fluorine atoms. For the lowest energy structures with an even number of electrons, a singlet multiplicity has been assumed, whereas doublet and quartet multiplicities were considered for clusters with an odd number of electrons. In the later case, clusters with lower multiplicity were energetically more favorable. The structures of the 2D and 3D charged clusters are shown in Figure 4. It has been previously reported that the fully planar D_{3h}B_{12}H_{6 }cluster undergoes structural distortions if charged with one electron, although the quasiplanarity is preserved [6]. In general, all the charged 2D B_{12}X_{6 }(X = H, F) clusters are quasiplanar rather than fully planar, as can be seen in Figure 4. In Figure 5, we have plotted the energy difference between 2D and 3D [B_{12}X_{6}]^{q }(X = H, F) structures as a function of the cluster charge state q. We have found that the addition of one or two electrons to fully planar B_{12}H_{6 }and B_{12}F_{6 }clusters (or the removal of one electron from them) makes those structures even less energetically favorable with respect to the corresponding 3D isomers. This is, however, less pronounced for B_{12}H_{6 }than for B_{12}F_{6 }as shown in Figure 5. Finally, it should be noted that the quasiplanar B_{12}F_{6}^{2+ }cluster is much more stable than its 3D isomer. Finally, all structures and energies are provided in Additional file 1.
Figure 4. The structures of charged 3D and 2D clusters of B_{12}H_{6 }and B_{12}F_{6}. The symmetry of each cluster is provided.
Figure 5. Energy difference between quasiplanar and icosahedral [B_{12}X_{6}]^{q }(X = H, F) clusters as a function of cluster charge.
Additional file 1. Total electronic energies of the boron structures. Electronic supplementary material Figure S1 shows a fully planar boronbased nanostructure, B_{504}H_{36}, which was obtained starting from planar B_{12}H_{6 }building blocks. Table S1 recollects total energies and Cartesian coordinates of the optimized structures shown in Figures 1 and 2.
Format: DOC Size: 3.9MB Download file
This file can be viewed with: Microsoft Word Viewer
Conclusions
Our density functional theory and quantum Monte Carlo results show that the B_{12}H_{n }and B_{12}F_{n }(n = 0 to 4) quasiplanar structures are energetically more favorable than the corresponding icosahedral clusters and that the fully planar B_{12}F_{6 }cluster is more stable than the 3D counterpart. We have also shown that negative or positive charge further stabilizes the 3D over the 2D B_{12}X_{6 }(X = H, F) clusters (except for B_{12}X_{6}^{2+}, where the opposite is observed). Our findings are potentially useful for designing larger boronbased nanostructures starting from quasiplanar or fully planar building blocks.
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
NGS conceived the study, did the calculations, and drafted the manuscript. CJT conceived the study and revised the manuscript. All authors read and approved the final manuscript.
Acknowledgements
The authors would like to acknowledge the support given by the Robert A. Welch Foundation (grant J1675), the NSF CREST CRCN center (grant HRD1137732), and the Texas Southern University High Performance Computing Center (http://hpcc.tsu.edu/; grant PHY1126251).
References

Pitochelli AR, Hawthorne FM: The isolation of the icosahedral B_{12}H_{12}^{2 }Ion.

Peryshkov DV, Popov AA, Strauss SH: Direct perfluorination of K_{2}B_{12}H_{12 }in acetonitrile occurs at the gas bubble  solution interface and is inhibited by HF. Experimental and DFT study of inhibition by protic acids and soft, polarizable anions.
J Am Chem Soc 2009, 131:1839318403. PubMed Abstract  Publisher Full Text

Atiş M, Özdoğan C, Güvenç ZB: Structure and energetic of B_{n }(n = 212) clusters: electronic structure calculations.
Int J Quantum Chem 2007, 107:729744. Publisher Full Text

Zhai HJ, Kiran B, Li J, Wang LS: Hydrocarbon analogues of boron clusters  planarity aromaticity and antiaromaticity.
Nat Mater 2003, 2:827833. PubMed Abstract  Publisher Full Text

Alexandrova A, Koyle E, Boldyrev A: Theoretical study of hydrogenation of the doubly aromatic B_{7}^{ }cluster.
J Mol Model 2006, 12:569576. PubMed Abstract  Publisher Full Text

Gonzalez Szwacki N, Weber V, Tymczak CJ: Aromatic borozene.
Nanoscale Res Lett 2009, 4:10851089. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Sahu S, Shukla A: Probing aromaticity of borozene through optical and dielectric response: a theoretical study.
Nanoscale Res Lett 2010, 5:714719. PubMed Abstract  Publisher Full Text  PubMed Central Full Text

Forte G, La Magna A, Deretzis I, Pucci R: Ab initio prediction of boron compounds arising from borozene: structural and electronic properties.
Nanoscale Res Lett 2010, 5:158163. Publisher Full Text

Chen Q, Li SD: πAromatic B_{16}H_{6}: a neutral boron hydride analogue of naphthalene.
J Clust Sci 2011, 22:513523. Publisher Full Text

Böyükata M, Güvenç ZB: DFT study of Al doped cage B_{12}H_{n }clusters.
Int J Hydrogen Energ 2011, 36:83928402. Publisher Full Text

Bai H, Li SD: Hydrogenation of B_{12}^{0/}: a planartoicosahedral structural transition in B_{12}H_{n}^{0/ }(n = 16) boron hydride clusters.
J Clust Sci 2011, 22:525535. Publisher Full Text

Ohishi Y, Kimura K, Yamaguchi M, Uchida N, Kanayama T: Formation of hydrogenated boron clusters in an external quadrupole static attraction ion trap.
J Chem Phys 2008, 128:124304124307. PubMed Abstract  Publisher Full Text

Ohishi Y, Kimura K, Yamaguchi M, Uchida N, Kanayama T: Energy barrier of structure transition from icosahedral B_{12}H_{6}^{+ }to planar B_{12}H_{5}^{+ }and B_{12}H_{4}^{+ }clusters.

Ohishi Y, Kimura K, Yamaguchi M, Uchida N, Kanayama T: Synthesis and formation mechanism of hydrogenated boron clusters B_{12}H_{n }with controlled hydrogen content.
J Chem Phys 2010, 133:074305074307. PubMed Abstract  Publisher Full Text

Pardoe JAJ, Norman NC, Timms PL, Parsons S, Mackie I, Pulham CR, Rankin DWH: The surprising structures of B_{8}F_{12 }and B_{10}F_{12}.
Angew Chem Int Edit 2003, 42:571573. Publisher Full Text

Bock N, Challacombe M, Gan CK, Henkelman G, Nemeth K, Niklasson AMN, Odell A, Schwegler E, Tymczak CJ, Weber V: FreeON. [http://freeon.org/] webcite

Schmidt MW, Baldridge KK, Boatz JA, Elbert ST, Gordon MS, Jensen JH, Koseki S, Matsunaga N, Nguyen KA, Su S, Windus TL, Dupuis M, Montgomery JA: GAMESSUS, version 1 OCT 2010 (R1). [http://www.msg.chem.iastate.edu/gamess/] webcite

Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA Jr, Vreven T, Kudin KN, Burant JC, et al.: Gaussian 03, revision E.01. [http://www.gaussian.com/] webcite

Wagner LK, Bajdich M, Mitas L: QWalk, version 0.95.0. [http://www.qwalk.org/] webcite

CRENBL ECP [http://bse.pnl.gov/] webcite