Comment on "Beyond-classical computation in quantum simulation"

技術要約：「量子シミュレーションにおける古典を超えた計算」

問題提起
本コメントは、最近のベンチマーク研究 [1]（King et al., Science 388, 199-204, 2025）を取り上げている。当該研究では、スピングラス・ハミルトニアンにおける量子アニーリング動力学をシミュレートするために、様々な古典的数値手法と超伝導量子プロセッサユニット（QPU）を比較した。元の研究では、ニューラルネットワークに基づく変分アンザッツであるニューラル量子状態（NQS）が、特定のインスタンス（具体的には、$L=8$ の二次元正方格子、アニーリング時間 $t_a = 7$ ns の場合）において、QPUと同等の精度に達することに失敗したと結論付けている。本コメントの著者らは、報告されたNQSの失敗は、手法自体の根本的な限界ではなく、統計的サンプリングの制限によるアーティファクト（偽像）である可能性があると主張している。

手法
著者らは、元の研究で使用された特定の無秩序実現（disorder realization）とパラメータを再検討し、NQSの時刻依存変分原理（TDVP）進化に焦点を当てている。彼らの分析は以下の手順で行われる：

1. 相関の再評価： 元の論文の集約された相関誤差指標（$\epsilon_c$）のみに頼るのではなく、距離分解された二点相関の大きさ $C(d) = \langle |c_{i,j}| \rangle$ （ここで $c_{i,j} = \langle \sigma^z_i \sigma^z_j \rangle$）を分析する。
2. 統計ノイズ分析： 相関の推定におけるモンテカルロ（MC）サンプリングノイズの影響を調査する。彼らは、真の値が小さい長距離相関において、MCノイズが $\sigma_{MC} \approx 1/\sqrt{N_{MC}}$ とスケールすることを導出した。元のサンプリング予算 $N_{MC} = 10^6$ では、ノイズフロアは約 $10^{-3}$ となる。
3. 誤差分解： 著者らは、サンプリング予算 $N_{MC}$ の関数としての相関誤差 $\epsilon_c(N_{MC})$ に関する解析モデルを提案する：
   $$\epsilon_c(N_{MC}) = \sqrt{\epsilon_{sys} + \frac{\epsilon_{stat}}{N_{MC}}}$$
   ここで、$\epsilon_{sys}$ は系統的かつ手法固有の誤差を表し、$\epsilon_{stat}$ は統計的誤差の寄与を表す。
4. フィッティングと外挿： データを上記のモデルに適合させることで、サンプリングの無限極限（$N_{MC} \to \infty$）へと外挿し、NQS手法に固有の系統誤差を分離する。

主な結果

• プラトー挙動： 著者らは、距離分解された相関において、長距離での $C(d)$ が $10^{-3}$ 付近でプラトー（停滞）することに注目している。彼らは、このプラトーが $N_{MC} = 10^6$ における期待されるMCノイズの大きさと正確に一致していることを示しており、これは「長距離相関を捉えられない」という現象が、変分アンザッツの失敗ではなく、信号が統計ノイズに埋もれていることによるものであることを示唆している。
• 収束分析： 距離クラスごとにグループ化して相関を分析すると、短距離の相関は参照値である行列積状態（MPS）の基底状態へ急速に収束する。対照的に、長距離の相関は統計ノイズとして予想される $1/\sqrt{N_{MC}}$ のスケーリングに従っており、元のサンプリング予算では収束していないことが確認された。
• 修正された精度： 誤差モデルを用いて外挿を行うと、NQS手法の系統誤差（$\epsilon_{sys}$）は、元の研究で報告されたQPUの誤差よりも低いことが判明した。さらに、外挿されたNQSの誤差は、有限のサンプリングに基づいた元の論文のNQS誤差よりも低い。
• 自己相関の考慮事項： 著者らは、より長いアニーリング時間（例：$t_a = 7$ ns）の場合、最終状態が非常に鋭いボルン分布（Born distribution）を発達させ、マルコフ連鎖MCサンプリングにおいて長い自己相関時間を引き起こすことに言及している。これはサンプリングを困難にするが、相関の大きさ自体は増大するため、固定された $N_{MC}$ における相対誤差は減少する。

意義と主張
本論文は、NQSが元のベンチマークにおいて見せた見かけ上の性能不足が、サンプリング予算およびそれに伴う統計ノイズ（特に消失しつつある長距離相関に対して）によって大きく影響されていたことを示すことで、King et al. の知見に文脈を与えるものであると主張している。

• 競争力のある性能： 著者らは、NQSが量子アニーリング・プロトコルに対して競争力のある結果を提供できると結論付けている。分析された特定のインスタンスにおいて、サンプリングノイズを考慮したNQSの手法は、QPUを上回る相関誤差を実現している。
• 限界と範囲： 著者らは、自身の分析が最も単純な格子幾何学（$L=8$）および $t_a = 7$ ns までのアニーリング時間に限定されていることを認め、控えめな範囲設定を維持している。彼らは、元の研究がより複雑な幾何学構造やより長い時間（$t_a = 20$ ns）、および異なる結合分布（バイモーダル vs 一様）を扱っていたことを指摘している。
• 今後の方向性： 普遍的な優位性を主張するのではなく、より広い範囲におけるNQSの強みと限界を理解するためのさらなる数値実験を促している。その際、ダイヤモンド構造のより大きな格子への成功例 [4] を前例として挙げている。

要約すると、本論文は、NQSは検討された問題に対して本質的に劣っているのではなく、その真の系統的な精度を明らかにするためには、モンテカルロノイズと自己相関時間を慎重に考慮する必要があると主張している。