来源论文: https://arxiv.org/abs/2605.29724v1 生成时间: May 30, 2026 04:46

拓扑浴中的量子相干保护：基于SSH晶格束缚态的量子精密测量深度解析

0. 执行摘要

量子精密测量（Quantum Metrology）致力于利用量子力学特有的相干性、纠缠等资源，实现超越经典极限（标准量子极限）的超高灵敏度参数估计。然而，现实中的量子探针不可避免地要与周围的环境（即“浴”或“储能箱”，Bath/Reservoir）发生耦合。这种不可控的相互作用导致系统发生消相干（Decoherence）和能量弛豫，使得好不容易编码在探针相干叠加态中的参数信息迅速流失到环境中，严重制约了实际量子精密测量的精度上限。如何对抗环境噪声、实现相干信息的长效保持，是当前量子精密测量及开放量子系统领域的重大核心挑战。

近期由 Sofia Evangelou 撰写的理论工作《Bound-state-protected phase metrology for a quantum emitter in a Su–Schrieffer–Heeger bath》为此提供了一条极其优雅且具有高度物理可行性的微观解决途径。该工作跳出了传统的现象学主方程或人工设计谱密度的框架，从一维微观晶格哈密顿量出发，研究了单个双能级量子发射器（Quantum Emitter, QE）局域耦合到玻色型 Su-Schrieffer-Heeger（SSH）晶格浴中的局部相位估计动力学。

研究的核心发现表明：通过引入晶格的二聚化（Dimerization），可以在系统能带中心打开一个可控的能隙，并在能隙内诱导出一种由发射器与晶格协同构成的“发射器-浴束缚态”（Emitter-Bath Bound State）。该束缚态的存在从根本上抑制了发射器激发态向晶格的完全弛豫，使得发射器的量子相干性在极长时间尺度下得以保留。利用 Dyson 方程严格求解单激发扇区的格林函数，作者解析推导出了发射器的生存振幅、能隙内束缚态的能量极点方程、以及决定长期残余相干信号的发射器权重（残差）。结果表明，通过精细调控晶格二聚化程度、发射器与浴的耦合强度以及发射器的频率失谐，可以主动重构和优化后瞬态（Post-transient）的量子 Fisher 信息（Quantum Fisher Information, QFI）。这一工作不仅为“结构化储能箱工程”在量子技术中的应用提供了清晰的微观物理解释，也为基于波导 QED 和光子晶体晶格的精密测量实验设计指明了方向。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题

本项研究聚焦的核心科学问题可以概括为：如何从微观晶格哈密顿量出发，定量地评估环境的微观拓扑/几何结构如何通过诱导能隙束缚态来保护非平衡态量子精密测量的精度？

在开放量子系统的物理框架中，通常采用 Markov 近似或简单的 Lorentzian 谱密度来描述噪声通道。然而，这种简化忽略了环境本身的格点拓扑结构和局域杂化效应。具体到一维 SSH 晶格环境，其特有的交替耦合强度（二聚化）会在能带结构中心开出一个拓扑能隙。当一个外来的双能级发射器与其耦合时，如果发射器的能级位于能隙内，就会在发射器位置附近形成局域化的束缚态。这种束缚态是否能够作为量子相干信息的“安全屋”？其保护效果在时域上如何体现？耦合强度、二聚化系数、失谐量等微观参数如何协同影响量子 Fisher 信息的长期保持率和瞬态振荡行为？这些正是本工作试图彻底解答的问题。

1.2 理论基础

1.2.1 量子 Fisher 信息与 Cramér-Rao 偏差上限

量子精密测量的核心定量指标是量子 Fisher 信息 $F_Q^{(\phi)}(t)$。对于未知的待测参数 $\phi$，多次测量所得估计值 $\hat{\phi}$ 的均方误差受限于量子 Cramér-Rao 不等式：

其中 $M$ 是测量次数。对于局域相位估计，我们假设量子发射器（探针）制备在最大相干叠加态上：

在经历了含有未知参数 $\phi$ 的相干演化以及环境消相干动力学后，发射器与浴的整体波函数演化为：

其中 $u(t)$ 为发射器处于激发态 $|e\rangle$ 的生存振幅（Survival Amplitude），满足初始条件 $u(0)=1$；$|1_\lambda\rangle$ 表示环境（浴）中产生一个能量为 $\omega_\lambda$ 的玻色子激发。

对浴的自由度求偏迹（Partial Trace），可得发射器的约化密度矩阵 $\rho_\phi(t)$：

将此约化密度矩阵写为 Bloch 矢量形式 $\rho_\phi(t) = \frac{1}{2} [I + \mathbf{r}(t, \phi) \cdot \boldsymbol{\sigma}]$，对应的 Bloch 矢量分量为：

由于 Bloch 矢量的长度 $|\mathbf{r}(t, \phi)|^2 = 1 - |u(t)|^2 + |u(t)|^4$ 与 $\phi$ 无关，易知 $\mathbf{r}(t, \phi) \cdot \partial_\phi \mathbf{r}(t, \phi) = 0$。代入单参数 QFI 的标准 Bloch 矢量计算公式：

