当前位置：首页 > news >正文

离散模型解析嵌入式束缚态与法诺共振：从原理到光子器件设计

news 2026/5/26 20:53:18

1. 离散模型解析共振：从法诺共振到嵌入式束缚态

在光子晶体、光学传感和滤波器的设计里，我们常常需要精确操控光。比如，如何让特定频率的光几乎完全透过，而让相邻频率的光被强烈反射？这背后往往是一种被称为“共振”的物理现象在起作用。传统上，我们理解共振，就像理解一个钟摆：给它一个接近其自然摆动频率的推力，它就会越摆越高。在光学系统里，这个“钟摆”可能是光学微腔中的一个特定模式，当入射光的频率与之匹配时，能量会在腔内剧烈振荡，导致透射或反射谱上出现尖锐的峰或谷。

然而，有一类更微妙、也更强大的共振，其根源并非一个孤立的、稳定的“钟摆”，而是一种被称为“嵌入式束缚态”的特殊状态。想象一下，在一条繁忙的高速公路（连续谱）中间，存在一个完美的、静止的停车位（束缚态）。这个车位本身是高速公路的一部分，但车流却奇迹般地绕开了它。这种状态极其脆弱，只要道路标线稍有歪斜（结构不对称），这个车位就会立刻被车流淹没，但同时，车流会在原车位附近形成一个复杂的拥堵和疏通车流模式——这就是共振。本文要探讨的，正是如何用一个简洁的“玩具”模型——离散质量-弹簧链模型，来透彻地理解这种由嵌入式束缚态破坏所引发的、高度可调的共振现象，特别是著名的法诺共振线型。

2. 核心概念：从连续谱中的束缚态到共振

2.1 什么是嵌入式束缚态？

在典型的量子力学或波动力学系统中，束缚态的能量（或频率）通常位于连续谱之外。例如，原子中的电子能级是分立的，低于电离阈值。然而，在某些特殊的对称结构中，可以存在其能量（频率）恰好落在系统允许的传播模式所对应的连续谱范围内的束缚态。这就是“嵌入式束缚态”。

它的关键特性是“非鲁棒性”。这意味着该束缚态的存在强烈依赖于系统参数的精确取值（如特定的对称性、特定的入射角）。一旦这些参数发生哪怕微小的扰动，破坏了原有的对称性，这个束缚态就不再是严格的束缚态，其能量会获得一个虚部，变成一个“泄漏模”或“准导模”，并与连续谱中的传播模式发生强耦合。这个耦合过程，就在透射谱中表现为一个尖锐的共振特征。

2.2 法诺共振：非对称线型的典范

法诺共振是一种特殊的共振线型，其透射率曲线呈非对称的峰-谷形状，而非简单的洛伦兹型峰。它源于两种路径的干涉：一是光直接通过系统的背景路径，二是光通过共振态（即被破坏的嵌入式束缚态）的间接路径。这两种路径之间的相消或相长干涉，导致了在某些频率下透射几乎为零（谷），而在邻近频率下透射接近百分之百（峰）。

经典的Fano公式描述了这一线型：T(ε) ∝ |q + ε|² / (1 + ε²)，其中ε是归一化失谐频率，q是法诺参数，决定了线型的非对称程度。q很大时接近对称的洛伦兹峰；q很小时表现为尖锐的谷；q为中间值时呈现典型的峰谷结构。

2.3 离散模型的价值

面对复杂的连续系统（如光子晶体平板），完全解析求解麦克斯韦方程组是极其困难的。离散模型的价值在于，它用简单的差分方程（类比于质量-弹簧系统）抓住了连续系统的核心物理，同时保留了数学上可解析处理的优势。通过构建一个能同时支持传播模式和倏逝模式的离散格点模型，我们可以清晰地追踪嵌入式束缚态如何形成、又如何被破坏，并精确计算出共振线型随系统参数（如不对称性）的变化，为器件设计提供直观的“调控旋钮”。

3. 双链离散模型：构建与机理

为了模拟光子晶体平板中平面波与导模的相互作用，我们需要一个能同时支持传播波和倏逝波的“环境”。在连续系统中，这对应于在某个频率区间内，系统既存在辐射模式（传播衍射级次），也存在指数衰减的模式（倏逝衍射级次）。在离散模型中，最简化的实现方式是使用两条平行的、无限长的质量-弹簧链。

3.1 模型设定与运动方程

