Extending the FDTD GVADE method nonlinear polarization vector to include anisotropy

TL;DR

This paper extends the FDTD GVADE method to incorporate anisotropic nonlinear polarization, enabling more accurate simulations of anisotropic media at optical frequencies.

Contribution

The paper introduces a theoretical extension of the FDTD GVADE method to include anisotropic nonlinear polarization in 3D Cartesian coordinates.

Findings

Simulation results show improved accuracy when including anisotropy.

Comparison with isotropic models highlights the significance of anisotropic effects.

Application to fused silica and carbon disulfide demonstrates method versatility.

Abstract

In this paper the finite-difference time-domain general vector auxiliary differential equation method [Greene, J. H. and A. Taflove, Opt. Express 14, 8305 (2006)], 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.