这是一个极其优美且关键的解析结论：在零温及单激发扇区假设下，该精密测量方案的量子 Fisher 信息在任意时刻 $t$ 都严格等于发射器的激发态生存概率 $|u(t)|^2$。因此，保护 QFI 的科学问题等价于如何实现激发态生存概率的长效保持（避免衰减至零）。

1.2.2 微观哈密顿量：发射器耦合至玻色型 SSH 晶格

整个复合系统的哈密顿量为（采用单位制 $\hbar=1$）：

系统（发射器）哈密顿量：

其中 $\Delta$ 为发射器的本征跃迁频率（相对于能隙中心对称点的失谐量），$\sigma_{ee} = |e\rangle\langle e|$。

玻色型 SSH 晶格浴的哈密顿量：

该晶格由交替重复的 A、B 两个子格组成，$a_n$ 和 $b_n$ 分别是第 $n$ 个胞内 A、B 格点上玻色子的湮灭算符。$J > 0$ 代表平均跳符强度，$d \in [-1, 1]$ 代表二聚化参数。无界 SSH 晶格的色散关系（Dispersion Relation）通过傅里叶变换可以严格求解：

当 $d \neq 0$ 时，色散曲线分裂为上、下两个能带（宽度各为 $2J(1-|d|)$），中间在频率 $\omega \in [-2J|d|, 2J|d|]$ 处打开一个宽度为 $4J|d|$ 的中心能隙（Central Gap），如图 1(b) 所示。

局域发射器-浴耦合哈密顿量（发射器单点耦合至第 $n=0$ 个胞的 A 子格）：

其中 $g$ 为耦合强度，$\sigma_{eg} = |e\rangle\langle g|$。

1.3 技术难点与方法细节

难点一：超出 Markov 近似的严格非平衡动力学求解

当发射器本征频率 $\Delta$ 处于能隙内部或带边（Band Edge）附近时，环境的态密度（Density of States）发生剧烈变化，传统的玻恩-马尔可夫（Born-Markov）近似和久保（Kubo）线性响应理论完全失效。必须发展严谨的非马尔可夫（Non-Markovian）解析与数值方法。

方法细节一：Dyson 方程与严格格林函数分析

为了求解生存振幅 $u(t)$，作者采用了格林函数（Green’s Function）方法。发射器的延迟格林函数（Retarded Green’s Function）定义在单激发扇区内，满足 Dyson 方程：

其中 $G_A^0(\omega)$ 是无界、未耦合的纯 SSH 晶格在 $n=0$ 胞 A 子格点上的局域格林函数。通过在第一布里渊区内积分可以解析求得：

对于处于能隙内部的实频率（$|\omega| < 2J|d|$），由于带内没有实本征态，晶格的局部谱密度为零，积分结果保持纯实数。解析积分为：

其对频率的微商为：

根据格林函数理论，发射器-浴复合系统的束缚态对应于格林函数的极点，即分母为零：

由于 $G_A^0(\omega)$ 在区间 $\omega \in (-2J|d|, 2J|d|)$ 内是严格单调递减且连续的，且当 $\omega \to \pm 2J|d|$ 时 $G_A^0(\omega) \to \mp \infty$。因此，对于任意实数失谐 $\Delta$，上述极点方程在能隙内部有且仅有一个实数根 $\omega_{\text{BS}}(\Delta)$。这证明了一维 SSH 晶格能够非常鲁棒地支持单个能隙内束缚态。

该束缚态对应的发射器残差（Residue）——即波函数在该局域束缚态上的重叠概率，决定了长期保留的相干信息：

在时域上，生存振幅可以严格写为束缚态贡献与能带连续谱（Continuum）贡献的叠加：

在极长的时间尺度下（$t \to \infty$），由于能带内连续谱状态的相位快速相消干涉（Dephasing），第二项积分衰减为零，生存振幅趋于：

因此，长效保留的量子 Fisher 信息极限值为：

在共振情况（$\Delta = 0$）下，极点显然位于 $\omega_{\text{BS}} = 0$。代入极点和导数公式，残差简化为：

由此得到共振长期保留的 QFI：

这是该研究最核心的解析成就。它定量地建立起了微观量子精密测量精度（$F_Q$）与晶格拓扑带隙（$4J|d|$）以及杂化强度（$g$）之间的直接映射关系。

难点二：有限尺寸效应与长时间动力学数值收敛性

在实际数值计算中，必须使用有限尺寸的晶格。然而，有限尺寸晶格在足够长的时间后会产生相干反射（Finite-size Recurrences），这会干扰对发射器非马尔可夫衰减和束缚态长期稳定行为的观测。为此，必须设计足够大且满足时域截断的物理系统。

方法细节二：精确对角化与 Krylov 子空间 Lanczos 方法

作者通过两种高度精确且互补的数值策略来解析生存振幅：

