Error bounds for the truncated Baker--Campbell--Hausdorff and Zassenhaus formulas in unitary problems

技术摘要：酉算符问题中截断 Baker–Campbell–Hausdorff 与 Zassenhaus 公式之误差界

问题陈述
Baker–Campbell–Hausdorff (BCH) 公式与 Zassenhaus 公式是李理论中的基本工具，分别将非对易算符的指数乘积的对数，以及指数和的指数，表示为嵌套交换子的无穷级数。在实际应用中，特别是在微分方程的数值分析和量子模拟领域，这些无穷级数必须进行截断。虽然这些级数的收敛域已有深入研究，但针对截断近似（特别是针对嵌套交换子范数）的严格误差界在文献中长期以来是缺失的，仅在平凡的 Lie–Trotter 情况中有所涉及。这种误差界的缺失阻碍了对量子算法中门复杂度（gate complexity）和误差缩放（error scaling）的精确量化。

方法论
作者提出了一种通用策略，用于推导当算符为斜伴随（skew-adjoint）时（这一条件在酉量子演化问题中得到满足）的 BCH 和 Zassenhaus 公式的严格误差界。该方法依赖于以下步骤：

1. 微分方程表述： 通过定义两个随时间变化的酉算符 $F_1(t)$ 和 $F_2(t)$ 来重构问题，其中 $F_1(t)$ 代表精确的指数乘积（或和），而 $F_2(t)$ 代表截断近似。两者均满足形式为 $\frac{dF}{dt} = M(t)F$ 的初值问题。
2. 生成元分析： 作者分析了生成元算符 $M_1(t)$ 与 $M_2(t)$ 之间的差异。对于 BCH 情况，$M_1(t)$ 由乘积 $e^{tA}e^{tB}$ 导出，而 $M_2(t)$ 则由截断级数 $e^{Z^{[m]}(tA, tB)}$ 的指数导出。
3. 积分余项估计： 利用引理 2.1（该引理通过两个酉算符之差的生成元之差的积分来界定两个酉算符的差异），作者将误差估计简化为对 $\|M_1(t) - M_2(t)\|$ 的界定。
4. 带有余项的泰勒展开： 生成元按参数 $t$ 的幂次展开至截断阶数 $m$，并利用积分形式的余项（引理 2.2）。这使得作者能够分离出领先误差项，这些项被表示为嵌套交换子的线性组合。
5. 递归扩展： 该方法被推广到处理 $N$ 个算符（第 4 节）和对称 BCH 公式（第 5 节），通过定义生成元项的递归序列来实现。对于 Zassenhaus 公式（第 6 节），类似的递归程序也被应用于指数乘积。

主要贡献与结果

• BCH 的显式误差界： 本文推导了截断 BCH 级数 $Z^{[m]}$ 的显式误差界。对于在两个项后截断（$m=2$）的情况，作者给出了界限：
  $$\|e^{tA}e^{tB} - e^{Z^{[2]}(tA, tB)}\| \leq \frac{t^3}{4}\|[A, [A, B]]\| + \frac{t^3}{12}\|[B, [A, B]]\|$$
  该结果恢复了 $m=1$ 时已知的 Lie–Trotter 界，并为二阶近似提供了一个新的、显式的估计。
• Zassenhaus 的显式误差界： 作者提出了据其所知首个关于截断 Zassenhaus 公式的严格误差界。对于在第一个修正项后截断（$q=2$）的情况，界限为：
  $$\|e^{t(A+B)} - e^{tA}e^{tB}e^{t^2 C_2}\| \leq \frac{t^3}{3}\|[B, [A, B]]\| + \frac{t^3}{6}\|[A, [A, B]]\|$$
  作者指出这些系数是最优的。
• 向 $N$ 个算符及对称情况的推广： 该策略被扩展到涉及 $N$ 个斜伴随算符的 BCH 公式以及对称 BCH 公式（与回文分裂法相关）。定理 4.1 和定理 5.1 为这些情况提供了通用的误差界递归公式。
• 交换子缩放性质： 一个核心发现是，推导出的误差界是关于参数 $t$ 的多项式，其系数是嵌套交换子范数的线性组合。这体现了“交换子缩放性质”，当算符近乎对易时，该性质非常有利。
• 计算实现： 作者利用自由李代数的经典 Hall 基，提供了直到 $m=3$（以及通过提供的网站获取 $m=4$）的显式表达式。他们通过使用随机 $20 \times 20$ 矩阵的数值实验证明，推导出的界是紧致的，即使在 $t$ 相对较大时，对实际误差的过高估计也仅在 3 倍以内。

意义与主张
本文声称解决了文献中关于超越平凡一阶情况的截断 BCH 和 Zassenhaus 展开的严格误差分析的空白。这些结果的意义在于其在以下领域的应用：

• 量子模拟： 提供直接影响量子算法中乘积公式的门复杂度与深度优化的理论界。
• 误差量化： 提供依赖于嵌套交换子的界，这使得在算符近乎对易时能够获得更紧凑的估计。
• 严格分析： 建立了一个系统性的框架，用于生成任意截断阶数 $m$ 和 $N$ 个算符的误差界，并可通过计算机代数系统实现。

作者保持了谦逊的语气，指出虽然这些界是严格且显式的，但其计算复杂度随截断阶数 $m$ 指数级增长。他们强调，这些结果是针对有限维酉情况（斜伴随算符）特有的，但暗示该处理方法可以推广到这些公式的其他变体。