Extending the FDTD GVADE method nonlinear polarization vector to include anisotropy.

In this paper the finite-difference time-domain general vector auxiliary differential equation method [Opt. Express14, 8305 (2006)10.1364/OE.14.008305], nonlinear polarization vector, the nonlinear electric dipole moment per unit volume, is extended to include anisotropy, in nonlinear isotropic media at optical frequencies. The theory is presented for extending the numerical method in 3D Cartesian coordinates, and then example simulation results are presented for two isotropic media. First, the simplified 2D transverse magnetic case is revisited for the fused silica example introduced in the 2006 paper, including the anisotropic part of the nonlinear polarization vector in the simulation; the simulation results including the anisotropic part of the polarization vector were compared with the purely isotropic polarization vector simulation results. Second, a simplified 2D transverse magnetic example was simulated in carbon disulfide, with its strong molecular re-orientation-induced polarization anisotropy.