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

This paper presents a numerically efficient framework based on a Dyson scattering equation and quasinormal modes to calculate the scattered electromagnetic Green's function for systems of multiple coupled optical resonators with finite retardation, avoiding complex nested integrals while maintaining high accuracy.

The Gist

The Big Picture: The "Echo Chamber" Problem

Imagine you are in a large hall filled with different rooms (resonators), and you want to know exactly how sound travels from one room to another. In the world of light and quantum physics, these "rooms" are tiny optical cavities (like metal nanoparticles), and the "sound" is light.

To design future quantum computers or super-sensitive sensors, scientists need to know exactly how light behaves when it bounces between these cavities. This behavior is described by something called the Green's function. Think of the Green's function as the "ultimate map" that tells you: If I drop a pebble (a photon) here, exactly where and when will the ripples appear everywhere else?

The Problem:
Calculating this map directly is incredibly hard. It's like trying to predict the path of every single water ripple in a stormy ocean by simulating every drop of water. For complex systems with multiple cavities, the math gets so heavy that computers crash, or the calculations take forever.

The Old Way (The "Quasinormal Mode" Problem):
Scientists previously used a method called Quasinormal Modes (QNMs). Imagine these as the "natural notes" a cavity sings when you tap it.

• The Catch: These notes are "leaky." Because the cavities are open to the outside world, the sound fades away. Mathematically, this causes a problem: if you try to use these notes to describe the sound far away from the cavity, the math blows up (it goes to infinity). To fix this, you needed to use thousands of these "notes" to get a decent answer, which is computationally impossible for complex quantum systems.

The New Solution: The "Relay Race" Approach

The authors of this paper have invented a clever new framework to solve this. Instead of trying to calculate the whole ocean at once, they treat the system like a relay race.

Here is how their method works, step-by-step:

1. The Solo Run (Single Cavity)

First, they calculate the "natural note" (the QNM) for just one cavity in isolation. They know exactly how that single room sings.

• Analogy: You know exactly how a single bell rings when struck alone.

2. The Regularized Signal (Fixing the "Infinity" Problem)

The authors realized that while the "natural note" gets messy far away, you can create a "cleaned-up" version of that note that travels nicely through space without blowing up.

• Analogy: Instead of shouting the raw, echoing sound of the bell, you record it and play it back through a high-quality speaker that projects the sound clearly across the room without distortion.

3. The Relay Race (The Dyson Equation)

This is the core innovation. Instead of calculating the interaction between 10 cavities all at once (which is a nightmare), they build the answer step-by-step:

• Step 1: Calculate the map for Cavity 1.
• Step 2: Add Cavity 2. Use the map from Step 1 as the "input" to see how Cavity 2 scatters the light coming from Cavity 1.
• Step 3: Add Cavity 3. Use the map from Step 2 as the input.
• Analogy: Imagine a game of "Telephone." Person A whispers a message to Person B. Person B doesn't try to remember the whole history of the conversation; they just take what they heard from A and pass it to C. By doing this in a chain, you get the final message without needing a supercomputer to track every whisper simultaneously.

4. Handling the "Time Delay" (Retardation)

Light doesn't travel instantly; it takes time to cross the gap between cavities. This is called retardation.

• The Innovation: Many old methods ignored this time delay or approximated it poorly, leading to errors when cavities were far apart. This new method explicitly includes the time it takes for light to travel between the "relay stations."
• Analogy: If you are shouting across a canyon, you have to account for the time it takes for your voice to bounce back. This method ensures the "echo" arrives at the exact right moment, not too early and not too late.

Why This Matters (The "So What?")

The authors tested their theory using two metal "dumbbells" (dimers) with tiny gaps, placing light sources inside them. They compared their new "relay race" math against the "brute force" computer simulation (the gold standard).

The Result: Their new method matched the brute force simulation perfectly, even when the cavities were far apart.

Why is this a big deal?

1. Speed: It turns a calculation that might take days or require a supercomputer into something that takes seconds on a laptop.
2. Scalability: Because the math breaks down into simple "two-cavity" steps, it can easily be extended to 10, 20, or 100 cavities.
3. Quantum Tech: This is crucial for building quantum networks. If you want to send quantum information from one node to another, you need to know exactly how they talk to each other. This method gives engineers a precise, fast, and easy-to-use blueprint for designing these devices.

Summary Metaphor

Imagine you are trying to predict the traffic flow in a massive city with thousands of intersections.

• The Old Way: Simulate every single car, every driver's decision, and every traffic light change for the whole city at once. It's impossible.
• The New Way: Figure out the traffic rules for one intersection. Then, figure out how Intersection A talks to Intersection B. Then, how B talks to C. By chaining these simple interactions together, you can predict the traffic flow for the whole city instantly, without ever needing to simulate a single car.

This paper provides the "traffic rules" for light in complex quantum systems, making the design of future quantum technologies much faster and more reliable.

Technical Summary

1. Problem Statement

The electromagnetic Green's function is fundamental for modeling light-matter interactions in modern quantum photonic devices (e.g., nanolasers, quantum networks, and sensors). However, calculating the full Green's function for systems with multiple spatially separated open cavities is computationally challenging.

