Abstract
Structure and electronic properties of two unusual boron clusters obtained by fusion of borozene rings have been studied by means of first principles calculations based on the generalizedgradient approximation of the density functional theory. Moreover, a semiempirical tightbinding model has been appropriately calibrated for transport calculations on these clusters. Results show that the pure boron clusters are topologically planar and characterized by (3c–2e) bonds, which can explain, together with the aromaticity (estimated by means of NICS), the remarkable cohesive energy values obtained. Such feature makes these systems competitive with the most stable boron clusters to date. The energy gap values indicate that these clusters possess a semiconducting character, while when the larger system is considered, zerovalues of the density of states are found exclusively within the HOMO–LUMO gap. Electron transport calculations within the Landauer formalism confirm these indications, showing semiconductorlike low bias differential conductance for these structures. Differences and similarities with carbon clusters are highlighted in the discussion.
Keywords:
Boron clusters; Borozene; DFT; NICS; TransportIntroduction
Boron is the first element in group IIIA of the periodic table, presents the external electronic configuration s^{2}p^{1} and possesses a variety of compounds second only to carbon. In several boron compounds [1,2] the existence of multicenter bonds has been discovered, which arise from the electron deficiency of this element. Moreover, boron has a special place among the elements of the periodic table because of the wide variety of crystalline structure forms, i.e., polymorphism, which include nanotubes [3], nanoribbons [4], and nanoclusters [5,6].
The interest in boronbased nanostructures has recently increased due to new studies of both the synthesis of single walled boron nanotubes (SWBNTs) and the prediction of ballistic conduction in SWBNTs [3,6]; these findings, together with the properties that all boron nanotubes (a) are predicted to be metallic [7] and (b) are superconducting at low temperatures [8,9], promoted their prospective applications in fabrication of novel electronic devices. The most stable boron structure is the αrhombohedral bulk where boron icosahedra are centered on the edges of a rhombohedral unit cell [10].
Unlike the bulk boron compounds, boron clusters B_{n} (n < 20) are quasiplanar, or even planar, with a symmetrical bond distribution, aromatic [1113] and their existence is confirmed by the experiment [14]. From the Aufbau principle postulated by Boustani [5] it follows that these quasiplanar isomers are more stable than their icosahedral counterparts. Recently Szwacki et al. [15] have predicted the existence of a planar and aromatic boron compound, named borozene, which has strong similitudes with benzene.
Motivated by these findings we present a work regarding a first principles study, within the generalized gradient approximation (GGA), in terms of structural and electronic properties of two boron compounds, B_{60}H_{12} and B_{228}H_{24}, in which the molecule of borozene can be considered as the building block, as the benzene ring represents the embryo of compounds such Coronene, Coronene 19 etc. In general, we will refer to these compounds as boron clusters whose external dangling bonds are saturated by hydrogen atoms; they are obtained by fusing together the outer boron pairs of borozene molecules bonded to a hydrogen atom, see Fig. 1. We have organized the rest of the paper as follows: the computational methods adopted are presented in “Computational Method” section. In “Results” section we present and discuss our results in comparison also with carbon compounds, and, finally, we give a summary in “Conclusions” section.
Figure 1. Clusters B60H12 (left) and B228H24 (right) obtained after geometry optimization. Highlighted in red are the molecule of borozene (left) and the cluster B60H12 (right). In yellow are highlighted the boron atoms in contact with the gold leads in the transport calculations
Computational Method
The molecule B_{60}H_{12} here considered was built by fusing six borozene rings, for this reason it can be considered as the boron counterpart of coronene, whereas the structure of B_{228}H_{24} was obtained by surrounding B_{60}H_{12} with one series of borozene rings, therefore this cluster is constituted by a total of 24 borozene rings. A first optimization energy procedure was performed in the framework of the molecular mechanics approximation applying the CVFF Force Field [16,17] which is enclosed in the Materials Studio package [18]. The geometries obtained were fully optimized at a B3LYP/STO3G [1923] and B3LYP/6311G [24,25] level by using the quadratically convergent Self Consistent Field procedure [26].
In detail, due to the large size, the optimizations of B_{228}H_{24} were carried out by means of the minimal basis set STO3G whereas the more extended Pople basis set 6311G was used in the optimizations of B_{60}H_{12}. In order to estimate the degree of aromaticity, the calculation of Nuclear Independent Chemical Shifts [27] on the plane of the aromatic system (NICS0) was computed using the Gaussian 03 package [28], applying the GIAO method [29,30]. To obtain the contour plot of NICS, ghost atoms were placed on the plane of the molecule with a step size of about 1 Å.
Finally, electronic transport has been evaluated in the framework of the Nonequilibrium Green’s Function theory using a Landauer expression for the calculation of the current–voltage (I–V) characteristics [31]. Consistently to the electron structure findings we assume that only p_{z} orbitals contribute to the low bias transport along the molecules and that the Fermi energy is at the center of HOMO–LUMO gap. The molecular device configuration considered consists of two vertical gold leads in contact with the horizontal molecules forming ideal bonds with the boron atoms indicated in Fig. 1.
Results
Structural Properties
The analysis of the smaller cluster, henceforth named B6, was performed by using both basis sets mentioned above in order to make a consistent comparison with B_{228}H_{24}, henceforth named B24, analyzed only with the minimal basis set. We point out that the results obtained in the two cases are qualitatively equivalent, for this reason, unless specified, the data shown below are referred to the minimal basis set. As far as the boron–boron bond length is concerned, a shorter value has been found in the STO3G optimized structure, in particular, taking into account cluster B6, the average value of this parameter is calculated to be, respectively, 1.649 and 1.638 Å for 6311G and STO3G basis set, while a bond length average of 1.629 Å is obtained for the cluster B24.
The decrease of the bond length average by increasing the size of the cluster is also seen in the Coronene 19, i.e., a molecule of Coronene surrounded by a series of benzene rings, where, by using the same level of calculation, a 0.021 Å decrease of the same parameter is found with respect to the Coronene. It is also interesting to note that the reduction of the bond length takes place in particular in the inner bonds which tend to have the same value. The cohesive energies, evaluated in the minimal basis set for both clusters, were of 6.437 eV for B6 and 6.449 eV for B24. These values were calculated from
where
In the expressions above n is the number of boron atoms and the value of the B–H bond energy is calculated in the first approximation as 1/3 of the binding energy value of BH_{3}. Structural parameters evaluated are competitive, in terms of stability, with the more stable flat twodimensional structures considered up to date [32,33]. It is well known that Boron has a variety of compounds containing multicenter bonds, in particular the threecenter, twoelectron (3c, 2e) bond is present in molecules such diborane [34], boron clusters [32] and boron sheets [35].
Previous works have shown that (3c, 2e) bonds preclude the formation of boron rings in boron clusters [6,36], whereas, on the other hand, more recently the threecenter bonding has been proposed to explain the stability of boron fullerenes [34,37]. This peculiar feature is also seen for B6 and B24, while it is not present in Coronene and its larger clusters such Coronene 19, Coronene 37 and Coronene 61. Hence it is logical to assume that, as for boron fullerenes, it plays a pivotal role in maintaining a twodimensional stable planar structure.
The presence of (3c, 2e) bonds can be evaluated by means of the Mayer Bond Order indices, calculated from the canonical MOs in the canonical AO basis [34,3840], which, for closedshell species with 3 center bonds involving the atoms A, B and C, can be expressed as follows:
where P is the total density matrix and S is the overlap matrix. The bond order indices of threecenter bonds are positive with a theoretical maximum of ≈0.296. In Fig. 2, we report the more relevant values of I for both clusters, in general one can affirm that each boron is involved at least in two different threecenter bonds, i.e., each boron is directly linked at least with four boron atoms.
Figure 2. A section of the cluster B228H24; labels from 1 to 17 are also referred to B60H12 in black (red) are reported the principal threecenter–twoelectron Mayer bond order indices of cluster B6 (B24)
Electronic Properties
It has been suggested that the anomalous stability of the boron planar clusters depends on the aromaticity which arises from the delocalization of πelectrons and involves unoccupied 2p_{z} atomic orbitals [11,12,14]. As it will be discussed below, B6 and B24 show these features; before analyzing in detail we underline that a more extended electronic delocalization gives rise to a smaller GAP in carbon clusters [41], described as the HOMO–LUMO energy difference. In accordance with these calculations, the GAP values obtained for B6 and B24 are 1.33 and 1.17 eV, respectively; their density of states (DOS) as a function of energy (eV) is shown in Fig. 3. At energy E the density of states is written as
where the summation index i goes over all energy levels and δ is the Delta function.
Figure 3. Above Density of states and HOMO (in the insets) for cluster B6 (black line) and B24 (red line).Below Density of states and HOMO (in the inset) for Coronene
From Fig. 3 we note that: (a) both curves show a similar profile, cluster B24 has a larger density of states, while differently from cluster B6, it shows a zero value of the DOS only within the HOMO–LUMO gap; (b), the composition of the molecular orbitals, calculated by means of Mulliken Population Analysis, reported in Table 1, clearly highlights how this contributes to the HOMO of both clusters, shown in the insets of Fig. 3 and to their nearer molecular orbitals, are mainly due to the p_{z} atomic orbitals, confirming the stabilizing effect of πdelocalization. This result is in agreement with the one evidenced in Coronene, whose HOMO and DOS are reported, in Fig. 3; (c), as already observed for the set of Carbon clusters previously studied [41], the peaks near the HOMO energy can be joined by an almost straight line, reproducing the linear dependence shown by the infinite system near the Fermi level.
Table 1. Energies and percentual contributes of p_{z} orbitals to the composition of HOMO, LUMO and their nearest MOs
Now we turn to aromaticity which, as already mentioned, is considered as the basis of stability for boron planar clusters. Szwacki et al. [15] have discussed about the regions of aromaticity of borozene and since this molecule can be indicated as the embryo of B6 and B24, we find it necessary to investigate this aspect. Figure 4 shows the plot of the nucleus independent chemical shift (NICS), which represents the magnetic criterion to evaluate the ring current for cluster B24.
Figure 4. In plane contour plot of NICS (X,Y) for B24. The step size of the ghost atoms is about 1 Å
Negative value of NICS arise when diatropic ring current dominates, meaning that the system considered is aromatic, on the other hand a paratropic current gives rise to a positive value of NICS, therefore the corresponding system is antiaromatic. From Fig. 4 it is evident that inner bonds give rise to a paratropic current inside the round areas which can be considered as the expansion of the inner triangle antiaromatic area found in borozene [15], whereas a flow of diatropic current is homogeneously present in the rest of the cluster. In Fig. 5 the low bias differential conductance of the two clusters ideally contacted with two gold leads is shown.
Figure 5. Red line (black line): low bias differential conductance of the B6 (B24) molecule ideally contacted with gold leads in the two terminal configurations
The symmetry of the plots reflects the assumed symmetry in the device configuration. A semiconductorlike behavior is evidenced in both structures. However, the zero bias differential conductance is one order of magnitude higher for the B24 cluster, this is due not only to the lower HOMO–LUMO gap but also to the larger value of the DOS. Delocalized (along the cluster) p_{z} molecular orbitals allow an efficient charge transport through the cluster for larger bias (of the order of the gap) and a diodelike characteristic can be observed.
Conclusions
In this work first principles and semiempirical calculations were carried out to investigate both structural and electronic properties of two clusters obtained by condensation of 6 and 24 borozene molecules, considered as the analogs of Coronene and Coronene 19. Calculations predict a planar geometry for the pure clusters. Both (3c–2e) bonds and wide regions of aromaticity contribute to this stabilization, with cohesive energy values that are comparable with the most stable boron clusters considered to date.
Due to the high connectivity among boron atoms, we hypothesize that planar geometry could be compromised when impurities are introduced. Density of states spectra evidence a small gap, which decreases by increasing the cluster size, suggesting, at variance with carbon clusters, a semiconducting character in small sized clusters; furthermore the population analysis shows that the main contribution to the molecular orbitals near the GAP is due to πbonds which derive from p_{z} orbitals. Calculated low bias differential conductance for these structures confirms this semiconductorlike character.
Acknowledgment
G. Forte wishes to thank the Consorzio Interuniversitario Cineca for the computational support.
References

