Structure, Stability, and Electronic Properties of Boron Suboxide: A Density Functional Theory Study

Abstract

Boron suboxide (B 6 O) is a boron-rich compound derived from the αrhombohedral boron lattice with extreme hardness and unusual semiconducting properties. In this work, density functional theory (DFT) was used to show that unit cell volume, mechanical strength, band gaps, and thermodynamic stabilities of B 6 O were influenced by the interstitial elements and point defects at the icosahedral sites. While the hexagonal unit cell volume (HUCV) varies with interstitial occupancy, it is the icosahedral defect that weakens the intrinsic bulk modulus of B 6 O. Using the hybrid HSE functional, we confirmed that the perfect B 6 O bulk is a p-type semiconductor with a direct band gap of 2.8 eV. Furthermore, by screening α-boron compounds systematically, we found that a simple octet rule may offer a consistent explanation for the variations in the computed electronic structures. The formation free energies calculated over a wide range of temperatures (0−2500 K) and pressures (0−80 GPa) predict that formations of interstitial defects become favorable only at higher temperatures (ca. 1800 K) in bulk B 6 O lattices. 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. More complex B diffusion paths involving the reorganization of icosahedral boron atoms incur higher barriers (>1 eV). In contrast, the diffusion of interstitial O atoms is facile with a barrier of 0.4 eV. Lastly, successive O insertions into the α-B lattice were performed using DFT to generate a basic understanding of the oxidation process. These calculations provide fundamental atomistic insights into the growth of B 6 O crystals and control of their point defects.

Predictions of Boron Phase Stability Using an Efficient Bayesian Machine Learning Interatomic Potential

“…For defect-free systems, vibrational entropy plays a central role in dictating the free energies of solid-state boron crystals . Here, the vibrational entropies and zero-point energies were extracted from the MD simulation trajectories using Dump2phonon and its auxiliary post-processing tool Phana. , Unlike Phonopy, where atom displacements are generated manually, Dump2phonon constructs a dynamical matrix from the actual displacements occurring during MD simulations.…”