Modal Expansion Methods (MEMs)

• MEMs are spectral representations that express unknown fields as a superposition of modes derived from problem-specific operators.
• They rely on an expansion ansatz, spectral problem formulation, and coefficient-extraction rules using orthogonal or bi-orthogonal modal bases.
• MEMs enable reduced-order modeling and fast time-domain reconstruction for applications in elastic scattering, PEC structures, and open resonators.

Searching arXiv for the specified MEM-related papers and closely related terminology. arxiv_search(query="(Chen et al., 2022) OR (Lian, 2018) OR (Vial et al., 2013) modal expansion methods quasimodal expansion", max_results=10, sort_by="relevance") Retrieving bibliographic matches for the elastic, electromagnetic-power-based, and quasimodal expansion papers. arxiv_search(query="id:(Chen et al., 2022)", max_results=5, sort_by="relevance") arxiv_search(query="id:(Lian, 2018)", max_results=5, sort_by="relevance") arxiv_search(query="id:(Vial et al., 2013)", max_results=5, sort_by="relevance") Modal Expansion Methods (MEMs) are spectral representations in which an unknown field or current is written as a superposition of modes associated with an operator adapted to the governing wave problem. In the literature considered here, MEMs appear in elastodynamic scattering by metamaterial quasiparticles, in electromagnetic scattering by perfect electric conductors (PECs), and in open electromagnetic resonators treated by perfectly matched layers (PMLs). Across these settings, the common structure is an expansion ansatz, a spectral problem that defines the modal basis, and a coefficient-extraction rule; when the operator is non-self-adjoint, the expansion is built on a bi-orthogonal direct–adjoint pair rather than on a single orthogonal basis (Chen et al., 2022, Lian, 2018, Vial et al., 2013).

1. Conceptual framework

MEMs organize wave solutions into modal contributions that are intended to separate geometry, material response, and excitation dependence. In the PEC formulation, the total induced surface-current density is expanded as

$$\mathbf J(\mathbf r)=\sum_n a_n\,\mathbf J_n(\mathbf r),$$

or, in matrix form, $\mathbf J=\mathbf \Phi\,\mathbf a$, where $\mathbf \Phi$ is the modal-shape matrix. In the elastodynamic formulation, the scattered field is represented through a layer-potential density that is itself expanded in a modal basis. In the quasimodal formulation for open structures, the diffracted field is expanded over discrete quasimodes together with radiation-mode contributions (Lian, 2018, Chen et al., 2022, Vial et al., 2013).

Setting	Expanded quantity	Spectral structure
Elastic scattering by a metamaterial quasiparticle	Boundary density and scattered field	Static Neumann–Poincaré basis and polariton poles
PEC scattering	Induced surface current	Characteristic modes of the MoM impedance matrix
Open 2-D electromagnetic structures	Diffracted field	Direct–adjoint quasimodes with PMLs

This suggests that MEMs are less a single algorithm than a class of operator-based reductions. The precise meaning of “mode” depends on the problem: an eigenfunction of a static Neumann–Poincaré operator, a characteristic current of a Hermitian generalized eigenproblem, or a leaky quasimode of a non-self-adjoint PML-transformed operator. What remains invariant is the idea that the physically observed solution is reconstructed from spectrally distinguished building blocks.

2. Operator formulations and modal bases

For time-domain elastic scattering, the construction begins after a Fourier transform in $t$. The scattered field $u^{(s)}(\omega,x)$ satisfies, for $x\notin D$,

$\mathcal L_{\lambda,\mu}u+\omega^2u=0,$

together with transmission and radiation conditions at $\partial D$. It is represented by a single-layer potential

$$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$

where $\Gamma^\omega$ is the Kupradze matrix-fundamental solution. The unknown density solves a boundary system $\mathbf J=\mathbf \Phi\,\mathbf a$0, where $\mathbf J=\mathbf \Phi\,\mathbf a$1 involves the frequency-dependent Neumann–Poincaré operator $\mathbf J=\mathbf \Phi\,\mathbf a$2. In the static limit, one introduces

$\mathbf J=\mathbf \Phi\,\mathbf a$3

and, on a sphere $\mathbf J=\mathbf \Phi\,\mathbf a$4, $\mathbf J=\mathbf \Phi\,\mathbf a$5 is compact and self-adjoint with eigenpairs indexed by the three families $\mathbf J=\mathbf \Phi\,\mathbf a$6 (Chen et al., 2022).

For PEC objects, the operator framework is the Method of Moments (MoM) impedance matrix