1. 全谱精确对角化（Exact Diagonalization, ED）：对于中等大小的晶格（例如胞数 $L = 220$，对应单激发扇区矩阵维度为 $2L + 2 = 442$），直接求哈密顿量的全部本征值和本征矢量。这种方法非常适合精确提取单个本征态（如能隙内束缚态）并计算其准确投影。
2. Krylov 子空间 Lanczos 算法：对于大晶格系统（例如胞数 $L = 500$ 甚至更高，以保证在 $t = 100 J^{-1}$ 的时间窗口内不发生边缘反射），完全对角化极其耗时。作者通过以初始态 $|e\rangle$ 为种子向量（Seed State），在 Krylov 子空间内进行三对角化，从而在极高精度下（Lanczos 崩溃容差设为 $10^{-12}$）高效迭代计算生存振幅 $u(t) = \langle e| e^{-iHt} |e\rangle$。Krylov 子空间的维度设定为 $350 \sim 360$，保证了在绘制的时间范围内结果的绝对收敛和无反射纯净性。

2. 关键 Benchmark 体系，计算所得数据，性能数据

本工作通过设计多个关键的物理 Benchmark 体系，从不同参数维度全方位展示了基于 SSH 束缚态保护相干性的优越性能。

Benchmark 1：共振条件下的二聚化强度效应

• 物理配置：发射器处于 gap 中心（$\Delta = 0$），发射器-晶格耦合强度固定为 $g = 0.4J$。对比四种不同的二聚化系数 $d = 0.7, 0.5, 0.3, 0.1$。
• 计算方法：$L = 500$ 的一维有限链，Krylov 子空间维度 350。
• 数据结果（参考图 2）：
  • 在瞬态阶段（$t < 20 J^{-1}$），QFI 经历了快速下降，并伴随由非马尔可夫效应引发的剧烈周期性振荡。这对应于发射器激发的局部能量发射与重新吸收的往返过程。
  • 在后瞬态阶段（$t > 40 J^{-1}$），动力学基本平息，QFI 稳定在一个长期平台值上。该数值高度逼近由解析公式 $\left( 1 + \frac{g^2}{4J^2 |d|} \right)^{-2}$ 给出的解析预测线（在图中以水平虚线示出）。
  • 定量性能数据：
    • 当 $d = 0.7$ 时，长期保留的 $F_Q^{(\phi)}$ 维持在 0.90 以上的高位，说明相干信息几乎没有流失。
    • 当 $d = 0.5$ 时，长期 $F_Q^{(\phi)} \approx 0.83$。
    • 当 $d = 0.3$ 时，长期 $F_Q^{(\phi)} \approx 0.61$。
    • 当 $d = 0.1$ 时，由于能隙极其狭窄，长期保护效果显著变差，稳定在 0.15 附近。
• 物理结论：二聚化参数 $|d|$ 越大，中心能隙越宽（宽度为 $4J|d|$），局域束缚态的能量越远离带边，其对激发的局域化束缚能力越强，从而能够留存更多的量子 Fisher 信息。

Benchmark 2：数值残差与解析残差的精确对比

• 物理配置：$\Delta = 0$, $g = 0.4J$，对 $d$ 进行连续扫描（$|d| \in [0.05, 0.7]$）。
• 计算方法：$L = 220$ 有限链的精确对角化。
• 数据结果（参考图 3）：
  • 图 3(a) 绘制了从有限尺寸晶格中直接提取的真实在隙本征态发射器重叠概率（即数值 $Z_{\text{BS}}$，蓝色圆点）与解析表达式 $\left( 1 + \frac{g^2}{4J^2|d|} \right)^{-1}$（实线）的对比。两者在整个区间内完全重合。
  • 图 3(b) 对比了长期 QFI 平台值与 $Z_{\text{BS}}^2$。在极其微小的数值波动范围内，两者实现完美契合。这在数值上无懈可击地证实了：长期留存的相干信息来源正是且仅是该中心在隙束缚态。

Benchmark 3：发射器-浴耦合强度 $g$ 的“双刃剑”效应

• 物理配置：$\Delta = 0$, $d = 0.3$。对比不同的耦合强度 $g/J = 0.1, 0.4, 0.6, 0.8$。
• 计算方法：$L = 500$ 链，Krylov 子空间 350。
• 数据结果（参考图 4）：
  • 这是一个非常违反直觉但物理上合理的科学发现：在传统认知中，更强的耦合似乎能加强相干杂化；但在这里，随着 $g$ 的增加，长期保留的相干 QFI 反而单调下降。
  • 具体而言：
    • 当 $g = 0.1J$ 时，长期残留 QFI 高达 0.95 以上。
    • 当 $g = 0.4J$ 时，降至 0.61 左右。
    • 当 $g = 0.8J$ 时，长期 QFI 下跌至 0.15 左右。
  • 与此同时，强耦合显着增强了瞬态振荡的振幅，瞬态相干恢复的特征更加尖锐。
