有损高斯玻色采样能否突破平方根区间被经典模拟？

平方根区间限制的是一种方法

场景是标准的有损高斯玻色采样。N 个具有共同固定压缩幅度 r 的压缩真空态，加上 M−N 个真空模式，经过光子存活概率为 η 的均匀纯损耗和一个无源干涉仪 V，再进入光子数分辨探测器。目标 P 是全部 M 个输出模式上的完整光子计数分布，精度用总变差距离 TV(P,Q) 衡量，即计数概率绝对差之和的一半。平均存活光子数为 N_out = Nη·sinh²r。

一种广泛使用的经典方法把每个有损输入替换为最接近的经典高斯态，其余部分精确处理。在固定压缩和固定精度下，只要 N_out 增长不快于 N 的平方根，该方法就保持准确。论文的第一个结果表明，这不是选择高斯参考态造成的假象。定理 1 构造了显式干涉仪——把成对等压缩输入送入分束器，其余输出留在真空——在该族上，任何具有非负 Glauber–Sudarshan P 表示的状态，无论是否高斯、无论模式间如何关联，其总变差误差至少为 1 − exp(−c_r Nη²)。在该族上，P 经典集合内的最优误差为 Θ_r(min{1, Nη²})。

机制来自经典计数概率上的一个符号约束。对强度为 I 和 J 的乘积相干态，配对概率满足 q₁₁ − q₂₀ − q₀₂ = −½e^(−I−J)(I−J)² ≤ 0，任何 P 经典混合态都继承这一不等式。有损配对违反它，p₁₁ − p₂₀ − p₀₂ = n₀η²/[1 + n₀(2η − η²)]² > 0，量级为 η²。把该约束跨配对求和并应用 Cauchy–Schwarz 不等式，就把局部违反放大为全局下界。

适用范围与这个界同样重要。定理 1 是关于一种方法的陈述，不是复杂度下界：它不覆盖高斯但非经典的参考态，也不覆盖张量网络方法保留的非经典核心，更没有说明分布本身的困难程度。它排除的是任何先替换输入为经典光、再进行探测的算法。

近似探测为何能越过障碍

由于这一障碍是近似类别的性质，而不是经典计算的性质，保持输入精确、近似探测的方法不受它约束。论文把光子数分辨探测改写为一个虚拟吸收过程：每个输出模式耦合到一个完美吸收体，吸收体一次移除一个光子并记录该光子来自哪个模式。一直运行到所有光子被吸收，模式 i 中记录的移除次数就恰好是它的光子数，因此记录精确复现 PNR 分布。虚拟时间只用于给点击排序，最后被积分掉；实验中不施加任何操作。

在这一过程中，条件态遵循两条规则。没有点击发生时，它按 e^(−t n̂/2) 衰减，其中 n̂ 是总光子数算符——即尚未移除任何光子的条件态。模式 i 发生点击时，它变为 b_i ρ b_i†/λ_i，即一次光子减除，λ_i = ⟨b_i† b_i⟩ 是该模式的速率。两条规则都不是幺正操作。对高斯态两者都有闭式表达：记 Σ = I/2 + K，无点击演化给出 K̇ = −K − K²，一次点击把 K 更新为 K + Δ_i，其中 Δ_i = KΠ_iK/λ_i，λ_i = ½Tr(Π_iK)。下一次点击的等待时间也有闭式表达，S_K(t) = det[I + (1 − e^(−t))K]^(−1/2)，其极限 det(I + K)^(−1/2) 是当前态的真空概率，也就是没有光子剩余的概率。

采样器只存储 K。从 K₀ = ηO diag(p_r I_N, −d_r I_N)Oᵀ 出发，其中 O 是 V 的实化，p_r、d_r 由 r 确定，它抽取记录是否已经结束，否则抽取下一次点击的时间，按速率 λ_i 的比例选择被点击的输出，并施加更新。每个输出上的点击构成样本。无点击步骤是精确的；唯一的近似是点击之后，把已不再高斯的精确态替换为协方差同为 I/2 + K + Δ_i 的中心高斯态。不需要计算 hafnian，光子数和协方差都没有上限。

误差按 Nη³ 而非 Nη² 缩放

定理 2 表明，对每个固定的 r > 0、每个 M ≥ N 和每个无源等距 V，都存在一个多项式期望算术代价的随机经典算法，其输出分布 Q 满足 TV(P, Q) ≤ min{1, C₃(r)Nη³}，0 ≤ η ≤ 1，其中 C₃(r) 与 M、N、V 无关。该定理适用于上述压缩和真空输入，对干涉仪深度或连通性没有限制，且包含输出碰撞。

三次方标度有一个简单读法。低透射率下每个输出约含 O(η) 个光子，因此一次减除使态改变 Δ_i = O(η)，在典型记录上高斯替换每次点击的代价为 O(η²)。点击的总速率约为 Nη，乘积就是 Nη³。证明还覆盖了协方差变化很大的稀有记录，使用位移 Y_t = K_t + b e^(−t)P，其中 b = η/(1−η)、P = OOᵀ，以及非负泛函 Φ_t = Tr(Y_t K_t²)。不等式 Y_t ⪰ 0、0 ⪯ Δ_i ⪯ 2Y_t 和 Σ_i λ_iΔ_i = K_t² 依赖压缩输入的成对正交结构，而不是任何输出占据数的下界。