$\mathbf J=\mathbf \Phi\,\mathbf a$7

The associated power-based modal basis is commonly chosen through Characteristic-Mode Theory (CMT), which solves

$\mathbf J=\mathbf \Phi\,\mathbf a$8

These mode currents satisfy

$\mathbf J=\mathbf \Phi\,\mathbf a$9

With the normalization $\mathbf \Phi$0 and $\mathbf \Phi$1, the basis simultaneously diagonalizes radiated-power and stored-power bilinear forms in the sense used by the electromagnetic-power-based MEM (Lian, 2018).

For open two-dimensional resonators, the spectral problem is modified by PMLs. The direct problem is written as

$\mathbf \Phi$2

and the adjoint problem as

$\mathbf \Phi$3

Because the PML coordinate stretch is non-real, $\mathbf \Phi$4, so classical modal expansion fails. The spectrum is then complex and contains discrete quasi-normal modes $\mathbf \Phi$5 with $\mathbf \Phi$6, together with a discretized approximation of the continuous radiation spectrum, the so-called PML or Bérenger modes (Vial et al., 2013).

3. Modal coefficients, orthogonality, and poles

In the elastic construction, the density is expanded in the static Neumann–Poincaré basis,

$\mathbf \Phi$7

and orthogonality of the basis yields scalar equations for the modal coefficients

$\mathbf \Phi$8

The scattered field then takes the form

$\mathbf \Phi$9

This construction already exhibits the characteristic MEM pattern: modal numerators encode excitation overlap, while modal denominators encode spectral response (Chen et al., 2022).

In the PEC formulation, the complex power delivered by the total electric field to the current is

$t$0

with

$t$1

If the modes are normalized by $t$2, then the projection of any excitation-solution onto mode $t$3 is

$t$4

In the algorithmic workflow, after solving $t$5, one expands the physical solution as

$t$6

Here, orthogonality is not merely algebraic convenience; it ties each coefficient directly to radiated-power weighting (Lian, 2018).

In QMEM, the direct and adjoint eigenmodes satisfy the bi-orthogonality relation

$t$7

with

$t$8

For the scattering problem $t$9, the diffracted field is expanded as

$u^{(s)}(\omega,x)$0

and the discrete coefficients are

$u^{(s)}(\omega,x)$1

Each resonance $u^{(s)}(\omega,x)$2 therefore appears as a simple pole of $u^{(s)}(\omega,x)$3 (Vial et al., 2013).

4. Classification, subspaces, and common misconceptions

The electromagnetic-power-based MEM gives a detailed modal classification. For a mode $u^{(s)}(\omega,x)$4, one defines

$u^{(s)}(\omega,x)$5

A mode is radiative if $u^{(s)}(\omega,x)$6 and non-radiative if $u^{(s)}(\omega,x)$7, equivalently $u^{(s)}(\omega,x)$8. It is capacitive if $u^{(s)}(\omega,x)$9, resonant if $x\notin D$0, and inductive if $x\notin D$1. The distinction between intrinsic and non-intrinsic resonance is sharper: intrinsic resonance means $x\notin D$2, whereas non-intrinsic resonance means $x\notin D$3 but $x\notin D$4. The resulting subclasses include non-radiative intrinsically resonant, radiative intrinsically resonant, radiative non-intrinsically resonant, radiative intrinsically capacitive, radiative non-intrinsically capacitive, radiative intrinsically inductive, and radiative non-intrinsically inductive modes (Lian, 2018).

A central structural result is that only some of these classes form linear spaces. The set of all non-radiative modes is the null-space of $x\notin D$5 and therefore a linear subspace. The set of all intrinsically resonant modes is the null-space of $x\notin D$6 and is likewise a linear subspace. Their intersection, the set of non-radiative intrinsically resonant modes, is also a linear subspace. By contrast, neither the set of all resonant modes defined by $x\notin D$7 nor the sets of all capacitive or all inductive modes are closed under addition. The same source also states that, by including the mode $x\notin D$8 into the intrinsically capacitive mode set and the intrinsically inductive mode set, these two modal sets become linear spaces respectively (Lian, 2018).

This directly addresses a recurrent misconception: a physically meaningful modal label does not automatically define a vector space. Another recurrent misconception appears in open resonators, where it might be tempting to use ordinary orthogonal modes despite leakage; the PML-transformed operator is explicitly non-self-adjoint, and the quasimodal expansion therefore requires the adjoint problem and bi-orthogonality rather than a classical orthogonal decomposition (Vial et al., 2013).

