where ω$_{L}$ is the frequency of the longitudinal phonon mode, and ω$_{T}$ is the frequency of the transverse phonon mode. Since there is no dynamic ME coupling, the dispersion relation for bulk polaritons can be calculated from the macroscopic electromagnetic wave equation For the TE modes, the constitutive connecting ⃗ B with ⃗ H are defined by ⃗ µ = 1 + 4 π⃗ χ$_{m}$ . Using these relations, the solution of the wave equations for the bulk modes is of the form and an implicit expression for the bulk mode frequency is: Equation (26) has two zeros, one is for ϵ$_{x}$ = 0 and the other for f ( µ ) = µ$_{y}$µ$_{z}$ − µ 2 yz = 0. It also has two resonance poles from ϵ$_{x}$ and another from the condition µ$_{y}$ = 0. As a result, there are three bands of bulk polariton modes. The dispersion relation for surface modes is obtained by assuming solutions in the form: and where β and β$_{o}$ are positive real attenuation constants for the sample and vacuum, re- spectively. An implicit relation for the attenuation factor β of the medium is derived by substituting Eq. (27) into the wave equation Eq. (24): An explicit relation for the attenuation constant β$_{o}$ (in vacuum) is given by 7 β 2 o = k 2 y − ( ω c ) 2 . (30) µ$_{z}$β 2 = µ$_{y}$k 2 y − ϵ$_{x}$ ( ω c ) $^{2}$( µ$_{y}$µ$_{z}$ − µ 2 yz ) (29) ⃗ H ∼ e − β$_{o}$$^{z}$ei ( k$_{y}$y − ωt ) for z > 0 (28) ⃗ H ∼ e $^{βz}$ei ( k$_{y}$y − ωt ) for z < 0 (27) µ$_{y}$k 2 y = ϵ$_{x}$ ( ω c ) $^{2}$( µ$_{y}$µ$_{z}$ − µ$_{yz}$ $^{2}$) . (26) ⃗ H ∼ e i ( k$_{y}$y + k$_{z}$z − ωt ) (25) ∇ 2 ⃗ H − ∇ ( ∇ · ⃗ H ) − ϵ c $_{2}$⃗ µ · ∂ 2 ⃗ H/∂t 2 = 0 . (24) IV. RESULTS AND DISCUSSION