Green's function expansion for multiple coupled optical resonators with finite retardation using quasinormal modes

Basic Information

• Paper ID: 2510.12511
• Title: Green's function expansion for multiple coupled optical resonators with finite retardation using quasinormal modes
• Authors: Robert Fuchs¹, Juanjuan Ren², Stephen Hughes², Marten Richter¹
• Institutions: ¹Technische Universität Berlin, ²Queen's University
• Classification: cond-mat.mes-hall physics.optics quant-ph
• Publication Date: October 15, 2025
• Paper Link: https://arxiv.org/abs/2510.12511

Abstract

The electromagnetic Green's function is a key element in theoretical studies of modern photonic quantum devices, but is typically difficult or impossible to compute directly. This paper presents a numerically efficient framework for calculating scattering electromagnetic Green's functions for multi-cavity systems with spatially separated open cavities (arbitrary shape, dispersion, and loss) and finite retardation time. The framework is based on the Dyson scattering equation and enables construction of Green's functions from quasinormal modes (QNMs) of individual resonators within few-mode approximation and finite iteration steps, without requiring nested integrals. The method demonstrates excellent agreement with complete numerical Green's functions in the example of two coupled dipoles located in the gap of a metallic dimer serving as a QNM cavity, and can be easily extended to arbitrarily large separations and multiple cavities.

Research Background and Motivation

Problem Background

1. Core Issue: The electromagnetic Green's function is crucial in modern photonic quantum device theory, but direct computation of Green's functions for multi-cavity coupled systems is extremely difficult or infeasible.
2. Importance: Green's functions are widely applied to:
   • Purcell enhancement calculations in optical cavities
   • Quantization of electromagnetic fields
   • Derivation of coupling elements between photons and quantum emitters
   • Quantum optics and quantum dynamics applications
3. Limitations of Existing Methods:
   • Single-cavity QNM method: For positions far from resonators, QNMs diverge due to complex eigenfrequencies, requiring numerous modes for good approximation
   • Coupled quasinormal theory (CQT): Relies on divergent QNMs, insufficiently accurate for large separation distances
   • Multi-cavity exact QNM computation: Often infeasible for multiple spatially separated coupled structures
4. Research Motivation: Need for a method that can:
   • Handle spatially separated open cavities with arbitrary shape, dispersion, and loss
   • Include finite retardation effects
   • Efficiently compute within few-mode approximation
   • Apply to multi-cavity systems with large spatial separations

Core Contributions

1. Proposed a Dyson scattering equation-based framework for multi-cavity Green's function computation that constructs multi-cavity system Green's functions from single-cavity QNMs
2. Developed an iterative construction method that progressively adds cavities through finite iteration steps, avoiding computational complexity of nested integrals
3. Implemented regularized QNM field applications, enabling the method to handle systems with large spatial separations and significant retardation effects
4. Verified method accuracy by demonstrating excellent agreement with complete numerical Green's functions in metallic dimer coupled systems
5. Provided a scalable theoretical framework that naturally decomposes multi-cavity scattering into products of two-cavity scattering processes

Methodology Details

Task Definition

Input: N spatially separated open cavity systems, each with arbitrary shape, dispersion, and loss characteristics Output: Complete electromagnetic Green's function G^(N)(r,r',ω) of the multi-cavity system Constraints: Few-mode approximation with finite retardation effects, avoiding computation of complex nested integrals

Theoretical Framework

1. System Description

The complete permittivity of the multi-cavity system is expressed as:

ε(r,ω) = ε_back(r,ω) + Σ(i=1 to N) V_i(r,ω)

where V_i(r,ω) is the permittivity perturbation caused by the i-th cavity.

2. Single-Cavity Quasinormal Modes

Single-cavity QNM f̃_iμ(r) satisfies the source-free Helmholtz equation:

∇×∇×f̃_iμ(r) - (ω̃²_iμ/c²)ε_i(r,ω̃_iμ)f̃_iμ(r) = 0

with complex eigenfrequencies ω̃_iμ = ω_iμ - iγ_iμ, γ_iμ > 0.

