We use the first-principles approach to clarify the thermodynamic stability as a function of pressure and temperature of three different α-rhombohedral-boron-like boron subnitrides, with the compositions of B 6 N, B 13 N 2 , and B 38 N 6 , proposed in the literature. We find that, out of these subnitrides with the structural units of B 12 (N-N), B 12 (NBN), and [B 12 (N − N)] 0.33 [B 12 (NBN)] 0.67 , respectively, only B 38 N 6 , represented by [B 12 (N − N)] 0.33 [B 12 (NBN)] 0.67 , is thermodynamically stable. Beyond a pressure of about 7.5 GPa depending on the temperature, also B 38 N 6 becomes unstable, and decomposes into cubic boron nitride and α-tetragonalboron-like boron subnitride B 50 N 2 . The thermodynamic stability of boron subnitrides and relevant competing phases is determined by the Gibbs free energy, in which the contributions from the lattice vibrations and the configurational disorder are obtained within the quasiharmonic and the mean-field approximations, respectively. We calculate lattice parameters, elastic constants, phonon and electronic density of states, and demonstrate that [B 12 (N − N)] 0.33 [B 12 (NBN)] 0.67 is both mechanically and dynamically stable, and is an electrical semiconductor. The simulated x-ray powder-diffraction pattern as well as the calculated lattice parameters of [B 12 (N − N)] 0.33 [B 12 (NBN)] 0.67 are found to be in good agreement with those of the experimentally synthesized boron subnitrides reported in the literature, verifying that B 38 N 6 is the stable composition of α-rhombohedral-boron-like boron subnitride.