2025-11-23T19:31:17.813110

A Monte Carlo approach to bound Trotter error

Blunt, Ivanov, Bay-Smidt

Trotter product formulas are a natural and powerful approach to perform quantum simulation. However, the error analysis of product formulas is challenging, and their cost is often overestimated. It is established that Trotter error can be bounded in terms of spectral norms of nested commutators of the Hamiltonian partitions [Childs et al., Phys. Rev. X 11, 011020], but evaluating these expressions is challenging, often achieved by repeated application of the triangle inequality, significantly loosening the bound. Here, we show that the spectral norm of an operator can be upper bounded by the spectral norm of an equivalent sign-problem-free operator, which can be calculated efficiently to large system sizes using projector Monte Carlo simulation. For a range of Hamiltonians and considering second-order formulas, we demonstrate that this Monte Carlo-based bound is often extremely tight, and even exact in some instances. For the uniform electron gas we reduce the cost of performing Trotterization from the literature by an order of magnitude. For the Pariser-Parr-Pople model for linear acene molecules, which has $\mathcal{O}(N^2)$ long-range interaction terms, we show that it suffices to use $\mathcal{O}(N^{0.57})$ Trotter steps and circuit depth $\mathcal{O}(N^{1.57})$ to implement Hamiltonian simulation. We hope that this approach will lead to a better understanding of the potential accuracy of Trotterization in a range of important applications.

academic

トロッター誤差を制限するためのモンテカルロアプローチ

基本情報

論文ID: 2510.11621
タイトル: A Monte Carlo approach to bound Trotter error
著者: Nick S. Blunt, Aleksei V. Ivanov, Andreas Juul Bay-Smidt
分類: quant-ph physics.chem-ph
発表日: 2025年10月14日
論文リンク: https://arxiv.org/abs/2510.11621

概要

トロッター積公式は量子シミュレーションを実行するための自然で強力な方法である。しかし、積公式の誤差分析は困難であり、その成本はしばしば過大評価される。トロッター誤差はハミルトニアン分割のネストされた交換子のスペクトルノルムで制限できることが知られているが、これらの式を評価することは困難であり、通常は三角不等式を繰り返し適用することで達成され、制限を大幅に緩和する。本論文は、演算子のスペクトルノルムが等価な符号なし問題演算子のスペクトルノルムで上界できることを示し、後者は投影モンテカルロシミュレーションを使用して大規模システムまで効率的に計算できる。一連のハミルトニアンと二次公式に対して、著者らはこのモンテカルロベースの制限が通常極めて厳密であり、場合によっては正確であることを証明した。

研究背景と動機

問題定義

中核的問題：量子シミュレーションにおけるトロッター積公式の誤差分析の困難さにより、その成本が系統的に過大評価される
重要性：量子シミュレーションは量子計算の最も有望なアプリケーションの一つであり、トロッター化は量子シミュレーションを実現するための単純で自然な方法である
既存の制限：
- 既存のトロッター誤差制限はネストされた交換子のスペクトルノルムに基づいているが、これらの式を計算することは計算上困難である
- 通常、三角不等式を使用して計算を簡略化するが、これは制限を大幅に緩和する
- 耐性量子計算のリソース推定研究では、トロッター化のパフォーマンスが大幅に過小評価されている

研究動機

著者らは、トロッター化のパフォーマンスが過小評価される理由が2つあると考えている：

トロッター回路は効率的に並列化でき、回路深度と体積を削減できる
トロッター誤差分析は非常に困難であり、積公式の誤差はしばしば数桁過大評価される

中核的貢献

理論的貢献：任意の演算子のスペクトルノルムが等価な符号なし問題演算子のスペクトルノルムで上界できることを証明： $\|A\| \leq \|\text{abs}(A)\|$
方法的革新：投影モンテカルロ法（FCIQMC）を使用して符号なし問題演算子のスペクトルノルムを効率的に計算する方法を提案
実用的応用：
- 均一電子ガスに対して、文献中のトロッター化成本を1桁削減
- パリサー・パー・ポップルモデルの線形アントラセンに対して、 $O(N^{0.57})$ 個のトロッター段階のみが必要であることを証明
ベンチマーク結果：複数の重要なハミルトニアンに対して正確なトロッター誤差ノルムベンチマークを提供

方法の詳細

中核的理論基礎

スペクトルノルム上界定理：行列 $A$ に対して、 $\text{abs}(A)$ を要素が $|A_{ij}|$ である行列として定義すると： $\|A\| \leq \|\text{abs}(A)\|$

証明の要点：

スペクトルノルムの定義を利用： $\|A\| = \max_{\|v\|_2=1} \|Av\|_2$
コーシー・シュワルツ不等式を通じて不等式関係を確立
重要な点は $\|\text{abs}(A)\|$ がペロン・フロベニウス定理を通じて効率的に計算できることである

トロッター誤差制限

ハミルトニアン $H = T + V$ に対して、二次トロッター公式： $S_2(t) = e^{-iVt/2}e^{-iTt}e^{-iVt/2}$

誤差制限は： $\|S_2(t) - U(t)\| \leq W t^3$

ここで： $W_{VTV} = \frac{1}{12}\|[[V,T],T]\| + \frac{1}{24}\|[[V,T],V]\|$

モンテカルロ計算方法

FCIQMC アルゴリズム：

状態更新： $C_i(\tau + \Delta\tau) = C_i(\tau) - \Delta\tau \sum_j (A_{ij} - S\delta_{ij})C_j(\tau)$
符号なし問題処理： $A \rightarrow -\text{abs}(A)$ を設定して符号問題を排除
固有値推定：シフト推定器と混合推定器を使用