我们的模型如图1所示（概念图）：两条无限长的链（n=0, 1），每条链由一系列质量为1的珠子（m = …, -1, 0, 1, …）组成，珠子间用张力为1的弹簧连接。此外，每条链上的珠子还通过弹簧与一个固定的锚点连接。所有珠子被限制在垂直于链平面的方向运动（出平面振动）。

关键扰动参数：

环境不对称性 (δ)：通过改变两条链与各自锚点之间的弹簧张力引入。具体地，对于 m ≠ 0 的珠子，上链（n=0）锚点弹簧张力增加 δ，下链（n=1）锚点弹簧张力减少 δ。当 δ=0 时，环境完全对称。
散射体不对称性：位于 m=0 位置的两个珠子充当散射体，其质量 M₀ 和 M₁ 可以不同于标准质量1。当 M₀ = M₁ 时，散射体对称。

设 U_mn(t) 为位置 (m,n) 处珠子的位移。在谐波振动假设 U_mn(t) = u_mn e^{-iωt} 下，我们可以得到离散化的运动方程。对于 m ≠ 0 的均匀链部分，方程为：[ω² - 4 + (-1)^n δ] u_mn + u_{m+1,n} + u_{m-1,n} + u_{m,ñ} = 0其中 ñ = 1-n，表示另一条链上的对应位置。对于 m=0 处的散射体，方程为：(M_n ω² - 4) u_{0n} + u_{1n} + u_{-1n} + u_{0ñ} = 0

这个方程组描述了一个耦合的双链波动系统。参数 δ 控制了链间耦合的对称性，而 M₀ 和 M₁ 定义了局部缺陷（散射体）的特性。

3.2 均匀链的解：传播模与倏逝模

在远离散射体（m ≠ 0）的均匀区域，我们寻找形如u_mn = q_n h_m的解。代入运动方程后，可以发现存在两类解，对应两种“谐波”：

传播谐波：其空间依赖为h_m = e^{iκm}，其中 κ 为实数。这对应于能量可以沿链无限传播的模式，类比于光子晶体平板中的辐射模式。
倏逝谐波：其空间依赖为h_m = e^{-β|m|}，其中 β > 0。这对应于指数局域的模式，能量在空间上衰减，类比于平板中的导模衰减尾场。

通过分析特征方程可以证明，在一个特定的频率区间内（例如，当 δ=0 时，对应 1 < ω² < 3），系统恰好支持一个传播谐波和一个倏逝谐波。这正是模拟嵌入式束缚态所需的关键环境：一个模式可以携带能量至无穷远（连续谱），而另一个模式则被局域。在对称情况（δ=0）下，传播谐波中两条链同相运动（q₀/q₁ = 1），而倏逝谐波中两条链反相运动（q₀/q₁ = -1）。

3.3 嵌入式束缚态（陷获模）的存在条件

一个“陷获模”是指在没有任何入射波的情况下，系统自身支持的一个非零、且能量完全局域在散射体附近的稳态解。数学上，这要求当入射波振幅 I^p = 0 时，方程组存在非零解，且该解在 |m| → ∞ 时指数衰减（即只包含倏逝谐波部分）。

将只包含倏逝谐波的形式u_mn = B^e q_n^e e^{-β|m|}代入 m=0 处的边界条件（运动方程），我们得到关于 B^e 的齐次方程组。要使 B^e 有非零解，其系数矩阵的行列式必须为零。这给出了一个关于 ω 和 δ 的条件方程。对于给定的散射体质量 M₀ 和 M₁，该方程决定了唯一一对参数 (δ₀, ω₀)，使得系统在该参数下支持一个陷获模。重要的是，这个频率 ω₀ 位于系统连续谱的内部，因此它是一个“嵌入式”束缚态。

例如，当散射体对称（M₀ = M₁ = 2）且环境对称（δ₀ = 0）时，可以解析求得陷获模频率 ω₀ = √(21/10) ≈ 1.449。这个频率确实位于均匀链的传播带（1 < ω² < 7）内。

注意：嵌入式束缚态的存在是共振的前提，但它本身不对应于散射实验中的峰或谷。它是一个“幽灵”状态，只有在参数精确匹配时才存在。一旦参数偏离 (δ₀, ω₀)，这个状态就“溶解”到连续谱中，其物理效应通过共振显现出来。

4. 散射问题与传输异常的计算

现在，我们考虑散射问题：从左侧（m = -∞）入射一个传播谐波I^p q_n^p e^{iκm}。这个入射波会被 m=0 处的散射体散射，产生向左的反射波A^p q_n^p e^{-iκm}和向右的透射波B^p q_n^p e^{iκm}。此外，在散射体附近还会激发起指数衰减的倏逝波场A^e q_n^e e^{βm}(m<0) 和B^e q_n^e e^{-βm}(m>0)。

