⚛️ 量子物理

Realistic simulation of universal multiport interferometers with losses within a Gaussian formalism

本文提出了一种高斯形式化方法，用于逼真地模拟具有非均匀光学损耗的通用多端口干涉仪，从而实现对噪声 NISQ 光子平台上的量子算法结果进行准确的表征与预测。

原作者： Sebastiano Corli, Dario Chemoli, Samuele Altilia, Edoardo Suerra, Dario Tamascelli, Simone Cialdi, Enrico Prati

发布于 2026-10-14

📖 1 分钟阅读🧠 深度阅读

原作者： Sebastiano Corli, Dario Chemoli, Samuele Altilia, Edoardo Suerra, Dario Tamascelli, Simone Cialdi, Enrico Prati

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在构建能够解决当今计算机无法处理的问题之机器的竞赛中，光已经成为了一个领先的竞争者。与其他依赖微型电子电路或捕获原子的方法不同，这些光子计算机使用被称为光子的单个光粒子来携带信息。其优势在于光的自然特性：它传播速度快，产生的热量少，并且可以通过镜片和透镜进行操纵以执行复杂的计算。然而，这条路径面临着一个单一且持久的障碍：损耗。当光子穿过计算机芯片中脆弱的玻璃通道时，一些光子不可避免地会消失，被材料吸收或散射到空气中。这种消失不仅仅是一个小故障；它是阻碍这些机器发挥其全部潜力的主要障碍。为了了解我们距离开发出可用设备还有多远，科学家必须能够利用经典计算机对这些机器进行模拟，但当系统不完美时，这样做是非常困难的。如果模拟假设光表现得完美无瑕，它预测的结果可能是真实的、有缺陷的机器永远无法产生的。

一组研究人员开发了一种新的建模方式，用于模拟这些不完美的机器，为光在系统中泄漏的行为创建了一个现实的蓝图。他们的工作专注于一种特定的光学计算机，这种计算机使用由微型干涉仪（一种可以分裂并重新组合光束的装置）组成的网格来处理信息。在一个理想的世界里，这些设备会将每一颗光子精确地引导到它需要去的地方。但在现实中，玻璃以及它们之间的连接都是不完美的，这导致每一阶段都会有一部分光发生损耗。研究人员意识到，试图追踪每一颗丢失的光子及其逃逸到的环境，会使计算机模拟的计算量过于庞大而无法运行。相反，他们设计了一种方法，将所有这些损耗压缩成一种易于处理的数学形式。通过将整个光路径网络视为一个与简化环境相互作用的单一系统，他们可以在不需要模拟无限种丢失能量的可能性时，计算出光的最终状态。

他们发现的核心技术是一种允许他们将复杂的、混乱的光路径网络描述为一个统一的、描述光如何变化的单一映射的方法。他们证明了即使当光的损耗是不均匀的——即某些路径比其他路径更脏或更受损——整体效应也可以通过一组特定的数字来捕捉，这些数字描述了系统的“模糊性”或不确定性。这种方法适用于两种不同的构建此类计算机的方式：一种是光同时流经许多并行导线，另一种是光在回路中循环往复、来回反弹。研究人员通过在经典计算机上运行数千次模拟来测试他们的方法，模拟了这些光学设备在各种条件下的行为。他们发现，虽然光的行为确实受到机器特定布局的影响，但新模型可以准确预测实验结果，包括在末端检测到特定数量光子的概率。

该研究的一个重要发现挑战了科学家们长期以来使用的一种快捷方式。多年来，许多研究人员为了简化计算，一直假设光学损耗是均匀的，即损耗在机器的每个部分都以完全相同的速率发生。新的模拟表明，这种假设通常是不正确的。在一种常见的设计中，即光通过三角形网格组件传输时，损耗是不均匀的；损耗会根据光在网络中的位置而变化。研究人员表明，假设均匀损耗会导致对机器输出的预测出现错误，尤其是在机器规模变大时。相比之下，他们的新模型考虑了损耗的具体且不均匀的性质，能够提供关于机器实际行为的更准确图像。这至关重要，因为如果科学家想要建立纠错系统来修复这些错误，他们首先需要知道这些错误究竟是什么样子的。

