A fast Huygens sweeping method for capturing paraxial multi-color optical self-focusing in nematic liquid crystals

We propose a numerically efficient algorithm for simulating the multi-color optical self-focusing phenomena in nematic liquid crystals. The propagation of the nematicon is modeled by a parabolic wave equation coupled with a nonlinear elliptic partial differential equation governing the angle between the crystal and the direction of propagation. Numerically, the paraxial parabolic wave equation is solved by a fast Huygens sweeping method, while the nonlinear elliptic PDE is handled by the alternating direction explicit (ADE) method. The overall algorithm is shown to be numerically efficient for computing high frequency beam propagations. - Highlights: • A simple strategy can already improve the fast Huygens sweeping method for high frequency asymptotic solutions. • The Alternating Direction Explicit method can significant improve the computational efficiency. • Numerical solutions match the Fredericks transition.