二阶 Strang 与高阶 Suzuki

• 命题框：S₂(t) = Πₖ₌₁ᴸ e^{−iαₖHₖδ/2} · Πₖ₌ₗ¹ e^{−iαₖHₖδ/2}（第二半为逆序；相邻两个 e^{−iα_LH_Lδ/2} 合并为整步 e^{−iα_LH_Lδ}）
• 推导/演算步骤：

1. 回文性质：S₂(−δ) = S₂(δ)⁻¹ ⇒ log S₂ 是 δ 的奇函数 → 误差指数只含 δ³、δ⁵、…

2. Σ₂ 相消：正序 log F(u) = −iHu − (u²/2)Σ₂ + O(u³)，逆序 log F̃(u) = −iHu + (u²/2)Σ₂ + O(u³)（每个新配对的交换子项反号）→ S₂ = F(δ/2)·F̃(δ/2) 合并时 ±(δ²/8)Σ₂ 精确相消

3. 误差：‖e^{−iHt} − S₂(t)‖ = O(Λ₃t³/r²)；由回文性，下一项已是 t⁵/r⁴

4. Suzuki 递推：S₂ₖ(t) = [S₂ₖ₋₂(pₖt)]²·S₂ₖ₋₂((1−4pₖ)t)·[S₂ₖ₋₂(pₖt)]²，pₖ = 1/(4−4^{1/(2k−1)})

5. 演算（k = 2）：p₂ = 1/(4−4^{1/3}) = 1/2.413 ≈ 0.414；五块时间参数 p₂+p₂+(1−4p₂)+p₂+p₂ = 4×0.414+(−0.657) ≈ 1.000 ✓；S₄ 每步含 5 个 S₂ 块 = 2L×5 个指数因子

• 结果框：回文结构不是装饰——偶次项相消正是二阶精度的来源；代价 2L·5ᵖ⁻¹ 随阶数指数增长，实际很少超过四阶