Simulating closed- and open-system dynamics of Tavis--Cummings lattices on fault-tolerant quantum computers

技术摘要：在容错量子计算机上模拟 Tavis–Cummings 点阵的封闭与开放系统动力学

问题陈述
Tavis–Cummings 点阵（TCL）模拟了每个包含多个双能级发射器的耦合腔阵列。这些系统承载着极化激元多体动力学，其实时相关性的经典模拟计算成本极高，尤其是在二维情况下。虽然对于具有相同发射器的封闭系统，置换对称性和激发守恒可以减少希尔伯特空间，但当腔内发射器具有不同的弛豫或退相干速率时，这种对称性在“原子分辨”场景下的开放系统中会失效。在这种情况下，状态空间仍然庞大到无法处理（例如，一个包含每腔四个发射器、光子截断值 $n_{max}=15$ 的 $4 \times 4$ 点阵，即使经过对称性简化，单个状态向量也需要约 23 TB 的存储空间）。此外，现有的经典工具（如量子蒙特卡洛或张量网络）在处理二维实时动力学或置换对称性被破坏的情况时表现不佳。本文旨在解决在容错量子计算机上模拟此类原子分辨 TCL 的逻辑资源需求问题，目标涵盖了针对光子损耗、发射器弛豫和退相干的封闭系统谱函数以及开放系统瞬态统计。

方法论
作者构建了一个基于**阶梯算符块编码（LOBE）和广义量子信号处理（GQSP）**的容错模拟框架。

1. 哈密顿量块编码：

   • 光子跳跃（Photon Hopping）： 跳跃项 $H_{hop}$ 通过在现有系统寄存器中“原地”更新源和目标光子占据数来进行编码。这与需要第二个光子寄存器来存储更新配置的稀疏访问实现方式不同。跳跃系数 $\sqrt{n_j(n_i+1)}$ 被分解为源函数和目标函数，通过从共享查找表中加载，并利用可逆模运算移位和对角算符进行应用。
   • 局部交换与对角项： 光-物质交换项（$H_{int}$）和对角项（$H_{diag}$）使用共享的角度表进行编码，该角度取决于光子占据数，从而避免了需要 $LN$ 个独立多路复用旋转的情况。
   • 组合： 跳跃块和局部块通过线性组合酉算符（LCU）结合，以形成完整哈密顿量 $H_{TCL}$ 的块编码。作者证明，移除均匀极化激元势能和单位偏移可以显著降低归一化因子 $\alpha$，从而减少所需的行走步数。

2. 封闭系统动力学：

   • 实时演化 $e^{-iH\tau}$ 使用 GQSP 进行近似。其成本由将指数函数近似到目标精度所需的多项式阶数决定，该阶数取决于归一化演化时间 $\alpha |\tau|$。
   • 目标观测量包括迟滞腔谱函数，通过涉及光子产生和湮灭算符的干涉相关测量进行估计。

3. 开放系统动力学：

   • Lindblad 生成器使用**采样一阶幺正扩张（sampled first-order unitary dilation）**进行模拟。算法不是相干地在所有 $J = L + 2LN$ 个跳跃算符中进行选择，而是根据其概率权重在每个时间步长内对单个跳扑算符进行采样。
   • 采样的生成器被扩张为一个作用于扩大系统（增加一个扩张量子比特）上的厄米特哈密顿量。这避免了大型相干跳跃索引寄存器的开销。
   • 目标观测量包括瞬态光子统计（平均占据数 $\bar{n}_i$ 和二阶相关函数 $g^{(2)}_i(0)$），通过在采样演化后进行计算基采样来估计。

主要贡献与结果

• 原地光子更新： 作者引入了一种电路构造，可以直接更新系统寄存器中的光子占据数。与“配对双寄存器”实现（使用第二个寄存器存储更新后的配置）相比，原地更新法消除了 $L \lceil \log_2(n_{max}+1) \rceil$ 个逻辑量子比特。对于一个 $16 \times 16$ 的点阵（$L=256$），这节省了 1,024 个量子比特。
• 门计数减少： 在三个均匀系数示例（$L=16, 64, 256$）中，假设匹配了共享准备、精度和工作空间约定，原地控制行走比最便宜的配对双寄存器实现在每步使用的 Toffoli 门数减少了 4.75% 至 8.11%。
• 资源评估：
  • 封闭系统： 本文提供了估计迟滞谱函数的 Toffoli 门计数和逻辑量子比特的具体界限。对于特定的实例（$L=16, N=1, n_{max}=3$），移除均匀势能可以将归一化因子 $\alpha$ 降低约两倍，从而显著降低模拟成本。
  • 开放系统： 作者推导了估计瞬态统计所需的时间步数和总门复杂度。他们指出，一阶采样法需要大量的步数（$N_{step} \propto T^2/\epsilon$）来控制采样误差。对于特定示例，估计成本约为 $1.9 \times 10^{20}$ 个 T 门（不包括状态准备），这主要由一阶时间步长决定。
• 验证： 作者通过小规模实例的门级模拟和独立的矩阵计算验证了其算符恒等式和电路分解。他们确认了采样扩张在时间步长方面具有一阶收敛性，而精确哈密顿层与耗散层的对称分裂则具有二阶收敛性。

意义与主张
本文声称提供了一个评估模拟原子分辨 TCL 容错成本的具体参考，其中对状态准备、测量和逻辑门计数进行了明确的核算。

• 关于优势的适度主张： 作者明确指出，其分析仅估算了量子成本，并未建立量子优势。他们承认，由于局域性，乘积公式（Trotter-Suzuki）在特定实例下可能提供更具实际意义的低成本，尽管信号处理方法（GQSP）在精度缩放方面表现更好。
• 开放系统的局限性： 本文强调，用于开放系统的第一阶采样扩张与高阶或对称分裂方法相比，在时间步长数量上产生了巨大的开销。虽然采样法通过简化电路（去除了相干跳跃路由）简化了过程，但作者指出，确定净资源节省情况需要编译精确的局部通道电路，这属于未来的研究工作。
• 测量开销： 研究强调，测量重复次数（샷噪声/shot noise）是总成本的一个主要因素，特别是在平均光子占据数较低时估计 $g^{(2)}$。作者建议，方差感知估计器或振幅估计可以减少这些计数，但并未在主要成本分析中实现它们。

总之，这项工作为模拟原子分辨 TCL 的逻辑资源分析提供了严谨的框架，证明了原地更新在量子比特数和门复杂度方面比稀疏访问方法具有切实的优势，同时也为高精度模拟开放系统动力学所需的巨大资源提供了基准。