Boron-rich compounds (IBCs) such as boron suboxide (B6O) and boron subphosphide (B12P2) are derived from the α-rhombohedral boron lattice with extreme hardness and unusual electrical properties. These compounds display an extraordinary self-healing ability to repair the lattice defects generated from exposure to high-energy irradiation. This unique ability will enable these boron-rich compounds to be deployed in nuclear battery devices for space or deep-sea exploration applications when coupled with their high hole mobility.Synthesis and control of the stoichiometry of IBCs are the key technical challenge. For instance, the electrical properties of B6O are susceptible to specific B:O ratios. Thus far, the highest quality B6O crystals were obtained from high-pressure, high-temperature (HP-HT) experiments using mixed boron and boron oxide powders. Nevertheless, a mature synthesis protocol yielding consistent, high-quality B6O crystals has not been established.In this talk, density functional theory (DFT) was used to understand how the structural, electronic, and thermodynamic stabilities of B6O (and other IBCs) are influenced by the interstitial elements and point defects at the icosahedral sites. Using the hybrid HSE functional, we confirmed that the perfect B6O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening the α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. Then, the thermodynamic calculations based on the DFT method predict that formations of interstitial defects become favorable only at higher temperatures.To understand the self-healing property, the nudged elastic band (NEB) method was employed to identify the minimum energy pathways for the diffusions of dislocated B and O atoms. The diffusion of icosahedral B atoms has an energy barrier of 0.16 eV, confirming that the boron vacancies can be repaired with little kinetic limitation. Nevertheless, a substantially perturbed lattice requires more complex structural reorganization and thus incurs higher barriers (> 1 eV).Lastly, reactive molecular dynamics simulations with a frozen core were performed to gain insights into the initial crystal nucleation process using a newly developed ReaxFF force field. For the first time, we can gain a molecular perspective on the formations of the boron icosahedra (B12) motif and the energetic profile. Such knowledge will aid the experimental design of IBC synthesis.
Read full abstract