通过将总波场表达式在 m=0 处满足运动方程，并利用位移和力的连续性条件，我们可以求解出散射系数 A^p, A^e, B^p, B^e。对于这个无耗散系统，能量守恒要求 |A^p|² + |B^p|² = |I^p|²。我们关心的透射系数T 定义为透射能量与入射能量之比：T = |B^p|² / |I^p|²。

通过解析求解这个 2x2 的线性方程组，我们可以得到透射系数 T 作为频率 ω 和环境不对称参数 δ 的显式函数。计算结果表明，当系统参数 (δ, ω) 接近陷获模参数 (δ₀, ω₀) 时，透射系数 T 会在 ω₀ 附近表现出急剧的变化，即传输异常。

4.1 对称性破缺对共振的调控

图4展示了不同散射体不对称性（M₀）下，透射谱随频率 ω 和环境不对称性 δ 的变化。这里固定 M₀ + M₁ = 4，因此 M₀ 的大小直接衡量了散射体的不对称程度。

情况 (a): M₀ = 2.0 (对称散射体)。此时陷获模对应 δ₀ = 0。当环境也对称（δ=0，黑线）时，透射谱在 ω₀ 处没有尖峰，而是平滑的（实际上，此时系统存在真正的束缚态，不参与散射）。当环境不对称性 δ 略微偏离0（例如 ±0.1，黄/绿实线），一个尖锐的、以 ω₀ 为中心的共振峰立即出现。峰和谷对称地分布在 ω₀ 两侧。
情况 (b)-(d): M₀ > 2.0 (不对称散射体)。例如 M₀=2.2 时，陷获模存在于 δ₀ ≈ 0.43 处。此时，如果我们固定环境不对称性 δ 为0（黑线），系统远离其束缚态条件，共振很弱。只有当 δ 接近 0.43 时，尖锐的共振才会出现。更重要的是，随着 δ 从 δ₀ 偏离，共振峰的中心频率会发生移动（失谐），不再固定在 ω₀。例如，在图4(c)中，当 δ 从 δ₀ 增加或减少时，共振峰的中心分别向更高或更低的频率移动。

这个现象直观地展示了结构不对称性如何作为“调谐旋钮”：散射体的不对称性（M₀）决定了在哪个环境不对称性（δ₀）下存在嵌入式束缚态；而实际的环境不对称性（δ）相对于 δ₀ 的偏离量\tilde{δ} = δ - δ₀，则直接控制了共振的中心频率和宽度。在光子晶体平板的语境下，δ 类比于平行于平板的波矢分量（或入射角），而 M₀ 类比于平板结构的几何或材料不对称性参数。

5. 通用公式：统一描述共振线型

离散模型和连续的光子晶体平板系统，其共振行为可以用一个统一的解析公式来描述。这个公式基于对非鲁棒束缚态附近进行解析摄动推导而来，具有普适性。对于我们的离散模型，该公式给出了在(δ, ω)接近(δ₀, ω₀)时，透射系数 T 的近似表达式：

T(δ, ω) ≈ 1 / [1 + D²]其中D = [r₀ |\tilde{ω} + ℓ₁\tilde{δ} + r₂\tilde{δ}²|] / [t₀ |\tilde{ω} + ℓ₁\tilde{δ} + t₂\tilde{δ}²| (1 + α\tilde{ω})]

这里，\tilde{ω} = ω - ω₀,\tilde{δ} = δ - δ₀。公式中的参数有明确的物理意义：

r₀, t₀: 在束缚态点(\tilde{δ}, \tilde{ω}) = (0, 0)处的相对反射和透射振幅，满足r₀² + t₀² = 1。它们决定了共振的对比度。
ℓ₁:线性失谐系数。它控制了共振中心频率随\tilde{δ}的线性移动：ω_r ≈ ω₀ + ℓ₁\tilde{δ}。当系统结构对称时（δ₀=0），有 ℓ₁=0，共振峰以 ω₀ 对称。当对称性破缺时，ℓ₁ ≠ 0，导致共振峰失谐。
r₂, t₂:二次项系数。它们决定了共振峰的宽度（半高全宽）随\tilde{δ}的变化关系，宽度正比于\tilde{δ}²。峰和谷的位置分别由方程\tilde{ω} + ℓ₁\tilde{δ} + r₂\tilde{δ}² = 0和\tilde{ω} + ℓ₁\tilde{δ} + t₂\tilde{δ}² = 0给出，因此峰谷间距为(r₂ - t₂)\tilde{δ}²。
α: 控制背景透射斜率的参数。

