3. 数值离散化与截断误差来源

在量子系统仿真中，我们几乎从不直接求解连续的薛定谔方程或连续的算符本征问题，而是先把连续对象“翻译”成有限维、有限精度的离散对象，再交给计算机处理。这个翻译过程本身就会引入误差。与浮点舍入误差不同，离散化误差不会因为换用更高精度的浮点类型而消失，它来自我们对连续世界的有限近似，因此必须被识别、量化并主动控制。本节聚焦四类最常见的离散化截断来源：时间步长、空间网格、基组截断与能级截断，并给出可操作的误差来源清单与初步控制策略。

四类截断误差的物理来源

时间步长截断出现在含时演化中。以最常用的四阶龙格–库塔（RK4）为例，单步局部误差为 O(Δt⁵)，全局误差为 O(Δt⁴)；而一阶欧拉法全局误差仅为 O(Δt)。当哈密顿量含快速振荡项（如驱动频率 ω_d 接近系统能隙）时，步长必须同时分辨振荡周期，否则会出现相位漂移与布居数虚假转移。判断依据不是“步长足够小”，而是“步长能分辨系统最快的物理时间尺度”。

空间网格截断出现在位置表象求解中。有限差分把二阶导数近似为相邻格点的线性组合，其截断误差为 O(Δx²)；同时，计算域被限制在 [−L/2, L/2] 内，波函数在边界处被人为截断或吸收。网格过粗会低估动能、扭曲色散关系；域过小则会让本应向外传播的波包在边界反射回来，污染后续演化。对束缚态问题，域必须显著大于波函数的衰减长度；对散射问题，则需配合吸收边界或复吸收势。

基组截断出现在以有限个基函数展开波函数或算符时。无论用平面波、谐振子本征态还是高斯型基组，截断都意味着丢掉高能分量。其误差量级取决于被丢弃分量在目标态中的权重：若目标态在截断基上的残余范数平方为 ε²，则能量误差通常按 ε² 量级下降。基组截断的隐蔽之处在于，它可能让低能本征值看起来收敛，却让高激发态或含时响应严重失真。

能级截断出现在多能级系统或腔 QED 类模型中。把无穷维希尔伯特空间截到 N 个能级，等价于假设高能级布居可忽略。当驱动强度或耦合强度 g 增大时，高能级布居随之上升，截断误差迅速放大。经验上，截断能级数应使最高保留能级的布居长期保持在 10⁻⁴ 以下，并随耦合强度增大而增加。

建立误差来源清单

把上述四类来源落成一份可执行的清单，关键是让每一项都绑定一个“可扫描的离散化参数”和一个“可观测的误差指标”。下面给出一个可直接套用的示例结构（参数为教学示例，非实测结果）：

初步控制策略：收敛扫描而非盲目加密

控制离散化误差最有效的手段是收敛扫描：固定其他参数，只把一个离散化参数按 2 倍或 1.5 倍逐级加密，观察目标量是否按理论阶数趋于稳定。例如对 RK4 时间步长，若把 Δt 减半后误差应缩小约 16 倍（O(Δt⁴)），实测缩小倍数明显偏离，说明误差已被其他来源主导，此时继续减小 Δt 是浪费算力，应转而检查空间网格或基组。这种“误差主导项识别”比单纯堆参数更重要。

对空间网格，推荐同时做两件事：一是把 Δx 减半观察能量变化，二是把计算域 L 扩大 1.5 倍观察边界反射是否消失。两者必须分开扫描，否则无法判断误差来自差分近似还是域截断。对基组与能级截断，则应直接监控残余范数与最高能级布居这两个“泄漏指标”，它们比能量本征值更敏感，能在能量尚未明显偏移时就预警截断不足。

综上，离散化误差不是仿真中需要“尽量消除”的瑕疵，而是必须被显式建模、逐项量化的系统属性。把时间步长、空间网格、基组截断、能级截断四类来源分别绑定到可扫描参数与可观测指标，再通过收敛扫描识别误差主导项，就能在精度与算力之间做出有依据的取舍，为后续的收敛性验证与精度评估打下基础。