Further stability criteria can be sought from the phonon dispersion relations in the Brillouin zone, i.e. the phonon band structures. Following the method presented in the 'Computational framework' section, the obtained phonon band structures for three BP polymorphs are shown in Fig. . In each panel the bands run along the main lines of the cubic BZ (Fig. ) and hexagonal BZ (Figs. and).