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.

Multiresolution of the one dimensional free-particle propagator

High-order Symplectic FDTD Scheme for Solving a Time-Dependent Schrödinger Equation

Local Split Step Fourier Propagator for Wave Equation Depth Migration

The discrete Green's function method for wave packet expansion via the free Schrödinger equation

A Spectral Multiscale Method for Wave Propagation Analysis: Atomistic-Continuum Coupled Simulation

Numerical Quantum Propagation with Time-Dependent Hamiltonian

Asymptotic initial value representation of the solutions of semi-classical systems presenting smooth codimension one crossings

A high-order integral equation-based solver for the time-dependent Schrodinger equation

Multi-phase Computations of the Semiclassical Limit of the Schrödinger Equation and Related Problems: Whitham Vs Wigner

A Semi-Lagrangian Time Splitting Method for the Schrödinger Equation with Vector Potentials

Study of multi-dimensional problems arising in wave propagation using a hybrid scheme

A vector Helmholtz electromagnetic wave propagator for inhomogeneous media

A pseudo-spectral multiscale method: Interfacial conditions and coarse grid equations

Polynomial propagators for classical molecular dynamics

Multiscale Simulation Based on the Meshless Local Petrov-Galerkin (mlpg) Method

Multiscale analysis and numerical algorithm for the Schrödinger equations in heterogeneous media

Fast Multiresolution Methods for Density Functional Theory in Nuclear Physics

Symplectic Hamiltonian Finite Element Methods for Semilinear Wave Propagation

Solitary Wave Propagation and Interactions for the Klein–Gordon–Zakharov Equations in Plasma Physics

Efficient hierarchical Liouville space propagator to quantum dissipative dynamics

Unified gas-kinetic wave-particle methods III: Multiscale photon transport