这个公式的强大之处在于，它用少数几个由系统本身决定的参数（ℓ₁, r₀, t₀, r₂, t₂, α），就精确刻画了共振线型的核心特征：中心频率、宽度、峰谷位置和不对称程度。图5展示了改变 ℓ₁ 如何影响共振峰的位置，直观体现了不对称性导致的失谐效应。

5.1 与法诺公式的联系

经典的Fano公式可以看作是上述通用公式的一个特例。当满足以下条件时：

ℓ₁ = 0 （共振无线性失谐，结构对称）。
r₂ 和 t₂ 为实数（共振能达到0%和100%的透射极值）。
背景平坦（α=0）。

通用公式将退化为Fano形式。此时，共振宽度 γ 与\tilde{δ}²成正比，即γ ∝ \tilde{δ}²，而法诺参数q = sgn(t₂) * (r₀/t₀)。这建立了离散模型中扰动参数\tilde{δ}与Fano共振特征宽度之间的直接联系，本质上是一种费米黄金定则的体现，将泄漏模复频率的虚部（衰减率）与共振宽度联系起来。

实操心得：在利用此类共振设计光学器件（如传感器）时，这个通用公式是强大的设计指南。如果你希望共振峰对某个物理量（如折射率变化、位移）敏感，就应该致力于让该物理量的变化主要影响 ℓ₁ 或 δ₀，从而产生显著的线性失谐。如果你更关注共振的锐度（高Q值），则需要关注 r₂ 和 t₂，并尽量使系统工作在\tilde{δ}较小的区域，因为宽度随\tilde{δ}²增大而展宽。

6. 非线性散射体的拓展：多稳态与连续共振

将模型从线性推广到非线性，可以揭示更丰富的物理现象，如光学双稳态，这在全光开关和逻辑门中具有应用潜力。在离散模型中，可以通过在散射体珠子（m=0处）的势能中加入一个四次项（如 Kerr非线性项β|U_{0n}|⁴）来引入非线性。

在线性系统中，对于一个固定的入射频率和强度，透射率是唯一确定的。但在非线性系统中，由于介质的折射率（或等效的恢复力）依赖于光强，系统的响应可能呈现多值性，即双稳态或多稳态：对于同一个入射频率和强度，可能存在多个稳定的透射状态。

在我们的双链非线性离散模型中，数值计算发现了比单链Anderson非线性模型更复杂的行为。不仅在与线性束缚态频率 ω₀ 相近的频率处会出现共振和非线性双稳态，在从 ω₀ 一直延伸到单传播谐波频率区间上限的连续频率范围内，都可能出现多个谐振态（见图6）。这似乎与双链结构提供了更丰富的模式空间有关。对于固定的入射光强和非线性系数，透射率-频率曲线可能出现多个分支，对应不同的场分布和透射率。

这种现象暗示，基于嵌入式束缚态的非线性光子结构可能提供更宽的频率调谐范围和更复杂的多稳态行为，为非线性光子器件设计提供了新的可能性。当然，这需要进一步研究来理解其内在机制，并探索连续介质系统（如非线性光子晶体平板）是否也存在类似行为。

7. 模型优势、应用与展望

7.1 离散模型的优势总结

物理直观：用质量、弹簧、张力等力学概念类比电磁参数，使得抽象的光学模式耦合过程变得可视化和可触摸。
解析可解：所有方程均为线性或可处理的非线性方程，关键参数（如束缚态频率、透射系数）均可获得解析或半解析解，避免了纯数值模拟的“黑箱”问题。
参数清晰：明确分离了环境不对称性（δ，类比入射角/波矢）和散射体不对称性（M₀，类比结构扰动）两个调控维度，清晰地揭示了它们对共振中心频率和宽度的独立与耦合影响。
普适性桥梁：其共振行为受前述通用公式支配，该公式同样适用于连续的平板系统，因此离散模型的结论可以定性地、甚至定量地指导真实光子器件的设计。

7.2 在光子器件设计中的潜在应用思路

高灵敏度传感器：利用共振峰对参数 δ（可对应环境折射率）的敏感移动。例如，将待测物折射率变化映射为 δ 的变化，通过追踪共振峰位移实现传感。由于嵌入式束缚态共振通常非常尖锐（高Q值），其灵敏度可以很高。
可调谐滤波器/开关：通过主动改变 M₀ 或 M₁（例如，利用电光、热光或载流子色散效应改变局部折射率），可以动态地移动或打开/关闭共振峰，实现滤波或开关功能。
非线性光学元件：基于第6节讨论的非线性效应，可以设计低阈值的光学双稳态开关或逻辑门。入射光强本身可以改变系统的有效参数，从而改变共振条件，实现光控光。
对称性工程：通过精心设计结构的对称性（M₀与M₁的关系，以及δ），可以精确控制法诺参数q，从而获得所需的线型，是完美吸收、临界耦合等特殊功能的基础。

