数字量子处理器能测量 Hubbard 准粒子散射与纠缠吗？

42 个量子比特上的两个准粒子对撞

这篇论文把量子处理器变成一台散射实验装置。研究对象是一维 Fermi-Hubbard 模型：费米子在相邻格点之间以振幅 J 跃迁，当两个费米子占据同一格点时还要付出在位相互作用能 U。在半满 Mott 背景上，相关的电荷激发是 doublon（同时占据两种自旋的格点）和 holon（空格点）。作者在 21 个格点链条的两侧各制备一个这样的激发，再让它们朝中心运动。

量子比特与格点的对应关系很具体：两个相邻量子比特构成一个 Hubbard 格点，每个量子比特对应一种自旋轨道，因此 21 个格点占用 42 个量子比特。作者采用带锯齿排序的 Jordan-Wigner 变换，使每个 Trotter 步只需最近邻门。为了让两个准粒子具有确定方向，演化被分成两段。在 0 < T ≤ 1 期间，指向较近链端的跃迁被关闭，每个激发因此形成只朝链中心传播的定向波包；在 1 < T ≤ 4 期间，所有跃迁恢复，两个波包对撞。作者报告，在 T ≤ 4 时，波包 93% 以上的权重集中在主导分量上。

由于初始 Neel 态并非哈密顿量的精确基态，每次测量都配一次参考运行：使用相同线路但不加激发，再减去参考值。带 Δ 前缀的量就是这种扣背景后的值。这一点影响对下文所有数字的解读：报告的密度和概率都是相对无激发基线的差值，而不是原始格点占据数。

相互作用强度决定透射还是反射

核心测量量是 doublon 透射概率 T，定义在 T = 3.2 时刻，即位于链中心另一侧（格点 i ≥ 11）的 doublon 数占总 doublon 数的比例。它是 doublon 的比例，不是总电荷的比例，而且只在碰撞后一个固定时刻取值，并非对整个演化过程积分。

在 U = 1 到 8 之间，实验显示出一个转变。U = 1 时，doublon 与 holon 波包互相穿过并保持原方向；U = 7 时，两者被反射回去。实测曲线与 Bethe-ansatz 两体结果 T(θ) = θ²/(1 + θ²) 相符，其中 θ = 2J(sin k₁ − sin k₂)/U，取最大相对群速度 c = 4J，且没有自由参数。若把 c 作为自由参数拟合同一形式，得到 c ≈ 3.6J；作者将其归因于波包携带围绕 k = π/2 的动量分布，因此平均相对速度低于最大值。

论文给出了强相互作用导致反射的机制。doublon 与 holon 若要交换位置，必须经过一个与初末态相差 U 的中间构型，因此交换振幅被压低到 J²/U 量级，两个准粒子表现得像硬核粒子。作者还指出，出射波包相对入射波包没有明显速度变化，他们将其解读为弹性碰撞，这与一维 Hubbard 模型的可积性一致。

电荷信号消失不等于对湮灭

在中间值 U = 3.5 时，碰撞后两个波包的电荷密度几乎消失。这是一个有歧义的信号：如果 doublon 与 holon 直接湮灭，也会出现同样现象。作者通过分别测量 doublon 密度 d_i = ⟨n_i↑ n_i↓⟩ 和 holon 密度 h_i = ⟨(1 − n_i↑)(1 − n_i↓)⟩ 来区分这两种可能。

质心轨迹给出了答案。U = 1 时，doublon 与 holon 的质心互相穿过；U = 7 时，两者折返回起始一侧；U = 3.5 时，两个质心都停留在链中心附近，同时两个分布都变宽。这正是透射与反射分量叠加 |ψ_out⟩ ∼ t|d_R h_L⟩ + r|d_L h_R⟩（其中 |t|² = T，|r|² = 1 − T）会产生的图像：当 T ≈ 1/2 时，两侧各含相当的 doublon 与 holon 权重，它们相反的电荷相互抵消。

