Perfectly matched layers (PMLs) were introduced as a better alternative to absorbing boundary conditions (ABCs) by Berenger[1]} to achieve a higher accuracy in domain truncation by eliminating undesired numerical reflections from boundaries, leading to exponential convergence of the numerical solution to the exact solution [2]}. PML implementation is done by stretching Cartesian coordinates such that the stretching is defined in a layer surrounding \(\Omega \) , as in [3]}, giving a larger computational domain \(\Omega _\text{PML}\) . In this regard, we assume the boundaries of the artificially truncated domain \(\Omega \)  are aligned with the coordinate axes.