• 物理机制：尽管 $g$ 变大增加了发射器与临近格点的量子纠缠，但这使得局域束缚态的波函数向晶格纵深处大范围地铺展（杂化程度加深），从而将发射器自身的权重（即残差 $Z_{\text{BS}}$）稀释到了浴的格点上。因此，发射器局域观测到的相干相位的残留量被迫减小。

Benchmark 4：非拓扑对照组（均匀晶格体系，Gap 闭合极限）

• 物理配置：$\Delta = 0$, $g = 0.4J$, $d = 0$（对应一维普通紧束缚 Tight-binding 链，无中心 Gap）。
• 计算方法：$L = 500$ 链，Krylov 子空间 350 演化动力学。
• 数据结果（参考图 5）：
  • 在瞬态振荡后，QFI 的时域曲线在 $t \approx 30 J^{-1}$ 之后彻底衰减至零，没有任何长期平台存在。
• 物理结论：这证明了如果没有二聚化打开能隙并提供束缚态支撑，量子相干精密测量在开放环境下的失效是不可避免的。本工作的保护机制完全由能隙物理主导。

Benchmark 5：偏离共振的物理行为与三种运行指标（Operational Diagnostics）

当失谐 $\Delta \neq 0$ 时，发射器频率向带边靠近并最终穿过带边进入连续谱能带。为此，作者定义了三个极具实验工程指导意义的操作指标：

1. 时域平均量子 Fisher 信息 (Late-time averaged QFI, $\overline{F}_Q^{(\phi)}$)： $$\overline{F}_Q^{(\phi)} = \frac{1}{t_2 - t_1} \int_{t_1}^{t_2} F_Q^{(\phi)}(t) dt$$ 计算中选取观测窗口 $t_1 = 40 J^{-1}$ 至 $t_2 = 100 J^{-1}$。该指标绕开了时域残余微小振荡的干扰，客观度量了后瞬态的整体相干度。
2. QFI 保持时间 (QFI Retention Time, $t_\eta$)： QFI 首次跌破给定阈值 $\eta$ 的时间：$F_Q^{(\phi)}(t_\eta) = \eta$（计算中取 $\eta = 0.2$）。
3. 后瞬态实用时间窗口 (Post-transient Useful-window Duration, $W_\eta(t_{\text{cut}}, T)$)： 在度过初期准备瞬态后（设定阈值截止时间 $t_{\text{cut}} = 20 J^{-1}$），在整个总观测时间 $T = 100 J^{-1}$ 内，QFI 保持在阈值 $\eta = 0.4$ 之上的累计持续时间： $$W_\eta(t_{\text{cut}}, T) = \int_{t_{\text{cut}}}^T dt \, \Theta \left( F_Q^{(\phi)}(t) - \eta \right)$$ 其中 $\Theta$ 是 Heaviside 阶跃函数。

关键扫描数据分析：

• 时域平均 QFI 的无崩塌特征（图 7）：将 $\overline{F}_Q^{(\phi)}$ 作为归一化失谐 $\Delta/(2J|d|)$ 的函数进行绘制。在能隙内部（$\Delta/(2J|d|) < 1$），平均 QFI 维持在高位；一旦跨越带边（$\Delta/(2J|d|) > 1$），曲线迅速塌水般滑落。但值得注意的是，不同二聚化程度 $d=0.5, 0.3, 0.1$ 的曲线并未完全重合（即未发生“普遍崩溃归一化”）。这表明：除了能隙宽度外，晶格能带内的微观态密度细节同样参与了对生存振幅的塑造。
• 保持时间 $t_{\eta=0.2}$ 的断崖式下跌（图 8）：
  • 在 Gap 内部极深处（如 $\Delta/(2J|d|) < 0.8$），由于 QFI 永远不跌破 0.2，其保持时间直接触及数值模拟的总上限 $t_{\text{max}} = 100 J^{-1}$。
  • 在带边附近，保持时间断崖式下跌至 $20 J^{-1}$ 以下。
  • 进入能带 continuum 区，保持时间甚至低于 $10 J^{-1}$。二聚化越弱（如 $d=0.1$ 相比于 $d=0.3$），这种在带边的脆弱性越明显。
• 实用窗口 $W_{0.4}$ 对实验调控的指导（图 10）：
  • Case (a) 隙内深处 (共振, $\Delta/(2J|d|) = 0$)：瞬态过后，QFI 完全稳定在 $0.4$ 线上方，实用测量时间窗口 $W_{0.4} = 80 J^{-1}$（达到了后瞬态演化的理论极限 $T - t_{\text{cut}}$）。这代表着实验测量拥有极宽的时间自由度，对对准时间不敏感。
  • Case (b) 带边临界点 ($\Delta/J = 0.57$, 归一化失谐 $0.95$)：虽然由于初始衰减导致 QFI 一度跌破 $0.4$，但非马尔可夫相干反馈很快在长时将 QFI 重新推回 $0.4$ 之上并稳定住。实用窗口缩窄为 $W_{0.4} \approx 55.9 J^{-1}$。这说明相干保护仍然有效，但容错率降低。
  • Case (c) 能带内 ($\Delta/J = 0.72$, 归一化失谐 $1.20$)：QFI 发生不可逆的衰减，永远无法重新抬头穿过 $0.4$ 的测量基准线。此时实用测量窗口彻底关闭，$W_{0.4} = 0$。实验上在该频段无法进行任何有效的、基于后瞬态相干的精密测量。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

