Stratton–Chu Representation: Integral Formulation in EM

• Stratton–Chu representation is a fundamental integral formula that expresses electromagnetic fields as surface or hybrid surface-volume integrals using Green’s functions.
• It converts volumetric Maxwell equations into boundary and hybrid integral equations, enabling efficient formulations for scattering and computational physics.
• The approach underpins advanced boundary integral methods and physical optics approximations, facilitating high-frequency and complex electromagnetic analyses.

The Stratton–Chu representation is a fundamental integral formula expressing electromagnetic fields as surface or mixed surface-volume integrals in terms of tangential field traces and Green’s functions. Originating from vector potential theory and Green’s theorems, it provides a key analytic and computational tool for the solution of frequency-domain Maxwell problems, particularly enabling the conversion of volumetric partial differential equations into boundary or hybrid boundary-volume integral equations. The representation underpins advanced formulations of scattering by penetrable and composite objects, supports robust boundary integral equation (BIE) methods, and serves as the foundation for efficient physical optics approximations in high-frequency and tightly focused field settings.

1. Mathematical Foundation and Representation Formulas

The Stratton–Chu representation is derived from Maxwell’s equations via vector Green identities. In a domain $\Omega \subset \mathbb{R}^3$ with Lipschitz boundary %%%%1%%%%, time-harmonic electromagnetic fields $(\mathbf{E}, \mathbf{H})$ satisfy

$$\nabla \times \mathbf{E} = i \omega \mu(\mathbf{x}) \mathbf{H}, \qquad \nabla \times \mathbf{H} = -i\omega \epsilon(\mathbf{x}) \mathbf{E},$$

with $\epsilon, \mu$ as permittivity and permeability (possibly variable in $\Omega$). The scalar and dyadic Green’s functions for wavenumber $\kappa$ are

$$G_\kappa(\mathbf{x}, \mathbf{y}) = \frac{e^{i\kappa|\mathbf{x} - \mathbf{y}|}}{4\pi|\mathbf{x} - \mathbf{y}|}, \qquad \mathbf{G}(\mathbf{x}, \mathbf{y}) = \left(I + \frac{1}{\kappa^2} \nabla_\mathbf{x} \otimes \nabla_\mathbf{x}\right) G_\kappa(\mathbf{x}, \mathbf{y})$$

(Labarca-Figueroa et al., 2024, Le et al., 20 Jan 2026).

For constant parameter domains, the classical Stratton–Chu formulas express the fields in terms of surface (layer) integrals involving tangential traces: $$\mathbf{E}(\mathbf{x}) = \nabla \times \int_\Gamma G_\kappa(\mathbf{x}, \mathbf{y})[\mathbf{n}(\mathbf{y}) \times \mathbf{E}(\mathbf{y})]\,dS_\mathbf{y} + i\omega\mu \int_\Gamma G_\kappa(\mathbf{x}, \mathbf{y})[\mathbf{n}(\mathbf{y}) \times \mathbf{H}(\mathbf{y})]\,dS_\mathbf{y} - \nabla \int_\Gamma G_\kappa(\mathbf{x}, \mathbf{y})[\mathbf{n}(\mathbf{y}) \cdot \mathbf{E}(\mathbf{y})]\,dS_\mathbf{y}$$ and a similar representation for $\mathbf{H}$ (Labarca-Figueroa et al., 2024, Le et al., 20 Jan 2026, Dumont et al., 2016).

For inhomogeneous interior media, the Stratton–Chu representation generalizes to include volume Newton potentials, yielding coupled surface-volume potentials and allowing for single-trace hybrid formulations (Labarca-Figueroa et al., 2024).

2. Surface and Volume Potentials: Classical and Generalized Forms

Classically, all integral terms in the Stratton–Chu representation are surface (boundary) integrals, with field values everywhere determined by boundary tangential traces (Cauchy data). This reduction is achieved using layer potentials—specifically, the Maxwell single-layer $\mathcal{T}_\kappa$ and double-layer $\mathcal{D}_\kappa$ operators: $$\mathcal{T}_\kappa[\beta](\mathbf{x}) = \frac{1}{\kappa^2} \text{curl curl } S_\kappa[\beta] + S_\kappa[\beta],\quad \mathcal{D}_\kappa[\alpha](\mathbf{x}) = S_\kappa[\mathbf{n} \times \alpha],$$ where $S_\kappa$ is the vector-valued single-layer potential (Labarca-Figueroa et al., 2024).

For domains with inhomogeneous coefficients, the interior Maxwell system is recast in terms of a constant-coefficient part plus volume sources dependent on contrasts $$p_e(\mathbf{x}) = 1 - \epsilon(\mathbf{x})/\epsilon_1$$, $p_m(\mathbf{x}) = 1 - \mu(\mathbf{x})/\mu_1$. The resulting representation includes both surface layer potentials (as above) and volume Newton-type potentials, e.g., for the electric field,

$$\mathbf{E}(\mathbf{x}) = -i\omega\mu_1 N_{\Omega, \kappa_1}[p_m \mathbf{H}] - \kappa_1^2 N_{\Omega, \kappa_1}[p_e \mathbf{E}] - N_{\Omega, \kappa_1}[\operatorname{div} \mathbf{E}] + \text{boundary layer potentials}$$

(Labarca-Figueroa et al., 2024).

A compact formulation is obtained using weighted volume and boundary operators, facilitating numerical schemes that combine both boundary and domain discretizations.

3. The Extinction Property and Composite Domains

