  This paper proposes a new boundary integral equation (BIE) methodology based
on the perfectly matched layer (PML) truncation technique for solving the
electromagnetic scattering problems in a multi-layered medium. Instead of using
the original PML stretched fields, artificial fields which are also equivalent
to the solutions in the physical region are introduced. This significantly
simplifies the study of the proposed methodology to derive the PML problem.
Then some PML transformed layer potentials and the associated boundary integral
operators (BIOs) are defined and the corresponding jump relations are shown.
Under the assumption that the fields vanish on the PML boundary, the solution
representations, as well as the related BIEs and regularization of the
hyper-singular operators, in terms of the current density functions on the
truncated interface, are derived. Numerical experiments are presented to
demonstrate the efficiency and accuracy of the method.
