---
title: Dressed Global Polarizability Matrix (dGPM)
url: https://www.emergentmind.com/topics/dressed-global-polarizability-matrix-dgpm
type: topic
---

# Dressed Global Polarizability Matrix (dGPM)

Searching arXiv for the primary dGPM paper and its GPM precursor to ground the article in the relevant literature.
The dressed Global Polarizability Matrix (dGPM) is a reduced-order electromagnetic framework for modeling complex meta-atoms in heterogeneous photonic environments, especially when scatterers touch, straddle, or are embedded within material interfaces. In this formulation, a complex scatterer is replaced by a compact operator acting on a small set of numerical electric and magnetic dipoles, while interface-mediated feedback is absorbed directly into the operator through learning in the stratified environment rather than being treated only during the subsequent many-body solve. The framework was introduced as an extension of the Global Polarizability Matrix (GPM), which had been formulated for dense ensembles of non-spherical particles in stratified media but calibrated in a homogeneous background [2606.21302; 1907.12823].

## 1. Definition and conceptual scope

The dGPM extends the original GPM to heterogeneous photonic environments in which meta-atoms interact strongly with nearby planar interfaces. In both GPM and dGPM, a complex scatterer is represented by a compact operator that maps a small basis of numerical electric and magnetic dipoles to effective induced moments, from which the scattered field is reconstructed through appropriate Green tensors [2606.21302]. The essential distinction is that the GPM is learned in a homogeneous background and can be freely translated and rotated in three dimensions, whereas the dGPM is learned directly in the stratified environment and therefore incorporates interface-mediated feedback into the operator itself [2606.21302].

This distinction is operational rather than merely terminological. In the GPM formalism, environmental effects are included during multiple-scattering simulation through the dyadic Green tensor of the stratified medium and by subtracting the background contribution already encoded in the learned operator [1907.12823]. In the dGPM formalism, by contrast, direct, reflected, transmitted, and evanescent contributions from nearby planar interfaces are absorbed into the learned operator during training, so that self-interactions mediated by the interface are transferred from the coupling operator into the dressed scatterer description [2606.21302].

A plausible implication is that the dGPM should be understood as a renormalized scatterer representation specialized to layered environments. The paper makes this intuition explicit by noting that one may think heuristically of a renormalization of the bare polarizability by interface-induced self-energy, although the dressed operator is learned directly rather than assembled from an explicit self-energy term [2606.21302].

## 2. Operator structure and electromagnetic representation

In the dGPM, each meta-atom is represented by $N_p$ numerical dipoles, each carrying six moment channels: three electric and three magnetic. The induced-moment vector at numerical dipole $n$ is
$$
\mathbf{P}_n = [p_x,p_y,p_z,m_x,m_y,m_z]^\top,
$$
and the corresponding exciting field vector is
$$
\mathbf{F}_n = [E_x,E_y,E_z,H_x,H_y,H_z]^\top
$$
[2606.21302].

The single-particle dGPM is a $6N_p \times 6N_p$ matrix $\tilde{\mathbf{A}}$ organized into $6\times 6$ blocks $\tilde{\mathbf{A}}_{m\to n}$:
$$
\mathbf{P}=\tilde{\mathbf{A}}\,\mathbf{F},
$$
with
$$
\tilde{\mathbf{A}}_{m\to n}\in\mathbb{C}^{6\times 6}
$$
[2606.21302]. Unlike the discrete dipole approximation, these numerical dipoles are not local physical dipoles obeying local constitutive relations. The operator is generally fully populated, and its off-diagonal blocks encode nonlocal internal couplings inside the scatterer [2606.21302]. This nonlocality is central to both GPM and dGPM: it is what allows a small number of numerical dipoles to reproduce the external electromagnetic response of complex, high-index, high-aspect-ratio meta-atoms [2606.21302; 1907.12823].

Higher-order multipolar physics is not introduced through an explicit quadrupolar or higher multipole basis. Instead, it is captured implicitly through the nonlocal coupling between numerical sites and through the layered Green tensors [2606.21302]. This suggests that the compactness of the representation does not derive from truncating the physics to dipolar order in the usual sense, but from compressing the scatterer response into an operator basis calibrated against full-wave data.

## 3. Learning in a stratified environment

