Abstract
A family of unusually stable boron cages was identified and examined using firstprinciples localdensity functional method. The structure of the fullerenes is similar to that of the B_{12}icosahedron and consists of six crossing doublerings. The energetically most stable fullerene is made up of 180 boron atoms. A connection between the fullerene family and its precursors, boron sheets, is made. We show that the most stable boron sheets are not necessarily precursors of very stable boron cages. Our finding is a step forward in the understanding of the structure of the recently produced boron nanotubes.
Keywords:
Boron clusters; Boron fullerenes and nanotubes; Boron sheets; Quantummechanical modelingIntroduction
The chemistry of boron resembles that of carbon in its ability to catenate and form molecular networks. Unlike carbon, bulk boron cannot be found in nature and all known boron allotropes where obtained in the laboratory. All of them are based on different arrangements of B_{12} icosahedrons. It is very natural to believe that not only carbon but also boron posesses molecular allotropes (fullerenes and nanotubes). Experimental and theoretical research on the chemistry of boron nanomaterials is developing rapidly. The existence of quasiplanar [1] and tubular [2] boron clusters was predicted by theory and confirmed more recently by experiment [3,4]. Up to now, however, very little is known about the properties of these novel boron nanostructures.
In this work, we will describe the properties of a family of boron nearly round cages which are built from crossing boron doublerings (DRs). The smallest members of the family, B_{12} and B_{80}, where previously studied using firstprinciples methods [5]. Here, we will show how to construct bigger cages with similar structural characteristics to those found in B_{12} and B_{80}. We will also show the connection between boron cages and nanotubes as well as their precursors—boron sheets.
Method
The calculations were performed within the density functional theory framework, using ultrasoft Vanderbilt pseudopotentials [6] and localdensity approximation for the Perdew–Burke–Ernzerhof exchangecorrelation potential [7]. Computations were done using the planewavebased QuantumESPRESSO package [8]. The optimized geometries of the structures were found by allowing the full relaxation of the atoms in the cell until the atomic forces were less than 10^{−3} and 10^{−4} Ry/Bohr for the atomic cages and sheets, respectively. A proper kpoint sampling for each system together with a 35 Ry cutoff for the planewave basis set have been used to ensure energy convergence to less than 1 meV/atom. To study properties of the fullerenes (nanotubes) the supercell geometry was taken to be a cubic (tetragonal) cell with sufficiently large lattice constant (constants) to avoid interactions between periodic replicas of the cluster. For infinitely long structures the supercell was optimized using variable cell optimization methods included in the program package.
Boron Fullerenes
The unusual stability of the recently proposed fullerene of 80 boron atoms [5] and its structural similarities to the B_{12} icosahedron motivates us to investigate larger cages with similar structural characteristics, that is, built from six crossing DRs but with larger diameters. It is known now that boron clusters with a number of atoms smaller than 20 are rather planar or quasiplanar and the B_{12} icosahedron is energetically less favorable than the quasiplanar convex structure of C_{3v} symmetry [1,3,9], however, we would like to add formally the B_{12} cage to the fullerene family as its smallest member. In Fig. 1, we have shown what we call a next member of the fullerene family, which is made up of 180 atoms. Similarly to B_{80}, this cage posses I_{h} symmetry but is more stable in energy than the B_{80}. As discussed previously [5] and also shown in Fig. 1b, we can clearly identify DRs as fragments of the B_{180} round cage. For B_{12}, B_{80}, and B_{180} the DRs have 10, 30, and 50 atoms, respectively. The only member of the family which is completely close is the B_{12} icosahedron and the surface of the cage consists of 20 triangles. Since the fullerenes are built from a fixed number of DRs, in larger cages the DRs cannot cover the whole sphere with a triangle network of atoms and the fullerenes will exhibit empty spaces—holes. All cages larger than B_{12} will have (at least) 12 holes. The holes in B_{80} are in the shape of pentagons, while the holes in B_{180} are rather circular (more precisely they are closer to decagons; see e.g. Fig. 1b). The size of the empty spaces increases with increasing cage diameter (see Fig. 3a–c).
Figure 1. Top (a) and side (b) views of the optimized B_{180} fullerene. It can be observed in (a) the almost perfect spherical shape of the cage. In (b) it is outlined a DR of 50 atoms. The black and white boron atoms represent additional 20 atoms with respect to the DR in B_{80}. Note that the black atoms are not shared by the DRs but each white atom is shared by two of them. The interatomic distances between neighboring boron atoms are shown in (c) using a fragment of the cage
The increase of the number of atoms by 20 between DRs belonging to two consecutive members of the fullerene family can be explained using the B_{80} and B_{180} cages as follows: each of the DRs in B_{80} is adjacent to 10 pentagons (holes), so if, in order to obtain the B_{180}, we add one atom to each side of the pentagons the number of atoms of every DR increases already by 10. In addition each DR requires still another 10 atoms in order to preserve its structure, so the total number of atoms for every DR increases by 20. We highlighted this in B_{180} in Fig. 1b, where in one of the DRs the additional 20 atoms (respect to a DR in B_{80}) described above were colored in black and white.
The next in size fullerene after B_{180} is made up of 300 atoms (it is built from DRs with 70 atoms). The optimized structure of this cage is shown in the left of Fig. 3c. The total number of atoms in B_{12}, B_{80}, and B_{180} cannot be described by one general formula since the number of atoms shared by the DRs varies from one cage to another. Note in Fig. 1b that the black atoms belong only to one DR, in contrast, the white atoms are shared by two crossing DRs. However, for B_{300}and all larger fullerenes the number of shared atoms is constant and the total number of atoms in the cage can be obtained using a simple formulaN(n) = N(B_{180}) + 120(n − 3), where withn(n ≥ 4) we label the fullerenes starting from the smallest cage. The number of atoms in each of the nfullerene DRs can be expressed by the formula N_{DR}(n) = N_{DR}(B_{12}) + 20(n − 1), where n ≥ 1.
In Fig. 2a, we have plotted the cohesive energy (E_{coh}) of the four fullerenes discussed above versus the number of atoms in the cage. The less stable of them is the B_{12} (E_{coh} = 5.04 eV/atom) and the most stable is the B_{180} (E_{coh} = 5.77 eV/atom), which is 10 meV/atom more stable than the B_{80}. The B_{300} cage has the same E_{coh} as the B_{80} fullerene to within 2 meV/atom. In the case of infinitely large cage E_{coh} cannot be of course calculated exactly, however, it can be approximated by the E_{coh} of an infinitely long stripe made up of boron atoms. Indeed, in a very large cage the six crossing DRs are almost isolated one from each other, and the atoms from the regions where the DRs cross give insignificant contribution to E_{coh} in comparison to the rest of the atoms. The E_{coh} = 5.69 eV/atom of the stripe will be then the lower limit for the energetic stability of large cages. In Fig. 2a, it is shown that the (−E_{coh}), after its minimum at N = 180 atoms, increases for B_{300} (by ∼10 meV/atom), and that tendency should prevail also for larger cages until the lower limit for the E_{coh} is reached. On the other hand it is interesting to note that the stability of DRs increases with increasing radius of the structures and the most stable is a DR with an infinite radius (stripe) [5]. Therefore, the E_{coh} for the boron stripe is not only the lower limit for the energy of large cages but also the upper limit for the energetic stability of DRs.
Figure 2. (a) Cohesive energy per atom as a function of the number of atoms N in the B_{N} cluster. The horizontal line corresponds to the E_{coh} = 5.69 eV/atom of the strip of boron atoms. The lines are a guide to the eye. (b) Radial distance, r(θ), of boron atoms belonging to B_{80} (red circles), B_{180} (blue circles), and B_{300} (black circles) cages, from the center of mass of each cage, as a function of the spherical angle θ. (The red, blue, and black lines correspond to the average values ofr(θ) for B_{80}, B_{180}, and B_{300} cages, respectively
The B_{180} is not only the most stable (in energy) cage from those studied but also possesses almost perfect spherical shape. In Fig. 2b, we have plotted the radial distances, r(θ), of boron atoms belonging to B_{80}, B_{180}, and B_{300} cages, from the center of each cage (center of mass), as a function of the spherical angle θ. The average values for r(θ) are R_{1}= 4.13,R_{2} = 6.85, and R_{3} = 9.39 Å for B_{80}, B_{180}, and B_{300}, respectively. The circles in Fig. 2b represent the positions of boron atoms. In the case of B_{80} there are 20 inner atoms and 60 outer atoms. The 60 outer atoms form a frame of the same structure as exhibited by the C_{60}fullerene and lie on a sphere. The radius of this sphere is slightly larger thanR_{1}. The 20 inner atoms are lying almost exactly at the centers of the hexagons of the B_{60}frame, at a radial distance of ∼0.4 Å from the larger sphere. In the case of the B_{180}fullerene the atoms are lying almost perfectly (to within 0.1 Å) on a sphere of radiusR_{2}. More complex picture is present in the case of the B_{300}fullerene. Here half of the atoms lie inside a sphere of radiusR_{3}and half of them outside of this sphere. The more distant lie ∼0.2 Å above or below the sphere surface. Should be pointed out that B_{80}and B_{300}exhibit braking symmetry distortions, which are however very small. A more detailed analysis has to be done to determine the nature of these structural distortions. In Table 1, we summarized our results forE_{coh}and interatomic distances between neighboring boron atoms for boron fullerenes and also boron sheets (BSs) which are studied in the next section.
Table 1. Point symmetries, cohesive energies, and interatomic distances,d_{BB}, between neighboring boron atoms for fullerenes and sheets
The structure of the cages influences also their electronic properties. For B_{180} the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are triply degenerate. (Similar result was reported previously for B_{80}[5].) For B_{300} the triple degeneracy both of HOMO and LUMO is slightly lifted (by less than 95 meV). This is most probably the result of the structural distortions mentioned above. The HOMO–LUMO energy gaps are 0.43 and 0.10 eV for B_{180} and B_{300}, respectively. We have also calculated the electronic properties of the boron strip which was found to be metallic.
Boron Sheets
To deepen the understanding of boron nanotubes and also boron cages it is important to know what the structure of the BS is. Several theoretical efforts [1015], using firstprinciples methods, have been done so far to understand the structure and properties of BSs. Most of these investigations have determined that the buckled triangle atomic lattice represents the most stable structure for the BS. However, it is intuitively understandable that the buckling of boron atoms is a response of the sheet to the internal stress imposed by the arrangement of the atoms in a triangle lattice. In principle it could not be discarded that there is an alternative arrangement in which the boron atoms may stay covalently bounded in a plane without buckling. Lau et al. [11] proposed that the BS can be formed by a network of triangle–square–triangle units of boron atoms. This structure, although planar, is less stable than the buckled triangle atomic lattice what was later confirmed by the authors [12]. Very recently it was proposed a new BS which resembles the carbon honeycomblike structure [16]. It is fully planar, possesses metallic properties, and is energetically more stable than the boron buckled TS. This structure is shown in Fig. 3a (bottom) and will be labelled by us as αBS.
Figure 3. Three members of the fullerene family are shown: (a) B_{80}, (b) B_{180}, and (c) B_{300}. Each fullerene is accompanied by its precursor sheet. In all cages and sheets are highlighted (in blue) the corresponding atomic motifs which are discussed in the text. The unit cells used for calculations of sheets are shown in red
Our investigation of possible candidates for the BS we restricted to those which have structural similarities with the fullerenes studied above. In Fig. 3, we have shown three cages, B_{80}, B_{180}, and B_{300}, and their corresponding sheets αβ, and γ, respectively. In the top of Fig. 3a–c, we have highlighted on each cage the characteristic atomic motif of the fullerene which also will appear on the sheet corresponding to it (see the sheets in the bottom of Fig. 3a, b and the sheet in the right of Fig. 3c). Let us forget for a moment about the holes in the cages, then the characteristic motif for the B_{80} is a cluster of 7 atoms with one central atom lying almost in plain defined by a hexagonal chain of 6 atoms. This cluster has C_{3v} symmetry. An isolated neutral cluster made up of 7 atoms has C_{2v} symmetry [9]. The next two cages, B_{180} and B_{300}, have motifs which are similar in shape and consist of quasiplanar structures of 12 and 18 atoms, respectively, and C_{3v} symmetry. These clusters have 9 and 12 peripheral atoms and 3 and 6 central atoms in B_{180} and B_{300} cages, respectively. The interatomic distances between neighboring boron atoms in the fragment of B_{180} are shown in Fig. 1c. As it was mentioned in the previous section the isolated neutral cluster of 12 atoms is a quasiplanar convex structure of C_{3v} symmetry. It was experimentally determined that this cluster has unusually large HOMO–LUMO gap of 2 eV [3]. It was also suggested that this cluster must be extremely stable electrically and should be chemically inert [9]. Perhaps these unusual characteristics are responsible also for the outstanding stability of the B_{180} fullerene.
All BSs studied in this work are fully planar and have metallic properties. We have found that the αBS which is a precursor of the B_{80} cage is the most stable sheet over all studied (in agreement with recent findings [16]). The point symmetry of that sheet is D_{6h}. Figure 3a (bottom) shows the unit cell used for calculations of the αBS. The BSs β and γ, corresponding to B_{180} and B_{300} fullerenes, respectively, have D_{3h} symmetry. The shape of the holes in BSs will be determined by the type of atomic motif we are using to build the sheet. In the case of the αBS the holes are hexagons and in the case of β and γ BSs the holes are distorted dodecagons and hexagonslike, respectively. It is important to observe that the BSs shown in Fig. 3 can also be seen as built from interwoven boron stripes. This observation may help to understand not only structural but also electronic properties of BSs.
TheE_{coh}for theαBS is 5.94 eV/atom and is bigger than theE_{coh}for theβandγBSs by 0.14 and 0.17 eV/atom, respectively. This is an interesting result since it means that the most stable structure for the BS does not necessarily have to be a precursor of very stable (in energy) boron cages.
Although there are some discrepancies between the E_{coh} values obtained in this work (see Table 1) and reported in the literature, the differences between E_{coh} values corresponding to different sheets match very well. Indeed the E_{coh} for the buckled TS is higher in energy than E_{coh} for the flat TS by 0.24 eV/atom and this value is close to 0.21 and 0.22 eV/atom reported in Refs. [16,12], respectively. Similarly, E_{coh} for the sheet α is higher in energy than E_{coh} for the buckled sheet by 0.08 eV/atom, what represents slightly smaller value than 0.11 eV/atom obtained in previous calculations [16].
Fullerenederived Nanotubes
Carbon nanotubes and fullerenes are closely related structures. Capped nanotubes are elongated fullerenelike cylindrical tubes which are closed at the rounded ends [2]. To look for similar connections between boron nanostructures we have investigated boron nanotubes derived from the B_{80} cages.
B_{80}derived tubule can be obtained, in simplest way, by bisecting a B_{80} molecule at the equator and joining the two resulting hemispheres with a cylindrical tube one monolayer thick and with the same diameter as the B_{80}. Almost all boron nanotubes studied theoretically so far are rolled up TSs of boron atoms [10,1315]. It was natural then to take a cylindrical tube made of a triangle lattice as a candidate to join the two cups obtained from B_{80}. We have investigated only the simplest case when the fullerene is divided exactly at the middle of one of the DRs. The two hemispheres were then joined with a (15, 0) tube of “zigzag” geometry. However, upon relaxation this cupped nanotube significantly deformed. Since the B_{80} cage can be seen as built from DRs we decided to design similarly the body of the nanotube as built of crossing stripes of double chains of boron atoms. Schematic representation of the resulting nanotube of 240 atoms is shown in Fig. 4a. The two B_{80} hemispheres are shown in red and the nanotube is shown in blue. After structural optimization (see Fig. 4c) the nanostructure preserved its tubular form and the cups remained almost unchanged. It turns out that the body of this nanotube has a similar structure as the αBS.
Figure 4. (a) Schematic representation of the B_{240} nanotube described in the text. Optimized structures of the capped B_{160}(b) and B_{240}(c) nanotubes and the infinite (15, 0) nanotube (d). The clusters in (b) and (d) are shown in side (left) and front (right) views. The caps are two hemispheres of the B_{80} molecule and the tubes are wrapped αBSs. For the infinite nanotube a unit cell of 80 atoms was used for calculations
The shortest nanotube is, of course, the B_{80} cage and its E_{coh} is 5.76 eV/atom [5]. We have also optimized the structures of two longer finite tubes – the B_{160} (see Fig. 4b) and the previously described B_{240} (see Fig. 4c). As expected the stability of the nanostructures increases with increasing lengths of the tube: E_{coh}(B_{160}) = 5.81 eV/atom and E_{coh}(B_{240}) = 5.84 eV/atom. The most stable is the infinite nanotube (see Fig. 4d) with E_{coh} = 5.87 eV/atom. The HOMO–LUMO energy gaps for the B_{160} and B_{240} clusters are 0.33 and 0.02 eV, respectively, and the infinite tube was found to be metallic.
Summary
We are predicting the existence of a family of very stable boron fullerenes. The cages have similar structure consisting of six interwoven boron DRs. The most stable fullerene is made up of 180 atoms and has almost perfect spherical shape. A recently proposed very stable BS of triangular and hexagonal motifs is a precursor of the B_{80}cage. However, it was shown that the most stable sheets are not necessarily the precursors of very stable boron cages. Finally, we have shown that the proposed fullerenes and novel boron nanotubes are closely related structures.
Acknowledgment
The Interdisciplinary Centre for Mathematical and Computational Modelling of Warsaw University is thanked for a generous amount of CPU time.
References

Phys. Rev. B. 1997, 55:16426.
COI number [1:CAS:528:DyaK2sXkt1yls7Y%3D]
Publisher Full Text 
Gindulyte A, Lipscomb WN, Massa L:
Inorg. Chem.. 1998, 37:6544.
COI number [1:CAS:528:DyaK1cXntl2qtr8%3D]
Publisher Full Text 
Zhai HJ, Kiran B, Li J, Wang LS:
Nat. Mater.. 2003, 2:827.
COI number [1:CAS:528:DC%2BD3sXptlemtrY%3D]
Publisher Full Text 
Ciuparu D, Klie RF, Zhu Y, Pfefferle L:
J. Phys. Chem. B. 2004, 108:3967.
COI number [1:CAS:528:DC%2BD2cXhvFWltbs%3D]
Publisher Full Text 
Gonzalez Szwacki N, Sadrzadeh A, Yakobson BI:
Phys. Rev. Lett.. 2007, 98:166804. Publisher Full Text

Phys. Rev. B. 1990, 41:7892. Publisher Full Text

Perdew J P, Burke K, Ernzerhof M:
Phys. Rev. Lett.. 1996, 77:3865.
COI number [1:CAS:528:DyaK28XmsVCgsbs%3D]
Publisher Full Text 
Alexandrovaa A N, Boldyreva AI, Zhaib HJ, Wangb LS:
Coord. Chem. Rev.. 2006, 250:2811. Publisher Full Text

Phys. Rev. B. 2006, 74:035413. Publisher Full Text

Lau K C, Pati R, Pandey R, Pineda AC:
Chem. Phys. Lett.. 2006, 418:549.
COI number [1:CAS:528:DC%2BD28XotlCgsw%3D%3D]
Publisher Full Text 
J. Phys. Chem. C. 2007, 111:2906.
COI number [1:CAS:528:DC%2BD2sXptVGrtg%3D%3D]
Publisher Full Text 
Kunstmann J, Quandt A, Boustani I:
Nanotechnology. 2007, 18:155703. Publisher Full Text

Cabria I, López MJ, Alonso JA:
Nanotechnology. 2006, 17:778.
COI number [1:CAS:528:DC%2BD28Xis1yhsb4%3D]
Publisher Full Text 
Evans MH, Joannopoulos JD, Pantelides ST:
Phys. Rev. B. 2005, 72:045434. Publisher Full Text

Phys. Rev. Lett.. 2007, 99:115501. Publisher Full Text