Friedrichs模型：量子系统与连续谱耦合的理论与应用

1. Friedrichs模型基础与物理背景

Friedrichs模型是量子光学和开放量子系统研究中一个极具代表性的理论框架，它描述了离散量子态与连续能谱之间的相互作用。这个模型最初由K.O. Friedrichs在1948年提出，用于研究量子场论中的共振现象。经过数十年的发展，它已成为分析量子系统与环境耦合的基础工具。

1.1 模型的核心物理图像

Friedrichs模型的核心在于描述一个或多个离散能级与连续能谱的耦合。这种耦合会导致离散态变成具有有限寿命的共振态，表现为量子系统的自发衰变。从物理上看，这个过程类似于一个激发态的原子通过发射光子衰变到基态，或者量子点中的激子通过辐射复合消失。

在N级Friedrichs模型中，系统哈密顿量可以表示为三个部分的组合：

• 离散系统的本征能级部分
• 连续谱的自由哈密顿量部分
• 离散态与连续态之间的耦合项

这种分解使我们能够清晰地看到系统各组成部分的贡献以及它们之间的相互作用机制。

1.2 因子化相互作用的优势

本文研究的Friedrichs模型采用了因子化相互作用假设，即耦合强度可以分解为两部分：与离散态相关的因子f_n和与连续态相关的因子g(ω)。这种假设带来了两个关键优势：

1. 解析可解性：因子化形式使得许多物理量（如自能Σ(E)）可以表示为相对简单的积分表达式，便于理论分析。

2. 普适性：这种形式涵盖了广泛的物理系统，包括原子-光子耦合、电子-声子相互作用等常见场景。

在实际应用中，因子化假设对应于系统与环境耦合的"旋转波近似"，即忽略高频振荡项，这在大多数量子光学系统中都是合理的近似。

2. 束缚态的存在条件与分类

束缚态是Friedrichs模型中最重要的物理现象之一，它们的存在会显著改变系统的动力学行为。根据能量位置的不同，束缚态可以分为连续谱外的束缚态和连续谱中的束缚态(BIC)两类。

2.1 连续谱外束缚态的判定准则

判断束缚态是否存在的核心方程是： K(E) = Σ^{-1}(E)

其中K(E)是离散态性质的函数，Σ^{-1}(E)则反映了连续谱的特性。这个方程的图解解法非常直观：在能量轴上寻找两条曲线的交点。

具体判定时需要考虑以下关键因素：

1. 离散能级ϵ_n在连续谱外的数量N_out
2. 函数K(E)在连续谱边界处的值
3. 自能函数Σ(E)在连续谱边界处的行为

特别值得注意的是，当连续谱上边界延伸到无穷大时，系统只能在连续谱下方形成束缚态。这种情况在光子晶体波导等系统中很常见。

2.2 连续谱中束缚态(BIC)的形成机制

连续谱中的束缚态(Bound States in the Continuum, BIC)是一种反直觉的量子现象，它们存在于连续能谱中却保持局域化不衰变。BIC的形成需要两个必要条件：

1. 在能量E处谱密度J(E)=0
2. 同时满足束缚态方程(9)

在物理实现上，BIC可以通过以下方式产生：

• 连续谱中存在点状能隙
• 特定离散态与连续谱的耦合为零(g(E)=0)
• 系统具有特殊的对称性

BIC在量子信息存储、高Q值光学腔设计等方面有重要应用价值。

3. 系统动力学与衰变行为

Friedrichs模型的动力学演化展示了开放量子系统的丰富物理现象。通过投影算子方法，我们可以解析地求解系统的含时演化。

3.1 散射态与Lippmann-Schwinger方程

当系统能量E落在连续谱内时，我们需要考虑散射态的解。这时有效哈密顿量变为非厄米的： H_eff^±(E) = Σ_nϵ_n|n⟩⟨n| + Σ^±(E)Σ_{n,n'}f_nf_{n'}^*|n⟩⟨n'|