The dressed operator is retrieved from full-wave data in the stratified medium using source and learning surfaces. The construction employs a source vector $\mathbf{s}\in\mathbb{C}^{6N_s}$, a stratified Green matrix $\mathbf{G}_{\mathrm{sd}}$ mapping the source surface to the numerical dipoles, and a stratified Green matrix $\mathbf{G}_{\mathrm{dl}}$ mapping the induced moments at the numerical dipoles to the learning surface [2606.21302]. The dGPM-predicted scattered field on the learning surface is
$$
\mathbf{w}_{\mathrm{dGPM}}=\mathbf{G}_{\mathrm{dl}}\,\tilde{\mathbf{A}}\,\mathbf{G}_{\mathrm{sd}}\,\mathbf{s},
$$
and the dressed operator is retrieved through the regularized least-squares problem
$$
\min_{\tilde{\mathbf{A}}}\left\lVert \mathbf{G}_{\mathrm{dl}}\,\tilde{\mathbf{A}}\,\mathbf{G}_{\mathrm{sd}}\,\mathbf{s}-\mathbf{w}_{\mathrm{ref}}\right\rVert^2
$$
[2606.21302].

The precomputation pipeline reported for a single meta-atom consists of three steps: generating the reference dataset with COMSOL using random electric and magnetic dipolar sources on a source surface, building the stratified Green matrices $\mathbf{G}_{\mathrm{sd}}$ and $\mathbf{G}_{\mathrm{dl}}$, and solving the inverse problem for $\tilde{\mathbf{A}}$ [2606.21302]. Typical runtime for $N_p\le 25$ is reported as approximately $40$–$64$ minutes with fine meshes ensuring negligible discretization error [2606.21302].

