A perfectly matched layer approach for the spectral split-step Padé method

Abstract

The split-step-Padé (SSP) method is widely used to model wave phenomena in various applications, including radio physics, optics and acoustics. In this method, the propagator of the one-way counterpart of the Helmholtz equation is computed through its Padé approximant and a finite-difference discretization of the transverse operator. This work develops and validates numerically a spectral counterpart of the SSP method. A key challenge in practical applications is inverting the transverse operator in the presence of perfectly matched layers (PMLs), which are commonly used to truncate the computational domain. Such inversion can be accomplished using Krylov subspace methods, which converge rapidly, provided that a suitable preconditioner is used. We also study the analytical properties of the spectral SSP marching scheme under periodicity conditions in the transverse variable. We validate the newly developed spectral SSP method numerically in two realistic test scenarios from radio physics and underwater acoustics.