Generalized Huang's Equation for Phonon Polariton in Polyatomic Polar Crystal

The original theory of phonon polariton is Huang's equation which is suitable for diatomic polar crystals only. We proposed a generalized Huang's equation without fitting parameters for phonon polariton in polyatomic polar crystals. We obtained the dispersions of phonon polariton in GaP (bulk), hBN (bulk and 2D), {\alpha}-MoO3 (bulk and 2D) and ZnTeMoO6 (2D), which agree with the experimental results in the literature and of ourselves. We also obtained the eigenstates of the phonon polariton. We found that the circular polarization of the ion vibration component of these eigenstates is nonzero in hBN flakes. The result is different from that of the phonon in hBN.