该团队还探索了这些损耗如何影响被称为高斯玻色采样（Gaussian boson sampling）的计算结果，这是一项旨在证明量子计算机可以超越经典计算机的任务。他们模拟了将挤压态光（一种不确定性在某一方向上降低的特殊状态）通过一个有损网络的处理过程。他们观察到，随着光在不完美的机器中传输，检测到光子的概率发生了剧烈偏移。在充满噪声的环境中，没有检测到光子的真空态比在完美环境中变得更加普遍。此外，他们发现特定的损耗模式会改变检测到特定数量光子的可能性，从而打破了完美机器中存在的对称性。通过将他们的现实模型与简化的均匀损耗假设进行对比，他们可以清楚地看到两者在哪里产生分歧。随着光路径数量的增加以及损耗程度的恶化，这种差异会变得越来越大，从而证实了简化视图对于大规模设备而言是不够的。

这项工作为光学量子计算的未来提供了重要的工具。它提供了一种方法，可以将关于特定芯片如何丢失光线的真实世界数据直接代入模拟，以观察量子算法在该特定硬件上的表现。这意味着工程师不再需要猜测机器的行为，而是可以在制造之前以高度真实的程度测试其设计。研究人员建议，该框架可以适配任何类型的连续变量光学架构，无论其使用的是空间路径还是基于时间的回路。通过准确表征这些损耗，科学界可以更接近于开发有效的纠错技术，而这对于构建容错量子计算机至关重要。这项研究并不声称已经解决了光损耗问题，但它提供了一张清晰且现实的地形图，指出了陷阱究竟在哪里以及如何绕过它们。这种清晰度是实现从理论上的光子量子计算承诺向实际应用转化的必要步骤。

标题：基于高斯形式化方法的具有损耗的通用多端口干涉仪的现实模拟

问题陈述
光子技术是展示量子优势的主要竞争者，但光子损耗是实现容错光学量子计算的主要障碍。准确表征这些损耗对于在噪声后端上对量子算法进行经典模拟至关重要。目前的模拟方法通常依赖于一个先验假设，即所有通道的光学损耗是均匀分布的。作者认为，这种假设可能会导致对特定硬件上算法结果的预测出现偏差。此外，虽然存在针对基于比特（qubit）设备的纠错码，但在连续变量（CV）量子计算中尚未建立成熟的纠正光子损耗的技术。现有的用于噪声环境下连续变量模拟的软件库与基于物质的比特（如 Qiskit, Cirq）相比仍显不足。通过为每一次相互作用引入外部环境模来模拟损耗的方法，会导致模式数量呈多项式级增长，使得在 Fock 空间中的模拟在计算上变得不可行。

方法论
作者提出了一种高斯形式化方法，用于表征和模拟多端口干涉仪中的光学损耗，该方法适用于空间分箱（spatial-bin）和时间分箱（time-bin）编码架构。该方法流程如下：

