Prediction of a new ground state of superhard compound B6O at ambient conditions.

Boron suboxide B6O, the hardest known oxide, has an Rm crystal structure (α-B6O) that can be described as an oxygen-stuffed structure of α-boron, or, equivalently, as a cubic close packing of B12 icosahedra with two oxygen atoms occupying all octahedral voids in it. Here we show a new ground state of this compound at ambient conditions, Cmcm-B6O (β-B6O), which in all quantum-mechanical treatments that we tested comes out to be slightly but consistently more stable. Increasing pressure and temperature further stabilizes it with respect to the known α-B6O structure. β-B6O also has a slightly higher hardness and may be synthesized using different experimental protocols. We suggest that β-B6O is present in mixture with α-B6O, and its presence accounts for previously unexplained bands in the experimental Raman spectrum.

Superhard materials are used in many applications, from cutting, grinding and drilling tools to wear-resistant coatings123. However, most superhard materials4, such as diamond5 and cubic-BN6, are synthesized at high pressure, which makes them expensive, but some (boron allotropes, B6O, B4C) are thermodynamically stable already at ambient conditions. The hardness of α-B6O7 was reported to be in the range 30–45 GPa89, making it the hardest known oxide91011.

Objects with icosahedral symmetry (Ih) bear a special fascination because of incompatibility of fivefold symmetry with crystalline periodicity. The discovery of multiply-twinned particles B6O, an icosahedral packing of B12 icosahedra with Ih symmetry, had aroused interest7. Here we report the prediction of a new phase of B6O, with space group Cmcm, which we name β-B6O. This structure is energetically almost degenerate with α-B6O (and slightly more stable), is predicted to have a higher hardness, and actually corresponds to twinned α-B6O structure.

Results

Discovery of β-B6O at ambient conditions

Our variable-composition evolutionary searches expectedly find B2O3 and B6O to be the only stable compounds in the B-O system. Interestingly, there are also several compounds very close to stability - B2O7 (this is a 2D-form of B2O3 intercalated with oxygen molecules) and oxygen-deficient versions of B6O with B6O-like structures and compositions between B and B6O. To our surprise, Cmcm-B6O (β-B6O, see Table 1 for structural parameters), instead of the well-known Rm-B6O (α-B6O)7121314, turned out to be the most stable structure at ambient pressure, as shown in Fig. 1; phonon calculations have confirmed its dynamical stability. Transmission electron microscopy recently confirmed its existence15. Structural parameters and some of the physical properties of β-B6O are shown in Table 1, in comparison with α-B6O and two related forms of pure boron.

Table 1

Calculated structural parameters, hardness, elastic moduli and band gap of α-B, Cmcm-B, α-B6O, and Cmcm-B6O phases at ambient pressure.

Phase	α-B	Cmcm-B	α-B6O	Cmcm-B6O
Space group	Rm	Cmcm	Rm	Cmcm
V, Å3/atom	7.248	7.262	7.387	7.384
Cell parameters	a = b = c = 5.050 Å, α = 58.04°	a = 4.883 Å, b = 8.852 Å, c = 8.064 Å α = β = γ = 90°	a = b = c = 5.153 Å, α = 63.10°	a = 5.393 Å, b = 8.777 Å, c = 8.736 Å α = β = γ = 90°
Atomic coordinates	B1(0.654, 0.010, 0.010) B2(0.630, 0.221, 0.221)	B1(0.000, 0.236, 0.568) B2(0.500, 0.937, 0.576) B3(0.823, 0.167,0.750) B4(0.797, 0.831,0.639) B5(0.682, 0.995,0.750)	B1(0.998, 0.998, 0.667) B2(0.676, 0.201, 0.201) O(0.622, 0.622, 0.622)	B1(0.000, 0.756, 0.588) B2(0.000, 0.549, 0.584) B3(0.165, 0.824, 0.750) B4(0.238, 0.155, 0.649) B5(0.334, 0.987, 0.750) O(0.000, 0.840, 0.439)
Hv(Chen), GPa	39	35	38	39
Hv(Lyakhov), GPa	33.0	32.7	31.6	31.7
B0, GPa	212(214a)	208	227(181b)	226
G, GPa	201	189	208	209
DFT band gap, eV	1.457	1.772	1.854c	1.805c

aTheory, reference 40.

bExperiment, reference 41.

cSee Supplementary Fig. S1 for detailed band structures.

Figure 1

Stability of Cmcm-B6O.

(a) Convex hull of the B-O system at ambient pressure. The solid (hollow) points represent the stable (metastable) structures. (b) Enthalpy difference between β-B6O and α-B6O, including zero-point energy.