A crucial property of the Stratton–Chu representation is “extinction”: for any subdomain $\Omega_k$, the potential constructed using its Cauchy data yields the physical field inside and vanishes identically outside $\overline{\Omega}_k$. Symbolically, for an operator $\mathbf{U}^{(k)}(\mathbf{x})$ formed from Cauchy data on $\Gamma_k$,

$$\mathbf{U}^{(k)}(\mathbf{x}) = \begin{cases} (\mathbf{E}(\mathbf{x}), -\mathbf{H}(\mathbf{x}))^T, & \mathbf{x} \in \Omega_k \ 0, & \mathbf{x} \notin \overline{\Omega}_k, \end{cases}$$

with tangential traces on unrelated boundaries vanishing (Le et al., 20 Jan 2026). This property is central in formulating multi-trace boundary integral equations for composite objects: “off-diagonal” contributions are precisely calibrated to cancel (via complement-region extinction), resulting in block systems composed of second-kind (identity plus compact) operators (Le et al., 20 Jan 2026).

4. Computational Implementations and Discretization

Discretization schemes leveraging the Stratton–Chu representation typically expand boundary unknowns in Rao–Wilton–Glisson (RWG) elements, suitable for $H^{-1/2}(\text{curl}_\Gamma)$, and Buffa–Christiansen elements for Petrov-Galerkin schemes, or use Nédélec edge elements for $H(\text{curl}, \Omega)$ in domain approaches (Labarca-Figueroa et al., 2024, Le et al., 20 Jan 2026). The presence of both surface and volume potentials in hybrid formulations allows for highly flexible numerical frameworks.

A major challenge is efficient evaluation of dense integral operators, which is addressed using $\mathcal{H}$- or $\mathcal{H}^2$-matrix compression for system matrix assembly. Special quadrature rules—e.g., Duffy transforms for singular integrals, high-order tensorized Gauss quadrature—are required for accurate evaluation of weakly and strongly singular layer and volume potentials (Labarca-Figueroa et al., 2024, Dumont et al., 2016). For time-dependent or broad-band sources, as in tightly focused laser pulses, the field expressions are evaluated over multiple frequencies, enabling full temporal reconstruction via inverse discrete Fourier transforms (Dumont et al., 2016).

Parallel strategies, such as embarrassingly parallel MPI algorithms assigning focal subdomains to each process, have demonstrated near-ideal strong scaling on high-performance computing platforms (Dumont et al., 2016).

5. Role in Boundary Integral Equation Formulations

The Stratton–Chu representation underpins a variety of BIE methods for frequency-domain electromagnetic scattering. In the classical Müller approach, it provides the layer potentials whose combination cancels hypersingularities and yields second-kind integral equations with favorable conditioning: $$[M]_{kk} = \tfrac12(\beta_k + \beta_0)I + \beta_k A^{(k)}_{kk} - \beta_0 A^{(0)}_{kk}, \quad [M]_{jk} = \tfrac12(\beta_k - \beta_0)I + \beta_k A^{(k)}_{jk} - \beta_0 A^{(0)}_{jk},\ j\neq k$$ (Le et al., 20 Jan 2026). The multi-trace extension introduces extinction-based off-diagonal corrections, ensuring that all blocks are identity plus compact and that the spectrum is well-clustered irrespective of mesh refinement and frequency.

For penetrable, inhomogeneous, or composite geometries, the Stratton–Chu-based coupled boundary and volume integral frameworks handle material variation directly through interior reference parameter selection and explicit volume contrast source terms (Labarca-Figueroa et al., 2024).

6. Applications in Physical Optics and High-Frequency Regimes

In high-frequency or geometrical optics, the Stratton–Chu representation is adapted (e.g., with the physical optics approximation, POA) to obtain rigorous surface-integral expressions for diffracted or reflected fields (Dumont et al., 2016). For perfectly conducting scatterers, boundary conditions and limiting procedures yield explicit formulas for fields at any observation point outside the scatterer, including regularization via the Hadamard finite part for hypersingular integrals: $$\lim_{\mathbf{r}' \to S} \iint_S \nabla\nabla G \cdot A\,dS = \text{Hadamard finite part} \int_S \nabla\nabla G \cdot A\,dS$$ (Dumont et al., 2016). In geometries such as parabolic mirrors, unexpected phase cancellations make the relevant integrals only weakly oscillatory, enabling efficient quadrature without specialized oscillatory integration (Dumont et al., 2016). These methods facilitate direct modeling of tightly focused laser pulses, predicting key features such as ellipticity, strong longitudinal field components, and record-level focal intensities.

7. Significance for Electromagnetic Scattering and Computational Physics

The Stratton–Chu representation provides a unifying mathematical construct allowing for the rigorous reformulation of Maxwell’s equations into integral form. Its extinction property permits decomposition and coupling of field representations in complex, multi-domain geometries. Modern formulations relying on this representation demonstrate robust conditioning, compatibility with high-order finite and boundary element spaces, and efficient numerical solvability even in low-frequency or high-contrast scenarios (Labarca-Figueroa et al., 2024, Le et al., 20 Jan 2026).

A plausible implication is that future advances may further generalize surface-volume hybrid representations, optimize discretizations for singular kernel integrals, and extend applicability to broad classes of frequency-domain and time-domain electromagnetic phenomena. The approach continues to be central in the theoretical analysis and practical computation of scattering, guiding the development of scalable, accurate solvers for complex electromagnetic environments.