3. Regularized QNM Field

To address QNM divergence in the far field, regularized QNM fields are introduced:

F̃_iμ(r,ω) = ∫_V_i d³r' Δε(r',ω)G_back(r,r',ω)·f̃_iμ(r')

4. Iterative Construction of Multi-Cavity Green's Function

The core innovation is iterative construction through a series of scattering equations:

G^(n)(r,r',ω) = G^(n-1)(r,r',ω) + ∫_V_n d³s Δε(s,ω)G^(n-1)(r,s,ω)·G^(n)(s,r',ω)

Termination Condition:

G^(k)(r,r',ω)|_(r,r'∈V_i) ≈ G_i(r,r',ω)

5. Surface Integral Representation

To enhance numerical efficiency, volume integrals are converted to surface integrals:

ω²/c² G^(n)(r,r')|_(r∉V_n) = χ_V_n(r') ω²/c² G^(n-1)(r,r') + [surface integral terms]

Technical Innovations

1. Decomposition Strategy: Multi-cavity scattering naturally decomposes into products of two-cavity scattering processes, avoiding nested integrals
2. Retardation Effect Handling: Finite retardation time is correctly included through phase factors e^(iωR_ij/c)
3. Pole Approximation: Pole approximation ω→ω_i applied to slowly varying envelope functions for improved computational efficiency
4. Surface Integral Optimization: Surface integrals used to avoid electric field discontinuity issues at medium boundaries

Experimental Setup

System Configuration

• Test System: Two metallic dimers serving as QNM cavities
• Material Model: Drude model, ε_j(ω) = 1 - ω_p²/(ω² + iωγ_pj)
• Parameter Settings:
  • Plasma frequency: ℏω_p = 8.2934 eV
  • Damping parameters: ℏγ_p1 = 0.0928 eV, ℏγ_p2 = 0.3ℏγ_p1
• QNM Eigenfrequencies:
  • Dimer 1: ℏω̃_1 = (1.6904 - i0.0652) eV (λ_1 ≈ 733.46 nm)
  • Dimer 2: ℏω̃_2 = (1.6482 - i0.0388) eV (λ_2 ≈ 752.24 nm)

Test Configuration

• Dipole Setup: Two identical z-polarized dipoles d_a = d_b = dê_z
• Positions: Located in the gaps of the two dimers respectively (r_a ∈ V_1, r_b ∈ V_2)
• Separation Distances: Two configurations with R_12 = 2020 nm and 760 nm

Evaluation Metrics

Dipole coupling strength:

g_ba(ω) = (1/ε_0) d_b · G(r_b, r_a, ω) · d_a

Comparison Methods

1. Complete Numerical Green's Function (G_full): Rigorous solution of 3D Maxwell equations
2. QNM Expansion (G_QNM): Proposed method, equations (12)+(14)
3. Pole Approximation (G_QNM^pole): Equations (12)+(15), with phase included in pole approximation
4. Single-Cavity Expansion (G_1): Traditional method considering only single-cavity QNM

Experimental Results

Main Results

1. Large Separation Distance Case (R_12 = 2020 nm)

• QNM expansion method shows excellent agreement with complete numerical solution
• Correct handling of retardation effects (e^(iωR_12/c)) is crucial for accuracy
• Pole approximation method (G_QNM^pole) deviates from numerical solution due to incomplete consideration of retardation effects
• Single-cavity expansion severely underestimates coupling strength, failing to capture inter-cavity QNM coupling effects

2. Moderate Separation Distance Case (R_12 = 760 nm)

• QNM expansion method remains in perfect agreement with numerical solution
• Pole approximation method performs better due to weaker retardation effects
• Single-cavity expansion fails both quantitatively and qualitatively

3. Key Numerical Results

For R_12 = 2020 nm: N_21 = (2.0694 - i0.1357) × 10^(-7) nm^(-2) For R_12 = 760 nm: N_21 = (5.3494 + i0.2984) × 10^(-7) nm^(-2)

Method Verification

• No Fitting Parameters: All calculations based on first principles without parameter adjustment
• Frequency Dependence: Correctly captures frequency-dependent characteristics of coupling strength
• Phase Effects: Retardation phase e^(iωR_12/c) is crucial for accuracy at large separations

Computational Efficiency

• Avoids Nested Integrals: Multi-cavity scattering decomposed into products of two-cavity processes
• Finite Iteration Steps: N cavities require N iterations
• Quadratic Scaling: Under single-mode approximation, number of scattering integrals scales quadratically with number of cavities

Related Work

Main Research Directions

1. Quasinormal Mode Theory: Sauvan et al.'s spontaneous emission theory for open cavities
2. Field Quantization: Electromagnetic field quantization methods based on QNMs
3. Coupled Cavity Theory: Coupled quasinormal mode theory (CQT) and its limitations
4. Green's Function Methods: Single-cavity QNM expansion and regularization techniques

Advantages of This Work

1. Retardation Effect Handling: Accurately handles large separation cases compared to CQT
2. Computational Efficiency: Avoids nested integral problems of traditional methods
3. Universality: Applicable to cavities with arbitrary shape, dispersion, and loss
4. Scalability: Naturally extends to multi-cavity systems

Conclusions and Discussion

Main Conclusions

1. Successfully developed multi-cavity Green's function computation framework based on Dyson scattering equation and single-cavity QNMs
2. Achieved high-precision computation with excellent agreement to numerical solutions at different separation distances
3. Solved retardation effect problem, enabling application to systems with large spatial separations
4. Provided scalable solution easily generalizable to arbitrary numbers of cavities

Limitations

1. Few-Mode Approximation: Method assumes each cavity is described by few dominant QNMs
2. Termination Condition: Depends on assumption that cavity Green's function can be well approximated by single-cavity QNMs
3. 3D Homogeneous Medium: Verification primarily focused on 3D cavities in uniform background medium
4. Computational Complexity: While avoiding nested integrals, still requires computing multiple surface integrals

Future Directions

1. Extension to Waveguide-Coupled Systems: Adaptation to different boundary conditions
2. Nonlinear Effects: Inclusion of nonlinear optical processes
3. Quantum Dynamics Applications: Specific applications in quantum optical systems
4. Algorithm Optimization: Further improvement of computational efficiency for large multi-cavity systems

In-Depth Evaluation

Strengths

1. Strong Theoretical Innovation: Iterative application of Dyson scattering equation is novel and elegant
2. High Practical Value: Addresses actual computational needs of multi-cavity quantum optical systems
3. Thorough Verification: Comparison with rigorous numerical solutions validates method accuracy
4. Clear Presentation: Rigorous theoretical derivation with clear physical interpretation

Weaknesses

1. Lack of Experimental Verification: Only theoretical and numerical verification, no experimental validation
2. Limited Application Scope: Testing primarily focused on specific cavity system types
3. Complex Systems: Applicability to strongly coupled or nonlinear systems remains unclear
4. Implementation Details: Some numerical implementation details may affect method reproducibility

Impact

1. Academic Contribution: Provides important tool for multi-cavity quantum optics theory
2. Technical Application: Potential applications in quantum networks, quantum computing, etc.
3. Methodological Significance: Demonstrates potential of QNM theory in complex systems
4. Reproducibility: Clear theoretical framework facilitates application and extension by other researchers

Applicable Scenarios

1. Quantum Optical Device Design: Theoretical analysis of multi-cavity coupled systems
2. Plasmonics: Electromagnetic response calculation of metallic nanostructures
3. Quantum Networks: Coupling analysis of spatially separated quantum nodes
4. Photonic Crystals: Study of defect cavity coupled systems

References

The paper cites 47 important references covering quasinormal mode theory, quantum optics, electromagnetic field quantization, and other related fields, providing solid theoretical foundation for this research.

Overall Assessment: This is a high-quality theoretical physics paper proposing an innovative method to solve the important problem of Green's function computation for multi-cavity optical systems. The theoretical derivation is rigorous, numerical verification is thorough, and the work has significant implications for quantum optics and photonics. The universality and scalability of the method provide broad application prospects.