For the ZnS-capped GaAs pillar at $\lambda=10.6\,\mu\mathrm{m}$, the reported training configuration used $N_p=15$, $N_s=30$ sources with six polarizations each, and approximately $N_\ell\simeq 1500$ learning points. This yielded mean relative error $\bar{E}\simeq 2.5\%$ with standard deviation $0.2\%$ on training and $\bar{E}\simeq 2.72\%$ with standard deviation $0.29\%$ on generalization [2606.21302]. The error metric is given explicitly by
$$
E(\mathbf{x},\mathbf{x}')=\frac{1}{N}\sum_{k=1}^{N}\frac{|x_k-x_k'|}{|x_k|+|x_k'|},
$$
with variance
$$
\sigma^2(\mathbf{x},\mathbf{x}')=\frac{1}{N-1}\sum_{k=1}^{N}(E_k-\bar{E})^2
$$
[2606.21302].

The 2019 GPM precursor used an inverse procedure with two-stage least-squares inversion based on SVD pseudoinverses, learning the operator from scattered fields sampled on a learning surface in a homogeneous background [1907.12823]. That earlier work did not explicitly coin the term dGPM, but it established the compact nonlocal operator paradigm that the 2026 work extends to interface-coupled training [1907.12823].

## 4. Multiple scattering, Green operators, and formal dressing

For a collection of meta-atoms, the many-body system is assembled through a block-diagonal operator
$$
\boldsymbol{\Alpha}^{\mathrm{d}}=\mathrm{diag}\big(\tilde{\mathbf{A}}^{(1)},\tilde{\mathbf{A}}^{(2)},\dots,\tilde{\mathbf{A}}^{(M)}\big),
$$
where each block may be either a dGPM for an interface-coupled scatterer or a GPM for a free-space scatterer [2606.21302]. The self-consistent multiple-scattering equation is
$$
\mathbf{P}=\boldsymbol{\Alpha}^{\mathrm{d}}\big(\mathbf{F}^{\mathrm{inc}}+\mathbf{G}\mathbf{P}\big),
$$
or equivalently
$$
\big(\mathbf{I}-\boldsymbol{\Alpha}^{\mathrm{d}}\mathbf{G}\big)\mathbf{P}
=
\boldsymbol{\Alpha}^{\mathrm{d}}\mathbf{F}^{\mathrm{inc}}
$$
[2606.21302].

The inter-particle coupling operator is decomposed as
$$
\mathbf{G}=\mathbf{G}^{(0)}+\mathbf{G}^{(\mathrm{int})},
$$
where $\mathbf{G}^{(0)}$ contains the direct contributions within a homogeneous layer and $\mathbf{G}^{(\mathrm{int})}$ contains reflected, transmitted, and evanescent couplings [2606.21302]. For planar interfaces, the interface-mediated Green tensor may be written in spectral form using Sommerfeld integrals over in-plane wavevectors and an $s/p$ (TE/TM) decomposition [2606.21302]. In practice, the source-to-dipole, dipole-to-learning-surface, and many-body Green operators are evaluated with the stratified-medium formalism of Ref. [24] in the 2026 paper, so that all direct, reflected, transmitted, and evanescent components are included [2606.21302].

The conceptual dressing relation is written as
$$
\boldsymbol{\Alpha}^{\mathrm{d}}=(\boldsymbol{\Alpha}_0^{-1}-\boldsymbol{\Sigma}_{\mathrm{int}})^{-1},
$$
where $\boldsymbol{\Sigma}_{\mathrm{int}}$ denotes the interface-mediated self-interaction [2606.21302]. The paper emphasizes, however, that the dressed operator is not computed from an explicit $\boldsymbol{\Sigma}_{\mathrm{int}}$; instead, the effect of $\boldsymbol{\Sigma}_{\mathrm{int}}$ is encoded implicitly through the stratified Green operators used during training [2606.21302].

This point connects directly to the 2019 GPM formalism. There, for a single particle in a stratified medium, the environmental contribution appears through
$$
\mathbf{\Psi}_{\mathrm{exc}}=\mathbf{\Psi}_b+\Delta\mathcal{G}\,A\,\mathbf{\Psi}_{\mathrm{exc}},
$$
with the equivalent dressed form
$$
A_{\mathrm{dressed}}=A(I-\Delta G\,A)^{-1}
$$
derived as an operational consequence of the self-consistent equation [1907.12823]. The 2026 dGPM may therefore be viewed as an explicit learning-based realization of dressing that had been implicit in the earlier GPM-plus-environment formulation.

## 5. Observables, benchmarks, and computational performance

Once the induced moments are known, scattered near fields at a point $\mathbf{r}$ are reconstructed by applying the stratified dyadic Green tensors:
$$
\mathbf{E}_{\mathrm{sca}}(\mathbf{r})
=
\sum_n \mathbf{G}^{EE}(\mathbf{r},\mathbf{r}_n)\mathbf{p}_n
+
\sum_n \mathbf{G}^{EM}(\mathbf{r},\mathbf{r}_n)\mathbf{m}_n,
$$
$$
\mathbf{H}_{\mathrm{sca}}(\mathbf{r})
=
\sum_n \mathbf{G}^{ME}(\mathbf{r},\mathbf{r}_n)\mathbf{p}_n
+
\sum_n \mathbf{G}^{MM}(\mathbf{r},\mathbf{r}_n)\mathbf{m}_n
$$
[2606.21302]. The framework reconstructs fields anywhere outside the learning surfaces; inside the learning surface, point-dipole singularities appear near numerical dipoles, but this does not affect reconstruction outside [2606.21302].

Far fields are obtained through a near-to-far transformation for layered media, with RETOP given as an example [2606.21302]. For the ZnS-capped GaAs pillar at $\lambda=10.6\,\mu\mathrm{m}$, angular radiation diagrams match COMSOL to within $0.68\%$ integrated relative error and reproduce features including critical-angle cutoffs and evanescent-to-propagating conversion [2606.21302]. Reflection and transmission cross sections $\sigma_R,\sigma_T$ are computed by normalizing radiated powers by the incident flux, while absorption for lossy meta-atoms is evaluated by energy balance from $\sigma_{\mathrm{ext}}-\sigma_{\mathrm{sca}}$, thereby avoiding volumetric integration [2606.21302].

The published benchmarks span single particles, dimers, and hybrid ensembles. For a single interface-straddling pillar, the reported training error is $\bar{E}_{\mathrm{train}}=2.47\%$ with standard deviation $0.27\%$, and generalization error is $\bar{E}_{\mathrm{gen}}=2.72\%$ with standard deviation $0.29\%$; near-field component errors are reported as low as $1.47\%$ across all six components, and the far-field integrated error is $0.68\%$ [2606.21302]. For a dimer at separation $d=2\,\mu\mathrm{m}$ and incidence angle $\theta=30^\circ$, the mean relative error in near fields is reported as $1.82\%$ for dominant $E_x$ and $2.43\%$ across six components, with $2.3\%$ far-field error [2606.21302]. Dimer cross-section benchmarks show mean relative errors below $1.8\%$ for separations $d\ge 0.8\,\mu\mathrm{m}$, rising to $0.9\%/2.2\%$ at $d=0.6\,\mu\mathrm{m}$ and $7.4\%/3.6\%$ at $d=0.4\,\mu\mathrm{m}$ [2606.21302].

The reported speed and memory advantages over full-wave simulation are substantial. For the dimer benchmark, COMSOL requires approximately $12$ minutes and $125$ MB per angle, while dGPM requires approximately $1$ s per angle after approximately $9$ s initialization with approximately $0.5$ MB, corresponding to approximately $700\times$ speed and approximately $250\times$ memory reduction [2606.21302]. In the precursor GPM paper, a dense ensemble of $16$ cylinders in a metallo-dielectric stack required approximately $78$ GB for the COMSOL reference and approximately $14.4$ MB for the MATLAB GPM solution [1907.12823]. These figures establish the method as a reduced-order surrogate for repeated scattering calculations rather than a direct replacement for full-wave training.

## 6. Comparison with T-matrix and other scattering formalisms

A central argument for dGPM is that conventional T-matrix methods rely on vector spherical wave expansions outside circumscribing spheres. For high-aspect-ratio meta-atoms such as micropillars, the circumscribing spheres can be much larger than the physical scatterers, and at practical separations the spheres may overlap, causing loss of convergence and preventing near-gap field resolution [2606.21302]. Standard T-matrix plus S-matrix combinations also cannot handle scatterers straddling interfaces [2606.21302].

The dGPM addresses these issues in three ways stated explicitly in the 2026 paper: it uses conformal learning surfaces, allowing close packing without circumscribing-sphere overlap; it learns directly in layered media, embedding reflection, transmission, and evanescent feedback in the operator; and it preserves reduced-order compactness and portability for translations parallel to the interface and rotations around the interface normal while remaining compatible with many-body Green-tensor multiple scattering [2606.21302].

The GPM literature had already positioned the operator approach relative to discrete dipole approximations, T-matrix methods, and BEM/FEM. The GPM uses a small number of non-locally coupled numerical dipoles with $36d^2$ degrees of freedom rather than local dipoles with $6d$ degrees of freedom, and thereby can achieve accurate scattering predictions for high-index, resonant, non-spherical particles with dramatically fewer unknowns than DDA [1907.12823]. In the 2019 benchmarks, a resonant Si cylinder represented by $N_d=10$ numerical dipoles achieved approximately $1.05\%$ error in scattering cross section with $0.27$ s runtime and $1.91$ MB, whereas DDSCAT with $107{,}400$ dipoles achieved $2.89\%$ error with $14.8$ s runtime and $168$ MB [1907.12823]. The 2019 paper also reports that T-matrix methods fail badly for two cylinders with strong near-field coupling and touching learning surfaces [1907.12823].

A plausible implication is that dGPM is best situated between full-wave solvers and traditional multipole formalisms: more specialized than general FEM or BEM, but more robust than circumscribing-sphere-based T-matrix schemes for interface-coupled, high-aspect-ratio, and densely packed structures.

## 7. Applications, assumptions, and extensions

The framework is intended for rapid parameter sweeps and large-scale design of ordered and disordered metasurfaces. The 2026 paper states that one may precompute operators once, possibly in shared libraries, assemble heterogeneous scenes mixing buried, suspended, and interfacial meta-atoms, perform angle and frequency sweeps, and evaluate scattering, reflection, transmission, absorption, and radiation patterns efficiently [2606.21302]. It further states that the operator-based solver can be integrated with automatic differentiation, with TorchGDM given as an example, to support inverse-design pipelines and GPU-accelerated optimization [2606.21302].

The formal assumptions are linear, time-harmonic electromagnetism, layered planar environments, and material dispersion as represented in the full-wave training stage [2606.21302]. Formal validity additionally requires non-overlap of learning surfaces across different meta-atoms [2606.21302]. Several limitations are identified explicitly. Extreme near-touching configurations with significant learning-surface overlap degrade accuracy; very high-order multipolar content may require increasing $N_p$ or refining dipole placement; strongly nonlocal or spatially dispersive materials are not explicitly modeled; and symmetry portability is reduced relative to GPM, since only translations parallel to the interface and rotations about the interface normal are retained [2606.21302].

The implementation scaling is also spelled out. Assembly of the coupling operator scales approximately as $\mathcal{O}(N^2)$ in the number of meta-atoms because of pairwise coupling, dense direct solves scale approximately as $\mathcal{O}(N^3)$, and memory scales approximately as $\mathcal{O}(N^2)$ [2606.21302]. In ordered arrays, identical Green blocks can be reused, making initialization cheaper and leaving the dense solve dominant from about $N\sim 150$; in disordered arrays, initialization dominates until about $N\sim 10^3$ [2606.21302]. In the reported MATLAB implementation without GPU acceleration, systems with hundreds to thousands of meta-atoms are tractable [2606.21302].

The extensions proposed in the 2026 paper include nonlinearities via augmented learned operators, time modulation, anisotropic substrates, multilayer stacks already supported by stratified Green tensors, and non-planar interfaces provided suitable Green operators are available [2606.21302]. Straightforward directions for improved scalability are also enumerated: iterative solvers, preconditioning, GPU implementations, and hierarchical low-rank compression of the coupling matrix $\mathbf{G}$ [2606.21302].

Taken together, these features define the dGPM as a compact operator formalism in which interface-mediated self-interactions are learned rather than appended, enabling accurate and efficient simulation of large collections of interface-coupled meta-atoms in layered media [2606.21302]. The broader significance lies in extending operator-based electromagnetic modeling from homogeneous backgrounds to heterogeneous photonic environments while preserving the compact scatterer-based representation established by the GPM methodology [2606.21302; 1907.12823].

Source: https://www.emergentmind.com/topics/dressed-global-polarizability-matrix-dgpm