Improved convergence radius of the Fer expansion for Hermitian generators

技术摘要：改进费尔展开（Fer Expansion）在厄米生成元下的收敛半径

问题陈述
由随时间变化的哈密顿量 $H(t)$ 控制的量子系统的时间演化，在形式上由时间排序指数（time-ordered exponential）描述。虽然戴森级数（Dyson series）提供了幂级数展开，但其有限阶截断通常无法保持酉性（unitarity）。为了解决这一问题，费尔展开将传播子表示为矩阵指数的无限乘积，$U(t) = \prod e^{F_j(t)}$，其中每个因子通过构造本身就是酉矩阵。与其它保持酉性的方法（如 Magnus 展开或 Trotter-Suzuki 公式）相比，费尔展开的一个显著优势是，一旦级数收敛，它具有双指数级的收敛速率。

然而，费尔展开的收敛性仅在有限半径内得到保证。前人的研究表明，对于一般的生成元，其时间积分的范数必须足够小（$\int \|A(s)\| ds \lesssim 0.8604$）。具体而言，对于厄米哈密顿量（即生成元为反厄米矩阵），已知的充分收敛半径为 $\rho_{\text{old}} = 2$。该值源自文献 [14]，其表示条件为 $\int_0^t \|H(s)\| ds < 2$。作者指出，尽管费尔展开具有更优越的收敛速度，但其收敛半径仍小于 Magnus 展开的半径 ($\pi$)。

方法论
本文采用了一种经典的固定点分析策略来推导更紧致的收敛充分条件。核心方法在于限制每一步迭代中有效哈密顿量的规模。

1. 递归结构： 费尔展开由一个递归定义，其中剩余哈密顿量 $H_{j+1}(t)$ 是通过线性映射 $R_{F_j}$ 作用于 $H_j(t)$ 生成的。该映射涉及酉共轭以及对对易子的积分。
2. 有限角度估计： 作者重点研究了限制对易子 $[e^{-\mu X}, Y]$ 的范数，其中 $X$ 是反厄米的。标准的初等估计通过将单位矩阵替换为自身（$c=1$），导致一个与弧长 $\mu\|X\|$ 成正比的界限。
3. 几何优化： 本文方法论上的关键创新在于观察到，在对易子界限中，单位矩阵可以被任何标量倍数 $c\mathbb{1}$ 替换。由于 $c\mathbb{1}$ 与所有算符都对易，因此对易子保持不变。作者通过优化 $c$ 的选择，使酉算符 $e^{-\mu X}$ 的谱（位于单位圆的一段弧上）与点 $c$ 之间的距离最小化。
   • 这被转化为一个平面几何问题：寻找包含弧 $|\theta| \le \mu\|X\|$ 的最小圆盘，其中心为 $c$。
   • 当 $\mu\|X\| \le \pi/2$ 时，最优中心为 $c = \cos(\mu\|X\|)$，从而得到一个与 $\sin(\mu\|X\|)$ 成正比的界限。
   • 当 $\mu\|X\| > \pi/2$ 时，该弧包含对跖点，最优半径为 1（以原点为中心）。
4. 标量递归： 这种几何优化产生了一个新的标量函数 $\psi(\|X\|)$，用于限制哈密顿量范数的增长。通过对该函数进行积分，作者定义了一个控制累积范数 $K_j(t) = \int_0^t \|H_j(s)\| ds$ 演化的映射 $\Psi(x)$。收敛半径由第一个非平凡正不动点 $\rho$ 决定，满足 $\Psi(\rho) = \rho$。

主要贡献与结果

• 改进的收敛半径： 作者推导出了厄米生成元下费尔展开的一个新的充分收敛半径。通过数值求解固定点方程，他们发现：
  $$\rho_{\text{new}} \approx 2.6058$$
  这比之前已知的 $\rho_{\text{old}} = 2$ 提升了约 30%。
• 保持收敛速率： 分析表明，在原点附近（小范数情况），新的映射 $\Psi(x)$ 在主阶项上与旧映射 $M(x)$ 一致。具体而言，$\Psi(x) = M(x) - x^4/48 + O(x^6)$。因此，渐近的双指数衰减率保持不变。这种改进仅扩展了保证这种快速收敛所需的初始数据（哈密顿量）的定义域。
• 显式误差界限： 本文提供了一个截断乘积的显式误差界限。如果初始累积范数 $K_1(t) < \rho_{\text{new}}$，则经过 $N$ 个因子后的误差由应用于 $K_1(t)$ 的 $\Psi$ 的 $N$ 次迭代给出，并呈双指数级衰减。

意义
本文声称，这种改进源于关于酉共轭谱的一个单一几何观察。通过优化对易子界限中的参考点，作者将线性的弧长估计替换为更紧致的几何界限（最小包围圆的半径）。

其意义在于扩大了费尔展开的适用范围。虽然其收敛速率已知优于其他方法，但严格的半径条件（$\rho=2$）限制了其在处理较强哈密顿量或较长时间间隔时的应用。新的界限（$\rho \approx 2.6058$）使得费尔展开能够更严谨地应用于更广泛的随时间变化的哈密顿量，同时保持其酉结构和快速收敛特性。作者指出，这一结果在不改变展开的基本结构或其渐近行为的前提下，改进了迄今为止已知的最佳数值。