低能量子散射的量子线路实时模拟

《The Journal of Physical Chemistry Letters》：Real-Time Simulation of Low-Energy Quantum Scattering with Quantum Circuits

　　

模拟原子与分子散射等化学动力学是量子计算的自然目标，因为常规（经典）硬件面临指数级代价。研究人员提出了一种乘积公式（product-formula）量子算法，通过对高斯波包（Gaussian wavepacket）进行实时演化来建模单通道氦–氦散射，且无需预先计算本征态。与伦纳德-琼斯（Lennard-Jones）势及离心势垒（centrifugal barrier）的相互作用以酉方式施加，采用闭式（closed-form）线路，并通过无误差量子线路编译直接处理演化。对波函数的追踪得到收敛的相移（phase shift），并由此计算偏转函数（deflection function）、散射振幅（scattering amplitude）和微分截面（differential cross section）。全同玻色子（identical-boson）积分截面与时间无关的Numerov方法和对数导数（log-derivative）基准一致到6.1%，残差局域在s波散射中。该构造可实现任意形如∑ir–mi的势，其门数（gate count）与单个1/r项相同，可高效扩展到多通道动力学，并为近似量子模拟提供基真（ground-truth）算符。

《The Journal of Physical Chemistry Letters》刊发的这项工作面向化学散射中经典计算受维数灾难限制的问题：多自由度量子散射的状态空间随体系规模指数增长，子空间方法和半经典近似虽能处理特定相互作用，但本征态预计算本身仍具指数复杂度；近期算法先在传统硬件求本征态再实时演化，虽降低测量开销，却未摆脱本征态预计算瓶颈。为此，研究人员开展无本征态预计算的实时波包演化研究，用二进制编码位置寄存器上的闭式线路实现逆幂律势，验证其与时间无关基准的一致性，并为多通道量子散射模拟提供可扩展框架。

关键技术方法上，研究人员采用二阶Trotter分解传播径向含时薛定谔方程；在二进制位置基底用默比乌斯反演将1/r展开为多线性形式，以多控制R_z相旋实现逆多项式势；动能项在居中量子傅里叶变换（centered quantum Fourier transform，cQFT）间作为动量多项式施加；通过空间快速傅里叶变换（fast Fourier transform，FFT）对有/无离心势的波函数快照作参考相减提取相移；以Qiskit进行无误差态矢量仿真，并用Numerov与log-derivative作传统硬件基准。

研究结果如下。

实时代数薛定谔演化：研究人员在二进制位置基x = Δ∑_j=0^n–12^jq_j、N = 2ⁿ格点、盒长L、分辨率Δ = L/N下，将动能写成居中量子傅里叶变换间的二阶级多项式，初始态取高斯波包ψ(r,t₀) = (1/2πδ²)^1/4e^{–(r–r₀)2/4δ2}，结论是无需本征态即可推进弹性碰撞演化。

逆幂律势的闭式构造：研究人员将1/x表示为子集S的乘积求和，系数由默比乌斯反演a_S = ∑_T?S(–1)^|S|–|T|/[Δ∑_j∈T2^j]给出；1/r映射为|S|重控制R_z(2παa_S)门。结论是1/r^m与1/r使用相同2ⁿ–1个多控制旋转，Lennard-Jones 12–6加离心项V_eff(r) = 4?[(σ/r)¹²–(σ/r)⁶] + l(l+1)/(2μr²)不增加额外深度，门数O(2ⁿ)，多通道仅加?log₂M?辅助寄存器。

实时Trotter演化：研究人员用U ≈ [e^–iTτ/2e^–iVτe^–iTτ/2]^k二阶乘积公式迭代；结论是只需动能、闭式势和初态即可完成散射传播，Walsh–Hadamard变换还可把多控制电路化成CNOT加无控制R_z。

氦–氦散射与相移提取：研究人员对l = 0–18共19个分波跑全V_eff与仅V_cf两组，用FFT得出射/入射振幅比S_LJ(l,k) = e^2iδ_LJ(l,E)；由相移算散射振幅f(θ,E) = (1/2ik)∑_l(2l+1)(e^2iδ_l–1)P_l(cosθ)、微分截面dσ/dΩ = |f|²、总截面σ_tot = (4π/k²)∑_l(2l+1)sin²δ_l。结论是相移、偏转函数、散射振幅、微分截面均与Numerov和log-derivative吻合。

全同粒子对称化与积分截面：研究人员对⁴He取玻色对称σ_bos = (8π/k²)∑_{l even}(2l+1)sin²δ_l；电路仿真85.14 a₀²、Numerov 80.28 a₀²、log-derivative 80.39 a₀²，可区分粒子为307.76、321.91、321.94 a₀²；时间无关法互差0.11%，电路对称化结果较Numerov高6.05%，残差主要在s波，剔除s波后一致到0.28%。

寄存器规模与动量分辨率：研究人员由n_min = ?log₂(L√(2μE_char)/π)?、Δ_p/k₀ = 2π/(L√(2μE))、ε_δ ? |dδ_l/dk|(π/L)、n_min = ?log₂(2/ε_k√(V_cut/E))?给出标度；He–He中V_cut = 2.27 E_h、ε_k = 0.145得n_min = 13.0。结论是s波相移对Δ_p最敏感，加大盒长L可使误差按1/L下降。

讨论部分指出，量子信号处理（quantum signal processing，QSP）虽有最优查询复杂度，但谱范数随强耦合多通道增大而不利；Trotter化在势能面交叉处不变且交换子范数趋零，非稀疏势能量面需O(2ⁿ)项量子化，故乘积公式仍有价值。闭式算符可作硬件无关基真，用以通过希尔伯特–施密特距离（Hilbert–Schmidt distance）和纠缠保真度基准化变分压缩、量子辅助编译与强化学习线路合成；二进制编码经二进制反射格雷码（binary-reflected Gray code，BRGC）证明可高效逼近。

研究结论部分译为：综上，研究人员展示了无本征态预计算的原子散射事件实时量子线路仿真；经空间FFT参考相减提取的相移复现了氦–氦散射的偏转函数、散射振幅、微分与积分截面，并与传统硬件上的Numerov和log-derivative独立计算一致；到BRGC构造的希尔伯特–施密特距离在数值精度内归零，给出无近似的算符级验证。当前实用部署受远低于现有门保真度的旋转角精度，以及相移提取所需的长演化时间限制，但闭式形式界定了这些代价并为高效近似指明路径。借助多通道扩展，该框架迈向传统计算方法不可及的耦合势能面量子模拟。