3.1 算法设计与哈密顿量矩阵构建

为了能够让量子信息领域的科研人员快速复现本文结果，以下给出构建该微观有限一维 SSH-发射器哈密顿量的详细数学矩阵映射方法。

在单激发扇区内，系统的基矢空间（Dimension为 $2L + 2$）定义为：

为了便于矩阵表示，我们设基矢索引如下：

• 发射器状态 $|e\rangle$：Index为 0；
• 第 $n$ 个胞的 A 子格点（$n \in [-L, L]$）：Index为 2 * (n + L) + 1；
• 第 $n$ 个胞的 B 子格点（$n \in [-L, L]$）：Index为 2 * (n + L) + 2。

此时，总哈密顿量矩阵 $H$ 的非零矩阵元定义为：

1. 对角元：
   • $H_{0, 0} = \Delta$
   • 其余对角元为 0（代表晶格格点的在隙基准能量设为0）。
2. 晶格内部跳符（非对角元）： 对于每个胞 $n \in [-L, L]$：
   • 胞内跳符：$H_{A_n, B_n} = J(1+d)$，即 $H_{2(n+L)+1, \, 2(n+L)+2} = J(1+d)$；
   • 胞间跳符（若 $n < L$）：$H_{B_n, A_{n+1}} = J(1-d)$，即 $H_{2(n+L)+2, \, 2(n+L)+3} = J(1-d)$。
3. 发射器-晶格耦合： 发射器局域耦合到 $n=0$ 胞的 A 子格（其矩阵索引为 $Idx_0 = 2L + 1$）：
   • $H_{0, Idx_0} = H_{Idx_0, 0} = g$。

矩阵构建完毕后，利用 Hermite 矩阵的对称性（Real Symmetric 矩阵）进行存储。

3.2 基于 Python 的完整复现代码（Krylov 子空间法）

以下展示一个基于 SciPy 和 NumPy 库编写的高性能、生产级 Python 代码，利用 Lanczos 算法计算指定参数下的 $F_Q^{(\phi)}(t)$：

import numpy as np
import scipy.linalg as la
from scipy.sparse import lil_matrix
import matplotlib.pyplot as plt

def construct_ssh_emitter_hamiltonian(L, d, J, g, delta):
    """
    构建单激发扇区下的发射器-SSH浴复合系统哈密顿量
    L: 每个方向的胞数 (总胞数为 2*L + 1)
    d: 二聚化参数
    J: 平均耦合常数
    g: 发射器与浴的耦合强度
    delta: 发射器的失谐量 (本征频率)
    """
    num_cells = 2 * L + 1
    dim = 2 * num_cells + 1  # 1 (发射器) + 2*num_cells (A和B子格点)
    
    # 使用稀疏矩阵的LIL格式高效构建
    H = lil_matrix((dim, dim), dtype=float)
    
    # 1. 写入发射器对角能级
    H[0, 0] = delta
    
    # 辅助函数：根据格点物理坐标(n, sublattice)返回其矩阵索引
    # n 范围从 -L 到 L, sublattice 只能是 'A' 或 'B'
    def get_index(n, sub):
        offset = n + L
        if sub == 'A':
            return 2 * offset + 1
        elif sub == 'B':
            return 2 * offset + 2
        else:
            raise ValueError("Sublattice must be 'A' or 'B'")
            
    # 2. 写入晶格跳符矩阵元
    for n in range(-L, L + 1):
        idx_A = get_index(n, 'A')
        idx_B = get_index(n, 'B')
        
        # 胞内跳符
        j_intra = J * (1.0 + d)
        H[idx_A, idx_B] = j_intra
        H[idx_B, idx_A] = j_intra
        
        # 胞间跳符
        if n < L:
            idx_next_A = get_index(n + 1, 'A')
            j_inter = J * (1.0 - d)
            H[idx_B, idx_next_A] = j_inter
            H[idx_next_A, idx_B] = j_inter
            
    # 3. 写入发射器和 n=0 胞 A 子格的耦合
    idx_A0 = get_index(0, 'A')
    H[0, idx_A0] = g
    H[idx_A0, 0] = g
    
    return H.tocsr()