5. Static limits, perturbation theory, and time-domain reconstruction

The elastic MEM is notable for connecting a time-harmonic expansion to a time-domain approximation with explicit control of small parameters. In the static validation step,

$x\notin D$9

and Taylor expansions of layer potentials show that $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$0 converges in norm to $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$1. Since $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$2 is diagonalizable in the $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$3-basis, the truncated modal series

$\mathcal L_{\lambda,\mu}u+\omega^2u=0,$4

converges to the exact $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$5 uniformly for small $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$6, and as $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$7 one recovers the static expansion

$\mathcal L_{\lambda,\mu}u+\omega^2u=0,$8

with $\mathcal L_{\lambda,\mu}u+\omega^2u=0,$9-convergence of the boundary densities (Chen et al., 2022).

For finite-mode perturbative approximation, one fixes $\partial D$0 and writes

$\partial D$1

with $\partial D$2 small in $\partial D$3 if $\partial D$4 is smooth on $\partial D$5. The operator admits the perturbative form

$\partial D$6

hence

$\partial D$7

Accordingly,

$\partial D$8

and

$\partial D$9

The source describes this as yielding a fast convergent finite-mode approximation for $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$0 (Chen et al., 2022).

The same framework identifies polariton resonances. A static polariton resonance occurs when $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$1. More generally, the full denominator $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$2 has simple zeros $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$3 solving

$$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$4

with

$$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$5

Each $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$6 is a simple pole of the scattering operator, and near $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$7 the corresponding modal coefficient has the resonant form $$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$8 (Chen et al., 2022).

Time-domain reconstruction proceeds by inverse Fourier transform of the low-frequency part,

$$u^{(s)}(x)=S_D^\omega[\phi](x),\qquad S_D^\omega[\phi](x)=\int_{\partial D}\Gamma^\omega(x,y)\,\phi(y)\,d\sigma(y),$$9

Because each modal term factorizes, the inversion yields convolution in time:

$\Gamma^\omega$0

By closing the contour into the lower half-plane one picks up the pole contributions and, for $\Gamma^\omega$1,

$\Gamma^\omega$2

The low-frequency scattered field can therefore be written as a sum of damped resonant exponentials plus a remainder. The sharp error estimate of Theorem 4.10 states that, if $\Gamma^\omega$3 retains modes $\Gamma^\omega$4 and frequencies $\Gamma^\omega$5, then for $\Gamma^\omega$6,

$\Gamma^\omega$7

for any large integer $\Gamma^\omega$8, with $\Gamma^\omega$9 a compact observation region. The constants depend on the $\mathbf J=\mathbf \Phi\,\mathbf a$00–Sobolev regularity of the incident pulse, the retained number of modes $\mathbf J=\mathbf \Phi\,\mathbf a$01, and the bandwidth $\mathbf J=\mathbf \Phi\,\mathbf a$02 (Chen et al., 2022).

The reduced-order viewpoint also appears in QMEM. There, the continuous sum/integral is replaced in practice by

$\mathbf J=\mathbf \Phi\,\mathbf a$03

Once the eigenpairs $\mathbf J=\mathbf \Phi\,\mathbf a$04 are computed by a single FEM spectral solve with PMLs, evaluating $\mathbf J=\mathbf \Phi\,\mathbf a$05 at many frequencies and angles reduces to inner products and algebraic formulas. The source further states that QMEM converges rapidly near resonances, since the nearest poles dominate (Vial et al., 2013).

6. Representative applications, performance, and limitations

The electromagnetic-power-based MEM is illustrated on canonical PEC scatterers. For a PEC sphere of radius $\mathbf J=\mathbf \Phi\,\mathbf a$06 with a mesh of approximately $\mathbf J=\mathbf \Phi\,\mathbf a$07 RWG elements, at $\mathbf J=\mathbf \Phi\,\mathbf a$08 the method identifies a $\mathbf J=\mathbf \Phi\,\mathbf a$09-fold degenerate set of non-radiative modes with $\mathbf J=\mathbf \Phi\,\mathbf a$10 and $\mathbf J=\mathbf \Phi\,\mathbf a$11 that exactly match the analytic $\mathbf J=\mathbf \Phi\,\mathbf a$12 internal-resonant eigenmodes, together with a second $\mathbf J=\mathbf \Phi\,\mathbf a$13-fold set of radiative intrinsically resonant modes with the same $\mathbf J=\mathbf \Phi\,\mathbf a$14 but $\mathbf J=\mathbf \Phi\,\mathbf a$15. For a PEC circular cylinder $\mathbf J=\mathbf \Phi\,\mathbf a$16, at $\mathbf J=\mathbf \Phi\,\mathbf a$17 there is one non-radiative mode, $\mathbf J=\mathbf \Phi\,\mathbf a$18, with $\mathbf J=\mathbf \Phi\,\mathbf a$19 and $\mathbf J=\mathbf \Phi\,\mathbf a$20, while at $\mathbf J=\mathbf \Phi\,\mathbf a$21 there is one radiative intrinsically resonant mode. The same study reports that away from resonant frequencies typically only a handful of radiative capacitive modes carry most of the scattered power, whereas near resonances one “turns on” one or two intrinsically resonant modes to capture the large reactive-field portions (Lian, 2018).

