This article introduces a boundary integral equation (BIE) for time-harmonic electromagnetic scattering by composite dielectric objects. The formulation extends the classical Müller equation to composite structures through the global multitrace (MT) method. The key ingredient enabling this extension is the use of the Stratton–Chu representation in complementary region, also known as the extinction property, which augments the off-diagonal blocks of the interior representation operator. The resulting block system is composed entirely of second-kind operators. A Petrov–Galerkin (mixed) discretization using Rao–Wilton–Glisson trial functions and Buffa–Christiansen (BC) test functions is employed, yielding linear systems that remain well conditioned on dense meshes and at low frequencies without the need for additional stabilization. This reduces computational costs associated with matrix-vector multiplications and iterative solving. Numerical experiments demonstrate the accuracy of the method in computing field traces and derived quantities.