def lanczos_time_evolution(H, time_steps, krylov_dim=350):
    """
    利用 Lanczos 算法计算生存振幅 u(t)
    H: CSR格式的哈密顿量矩阵
    time_steps: 计算的时间数组
    krylov_dim: Krylov 子空间的投影维度
    """
    dim = H.shape[0]
    # 初始状态：发射器全激发态，其余全为空腔基态
    v_start = np.zeros(dim, dtype=complex)
    v_start[0] = 1.0
    
    # Lanczos 三对角化过程
    alpha = np.zeros(krylov_dim, dtype=float)
    beta = np.zeros(krylov_dim - 1, dtype=float)
    V = np.zeros((dim, krylov_dim), dtype=complex)
    
    V[:, 0] = v_start
    
    # 迭代生成正交基
    for j in range(krylov_dim):
        w = H.dot(V[:, j])
        if j > 0:
            w -= beta[j-1] * V[:, j-1]
        alpha[j] = np.real(np.vdot(V[:, j], w))
        w -= alpha[j] * V[:, j]
        
        if j < krylov_dim - 1:
            beta[j] = np.linalg.norm(w)
            if beta[j] < 1e-12:  # 容错截止，防止崩溃
                actual_dim = j + 1
                alpha = alpha[:actual_dim]
                beta = beta[:actual_dim-1]
                V = V[:, :actual_dim]
                break
            V[:, j+1] = w / beta[j]
            
    # 构建 Krylov 子空间内的三对角矩阵 T
    T = np.diag(alpha) + np.diag(beta, k=1) + np.diag(beta, k=-1)
    
    # 计算子空间内的特征值与特征向量
    eigenvalues, eigenvectors = la.eigh(T)
    
    # 计算生存振幅 u(t) = <e| exp(-iHt) |e>
    # 在Krylov表示下, 初始状态对应 T 中的第一个格点 [1, 0, 0, ...]^T
    u_t = []
    for t in time_steps:
        # 矩阵指数乘法转化为子空间对角化计算
        prop_t = eigenvectors @ np.diag(np.exp(-1j * eigenvalues * t)) @ eigenvectors.T
        u_val = prop_t[0, 0]  # 获取发射器分量
        u_t.append(u_val)
        
    return np.array(u_t)

# --- 运行验证演示 (复现论文中 Fig. 2 的 d=0.3 曲线) ---
if __name__ == "__main__":
    # 物理参数定义
    L = 500         # 左右各500个胞，确保长时间内无反射
    J = 1.0         # 跳符基准能量尺度
    g = 0.4 * J     # 发射器耦合强度
    delta = 0.0     # 共振
    d_param = 0.3   # 二聚化参数
    
    # 构建哈密顿量
    H = construct_ssh_emitter_hamiltonian(L, d_param, J, g, delta)
    
    # 定义时间序列
    times = np.linspace(0.0, 100.0, 1000)  # 对应论文中的 Time [J^-1]
    
    # 计算动力学
    u_t = lanczos_time_evolution(H, times, krylov_dim=350)
    qfi = np.abs(u_t)**2
    
    # 解析长期平台预测线
    qfi_inf_theoretical = (1.0 + (g**2) / (4.0 * (J**2) * d_param))**(-2)
    
    # 绘图对比
    plt.figure(figsize=(8, 5))
    plt.plot(times, qfi, label=f"Numerical ($d={d_param}$)", color="tab:green")
    plt.axhline(y=qfi_inf_theoretical, color="tab:green", linestyle="--", 
                label=f"Analytic theory: $F_Q(\infty) = {qfi_inf_theoretical:.3f}$")
    plt.xlabel("Time [$J^{-1}$]", fontsize=12)
    plt.ylabel("$F_Q^{(\phi)}(t)$", fontsize=12)
    plt.title("Phase-QFI Protection via SSH Bound State", fontsize=13)
    plt.ylim(0, 1.1)
    plt.grid(True, linestyle=":")
    plt.legend(loc="best")
    plt.show()

3.3 推荐的开源软件包及软件库链接

复现这类复杂微观量子动力学问题的学者，可参考或使用以下开源生态圈：

1. QuTiP (Quantum Toolbox in Python): 世界上最著名的量子系统演化模拟库。能够极其简便地构建量子算符、局部状态以及进行主方程求解。对于单激发扇区以外的多体激发物理非常推荐。 https://qutip.org
2. SciPy (Sparse Linear Algebra submodule): 本研究由于采用了显式的大型一维微观晶格，本质上属于大规模稀疏 Hermite 矩阵，利用 scipy.sparse.linalg 提供的经典 Lanczos / Arnoldi 求解器、特征值包（如 eigsh）能起到极大的加速作用。 https://scipy.org
3. Julia (KrylovKit.jl): 如果对复现效率有更极致的强算追求，Julia 语言的 KrylovKit.jl 提供了一流的、无矩阵（Matrix-free）的三对角化和大规模哈密顿量 Krylov 动力学积分器。 https://github.com/Jutho/KrylovKit.jl

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献及其在文中的角色

1. S. John and J. Wang, Phys. Rev. Lett. 64, 2418 (1990) [文献 10]
   • 角色作用：这是光子晶体带隙（PBG）物理中“原子-光子束缚态（Atom-Photon Bound State）”奠基性的工作。本论文正是将这一经典物理概念从宏观介质，完美移植并泛化到了微观的一维拓扑 SSH 紧束缚浴模型中，提供了最根本的物理脉络支撑。