In order to further compare the stability of β-B6O and α-B6O, we calculated their enthalpies as a function of pressure, as shown in Fig. 1b. We found that the enthalpy of β-B6O is lower than that of α-B6O at ambient pressure, but the energy difference is only about 1.8 meV/formula within the GGA (and almost degenerate within the HSE06 hybrid functional). As pressure increases, β-B6O becomes progressively more favorable than α-B6O, indicating that β-B6O might be more easily synthesized under pressure. The enthalpies of the two structures are so close that it makes us think: will the two structures coexist? what is their relationship? how to synthesize β-B6O? In order to answer these questions, we perform a detailed comparison of their structure, Raman spectra and phonon densities of states (PDOS).

Comparison of crystal structures of α-B6O and β-B6O

β-B6O structure has hexagonal close packing of B12 icosahedra (ABAB… stacking), while α-B6O is based on the cubic close packing (ABCABC… stacking) of B12 icosahedra, as shown in Fig. 2d,c. As is the case of hcp and fcc metals, twinning of α-B6O can produce local β-B6O stackings. It may also be possible to obtain β-B6O-like stacking faults by deformation of α-B6O, through plane slips. Most properties of these two phases are very similar: e.g. predicted volume per formula (V(α-B6O = 51.71 Å3/formula; V(Cmcm-B6O) = 51.69 Å3/formula), hardness (Hv (α-B6O) = 38 GPa; Hv (Cmcm-B6O) = 39 GPa), elastic moduli (Table 1), DFT band gaps (α-B6O has a 1.85 eV direct band gap, while Cmcm-B6O has a 1.81 eV indirect band gap).

Figure 2

Crystal structures of (a) α-B, (b) Cmcm-B, (c) α-B6O, (d) Cmcm-B6O, and their local structures, (e,f) B12 icosahedra. Green (large) and red (small) spheres denote B and O atoms, respectively.

Another interesting aspect is that if we remove the oxygen atoms from α-B6O and Cmcm-B6O, they turn into α-B16 and Cmcm-B, respectively (Fig. 2a,b). α-B and Cmcm-B are energetically nearly degenerate structures of boron at low pressure17, while Cmcm-B is a newly predicted structure1819. As shown in Fig. 2b, displacements of layers I and II can transform this structure into α-B, and vice versa; Cmcm-B6O and α-B6O have a similar relationship (Fig. 2d). Furthermore, we found that the conversion of Cmcm-B6O and α-B6O will cause a deflection of B-O bond by 36°, as shown in their local structure (Fig. 2e,f).

However, it should be pointed out that it is not easy to change the stacking in covalent systems, though some examples are known – lonsdaleite (metastable “hexagonal diamond”) is formed in shocked cubic diamond. To obtain multiple polytypes, methods like physical vapor transport (PVT), also known as seeded sublimation growth, can be used: e.g., different polytypes of SiC were obtained using the PVT method20.

Comparison of Raman spectra

As mentioned above, β-B6O and α-B6O are energetically almost degenerate at zero pressure, structurally related and can coexist. To test the latter, we calculated their Raman spectra and compared with the experimental results121321. In Fig. 3, the topmost curve is the Raman spectrum reported by Solozhenko et al.12; below it are the Raman frequencies reported by Werheit and Kuhlmann21 and marked by vertical bars (|). The two curves below are our Raman spectra of β-B6O and α-B6O. The experimental data correspond to normal isotopic abundance, and so do our calculations. Correspondingly, the atomic mass of boron in the calculations adopt the weighted average value, i.e. 10.811, which is based on the isotopic abundance of boron. As one could expect, the Raman spectra of β-B6O and α-B6O are similar and match perfectly with the experimental spectra. For example, the Raman modes at 374, 499, 541, 737, 785, 833, 1088, 1119, 1141 cm−1 are consistent with Solozhenko’s results12. However, there are four Raman modes which are unique for β-B6O. The first two are computed to be at 516 and 533 cm−1 (marked by letter A in Fig. 3), and Werheit and Kuhlmann indeed observed these two modes at 519 and 534 cm−1. The third mode is predicted to be at 630 cm−1 (marked with letter B), and Werheit and Kuhlmann have indeed observed a Raman-active phonon at 627 cm−1. We note that α-B6O does not have Raman-active modes between 570 and 700 cm−1, thus the one observed by Werheit and Kuhlmann at 627 cm−1 cannot come from α-B6O, but β-B6O. The fourth mode has the theoretical Raman frequency of 896 cm−1 (marked with letter C), while α-B6O has no Raman-active modes between 850 and 1000 cm-1. Thus, this mode is also unique to the β-B6O structure, and again seen in experiments: Solozhenko et al. observed Raman active phonon at 889 cm−1, and Werheit and Kuhlmann observed a Raman-active phonon at 902 (909) cm−1. Moreover, Wang et al. also observed similar Raman spectra in their B6O samples (see Fig. 1 in ref. 13). This analysis clearly shows that experimental samples contain β-B6O.

Figure 3

Raman spectra of B6O experimental spectrum of Solozhenko et al.12, Raman mode frequencies experimentally observed by Werheit and Kuhlmann21, and our calculated Raman spectra of Cmcm-B6O and α-B6O.

Comparison of PDOSs and Gibbs free energy

Comparing phonon densities of states (PDOSs) of α-B6O and β-B6O at ambient pressure (Fig. 4), we once again see a great degree of similarity. In order to further confirm the stability of β-B6O, we have calculated the Gibbs free energy of β-B6O and α-B6O as a function of temperature, shown in Supplementary Fig. S2. We conclude that β-B6O remains more stable than α-B6O also when temperature is taken into account – and is even slightly stabilized by thermal effects.

Figure 4

Phonon density of states of α-B6O phase and Cmcm-B6O phase at ambient pressure.

For clarity, the PDOS of α-B6O was shifted.

Having found that β-B6O should be more stable than α-B6O, and demonstrated that the two phases actually coexist in experimental samples, we ask a question: why is the synthetic compound mostly α-B6O, instead of the more stable β-B6O? While there can be no definitive answer at this point, we suggest that this may be because of the use of α-rhombohedral-B (Rm)9 or β-rhombohedral-B (Rm)22 as a starting material (together with B2O3) for synthesis of rhombohedral-B6O. Hubert et al.23 used amorphous-B as the starting material, but they got not only B6O but also B6O twinning particles and some amorphous phases. And they observed a stacking of ABAB… around the planar defect, which is the same as β-B6O (We speculate it probably is β-B6O). To obtain β-orthorhombic-B6O, one would need to crystallize it from the melt (preferably at high pressure, to increase its thermodynamic advantage over α-B6O), or use other phases of boron as precursors. Using Cmcm-B would be ideal, but this phase (just 11 meV/atom higher in energy than α-B19) remains hypothetical, though likely to be eventually synthesized. Alternatively, the PVT method20 could be used.

In summary, to our big surprise, ab initio structure prediction calculations discovered a new ground state for the widely studied superhard compound B6O – our predicted β-B6O is more stable than experimentally known α-B6O. The two phases are polytypes and have nearly the same densities (β-B6O is slightly denser), energies (slightly lower for β-B6O), band gaps (slightly smaller and indirect, rather than direct, for β-B6O), hardnesses (β-B6O is slightly harder) and phonon densities of states, but have important differences in Raman spectra. By comparing calculated and experimental Raman spectra, we demonstrated that the experimental samples are actually a mixture of α-B6O and β-B6O. The discovery of β-B6O opens up new possibilities, in view of its greater stability, hardness and indirect band gap. Our findings also indicate possibilities of tuning the properties of B6O by obtaining phase-pure samples (probably not obtained to date), and the possibility of metastable oxygen-deficient compounds based on α-B6O or β-B6O - these can be obtained at high temperatures (where disordered oxygen vacancies will stabilize the structure) and low chemical potentials of oxygen.

Methods

We used the ab initio evolutionary algorithm USPEX24252627 to search for thermodynamically stable B-O compounds and their structures at ambient pressure. This methodology has shown its predictive power in many studies (e.g., ref. 17 and 27, 28, 29). All structures were relaxed; structure relaxations and total energy calculations were done using density functional theory (DFT) within the generalized gradient approximation (GGA)30 as implemented in the VASP code31, with the projector-augmented wave method32. We used plane-wave kinetic energy cutoff of 600 eV, and sampled the Brillouin zone with uniform Γ-centered meshes of is 2π*0.07 Å−1 resolution within structure search, and 2π*0.04 Å−1 for subsequent highly precise relaxations and properties calculations. In order to confirm the relative stability of α-B6O and β-B6O, we used local density approximation (LDA)33 and HSE06 hybrid functional34. Phonon spectra was computed by PHONOPY35 and VASP, and Raman spectra were calculated using the Fonari-Stauffer method36. Hardness was calculated with Chen model37 and Lyakhov-Oganov model38. Elastic moduli were computed using Voigt-Reus-Hill averaging39.

Additional Information

How to cite this article: Dong, H. et al. Prediction of a new ground state of superhard compound B6O at ambient conditions. Sci. Rep. 6, 31288; doi: 10.1038/srep31288 (2016).