The quasimodal approach is illustrated on open dielectric and periodic structures. For a triangular dielectric rod in vacuum, $\mathbf J=\mathbf \Phi\,\mathbf a$22 eigenfrequencies are computed in the complex plane; the two principal leaky modes have $\mathbf J=\mathbf \Phi\,\mathbf a$23 with $\mathbf J=\mathbf \Phi\,\mathbf a$24 and $\mathbf J=\mathbf \Phi\,\mathbf a$25 with $\mathbf J=\mathbf \Phi\,\mathbf a$26. The modal coupling maps $\mathbf J=\mathbf \Phi\,\mathbf a$27 show strong angular dependence, and QMEM reconstruction with as few as $\mathbf J=\mathbf \Phi\,\mathbf a$28–$\mathbf J=\mathbf \Phi\,\mathbf a$29 modes yields sub-percent accuracy with integrated relative error smaller than $\mathbf J=\mathbf \Phi\,\mathbf a$30. For LDOS, evaluation on a $\mathbf J=\mathbf \Phi\,\mathbf a$31 grid requires only the single spectral solve plus rapid modal summations instead of $\mathbf J=\mathbf \Phi\,\mathbf a$32 separate FEM source solves. For a mono-periodic Ge-slit grating on a ZnS substrate, $\mathbf J=\mathbf \Phi\,\mathbf a$33 eigenmodes are computed for $\mathbf J=\mathbf \Phi\,\mathbf a$34 and $\mathbf J=\mathbf \Phi\,\mathbf a$35; six leaky modes in the atmospheric-IR window $\mathbf J=\mathbf \Phi\,\mathbf a$36–$\mathbf J=\mathbf \Phi\,\mathbf a$37 have quality factors ranging from a few to $\mathbf J=\mathbf \Phi\,\mathbf a$38, and QMEM-reconstructed Fresnel amplitudes agree with full-wave FEM to absolute errors smaller than $\mathbf J=\mathbf \Phi\,\mathbf a$39 and integrated field error smaller than $\mathbf J=\mathbf \Phi\,\mathbf a$40 when all $\mathbf J=\mathbf \Phi\,\mathbf a$41 modes are used to represent the continuum well (Vial et al., 2013).

The elastic MEM is specialized to a small metamaterial quasiparticle $\mathbf J=\mathbf \Phi\,\mathbf a$42 of radius $\mathbf J=\mathbf \Phi\,\mathbf a$43 in $\mathbf J=\mathbf \Phi\,\mathbf a$44, and its final justification is explicitly tied to the subwavelength, low-frequency regime. The source states that the low-frequency part of the time-domain scattered field can be computed by a small number of damped resonant exponentials, the polariton modes, plus a rapidly decaying remainder. This indicates both the reach and the limitation of the construction: it provides explicit, sharp error controls in $\mathbf J=\mathbf \Phi\,\mathbf a$45, $\mathbf J=\mathbf \Phi\,\mathbf a$46, and $\mathbf J=\mathbf \Phi\,\mathbf a$47, but the approximation target is the low-frequency component $\mathbf J=\mathbf \Phi\,\mathbf a$48 rather than an unrestricted full-band representation (Chen et al., 2022).

Taken together, these formulations show that MEMs serve several distinct but related purposes: spectral interpretation, model reduction, resonance identification, and physically structured decomposition. The differences among them are not superficial. In closed or effectively self-adjoint settings, orthogonality with respect to a Hermitian form can be exploited directly; in open systems, bi-orthogonality replaces ordinary orthogonality; and in perturbative time-domain scattering, the modal picture is coupled to residue calculus and explicit truncation estimates. This suggests that “modal expansion” is best understood as a family of operator-specific representations rather than as a single universally transferable recipe.