First-principles calculation of the lattice thermal conductivities of α-, β-, and γ-Si$_3$N$_4$
Lattice thermal conductivities (LTCs) of α-, β-, and γ-Si$_3$N$_4$ single crystals are investigated from ab initio anharmonic lattice dynamics, within the single-mode relaxation-time approximation of the linearized phonon Boltzmann transport equation. At a temperature of 300 K, a κ$_{xx}$ of 70 and a κ$_{zz}$ of 98 (in units of Wm$^{-1}$K$^{-1}$) are obtained for α-Si$_3$N$_4$. For β-Si$_3$N$_4$, κ$_{xx}$ and κ$_{zz}$ are found to be 71 and 194, respectively, which are consistent with the reported experimental values of 69 and 180 for individual β-Si$_3$N$_4$ grains in a ceramic. The theoretical κ$_{xx}$ values of α- and β-Si$_3$N$_4$ are comparable, while the κ$_{zz}$ value of β-Si$_3$N$_4$ is almost twice that of α-Si$_3$N$_4$, which demonstrates the very large anisotropy in the LTC of the βphase. It is found that the large anisotropy in the LTC of β-Si$_3$N$_4$ was caused by the elongated Brillouin zone along the c* axis, where the acoustic phonons mostly contribute to LTC and have large group velocities even near the Brillouin zone boundary. The LTC of γ-Si$_3$N$_4$ is 81, which is as small as that of α-Si$_3$N$_4$, although γ-Si$_3$N$_4$ has much larger elastic constants. This means that elastic constants are not always a good indicator of LTC. We show that knowing the detailed distributions of both the group velocities and phonon lifetimes in the Brillouin zones is important for characterizing the LTC of the three phases.