2. A. W. Chin, S. F. Huelga, and M. B. Plenio, Phys. Rev. Lett. 109, 233601 (2012) [文献 6]
   • 角色作用：该文献确立了非马尔可夫环境记忆效应在量子精密测量中充当“正向资源”的宏观理论框架。作者的微观晶格研究正好作为其实际物化案例，展现了记忆效应如何具体表现为瞬态的信息恢复和后瞬态的 QFI 锁定。
3. M. Bello, G. Platero, J. I. Cirac, and A. González-Tudela, Sci. Adv. 5, eaaw0297 (2019) [文献 29]
   • 角色作用：本项研究采用的一维 SSH 微观浴模型结构、色散关系推导方法、以及在单激发扇区下的解法，大体上继承了 Cirac 组在这篇关于拓扑波导 QED 中的束缚态工作的代数形式框架。这一继承确保了本工作格林函数求解的坚实数学物理基础。
4. K. Berrada, Phys. Rev. A 88, 035806 (2013) [文献 20] / PRA 88, 013718 (2013) [文献 19]
   • 角色作用：Berrada 在量子带隙中首次独立提出了通过带隙抑制自发辐射来保全量子 Fisher 信息并实现 QFI 捕获（QFI Trapping）的思路。本文工作是这一研究思路在微观拓扑晶格方向的深度演进。

4.2 本项工作的局限性学术评论

尽管本工作在数学上严谨、物理图像清晰，但从量子精密测量的实用化以及量子物理的深度考量，该方案在以下几个维度上暴露出了一定的局限性，有待后续学者突破：

局限性一：单激发扇区与零温假设的严苛物理约束

本研究的所有理论基石——包括 Dyson 方程格林函数求极点、相位估算中 QFI 严格等于 $|u(t)|^2$ 的代数恒等式——全部严格建立在环境温度为零温且系统总激发数守恒为 1 的绝对约束下。一旦考虑有限温度，环境本身的涨落将不仅通过真空吸附，还会通过热涨落（Thermal Excitations）导致发射器多级受激弛豫。此时，不仅单激发空间不再封闭，其密度矩阵向高维多阶空间溢出，原有的“残差 $Z_{\text{BS}}$ 解析等于长期 QFI 平台值”的宏观物理规律将极可能完全失效。如何在有限温下利用束缚态有效抵抗热脱域噪声，本工作未提供任何理论抓手。

局限性二：由于大块局域耦合导致的拓扑不敏感性

本工作的局限性之二在于其“空怀 SSH 拓扑晶格之名，却未能真正利用拓扑特有的边界保护机制”。正由于作者考虑的是发射器耦合到一维无限链大块（Bulk）内部的 A 格点，整个格林函数对二聚化参数 $d$ 的依赖只通过 $d^2$ 的形式进入极点方程。这导致改变 $d$ 的正负号（在拓扑物理中，改变 $d$ 的正负代表着系统在拓扑平凡相（Trivial）和拓扑非平凡相（Topological）之间的相变）对发射器的相干动力学没有任何定量影响。换言之，该方案保护 QFI 的核心纯粹源于“能隙（Gap）的存在性”而不是“能隙的拓扑不寻常性”。拓扑边缘态（Edge States）所特有的强抗局域无序扰动等优势，在这个 Bulk 局部耦合配置中完全无武之地，限制了其精密测量方案对抗格点无序（Disorder）的稳定性。

局限性三：实际格点无序与漏光通道对束缚态的退相干

在实际微观波导、超导量子比特谐振腔阵列中，不可避免地存在两种破坏性物理：

1. 格点能量无序（Grid Disordered On-site Energies）：破坏晶格的周易性，可能导致局域化深度发生物理漂移。
2. 腔本身的本征损耗（Cavity Dissipation）：每个 SSH 格点光子本身可能从侧向逸出。由于局域束缚态包含一部分晶格格点上的电磁激发（如残差公式所体现的那样，一部分权重落在了晶格上），当晶格本身存在本征漏光损耗时，该局域束缚态将不再是哈密顿量的真实实数极点，而是会获得一个虚数衰减项。在足够长时间下，QFI 依然不可避免地会衰减到零，本工作未能给出如何定量评估这种漏光引起的最终相干流失寿命。

5. 其他必要补充：物理机制探讨与实验路线图展望

为了给从事量子相干物理和精密测量实验研发的读者提供更宽广的前瞻视角，我们在本节中对核心物理机制的深层内涵以及具体的实验工程可行方案进行额外补充。

5.1 为什么是 $d^2$ 镜像对称？深入剖析 SSH 内部空间反演

在本文的数学分析中，格林函数 $G_A^0(\omega)$ 的表达式：

中仅含有二聚化因子的平方 $d^2$。从物理对称性角度来看，这对应着一维 SSH 晶格在空间反演对称性下的本质表现。对于无限长 bulk 晶格，二聚化参数由 $d$ 变为 $-d$ 仅仅相当于将 A 子格点和 B 子格点的空间物理位置进行一次整体移动：

