2025-11-15T23:07:18.927821

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

Fuchs, Ren, Hughes et al.
The electromagnetic Green's function is a crucial ingredient for the theoretical study of modern photonic quantum devices, but is often difficult or even impossible to calculate directly. We present a numerically efficient framework for calculating the scattered electromagnetic Green's function of a multi-cavity system with spatially separated open cavities (with arbitrary shape, dispersion and loss) and finite retardation times. The framework is based on a Dyson scattering equation that enables the construction of the Green's function from the quasinormal modes of the individual resonators within a few-mode approximation and a finite number of iteration steps without requiring nested integrals. The approach shows excellent agreement with the full numerical Green's function for the example of two coupled dipoles located in the gaps of two metal dimers serving as quasinormal mode cavities, and is easily extended to arbitrarily large separations and multiple cavities.
academic

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

基本信息

  • 论文ID: 2510.12511
  • 标题: Green's function expansion for multiple coupled optical resonators with finite retardation using quasinormal modes
  • 作者: Robert Fuchs¹, Juanjuan Ren², Stephen Hughes², Marten Richter¹
  • 机构: ¹Technische Universität Berlin, ²Queen's University
  • 分类: cond-mat.mes-hall physics.optics quant-ph
  • 发表时间: October 15, 2025
  • 论文链接: https://arxiv.org/abs/2510.12511

摘要

电磁格林函数是现代光子量子器件理论研究的关键要素，但通常难以甚至无法直接计算。本文提出了一个数值高效的框架，用于计算具有空间分离开放腔体（任意形状、色散和损耗）和有限延迟时间的多腔系统的散射电磁格林函数。该框架基于Dyson散射方程，能够在少模近似和有限迭代步骤内从单个谐振器的准正常模式构建格林函数，无需嵌套积分。该方法在两个耦合偶极子（位于作为准正常模腔体的金属二聚体间隙中）的例子中与完整数值格林函数显示出优异的一致性，并可轻松扩展到任意大的分离距离和多个腔体。

研究背景与动机

问题背景

  1. 核心问题: 电磁格林函数在现代光子量子器件理论研究中至关重要，但对于多腔耦合系统，直接计算格林函数极其困难或不可行。
  2. 重要性: 格林函数广泛应用于：
    • 光学腔体的Purcell增强计算
    • 电磁场的量子化
    • 光子与量子发射器耦合元素的推导
    • 量子光学和量子动力学应用
  3. 现有方法局限性:
    • 单腔QNM方法: 对于远离谐振器的位置，QNM由于复特征频率而发散，需要大量模式才能获得良好近似
    • 耦合准正常模理论(CQT): 依赖发散的QNM，对于大分离距离不够准确
    • 多腔精确QNM计算: 对于多个空间分离的耦合结构往往不可行
  4. 研究动机: 需要一个能够：
    • 处理任意形状、色散和损耗的空间分离开放腔体
    • 包含有限延迟效应
    • 在少模近似下高效计算
    • 适用于大空间分离的多腔系统

核心贡献

  1. 提出了基于Dyson散射方程的多腔格林函数计算框架，能够从单腔准正常模式构建多腔系统的格林函数
  2. 开发了迭代构建方法，通过有限次迭代步骤逐步添加腔体，避免了嵌套积分的计算复杂性
  3. 实现了正则化QNM场的应用，使方法适用于具有大空间分离和显著延迟效应的系统
  4. 验证了方法的准确性，在金属二聚体耦合系统中与完整数值格林函数显示出优异一致性
  5. 提供了可扩展的理论框架，自然地将多腔散射分解为双腔散射过程的乘积

方法详解

任务定义