第二项检验在 T = 3.2 时比较总电荷 Σ_i |Δ⟨n_i⟩| 与总 doublon 数 Σ_i Δd_i。总电荷在 U ≈ 3.5 附近出现凹陷，但 doublon 数随 U 单调上升。复合会同时降低两者，因此 doublon 数的单调性排除了以湮灭解释电荷凹陷的可能。

纠缠由粒子数涨落推断而来

纠缠结论依赖的观测量与透射概率不同。作者从粒子数涨落推断纠缠，用其方程 (4) 定义的数熵 S_num 量化，并报告 S_num 与纠缠在 U/J = 4 时达到最大，也就是透射与反射权重相当的取值。

论文明确说明了这一结论能成立到什么程度。实测数据表明两个分量存在及其相对权重，但没有直接确认二者之间的相干性。把它读作相干叠加依赖两点：时间演化是幺正的，且观测到的 U 依赖符合方程 (2) 的预测。作者表示，一旦能够模拟更长时间，对碰撞后的波函数再做一次干涉实验即可确认这一点。

这正是最影响引用方式的分界。散射概率测量是在硬件上与 Bethe-ansatz S 矩阵的直接比较；纠缠陈述则是从粒子数涨落加上幺正演化假设推断出来的，并不是对碰撞后状态的层析重建。

doublon-doublon 对照没有相互作用依赖

对照实验让两个 doublon 对撞，而不是一个 doublon 与一个 holon。论文报告，这种碰撞既没有出现 doublon-holon 情形中的相互作用依赖，也没有出现粒子数涨落，作者称这符合不可区分粒子的预期。

正是这一对比让 doublon-holon 的结果可以被解释。相互作用依赖的转变和粒子数涨落，只在两个对撞对象是带相反电荷、可区分的电荷激发时出现；当两个对象完全相同时，同一硬件上并不出现这些现象。没有这个对照，doublon-holon 运行中随 U 变化的信号就可能被归因于器件行为，而不是准粒子散射。

误差处理决定吻合结论的分量

与理论的比较有多强，取决于其背后的误差预算。读出和 CZ 错误率超过校准阈值的量子比特与耦合器被作为硬性约束排除，再对剩余耦合器做 Bell 态保真度测试，选出这条 42 量子比特链。每个数据点使用 70,000 次采样，分配到 70 个 Pauli twirling 实例，并结合动态解耦与读出 twirling，从而把相干门误差和非对称读出误差近似转换为随机比特翻转通道。

后处理再利用问题的对称性。总粒子数守恒且初始态的总粒子数已知，因此作者拟合测得的总粒子数分布来估计均匀比特翻转概率 p，并用它恢复无噪声的期望值。反演与粒子-空穴变换的组合对称性 IC 也被利用，激发和对称化都定义在前 20 个格点上，使长度为偶数。作者报告，额外加入除以 γ_q 的消光步骤后，缓解后的观测量只有边际变化。

还有两项约束限定了模拟。Trotter 化被限制在 N = 20 步，双量子比特线路深度为 5N + 2，并在碰撞附近细化时间网格以降低晚期 Trotter 误差。作者与基于 MPS 的 TEBD 张量网络模拟比较，称吻合程度足以区分不同的动力学行为；这比逐点定量吻合是更弱的标准。

平台在散射矩阵之外还提供了什么

论文对自身进展的定位很具体：能够按需制备选定类型的准粒子，并编程控制它们的运动与相互作用，这使数字量子处理器成为探测准粒子散射所产生纠缠结构的工具。测得的透射概率被表述为首次在量子硬件上对 Hubbard 模型准粒子做此类测量，并与 Bethe-ansatz S 矩阵在无自由参数条件下相符。

因此，论文中最值得读者权衡的分界，是两项强度不同的结论。随 U 变化的透射曲线、碰撞的弹性特征以及 doublon-doublon 对照，构成一次与精确两体结果相比较的硬件测量。纠缠解释则是从电荷信号消失、质心分布变宽以及 U/J = 4 附近数熵达到最大推断出的叠加态，相干性仍留待未来的干涉实验确认。读者评估这项工作时，应把散射矩阵的吻合视为已确立的结果，把纠缠产生视为有充分动机的推断，其直接确认被作者本人推迟。

本文引用的文献