理论框架： 系统使用海森堡基底（正交分量 $\hat{q}$ 和 $\hat{p}$ ）而非 Fock 基底进行描述。高斯态的动力学由均值向量和协方差矩阵决定。光学损耗被建模为物理模与外部环境模之间的相互作用，表现为分束器操作。
不完美干涉仪的表征： 作者基于 Qi 等人 [52] 对不完美马赫-曾德尔干涉仪（MZI）的表征进行了扩展。他们通过引入输入端、输出端以及分束器内部的损耗因子来模拟一个有损干涉仪 $T_\eta$ 。
奇异值分解 (SVD)： 为了避免追踪大量外部环境模导致的计算爆炸，作者采用了 SVD 方法。由不完美干涉仪网格产生的非幺正算符 $A$ 被分解为 $A = U_2 \sqrt{H} U_1$ ，其中 $U_1$ 和 $U_2$ 是幺正矩阵， $\sqrt{H}$ 是包含耗散因子的对角矩阵（奇异值）。这把所有的损耗通道压缩到了 $M$ 个外部模中（其中 $M$ 是物理模的数量），将系统规模从多项式增长降低到可控的 $2M \times 2M$ 协方差矩阵。
协方差矩阵推导： 作者推导出了受这些损耗影响的高斯玻色采样（GBS）过程的协方差矩阵 $\sigma_S$ 的解析表达式。该矩阵是通过从全局辛变换中部分迹掉（partially tracing out）环境模来获得的。
数值模拟： 作者使用高斯后端（Python 库 PNN）实现了一个流水线来模拟 GBS。他们对比了：
- 理想的 GBS 输出。
- 使用推导出的非均匀损耗矩阵 $\sqrt{H}$ 的现实 GBS 输出。
- 假设均匀损耗（即 $\sqrt{H} \propto I$ ）的 GBS 输出。
  他们测试了 Clements 和 Reck 两种干涉仪网格分解方案。

核心贡献

解析形式化方法： 本文提供了一个通过 SVD 推导出的、用于描述执行不完美干涉仪的 GBS 过程协方差矩阵的闭式表达式（公式 24）。这使得在不借助 Fock 空间的情况下，能够在经典后端模拟噪声 GBS 和一般的检测型（dyne-based）算法。
高效的模式缩减： 通过利用 SVD，该方法将模拟损耗所需的辅助模数量从对电路深度的多项式依赖降低到了线性依赖（ $2M$ ），显著降低了计算成本。
验证均匀性假设： 本研究通过数值评估了假设光学损耗均匀性的有效性。研究表明，这种假设并非普遍成立，且高度取决于干涉仪的分解方案（Clements 与 Reck）以及模式数量。
特定设备表征： 该框架允许将有关特定硬件缺陷（例如，单个 MZI 的特定损耗率）的实验数据直接集成到模拟流水线中。

结果

分布偏移： 数值模拟显示，光学损耗会增加真空态的占据数，并允许检测到在理想挤压态输入中被禁止的奇数光子计数，从而与理想情况相比显著改变了概率分布。
分解方案依赖性： 研究发现，损耗矩阵 $\sqrt{H}$ $H$ 与单位矩阵的偏离程度（即非均匀性）在不同架构下表现不同：
- 在 Reck 方案（等效于基于环路的时间分箱编码）中，与单位矩阵的距离随模式数 $M$ 的增加而增加。作者将这一趋势拟合为幂律函数和对数函数，发现幂律拟合效果略好。
- 在 Clements 方案（空间编码）中，这种距离随 $M$ 的变化呈现振荡并收敛到一个有界范围内（在 $M=5$ 之后），这表明其对均匀性的偏离更为稳定。
Hafnian 散度： 在模拟矩阵的 Hafnian 采样时，作者观察到，相比于均匀损耗近似模型，来自现实非均匀损耗模型的模拟结果更快地偏离了理想情况。这种差异在光学损耗增加（透射率 $\eta$ 较低）以及模式数增加时变得更加显著。

意义与主张
作者声称这项工作为未来特定设备和架构的表征提供了“蓝图”。所提出的框架能够对噪声中等规模量子（NISQ）光学平台上的量子算法进行现实模拟。通过精确建模非均匀损耗，该方法旨在防止因假设均匀损耗分布而可能导致的对算法结果的错误预测。作者将这项工作定位为将现实噪声模型与误差缓解技术相结合，并可能开发未来连续变量计算纠错码的基础性一步。他们明确指出，本工作并非提出新的纠错码，而是一个用于更好地理解和表征噪声格局的模拟工具。

类似论文