Improved convergence radius of the Fer expansion for Hermitian generators

技術要約：エルミート生成子に対するFer展開の収束半径の改善

問題提起
時間依存のハミルトニアン $H(t)$ によって支配される量子系の時間発展は、形式的には時間順序付き指数関数によって記述される。ダイソン級数は冪級数展開を提供するが、その有限次までの打ち切りは一般にユニタリ性を保持できない。これに対処するため、Fer展開はプロパゲーターを、各因子が構成上ユニタリとなる行列指数の無限積 $U(t) = \prod e^{F_j(t)}$ として表現する。Magnus展開やTrotter-Suzuki公式などの他のユニタリ保存手法と比較した際のFer展開の大きな利点は、級数が収束一旦、二重指数関数的な収束率を持つことである。

しかし、Fer展開の収束は有限の半径内でのみ保証される。先行研究では、一般的な生成子に対して、そのノルムの時間積分が小さいこと（$\int \|A(s)\| ds \lesssim 0.8604$）が必要であると確立されていた。具体的には、エルミート・ハミルトニアン（生成子が反エルミートである場合）において、既知の収縮半径は $\rho_{\text{old}} = 2$ であった。これは、文献[14]で導出された条件 $\int_0^t \|H(s)\| ds < 2$ を表している。著者らは、この半径がMagnus展開の半径（$\pi$）よりも小さいことを指摘している。これは、Fer展開がより優れた収束速度を持つにもかかわらずである。

手法
本論文では、よりタイトな収束十分条件を導出するために、古典的な不動点解析戦略を採用している。手法の核心は、各ステップにおける有効ハミルトニアンの大きさを抑え込むことにある。

1. 再帰的構造： Fer展開は、前のステップの $H_j(t)$ に対して線形写像 $R_{F_j}$ が作用することで、残差ハミルトニアン $H_{j+1}(t)$ が生成される再帰によって定義される。この写像には、ユニタリ共役および交換子の積分が含まれる。
2. 有限角推定： 著者らは、反エルミートな $X$ に対して、交換子 $[e^{-\mu X}, Y]$ のノルムを抑え込むことに焦点を当てる。標準的な初等的な推定では、単位行列をそれ自身（$c=1$）で置き換えるため、弧長 $\mu\|X\|$ に比例する境界が得られる。
3. 幾何学的最適化： 手法の主要な革新は、交換子の境界において、単位行列を任意のスカラ倍 $c\mathbb{1}$ で置き換えられるという観察である。$c\mathbb{1}$ はすべての演算子と可換であるため、交換子は変化しない。著者らは、ユニタリ演算子 $e^{-\mu X}$ のスペクトル（単位円上の弧）と点 $c$ との距離を最小化するように $c$ を選択し、最適化を行う。
   • これは平面幾何学の問題として定式化される。すなわち、弧 $|\theta| \le \mu\|X\|$ を含む最小の円の半径を求める問題である。
   • $\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)$ を定義する。収束半径は、$\Psi(\rho) = \rho$ を満たす最初の非自明な正の不動点 $\rho$ によって決定される。

主な貢献と結果

• 改善された収束半径： 著者らは、エルミート生成子に対するFer展開の新しい収束十分半径を導出した。不動点方程式を数値的に解くことにより、以下の値を得た：
  $$\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$ 回の反復によって抑えられ、二重指数関数的に減衰する。

意義
この改善は、ユニタリ共役のスペクトルに関する単一の幾何学的観察に起因すると著者らは主張している。交換子の境界における参照点を最適化することで、線形な弧長推定を、よりタイトな幾何学的境界（最小包含円の半径）に置き換えている。

その意義は、Fer展開の適用範囲を拡大することにある。その収束率はすでに他の手法よりも優れていることが知られているが、厳格な半径条件（$\rho=2$）が、より強いハミルトニアンや長い時間間隔での使用を制限していた。新しい境界（$\rho \approx 2.6058$）により、ユニタリ構造と高速な収束特性を維持しつつ、より広いクラスの時間依存ハミルトニアンに対してFer展開を厳密に適用することが可能になる。著者らは、この結果が、展開の基本構造や漸近的挙動を変えることなく、現在知られている最高値を改善したものであると述べている。