4. 收敛性判据设计
在量子系统仿真中,“跑完”不等于“收敛”。收敛性判据的任务,是把“结果已经稳定到可以采信”这一模糊判断,转化为可计算、可复现、可自动执行的数值规则。本节围绕能量、保真度、可观测量期望值三类目标量,讨论绝对判据、相对判据、滑动窗口判据与统计置信判据的设计方法,并说明各自的适用边界。
一、先明确目标量的误差来源
不同目标量的收敛行为差异很大。能量 E 是标量,通常随变分参数或截断维度单调下降并趋于平台,适合用相邻迭代差判断;保真度 F = |⟨ψ_exact|ψ⟩|² 有界于 [0,1],接近 1 时绝对差会被压缩,必须用 1−F 或相对量刻画;可观测量期望值 ⟨O⟩ 可能正负交替、量级悬殊,且受统计采样噪声影响,单次差值不可靠。因此判据设计的第一步是区分“系统性截断误差”与“随机采样噪声”,前者靠加密参数观察平台,后者靠统计估计。
二、绝对判据与相对判据
绝对判据形式最简单:|E_n − E_{n−1}| < ε_abs。它适合能量这类有明确物理量纲、且我们关心绝对精度的场景,例如化学精度 ε_abs = 1.6×10⁻³ Ha。缺点是当目标量本身量级很大时,绝对阈值过松;量级很小时又过严。
相对判据 |E_n − E_{n−1}| / max(|E_n|, δ) < ε_rel 解决了量级依赖问题,分母加小量 δ 防止除零。对保真度更推荐用 1−F < ε,因为 F 接近 1 时 ΔF 的数值分辨率不足。实践中常采用“双判据与”策略:绝对判据保证物理精度,相对判据防止大数掩盖停滞,两者同时满足才判定收敛。
10⁻⁶~10⁻⁸;保真度取 1−F < 10⁻⁴;期望值相对判据取 10⁻³ 量级,并配合误差棒判断。三、滑动窗口判据:识别平台而非单步抖动
单步差值容易被偶然波动欺骗。滑动窗口判据维护最近 W 次迭代的目标量序列,计算窗口内的极差或线性拟合斜率:
窗口极差:R_W = max(E_{n−W+1..n}) − min(E_{n−W+1..n})
收敛条件:R_W < ε_window 且 |slope| < ε_slope
例如取 W = 8、ε_window = 10⁻⁵ Ha,要求连续 8 次能量波动不超过该值,可有效排除参数更新初期的振荡。窗口过小会误判,过大则浪费算力;一般取 W 为迭代周期的 2~3 倍。对期望值,窗口内还应同时检查标准差是否随采样数 1/√N 下降,以区分“真收敛”与“噪声被平均”。
四、统计置信判据
当仿真含有限采样(如量子线路多次测量、蒙特卡洛轨迹),目标量本身是随机变量。此时应给出置信区间而非点估计。设 M 次独立估计的样本均值 μ̂ 与标准误 SE = s/√M,则 μ̂ ± z_{1−α/2}·SE 为置信区间。收敛判据可写为:
半宽 H = z·s/√M
收敛条件:H < ε_stat 且 M ≥ M_min
取 α = 0.05 时 z ≈ 1.96。若要求 ⟨O⟩ 的 95% 置信半宽小于 0.01,而单次采样标准差 s ≈ 0.3,则需 M ≥ (1.96×0.3/0.01)² ≈ 3457 次。该判据的边界在于:它只控制统计误差,不控制系统性偏差;若截断维度不足,置信区间再窄也收敛到错误值。
五、判据的适用边界与组合策略
绝对判据适用于量纲清晰、目标精度已知的场景,但对大数或小数目标量不稳健;相对判据适合跨量级比较,却在目标量趋近零时失去意义,必须加保护项;滑动窗口判据能识别平台,但窗口长度与阈值需按迭代动力学调参,对非单调振荡序列可能误判;统计置信判据只处理随机误差,无法发现截断不足导致的系统偏差。
因此工程上推荐“分层组合”:先用滑动窗口确认目标量进入平台,再用绝对/相对判据给出精度门槛,最后用统计置信判据确认采样噪声已被压到门槛以下。一个可执行的检查顺序是:① 窗口极差达标;② 绝对与相对判据同时达标;③ 统计半宽达标且采样数足够;④ 主动加密一次截断维度或采样数,若结果变化仍小于门槛,才最终接受。第④步是区分“真收敛”与“假平台”的关键,也是判据设计中容易被忽略的一环。
综上,收敛判据不是单一阈值,而是由目标量特性决定的一组规则。能量重绝对精度与平台稳定性,保真度重 1−F 的绝对刻画,期望值重统计置信与噪声标度。把判据写进仿真循环的停止条件时,应同时记录窗口统计量、置信半宽与加密验证结果,使收敛结论可追溯、可复核。