Nat. Phys.. 2007, 6:44.
COI number [1:CAS:528:DC%2BD2sXot1On]

Science. 2007, 315:356.
COI number [1:CAS:528:DC%2BD2sXmt1KgsQ%3D%3D]; Bibcode number [2007Sci...315..356L]
PubMed Abstract  Publisher Full Text 
Ciuparu D, Klie FR, Zhu Y, Pfefferle LJ:
Phys. Chem. B. 2004, 108:3967.
COI number [1:CAS:528:DC%2BD2cXhvFWltbs%3D]
Publisher Full Text 
Xu TT, Zheng J, Wu N, Nicholls WA, Roth RJ, Dikin AD, Ruoff RS:
Nano Lett.. 2004, 4:963.
COI number [1:CAS:528:DC%2BD2cXis1Citrw%3D]; Bibcode number [2004NanoL...4..963X]
Publisher Full Text 
Phys. Rev. B. 1997, 55:16426.
COI number [1:CAS:528:DyaK2sXkt1yls7Y%3D]; Bibcode number [1997PhRvB..5516426B]
Publisher Full Text 
Comput. Lett.. 2005, 1(Special Issue: Clusters):259. Publisher Full Text

Europhys. Lett.. 1997, 39:527.
COI number [1:CAS:528:DyaK2sXmtV2hsL8%3D]; Bibcode number [1997EL.....39..527B]
Publisher Full Text 
Phys. Rev. Lett.. 2001, 86:2416.
COI number [1:CAS:528:DC%2BD3MXhvFSjt70%3D]; Bibcode number [2001PhRvL..86.2416K]
PubMed Abstract  Publisher Full Text 
Phys. Rev. Lett.. 2006, 96:057001.
COI number [1:STN:280:DC%2BD287gtF2msA%3D%3D]; Bibcode number [2006PhRvL..96e7001T]
PubMed Abstract  Publisher Full Text 
Fujimori M, Nakata T, Nakayama T, Nishibori E, Kimura K, Takata M, Sahata M:
Phys. Rev. Lett.. 1999, 82:4452.
COI number [1:CAS:528:DyaK1MXjtlelsbw%3D]; Bibcode number [1999PhRvL..82.4452F]
Publisher Full Text 
Int. J. Quantum Chem.. 1994, 52:1081.
COI number [1:CAS:528:DyaK2cXmvVKqs7o%3D]
Publisher Full Text 
J. Am. Chem. Soc.. 2005, 127:13324.
COI number [1:CAS:528:DC%2BD2MXpsVGhtrg%3D]
PubMed Abstract  Publisher Full Text 
Sergeeva AP, Zubarev DY, Zhai HJ, Boldyrev AI, Wang LS:
J. Am. Chem. Soc.. 2008, 130:7244.
COI number [1:CAS:528:DC%2BD1cXmtVGitbs%3D]
PubMed Abstract  Publisher Full Text 
Zhai HJ, Kiran B, Li J, Wang LS:
Nat. Mater.. 2003, 2:827.
COI number [1:CAS:528:DC%2BD3sXptlemtrY%3D]; Bibcode number [2003NatMa...2..827Z]
PubMed Abstract  Publisher Full Text 
Szwacki NG, Weber V, Tymczak CJ:
Nanoscale Res. Lett.. 2009, 4:1085.
Bibcode number [2009NRL.....4.1085S]
PubMed Abstract  Publisher Full Text  PubMed Central Full Text 
Hagler AT, Dauber P, Lifson S:
J. Am. Chem. Soc.. 1979, 101:5131.
COI number [1:CAS:528:DyaL3cXlvV2r]
Publisher Full Text 
Hagler AT, Dauber P, Lifson S:
J. Am. Chem. Soc.. 1979, 101:5122.
COI number [1:CAS:528:DyaE1MXmtFSks7Y%3D]
Publisher Full Text 
Accelrys Inc.: MS Modeling Getting Started. San Diego: Accelrys Inc.; 2003.