其中自能函数Σ^±(E)可以分解为实部(能移)和虚部(衰变率)： Σ^±(E) = Δ(E) ∓ iΓ(E)

这种分解直接反映了系统与环境的相互作用导致的能级移动和衰变加宽。

3.2 存活概率的演化特性

初始处于离散态的系统的存活概率p(t)展现了丰富的动力学行为：

1. 短时行为：通常表现为非指数衰变，这与量子芝诺效应相关
2. 长时极限：由束缚态主导，可能出现：
   • 完全衰变（无束缚态）
   • 稳定平台（存在BIC）
   • 持续振荡（多个束缚态共存）

特别有趣的是，当系统存在多个束缚态时，存活概率会表现出量子拍频现象，振荡频率由束缚态能级差决定。这种行为在量子信息存储中可用于设计相干控制协议。

4. Markovian极限与非厄米物理

当系统与环境的耦合较弱且环境关联时间很短时，系统会进入Markovian区域，此时动力学可以用能量无关的非厄米哈密顿量描述。

4.1 非厄米哈密顿量的性质

在Markovian极限下，有效哈密顿量变为： H_eff^± = Σ_nϵ_n|n⟩⟨n| ∓ iΓΣ_{n,n'}f_nf_{n'}^*|n⟩⟨n'|

这个哈密顿量有两个重要特征：

1. 本征值一般为复数，虚部代表衰变率
2. 左右本征矢不同，需要采用双正交基表示

这种非厄米性质导致了许多新奇现象，如异常点(Exceptional Points)附近的特殊行为。

4.2 动力学分类

根据非厄米哈密顿量的本征态结构，系统的衰变动力学可分为三类：

1. 简单指数衰变：所有本征态非简并
2. 多指数衰变：多个衰变通道竞争
3. 幂律修正的指数衰变：在异常点附近

特别地，当哈密顿量出现高阶异常点时，存活概率会表现出t^{N-1}e^{-γt}的独特行为，这在量子控制中可能有特殊应用。

5. 原子链-光子晶体波导耦合系统

作为Friedrichs模型的具体实现，我们考虑一个原子链耦合到光子晶体波导的系统。这个平台不仅展示了丰富的物理现象，还具有实验可行性。

5.1 系统建模与参数映射

该系统的紧束缚模型哈密顿量为： H = -λΣ_{μ=1}^{N-1}(|μ⟩⟨μ+1|+h.c.) - κΣ_{ν=1}^∞(|ν⟩⟨ν+1|+h.c.) + ξ(|1⟩⟨l|+h.c.)

通过Bloch态展开，可以将其映射到Friedrichs模型框架：

• 离散能级ϵ_n = -2λcos(πn/(N+1))
• 连续谱ω(k) = -2κcos(k)
• 耦合函数g(ω)包含波导耦合位置l的信息

这种映射使得我们可以直接应用前面发展的理论工具。

5.2 束缚态调控与设计

通过调节系统参数，我们可以精确控制束缚态的产生：

1. 耦合位置控制：改变原子链与波导的耦合点l可以调制谱密度J(ω)的零点位置
2. 能带工程：调节κ/λ比值可以移动离散能级相对于连续谱的位置
3. 对称性利用：特定参数下(如l=N-1, κ=λ)可实现全BIC状态

图3展示了不同参数下束缚态数量的相图，为实验设计提供了明确指导。

5.3 动力学行为实例分析

图4比较了三种不同耦合构型下的存活概率演化：

1. 完全衰变：l=1时无束缚态，系统完全弛豫
2. 部分局域：l=2时存在BIC，存活概率趋于非零常数
3. 持续振荡：l→∞时存在两个连续谱外束缚态，导致长时振荡

这些行为验证了理论预测，并展示了Friedrichs模型应用的广泛性。

6. 非厄米对称性与量子调控

原子链-波导系统为实现非厄米量子物理提供了理想平台，特别是对PT对称性的研究。

6.1 反PT对称哈密顿量的实现

在Markovian极限下，两原子系统(