代价是多项式的，但并不小。直接实现的期望代价为 O(MN + (1 + Nηn₀)(N³ + MN²))，内存为 O(MN + N²)，其中 S_K 的标量求逆计为基本操作，期望点击次数至多为 Nηn₀。误差界中出现因子 N 而非 M，来自被占据的协方差子空间。

Nη³ 换来多少存活光子

在固定压缩和固定精度下令 Nη³ = O(1)，就允许 η = c₀N^(−1/3)，其中 c₀ 是足够小的固定常数，因此平均存活光子数 N_out 正比于 N^(2/3)。这是相对平方根区间的主要改进，后者由 Nη² = O(1) 强制 N_out = O(√N)。

两个定理在同一目标族上形成分离。取 η = N^(−a)，1/3 < a < 1/2。此时 Nη² → ∞ 而 Nη³ → 0，于是每个经典光计数分布都与目标渐近最大程度地远离，而探测近似采样器保持准确。这是在同一族干涉仪上的直接分离，并不声称所有有损 GBS 实例在该损耗下都可经典模拟。

另一套配对构造走得更远。使用配对产生与放大，并对少数配对扇区精确处理，总变差距离被界定为 (3/4)√(N x⁴)，其中 x = (n₀ + 1)η，因此对固定 r 和精度，η 正比于 N^(−1/4) 时平均存活计数正比于 N^(3/4)。该构造适用于每个输入模式压缩量相同的情形。论文明确指出，有限多个精确扇区的界本身并不能确立 N^(3/4) 指数，因为在该区间核心配对数量随 N 增长。

探测近似在哪里失效

高斯重置并非一致地小。如果位移强度 |α_i|² 与核心占据数 n_c,i 相当或更小，重置可能产生量级为一的误差：在零位移处它把 S(u)|1⟩ 替换为压缩热态，迹距离为 3/4。在初始位移系综中，该区域概率量级为 √η，那里的探测速率量级为 η²，提示贡献量级为 η^(5/2)，但论文指出这些初始速率并不能界定整条轨迹的误差。

位移变体 Q′_G 把精确位移分解与高斯重置结合，保留采样得到的位移，并在每次点击后同时匹配均值和协方差。定理 2 适用于仅处理协方差的采样器 Q_G，不适用于该变体。条件矩更新中的负项 −v_i v_iᵀ/λ_i² 使相对跳跃界和三次方预算无法直接用于位移情形，论文也没有给出可能消除单噪声正交限制的双正交噪声变体的误差分析。

数值检验相应地范围很窄。它们使用最多三个模式中的一到两个压缩输入，r 取 {0.3, 0.5, 1}，η 取 {0.1, 0.2, 0.3}。在 η = 0.5 时，Q_G 在 r ≈ 1.1 以上比经典参考态更不准确。对 η = 0.3 的随机三模式干涉仪，参考态、Q_G、Q′_G 和 Q_pair 的经验误差在 r = 1 时为 (0.058, 0.058, 0.014, 0.015)，在 r = 1.5 时为 (0.062, 0.136, 0.04, 0.080)。对 Q′_G，两个经验距离都与采样噪声相当，无法分辨真实误差。这些少模式检验既不能确立一致误差界，也没有与张量网络方法比较运行时间。

采样器以均值偏低换取精度

一项数值检验显示了该近似在具体、可测量方式上的失效。对 r = 1、η = 0.3 的单模式，10⁶ 条记录的运行给出均值 0.3795 ± 0.0010，而精确值为 0.4143；超过 16 次点击的概率为 1.3 × 10⁻⁴，而目标为 2 × 10⁻⁶。采样均值偏小与理论结果一致，后者只给出期望计数的上界而非等式。论文指出，这些数值并不能确立一致可积性的失效。

公式本身独立于数值检验得到验证。约化协方差表示、求和跳跃恒等式、三次方生成元和局部代价不等式在精确有理算术下对 18 个高斯前缀和 54 个探测器跳跃进行了检验，其中 42 个跳跃不对易，复有理等距取 (N, M) = (1, 2)、(2, 3) 和 (2, 4)。在 18 个单模式 Fock 计算中，把维度从 40 提高到 56 使局部迹距离变化至多 1.4 × 10⁻¹⁶。解析界独立于这些检验得到证明，而这些检验并不能严格界定截断或舍入误差。

什么决定该方法是否适用

对任何评估有损 GBS 优势声明的人来说，实际问题在于正在做哪种近似，因为两种选择的天花板不同。如果方法把有损输入替换为经典光，相关阈值是 Nη² = O(1)，在显式干涉仪族上平均存活光子数最坏情况下不能超过 O(√N)，无论经典参考态是高斯还是非高斯且带关联。如果方法保持输入精确并近似探测，阈值变为 Nη³ = O(1)，允许 N_out ∝ N^(2/3)，对等压缩配对构造则为 N^(3/4)。

分离在 1/3 < a < 1/2 且 η = N^(−a) 的窗口中最清晰，此时经典光误差趋于最大，而探测近似误差消失。正是在这个窗口，而不是仅凭渐近指数，两种方法能在同一目标族上被区分开。

使用这些结果时应始终带上两条边界。下界是关于经典光替换方法的陈述，不是复杂度下界，因此不覆盖高斯但非经典的参考态或保留的非经典核心。上界则是对仅处理协方差的采样器 Q_G 证明的；在小规模检验中经验误差更低的位移变体 Q′_G 没有对应定理。因此论文中最具决策价值的比较不是指数本身，而是阈值与方法的配对：输入替换对应 Nη²，探测近似对应 Nη³，而 N^(3/4) 配对构造仅限于等压缩输入。

本文引用的文献