Green-Tensor QNM Calibration Methods

• The paper introduces a novel calibration procedure that employs normalized quasinormal modes to accurately reconstruct the Green tensor for both electromagnetic and gravitational systems.
• It utilizes systematic mode-solving, biorthogonal normalization, and iterative refinement to align QNM expansions with full-wave simulations and address near-field singularities.
• The approach is significant for advancing practical models in nanoplasmonics, cavity QED, and gravitational physics by enabling rigorous quantization and accurate observable predictions.

Green-tensor quasinormal-mode calibration refers to a suite of mathematical and computational techniques that utilize the system's quasinormal modes (QNMs)—the natural damped resonances of open, lossy, and/or radiative systems—to construct a calibrated (physically and mathematically accurate) expansion of the electromagnetic or gravitational Green tensor. This approach is central in nanoplasmonics, cavity quantum electrodynamics, and black hole perturbation theory for modeling light-matter interactions, spontaneous emission, energy transfer, and self-force effects across spatial and spectral domains. Calibration refers to the process by which the QNM expansion is regularized, normalized, and, if necessary, refined or matched to ensure the Green tensor recovers the underlying physical response everywhere of interest, including the treatment of singular near-field divergences, background response, and spectral completeness.

1. Fundamental Construction and Normalization of QNMs

