Numerical analysis of nonlinear electromagnetic waves in nematic liquid crystal cells

In the current work, the nonlinear problem of electromagnetic wave propagation in a Nematic Liquid Crystal (NLC) cell is solved numerically. The LC is sandwiched between two glass layers of finite thickness and a linearly polarized beam is obliquely incident to the cell. The dielectric properties of N-LCs depend on the tilt angle of the directors. When the excitation beam enters the cell, and providing the incident intensity is above the Freedericksz threshold, the directors reorient themselves changing the LC's relative permittivity tensor. In turn, this affects beam propagation throughout the crystal. The electromagnetic field is modeled by the time-harmonic Maxwell equations whereas the director field is governed by a nonlinear ordinary differential equation (ODE). Our solution method is iterative, consistently taking into account this interaction between the excitation beam and the director field. The Maxwell equations are solved employing the Mode-Matching Technique (MMT). The solution of the nonlinear differential equation for the director field is obtained with the aid of a finite difference (FD) scheme.