Cartesian perfectly matched layers for the Helmholtz equation with variable propagation speed: well-posedness and exponential convergence

We analyse the Cartesian perfectly matched layer (PML) approximation of the two-dimensional Helmholtz equation $-c^2\Delta u-{\omega }^{2}u=f$, in which the propagation speed $c$ is variable inside the region of interest and equal to a constant $c_{\infty }$ outside a disc $B_R$. Existing analyses of the Cartesian PML assume a homogeneous, or piecewise homogeneous, background. We assume of $c$ only a two-sided bound, so that discontinuous speeds -- layered media and compact inclusions -- are covered. This is possible because the speed enters the PML sesquilinear form only through its zeroth order coefficient, the principal part being built from the stretching profiles alone. Well-posedness of the physical problem at this regularity follows from a Lippmann-Schwinger equation, uniqueness from Rellich's lemma and unique continuation. A Morawetz multiplier yields a sharp non-trapping condition on $c$, namely $\nabla c(x)\cdot (x-x_0)<c(x)$ for some $x_0$, which reduces in the radial case to the Herglotz condition $\frac{d}{d\rho }(\rho /c(\rho ))>0$; its failure produces a trapped bicharacteristic. The variable speed is invisible to the essential spectrum: after factoring out $c^2$, the two operators differ by a compactly supported multiplication, so the spectral analysis of the homogeneous Cartesian layer applies verbatim. This yields well-posedness of the truncated problem for every $\omega >0$, with a stability constant independent of the layer width, and convergence to the physical solution in $H^1$ at a rate exponential in $\omega {\gamma }_M{\delta }_M$, where ${\gamma }_M$ and ${\delta }_M$ are the absorption strength and the layer thickness. Numerical experiments with a low-velocity inclusion confirm the exponential decay in both $\delta$ and $\omega$.