Central to Green-tensor QNM calibration is the expansion of the system’s dyadic Green tensor, $\mathbf{G}(\mathbf{r},\mathbf{r}';\omega)$, in terms of QNMs $\{ \tilde{\mathbf{f}}_\mu \}$, which are solutions to the source-free Maxwell equations (or their analogous wave equations in other fields) with complex eigenfrequencies $\tilde{\omega}_\mu = \omega_\mu - i\gamma_\mu$ and outgoing boundary conditions (Ge et al., 2013, Ge et al., 2015, Franke et al., 2020). The QNMs must be normalized via a biorthogonal, typically complex, volume-plus-surface inner product:

$$\langle\!\langle \tilde{\mathbf{f}}_\mu | \tilde{\mathbf{f}}_\nu \rangle\!\rangle = \lim_{V \rightarrow \infty} \left\{ \int_V \frac{1}{2\omega} \left. \frac{\partial [\varepsilon(\mathbf{r},\omega)\omega^2]}{\partial \omega} \right|_{\omega=\tilde{\omega}_\mu} \tilde{\mathbf{f}}_\mu(\mathbf{r}) \cdot \tilde{\mathbf{f}}_\nu(\mathbf{r}) d^3 r + i \frac{n_B c}{2 \tilde{\omega}_\mu} \oint_{\partial V} \tilde{\mathbf{f}}_\mu(\mathbf{r}) \cdot \tilde{\mathbf{f}}_\nu(\mathbf{r}) dS \right\} = \delta_{\mu\nu}$$

for electromagnetic fields in dispersive, open systems (Ge et al., 2013, Ge et al., 2015). This normalization fixes the residue structure of the QNM expansion and is essential for rendering the Green tensor physically consistent and ensuring convergence.

2. QNM Expansion and Calibration of the Green Tensor

The Green tensor is systematically constructed as an expansion over the normalized QNMs:

$$\mathbf{G}(\mathbf{r},\mathbf{r}';\omega) = \mathbf{G}_B(\mathbf{r},\mathbf{r}';\omega) + \sum_{\mu=1}^{M} \frac{\omega^2}{2\tilde{\omega}_\mu (\tilde{\omega}_\mu - \omega)} \tilde{\mathbf{f}}_\mu(\mathbf{r}) \otimes \tilde{\mathbf{f}}_\mu(\mathbf{r}')$$

where $\mathbf{G}_B$ is the analytic Green tensor of the homogeneous background (Ge et al., 2013, Ge et al., 2015). Calibration entails the following concrete workflow:

1. Mode-solving: Use frequency-domain or time-domain solvers (e.g., FDTD with PMLs) to obtain QNMs and complex resonances $\tilde{\omega}_\mu$.
2. Normalization: Compute the normalization integrals above.
3. Bulk separation: Compute $\mathbf{G}_B$ in analytic or simulated form.
4. Assembly: Superpose the QNM expansion and background as above.
5. Refinement: Compare the result to full-wave simulations on a test grid in $(\mathbf{r},\mathbf{r}',\omega)$; if residuals exceed acceptable tolerance, iteratively adjust $\tilde{\omega}_\mu$ and normalization until convergence.
6. Validation: Increase $M$ (number of QNMs) and validate by comparing key observables (e.g., Purcell factor, EELS spectra) with full simulations or analytic benchmarks (Ge et al., 2015).

This procedure ensures the calibrated QNM Green tensor accurately reproduces the full physical Maxwell Green tensor within the spectral and spatial region of interest.

3. Regularization and Quasi-static Corrections

In nanoplasmonic systems, field divergences in the deep near-field ($x \to 0$) require analytic regularization. As the local density of optical states (LDOS) near metal surfaces diverges as $1/x^3$, the calibrated Green tensor incorporates a quasi-static correction (Ge et al., 2013):

$$\mathbf{G}^{qs}(\mathbf{r},\mathbf{r}';\omega) = \mp\,\mathbf{G}^B(\mathbf{r},-\mathbf{r}';\omega) \frac{\varepsilon_{\text{MNP}}(\omega) - \varepsilon_B}{2[\varepsilon_{\text{MNP}}(\omega) + \varepsilon_B]}$$

where the sign $\mp$ depends on dipole polarization. This analytic term is simply added to the QNM+background expansion, providing a uniformly accurate and finite Green tensor from sub-nanometer separations to infinity. Such regularization is critical for modeling deep near-field effects, emitter-surface energy transfer, and ensuring numerical stability (Ge et al., 2013).

4. Application to Quantum Fluctuations and Loss Channels

The fluctuation-dissipation theorem, together with Green-tensor QNM calibration, underpins rigorous quantization in absorptive and open systems. The field commutators and zero-point fluctuations are preserved if and only if the volume (material loss) and surface (radiation loss) contributions are handled according to (Franke et al., 2020):

$$\text{Im}\,\mathbf{G}(\mathbf{r},\mathbf{r}') = \int_{\text{material}} d^3s\,\varepsilon_I(\mathbf{s})\,\mathbf{G}(\mathbf{r},\mathbf{s})\,\mathbf{G}^*(\mathbf{s},\mathbf{r}') + \frac{c^2}{2i\omega^2} \oint_S \{ \mathbf{C} - \mathbf{C}^\dagger \}$$

Calibration ensures that in the strict non-dissipative limit ($\gamma_n\to0$), the QNM commutators reduce to the canonical Hermitian quantization scheme. In the presence of pure radiation loss, the surface term dominates. This makes the scheme robust for quantum nanophotonics and open-cavity QED (Franke et al., 2020). A worked example (3D photonic-crystal cavity) shows quantitative recovery of the semiclassical Purcell factor and photon emission rates from the calibrated Green tensor (Franke et al., 2020).

5. Quasinormal-mode Calibration in Gravitational Physics

The same QNM expansion and calibration strategy underpins the construction of the retarded Green tensor for wave propagation on curved spacetimes, such as Schwarzschild. The spectral construction involves a sum over QNM residues and a branch-cut (BC) integral; the matching region between local Hadamard (quasilocal) and distant QNM+BC expansions is determined empirically, and smoothness is imposed at the matching surface (Casals et al., 2013, Dolan et al., 2011):

• QNM sum: Dominates the Green function for $\Delta t > \Delta t_m$, converging rapidly in overtone number and multipole moment.
• Matching: The Hadamard expansion is accurate near coincidence, while the QNM+BC sum is accurate in the distant past; in practice, their overlap enables sub-percent matching error. This is essential in self-force calculations, black hole perturbation theory, and gravitational wave modeling (Casals et al., 2013).

Tensorial generalizations (for gravitational perturbations) utilize spin-weighted harmonics, Regge-Wheeler/Zerilli or Teukolsky functions, and polarization bitensors in the QNM expansion. The calibrated Green tensor—matching residue normalization, excitation factors, and caustic structure—faithfully encapsulates the physical propagator across all relevant causal domains (Dolan et al., 2011).

6. Alternative Modal Expansions and Calibration Strategies

Alternative expansions, such as the GENOME approach, define modes via the inclusion permittivity as the eigenvalue at fixed real frequency, yielding stationary, decaying modes that are trivially normalized by an interior volume integral (Chen et al., 2017). The full Green tensor is expanded as:

$$\mathbf{G}(\mathbf{r},\mathbf{r}';\omega) = \mathbf{G}_0(\mathbf{r},\mathbf{r}';\omega) + \sum_m \frac{\mathbf{E}_m(\mathbf{r}) \otimes \mathbf{E}_m(\mathbf{r}')}{\Delta_m(\omega)}$$

with appropriate resonance denominators $\Delta_m(\omega)$. Completeness is ensured via the Lippmann-Schwinger formalism for sources and detectors both inside and outside the resonator. Calibration is performed by truncating the mode sum at $N$ modes such that the remaining weight is negligible, and validating the expansion against reference simulations or analytic theory (Chen et al., 2017).

This approach bypasses implementation complexities associated with outgoing (divergent) QNMs, surface normalization, and PML boundary layers, making it particularly attractive for rapid simulations and analytic modeling in nanophotonics. The method converges rapidly and attains arbitrarily accurate agreement when properly calibrated.

7. Summary of Methodological Steps

The methodology underlying Green-tensor quasinormal-mode calibration, as exemplified across electromagnetic and gravitational systems, is summarized in the following procedural steps:

Step	Core Procedure	Reference
Mode computation	Solve open-boundary eigenproblem (e.g., FDTD+PML, FEM) for QNMs and eigenfrequencies	(Ge et al., 2015)
Normalization	Biorthogonal volume+surface (or interior-only) normalization of QNMs	(Ge et al., 2015, Chen et al., 2017)
Green tensor assembly	QNM expansion $\mathbf{G} = \mathbf{G}_B + \sum_\mu (\cdots)$	(Ge et al., 2013, Chen et al., 2017)
Regularization/calibration	Add analytic quasi-static (near-field) corrections as necessary	(Ge et al., 2013)
Refinement/validation	Fit expansion to full Maxwell (or wave) Green tensor, refine parameters, increase $M$ as needed	(Ge et al., 2015)
Quantum/statistical calibration	Confirm preservation of field commutators, FD theorem, and correct quantization rules	(Franke et al., 2020)
Gravitational generalization	Use spectral sum (QNM+BC), match to local Hadamard form in overlap region	(Casals et al., 2013)

The calibrated expansion, once constructed, enables rapid and accurate evaluation of all physical observables derived from the Green tensor across classical and quantum regimes.

---

References:

• (Ge et al., 2013) Quasinormal mode approach to modelling light-emission and propagation in nanoplasmonics
• (Ge et al., 2015) Quasinormal mode theory and modelling of electron energy loss spectroscopy for plasmonic nanostructures
• (Franke et al., 2020) Fluctuation-dissipation theorem and fundamental photon commutation relations in lossy nanostructures using quasinormal modes
• (Casals et al., 2013) Self-Force and Green Function in Schwarzschild spacetime via Quasinormal Modes and Branch Cut
• (Dolan et al., 2011) Wave Propagation and Quasinormal Mode Excitation on Schwarzschild Spacetime
• (Chen et al., 2017) Generalizing normal mode expansion of electromagnetic Green's tensor to lossy resonators in open systems