输入: N个空间分离的开放腔体系统，每个腔体具有任意形状、色散和损耗特性 输出: 多腔系统的完整电磁格林函数G^(N)(r,r',ω) 约束: 在少模近似下，考虑有限延迟效应，避免计算复杂的嵌套积分

理论框架

1. 系统描述

多腔系统的完整介电常数表示为：

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

其中V_i(r,ω)是第i个腔体引起的介电常数扰动。

2. 单腔准正常模式

单腔QNM f̃_iμ(r)满足无源亥姆霍兹方程：

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

具有复特征频率ω̃_iμ = ω_iμ - iγ_iμ，γ_iμ > 0。

3. 正则化QNM场

为解决QNM在远场发散问题，引入正则化QNM场：

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

4. 多腔格林函数迭代构建

核心创新是通过一系列散射方程迭代构建：

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

终止条件：

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

5. 表面积分表示

为提高数值效率，将体积分转换为表面积分：

ω²/c² G^(n)(r,r')|_(r∉V_n) = χ_V_n(r') ω²/c² G^(n-1)(r,r') + [表面积分项]

技术创新点

  1. 分解策略: 多腔散射自然分解为双腔散射过程的乘积，避免嵌套积分
  2. 延迟效应处理: 通过相位因子e^(iωR_ij/c)正确包含有限延迟时间
  3. 极点近似: 对慢变包络函数应用极点近似ω→ω_i，提高计算效率
  4. 表面积分优化: 使用表面积分避免介质边界处电场不连续性问题

实验设置

系统配置

  • 测试系统: 两个金属二聚体作为QNM腔体
  • 材料模型: Drude模型，ε_j(ω) = 1 - ω_p²/(ω² + iωγ_pj)
  • 参数设置:
    • 等离子体频率：ℏω_p = 8.2934 eV
    • 阻尼参数：ℏγ_p1 = 0.0928 eV, ℏγ_p2 = 0.3ℏγ_p1
  • QNM特征频率:
    • 二聚体1：ℏω̃_1 = (1.6904 - i0.0652) eV (λ_1 ≈ 733.46 nm)
    • 二聚体2：ℏω̃_2 = (1.6482 - i0.0388) eV (λ_2 ≈ 752.24 nm)

测试配置

  • 偶极子设置: 两个相同的z极化偶极子d_a = d_b = dê_z
  • 位置: 分别位于两个二聚体间隙中(r_a ∈ V_1, r_b ∈ V_2)
  • 分离距离: R_12 = 2020 nm和760 nm两种配置

评价指标

偶极子耦合强度：

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

对比方法

  1. 完整数值格林函数 (G_full): 严格求解3D麦克斯韦方程
  2. QNM展开 (G_QNM): 本文方法，公式(12)+(14)
  3. 极点近似 (G_QNM^pole): 公式(12)+(15)，相位包含在极点近似中
  4. 单腔展开 (G_1): 仅考虑单腔QNM的传统方法

实验结果

主要结果

1. 大分离距离情况 (R_12 = 2020 nm)

  • QNM展开方法与完整数值解显示优异一致性
  • 正确的延迟效应处理(e^(iωR_12/c))对准确性至关重要
  • 极点近似方法(G_QNM^pole)由于未完全考虑延迟效应而偏离数值解
  • 单腔展开严重低估耦合强度，无法捕获腔间QNM耦合效应

2. 中等分离距离情况 (R_12 = 760 nm)

  • QNM展开方法依然与数值解完美匹配
  • 极点近似方法由于延迟效应较弱而表现更好
  • 单腔展开仍然在定量和定性上都失效

3. 关键数值结果

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

方法验证

  • 无拟合参数: 所有计算均基于第一性原理，无需调节参数
  • 频率依赖性: 正确捕获了耦合强度的频率依赖特性
  • 相位效应: 延迟相位e^(iωR_12/c)对大分离距离的准确性至关重要

计算效率

  • 避免嵌套积分: 多腔散射分解为双腔过程乘积
  • 有限迭代步骤: N个腔体需要N次迭代
  • 二次标度: 单模近似下，散射积分数量与腔体数量呈二次关系

相关工作

主要研究方向

  1. 准正常模理论: Sauvan等人的开放腔体自发辐射理论
  2. 场量子化: 基于QNM的电磁场量子化方法
  3. 耦合腔理论: 耦合准正常模理论(CQT)及其局限性
  4. 格林函数方法: 单腔QNM展开和正则化技术

本文优势

  1. 延迟效应处理: 相比CQT，能够准确处理大分离距离情况
  2. 计算效率: 避免了传统方法的嵌套积分问题
  3. 普适性: 适用于任意形状、色散和损耗的腔体
  4. 可扩展性: 自然扩展到多腔系统

结论与讨论

主要结论

  1. 成功开发了多腔格林函数计算框架，基于Dyson散射方程和单腔QNM
  2. 实现了高精度计算，在不同分离距离下与数值解优异一致
  3. 解决了延迟效应问题，使方法适用于大空间分离系统
  4. 提供了可扩展解决方案，易于推广到任意数量的腔体

局限性

  1. 少模近似: 方法基于每个腔体由少数主导QNM描述的假设
  2. 终止条件: 依赖于腔内格林函数可由单腔QNM良好近似的假设
  3. 3D均匀介质: 主要针对均匀背景介质中的3D腔体进行验证
  4. 计算复杂度: 虽然避免了嵌套积分，但仍需计算多个表面积分

未来方向

  1. 扩展到波导耦合系统: 适配不同边界条件
  2. 非线性效应: 包含非线性光学过程
  3. 量子动力学应用: 在量子光学系统中的具体应用
  4. 优化算法: 进一步提高大型多腔系统的计算效率

深度评价

优点

  1. 理论创新性强: Dyson散射方程的迭代应用是新颖且巧妙的
  2. 实用价值高: 解决了多腔量子光学系统的实际计算需求
  3. 验证充分: 与严格数值解的对比验证了方法的准确性
  4. 写作清晰: 理论推导严谨，物理图像明确

不足

  1. 实验验证缺乏: 仅有理论和数值验证，缺少实验验证
  2. 应用范围: 主要针对特定类型的腔体系统进行了测试
  3. 复杂系统: 对于强耦合或非线性系统的适用性尚不明确
  4. 计算细节: 某些数值实现细节可能影响方法的可重现性

影响力

  1. 学术贡献: 为多腔量子光学理论提供了重要工具
  2. 技术应用: 对量子网络、量子计算等领域具有潜在应用价值
  3. 方法论意义: 展示了QNM理论在复杂系统中的应用潜力
  4. 可复现性: 理论框架清晰，便于其他研究者应用和扩展

适用场景

  1. 量子光学器件设计: 多腔耦合系统的理论分析
  2. 等离激元学: 金属纳米结构的电磁响应计算
  3. 量子网络: 空间分离量子节点的耦合分析
  4. 光子晶体: 缺陷腔耦合系统的研究

参考文献

论文引用了47篇重要参考文献，涵盖了准正常模理论、量子光学、电磁场量子化等相关领域的关键工作，为本研究提供了坚实的理论基础。

总体评价: 这是一篇高质量的理论物理论文，提出了解决多腔光学系统格林函数计算这一重要问题的创新方法。理论推导严谨，数值验证充分，对量子光学和光子学领域具有重要意义。方法的普适性和可扩展性使其具有广泛的应用前景。