where V is the primitive-cell volume, ω is the phonon frequency, D is the dynamical matrix, and e is the eigenvector. The above equation can be approximated using the finite difference method. In our workflow, phonopy is used to compute the mode-dependent Grüneisen parameters on a regular mesh and along a high-symmetry path. The average Grüneisen parameters are obtained following Ref. 97 as