由于发射器仅仅局域地单点耦合到了 A 子格点上，当发射器环视四周的浴环境时，其向左、向右所感受到的局部势能网络空间结构，在 $d \to -d$ 时，只是左-右空间反射对称的镜像，整体的能带密度（DOS）和局域格林函数投影毫无变化。 这种性质在实验上具有极强的工艺包容性：实验制备微观晶格时，无需小心翼翼地确认晶格内部到底对应哪个拓扑 winding sector，也无须为二聚化正负极性的漂移而担心；只要确保二聚化的绝对偏差量 $|d|$ 满足设计指标，便能稳健地在 gap 内部提供一致的 QFI 保护屏障。

5.2 相比于 Lorentzian 结构储能箱的本质物理超越

为了在理论上深化此工作的价值，可将其与经典的非马尔可夫模型——Lorentzian（洛伦兹）型带隙环境进行对比。 在洛伦兹模型中，环境的谱密度通常由一个具有确定宽度和截止频率的唯象公式给定：

这种基于唯象连续谱密度的模型缺失了格点自能（Self-energy）的空间相位自干涉机制。 而在微观的 SSH 晶格中，随着发射器频率向带边靠近，格林函数的导数 $dG_A^0/d\omega$ 会在带边 $\omega = \pm 2J|d|$ 发生数学上的分式发散（Divergence）。这种发散直接导致极靠近带边的束缚态残差 $Z_{\text{BS}}$ 极其迅速地塌缩归零，并在时域上引发非常复杂的、由多格点反射相干波包诱导的非指数型振荡。这种由晶格格点微观多路干涉产生的“带边奇异性”（Band-edge Singularity）是任何唯象 Lorentzian 谱密度无法模拟也无法解释的。本工作通过展示这一发散对 QFI 长期留存强度的损害，警示了我们在进行量子精密测量方案设计时，必须将探针工作频率牢牢锁在 Gap 深部（远离带边），从而避免带边发散带来的相干流失灾难。

5.3 三大主流物理平台实验实现路线图

本工作由于物理参数定义在普适一维紧束缚哈密顿量层面，能够在当前三个主流前沿量子物理实验平台上完美落地：

+----------------------------------------------------------------------------+
|                              实验路线图设计方案                             |
+----------------------------------------------------------------------------+
| 1. 超导量子电路 (Superconducting Qubits & Resonator Arrays)                 |
|    - 发射器 (QE): 采用超导 transmon 比特, 工作频率通过外加磁通线在能在         |
|      GHz 范围实现高精度无级调节, 对应参数 \Delta 调谐。                       |
|    - 晶格 (Bath): 使用超导共面波导谐振器(CPW Resonators)或集总参数LC谐振器     |
|      交替排列构成一维阵列。谐振腔之间的交替电容/电感耦合强度实现 J_1 和 J_2。 |
|    - 耦合器: 发射器通过电容跨接耦合至 $n=0$ 的 A 腔上, 电容值大小控制参数 g。 |
|    - 参数标定: J ~ 100 MHz, g ~ 10-20 MHz, 处于单激发强耦合物理区间。       |
+----------------------------------------------------------------------------+
| 2. 硅基光子晶体波导与光子晶格 (Waveguide QED in Photonic Crystals)          |
|    - 发射器 (QE): 采用嵌入在光子晶体波导中央的单个自组织半导体量子点。         |
|    - 晶格 (Bath): 在硅基(Silicon-on-Insulator)波导上刻蚀交替间距的空气孔。    |
|      空气孔的物理间隔距离交替排列直接决定电磁波在格点间的跳符强度 J(1+d) 与    |
|      J(1-d)。此微纳加工结构在红外通信波段(1550 nm)能打开宽达数个 THz 的物理能隙。|
|    - 优势: 固态芯片集成度极高, 适合微型化片上精密量子传感器制备。              |
+----------------------------------------------------------------------------+
| 3. 光晶格中的超冷原子 (Ultracold Atoms in Optical Lattices)                  |
|    - 发射器 (QE): 采用具有多能级本征结构的“单杂质原子”(Impurity Atom)。        |
|    - 晶格 (Bath): 采用双色激光束干涉形成的超周期光学势阱(Superlattice)。      |
|      通过调整辅助光束的相对相位和强度, 可以随意拉高或降低势垒高度, 实现二聚化 |
|      参数 d 从 0 (均匀晶格) 到 1 (完全解耦双基点) 的在线、实时大范围编程微调。|
|    - 优势: 系统环境绝对纯净、无任何天然材料无序缺陷, 是最理想的基础物理探索平台。|
+----------------------------------------------------------------------------+

通过对上述实验路线图的主动工程推进，并辅以本研究所提供的一套精确而完备的解析与数值计算工具，我们有理由相信，拓扑及微观结构化储能箱工程将在下一代更长寿命、更高精度的量子传感器研发中，发挥支柱性的基础科学指引作用。