Abstract:The vector electric-field Helmholtz equation, containing cross-polarization terms, is factored to produce both pseudo-differential and exponential operator forms of a three-dimensional, one-way, vector, wave equation for propagation through inhomogeneous media. From this operator factorization we develop a high-order approximate, vector Helmholtz propagator that correctly handles forward-arc, high-angle scattering and diffraction from inhomogeneities at all resolved length scales, and seamlessly includes evanescent waves. A rational approximation/partial fraction decomposition of the exponential operator converts the propagator into a moderate number of large, sparse, linear solves whose results are summed together at each step to advance the electric field in space. We use a new AAA-Lawson rational interpolant for this approximation, rather than the more common Padé expansions that have appeared in the seismic and ocean acoustics literature previously. GMRES is used to solve these large systems. A direct-solve, free space propagation method proves to be an effective preconditioner for GMRES, but can also serve as a standalone propagator in homogeneous media. Scalar computational examples shown include plane wave diffraction by a circular aperture and Gaussian beam propagation through sine-product and homogeneous refractive index fields. The sine-product example compares its results to that of paraxial propagation through the same media, and demonstrates the substantial differences between these propagator paradigms when the scale of the inhomogeneities is of the order of the fundamental wavelength in the Helmholtz equation. We also examine the convergence of the homogeneous media beam results to fields generated by Clenshaw-Curtis evaluation of the first Rayleigh-Sommerfeld integral for the same initial conditions.

Analytic 4×4 Propagation Matrices for Homogeneous Media

Propagation Characteristic Analysis and Applications of Spatial Four-Waveguide System

Light Interaction with Multilayer Arbitrary Anisotropic Structure: an Explicit Analytical Solution and Application for Subwavelength Imaging

Analytical model of spatial four-waveguide coupler

Exact closed forms for the transmittance of electromagnetic waves in one-dimensional anisotropic periodic media

Analytical solution of Maxwell’s equations in lossy and optically active crystals

Explicit Analysis of Anisotropic Planar Waveguides by the Analytical Transfer-Matrix Method

Analysis Methods for Polarization State and Energy Transmission of Rays Propagating in Optical Systems

T-MATRIX ANALYSIS OF MULTIPLE SCATTERING FROM PARALLEL SEMI-CIRCULAR CHANNELS FILLED WITH CHIRAL MEDIA IN A CONDUCTING PLANE

A vector Helmholtz electromagnetic wave propagator for inhomogeneous media

Propagation Characteristics of the Two-Dimensional Off-Axial Hermite-cosh-Gaussian Beams Through Rectangular Apertured and Misaligned Optical Systems

A finite element beam propagation method for simulation of liquid crystal devices.

Full-vectorial Propagation Model and Modified Effective Mode Area of Four-Wave Mixing in Straight Waveguides.

Propagation of random electromagnetic beams through axially nonsymmetrical optical systems

Radiation Pressure on Inhomogeneous Active Chiral Structures with Propagation Matrix

Analytic Matrix Method for the Study of Propagation Characteristics of A Bent Planar Waveguide

Propagation Properties and Application of Radially Polarized Partially Coherent Twisted Vortex Hermite‐Gaussian Schell‐Model Beams

Unified approach to paraxial propagation in uniform media and media with linear or quadratic refractive index distribution

Multilayer optical calculations

Efficient simulation of bi-periodic, layered structures based on the T-matrix method