7.3 常见问题与概念辨析

问：嵌入式束缚态和普通的缺陷态或带边态有何区别？答：普通缺陷态的频率通常位于光子带隙中，不与任何传播模式频率重叠，因此是鲁棒的，微扰只会轻微移动其频率而不会使其“泄漏”。带边态位于带边，群速度为零，但仍在连续谱内。嵌入式束缚态则更深地嵌入连续谱内部，其存在完全依赖于精确的对称性，非鲁棒性是其最显著特征，也是产生尖锐共振的根源。
问：在实验中，如何判断观察到的共振是源于嵌入式束缚态，而不是其他机制（如法布里-珀罗共振）？答：一个关键特征是对对称性的极端敏感性。对于法布里-珀罗共振，其频率主要取决于腔长和有效折射率，对结构的细微不对称不敏感。而对于嵌入式束缚态共振，其出现和位置强烈依赖于结构的对称性参数（如入射角、结构左右/上下对称性）。通过微小地打破对称性，观察共振峰是否出现、是否发生显著频移，可以作为一个判断依据。此外，其线型通常符合非对称的法诺线型。
问：这个双链模型能否推广到更复杂的二维或三维光子晶体？答：核心思想是通用的。在更高维度的系统中，需要构建能同时支持至少一个传播模式（在连续谱内）和一个或多个指数局域模式的离散模型。这可能需要更多数量的耦合链或更复杂的格点结构来模拟多维度的布洛赫波矢。但描述共振线型的通用公式（第5节）由于其推导基于一般的摄动理论，对于满足能量守恒、存在非鲁棒束缚态的系统仍然适用。
问：模型中忽略了损耗（材料吸收、辐射损耗），实际器件中损耗的影响有多大？答：损耗是实际器件中必须考虑的因素。它会“抹平”共振特征，降低透射峰的峰值（达不到100%），抬高透射谷的谷底（达不到0%），并加宽共振线宽（降低Q值）。在通用公式中，可以引入一个复频率或复耦合系数来唯象地描述损耗。通常，只要损耗率远小于共振的固有宽度（由非鲁棒性决定，即Im(ℓ₂)˜δ²），共振的基本特征仍然清晰可辨。设计高Q值共振是抑制损耗影响的关键。

8. 从理论到设计：一个简化的设计流程参考

虽然实际的光子器件设计需要借助专业的电磁仿真软件（如 COMSOL, Lumerical FDTD, MEEP 等），但基于此离散模型的思想，我们可以梳理出一个概念性的设计流程：

目标定义：明确所需共振的中心频率、带宽（Q值）、线型（对称洛伦兹或非对称法诺）以及调谐方式（通过入射角、结构参数还是外部刺激）。
对称性分析：设计一个在目标频率附近支持嵌入式束缚态的基础对称结构。对于光子晶体平板，这通常意味着设计一个具有特定对称面（如垂直于平板的镜面对称）的周期性结构，使其在某个特定入射角下支持一个真正的导模。
引入可控的不对称性：确定一个易于调节的物理参数作为“扰动旋钮”。这可以是：
- 几何不对称：改变空气孔的大小、形状或位置。
- 材料不对称：在结构的不同区域使用不同的材料或通过外部场改变局部折射率（电光、热光效应）。
- 激励不对称：改变入射光的偏振或角度（对应模型中的 δ）。
参数扫描与优化：使用仿真软件，在基础对称结构上，系统地改变不对称性参数，计算其透射/反射谱。观察共振峰的出现、移动和线型变化。将仿真结果与通用公式T ≈ 1/(1+D²)进行拟合，可以提取出 ℓ₁, r₂, t₂ 等关键参数，量化不对称性的影响。
考虑非线性（如需要）：如果设计非线性器件，需在材料模型中引入三阶非线性极化率（χ⁽³⁾），并研究在不同入射光强下的响应，寻找双稳态区域和开关阈值。
容差与鲁棒性分析：评估制造误差（如尺寸偏差）对共振性能的影响。由于嵌入式束缚态对参数极其敏感，可能需要通过优化设计，使工作点位于对某些误差不敏感的区域，或引入反馈机制进行动态调谐。