Uniaxial crystal in 2D. Write a 2D TM mode FDTD simulation in the x-y plane for a uniaxial crystal with no = 1.5 and ne = 2.
(a) Create a slab of this crystal material that spans y and covers two wavelengths in x, and align the optic axis with y. Launch an z plane wave from x = 0 incident on the slab at frequency 1 THz and use 20 grid cells per wavelength. Observe the polarization and amplitude of the transmitted and reflected waves.
(b) Repeat (a), but rotate the crystal in the y-z plane through angles of 15, 30, 45, 60, and 90 degrees, using the rotation matrix in Equation (11.18). How do the transmitted polarization and amplitude change?