Phys. Rev. A. 1988, 38:3098.
COI number [1:CAS:528:DyaL1cXmtlOhsLo%3D]; Bibcode number [1988PhRvA..38.3098B]
PubMed Abstract  Publisher Full Text 
Phys. Rev. B. 1988, 37:785.
COI number [1:CAS:528:DyaL1cXktFWrtbw%3D]; Bibcode number [1988PhRvB..37..785L]
Publisher Full Text 
J. Chem. Phys.. 1993, 98:5648.
COI number [1:CAS:528:DyaK3sXisVWgtrw%3D]; Bibcode number [1993JChPh..98.5648B]
Publisher Full Text 
Hehre WJ, Stewart RF, Pople JA:
J. Chem. Phys.. 1969, 51:2657.
COI number [1:CAS:528:DyaF1MXltVGnsrY%3D]; Bibcode number [1969JChPh..51.2657H]
Publisher Full Text 
Collins JB, von Ragué Schleyer P, Binkley JS, Pople JA:
J. Chem. Phys.. 1976, 64:5142.
COI number [1:CAS:528:DyaE28XltlSrs78%3D]; Bibcode number [1976JChPh..64.5142C]
Publisher Full Text 
J. Chem. Phys.. 1980, 72:5639.
COI number [1:CAS:528:DyaL3cXksFCnu7c%3D]; Bibcode number [1980JChPh..72.5639M]
Publisher Full Text 
Krishnan R, Binkley JS, Seeger R, Pople JA:
J. Chem. Phys.. 1980, 72:650.
COI number [1:CAS:528:DyaL3cXpvFyitA%3D%3D]; Bibcode number [1980JChPh..72..650K]
Publisher Full Text 
Chem. Phys.. 1981, 61:385.
COI number [1:CAS:528:DyaL3MXmt1Gqtb8%3D]; Bibcode number [1981CP.....61..385B]
Publisher Full Text 
von Rague Schleyer P, Maerker C, Dransfeld A, Jiao H, van Eikema Hommes NJR:
J. Am. Chem. Soc.. 1996, 16:6317. Publisher Full Text