• Limitations of Existing Methods:
  • Direct Numerical Calculation: Often too computationally expensive for complex quantum dynamics where the Hilbert space scales exponentially with the number of modes.
  • Standard Quasinormal Mode (QNM) Expansion: While effective near a single resonator, QNMs diverge spatially at large distances due to their complex eigenfrequencies. This requires a large number of modes for accurate far-field approximations, rendering "few-mode" expansions impractical for multi-cavity systems.
  • Coupled QNM Theory (CQT): Existing theories for coupled cavities often rely on divergent QNMs and fail to accurately account for finite retardation effects (time delays) when cavities are widely separated.

2. Methodology

The authors propose a numerically efficient framework based on a Dyson scattering equation that constructs the multi-cavity Green's function iteratively from the QNMs of individual cavities.

Core Framework:

1. Iterative Scattering Approach:
   • The system consists of $N$ spatially separated cavities.
   • The total Green's function $G^{(N)}$ is built iteratively by adding one cavity at a time to the Green's function of the previous $N-1$ cavities ($G^{(N-1)}$).
   • This is governed by a scattering equation:
     $$G^{(n)}(r, r', \omega) = G^{(n-1)}(r, r', \omega) + \int_{V_n} d^3s \, \Delta\epsilon(s, \omega) G^{(n-1)}(r, s, \omega) \cdot G^{(n)}(s, r', \omega)$$
2. Termination Condition (Few-Mode Approximation):
   • To avoid infinite recursion, the term $G^{(n)}$ on the right-hand side (inside the $n$-th cavity) is approximated by the single-cavity Green's function of that specific cavity ($G_n$).
   • This assumption holds well when cavities are separated by distances on the order of the wavelength or larger, where inter-cavity influence is weak compared to intra-cavity fields.
3. Regularized QNM Fields:
   • To handle the spatial divergence of QNMs outside the cavity, the authors use regularized frequency-dependent fields ($\tilde{F}_{i\mu}$).
   • These are calculated using surface integrals (via the field-equivalence principle) rather than volume integrals, which avoids numerical issues related to field discontinuities at boundaries.
4. Inclusion of Retardation:
   • The framework explicitly includes the phase factor $e^{i\omega R/c}$ representing the finite time delay (retardation) for light traveling between cavities.
   • A "pole approximation" is applied to the slowly varying envelope of the regularized fields to simplify the calculation of coupling coefficients.

Key Mathematical Result:

The multi-cavity scattering process naturally decomposes into products of two-cavity scattering processes. For two cavities (1 and 2), the Green's function between them factorizes as:
$$G(r_b, r_a, \omega) \approx B_{21}(\omega) \tilde{f}_2(r_b) \tilde{f}_1(r_a)$$
where $B_{21}$ is a coupling coefficient derived from surface integrals involving the regularized fields and single-cavity QNMs. This decomposition avoids computationally expensive nested integrals.

3. Key Contributions

• Dyson Equation for Multi-Cavity Systems: Introduced a recursive Dyson scattering formalism that allows the construction of the full Green's function for $N$ cavities using only single-cavity QNMs.
• Handling Finite Retardation: Successfully integrated retardation effects into a few-mode QNM expansion, solving a known failure point of previous coupled QNM theories for large separations.
• Factorization of Scattering: Demonstrated that multi-cavity scattering terms decompose into products of two-cavity terms, significantly reducing computational complexity (scaling quadratically with the number of cavities in a single-mode approximation).
• Surface Integral Formulation: Utilized surface integrals for regularized fields to ensure numerical stability and accuracy, particularly for far-field interactions.

4. Results

The authors validated their framework using a system of two coupled metal dimers (acting as QNM cavities) with dipole emitters placed in the gaps.

• Validation Cases:
  • Large Separation ($R \approx 2.75\lambda$): The QNM expansion showed excellent agreement with full numerical solutions of Maxwell's equations. Crucially, including the frequency-dependent phase (retardation) was essential; omitting it led to significant errors.
  • Shorter Separation ($R \approx 1.04\lambda$): The method remained highly accurate.
  • Very Short Separation ($R \approx 0.18\lambda$): The method qualitatively captured the coupling shape but slightly underestimated the strength. The authors note this is expected as the "termination condition" (assuming the field inside a cavity is dominated only by its own QNMs) breaks down when cavities are in the near-field regime.
• Comparison with Single-Cavity Approximation:
  • A naive expansion using only the QNMs of the source cavity (ignoring the second cavity's influence) failed to capture the coupling strength and spectral shape, confirming the necessity of the multi-cavity expansion.
• Identical Dimers: The method was also tested on identical dimers, yielding excellent agreement with numerical results.

5. Significance

• Enabling Quantum Dynamics Simulations: By providing an accurate few-mode approximation, this framework makes it feasible to simulate complex quantum dynamics (e.g., entanglement generation, photon transport) in multi-resonator networks where the Hilbert space would otherwise be too large for direct numerical methods.
• Design Tool for Quantum Devices: It offers an intuitive and efficient tool for designing nanolasers, quantum sensors, and quantum networks involving spatially separated resonators.
• Bridging Theory and Computation: The approach bridges the gap between rigorous full-wave numerical simulations (which are slow) and analytical few-mode theories (which are often inaccurate for separated systems), offering a "best of both worlds" solution for open cavity quantum electrodynamics.

In summary, this work provides a robust, scalable, and physically accurate method to calculate electromagnetic interactions in complex, multi-resonator photonic systems, specifically addressing the challenges of spatial separation and finite retardation.