Gaussian 03, Revision D.02, Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA Jr, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V. Gaussian, Inc., Wallingford CT; 2004.

Phys. Rev.. 1962, 126:1028.
Bibcode number [1962PhRv..126.1028M]
Publisher Full Text 
Mol. Phys.. 1974, 27:789.
COI number [1:CAS:528:DyaE2cXkvVektrk%3D]; Bibcode number [1974MolPh..27..789D]
Publisher Full Text 
Nanotechnology. 2006, 17:5063.
COI number [1:CAS:528:DC%2BD28Xht1OjtLvP]; Bibcode number [2006Nanot..17.5063D]
Publisher Full Text 
Phys. Rev. Lett.. 2007, 99:115501.
Bibcode number [2007PhRvL..99k5501T]
PubMed Abstract  Publisher Full Text 
Phys. Rev. B. 2008, 77:041402.
Bibcode number [2008PhRvB..77d1402Y]
Publisher Full Text 
Chem. Phys. Lett.. 1990, 173:569.
COI number [1:CAS:528:DyaK3cXmtlygtLk%3D]; Bibcode number [1990CPL...173..569S]
Publisher Full Text 
Phys. Rev. B. 2006, 74:035413.
Bibcode number [2006PhRvB..74c5413K]
Publisher Full Text 
Chacko S, Kanhere GD, Boustani I:
Phys. Rev. B. 2003, 68:035414.
Bibcode number [2003PhRvB..68c5414C]
Publisher Full Text 
Szwacki NG, Sadrzadeh A, Yakobson BI:
Phys. Rev. Lett.. 2007, 98:166804.
Bibcode number [2007PhRvL..98p6804G]
PubMed Abstract  Publisher Full Text 
Chem. Phys. Lett.. 1992, 192:14.
COI number [1:CAS:528:DyaK38XktlCjtrs%3D]; Bibcode number [1992CPL...192...14K]
Publisher Full Text 
deGiambiagi MS, Giambiagi M, Fortes MD:
J. Mol. Struct. Theochem. 1997, 391:141.
COI number [1:CAS:528:DyaK2sXivFajt7g%3D]
Publisher Full Text 
J. Phys. Chem. A. 1997, 101:1738.
COI number [1:CAS:528:DyaK2sXhsFWnsLw%3D]
Publisher Full Text 
Forte G, Grassi A, Lombardo GM, La Magna A, Angilella GGN, Pucci R, Vilardi R:
Phys. Lett. A. 2008, 372:6168.
COI number [1:CAS:528:DC%2BD1cXhtFWnsLzI]; Bibcode number [2008PhLA..372.6168F]
Publisher Full Text