Full Uncertainty Quantification of Sign-Problem-Free Quantum Monte Carlo Methods and Nuclear Lattice Effective Field Theory Benchmarks

技術要約：符号問題のない量子モンテカルロ法の完全な不確かさ定量化と核格子有効場理論のベンチマーク

問題の所在
現代のab initio（第一原理）核理論は、カイラル有効場理論（EFT）と多体ソルバーを組み合わせることで、軽い原子核や無限核物質の観測量を再現することに大きな成功を収めてきた。しかし、これらの手法を中質量および重い原子核（$^{100}$Snまで）へと拡張するにあたっては、決定的なボトルネックに直面している。それは、系統的なアルゴリズム誤差に関する厳密な不確かさ定量化（UQ）の欠如である。このような制御がなければ、理論と実験の間の相違は、ソルバーの不備を補うために低エネルギー定数（LEC）を「再調整」することに帰せられてしまう。これは、理論的誤差の真の起源を隠蔽し、ab initio的手法の予測能力を損なうものである。

具体的には、近年の批判（文献[1]）は、変分境界、相関エネルギー、および核飽和メカニズムの問題を理由に、核格子有効場理論（NLEFT）のベンチマークの妥当性に異議を唱えている。これらの批判は、多体手法の真の失敗ではなく、異なるレギュレートされた格子理論（例：転送行列形式とハミルトニアン形式の混合）や、不整合な有限体積処理を比較していることに起因することが多い。

手法
著者らは、NLEFTの枠組み内で、特にLAT-OPT1作用を用いた、符号問題のない量子モンテカルロ（QMC）法、具体的には補助場量子モンテカルロ（AFQMC）を用いている。この研究は、偶ー偶原子核において正のモンテカルロ重みを維持するために必要な反ユニタリ対称性を保持しつつ、非摂動的なスピン軌道相互作用を組み込んだ、最小限の符号問題のない格子ハミルトニアンであるLAT-OPT1作用に焦点を当てている。

手法には、包括的な誤差予算分析と、ベンチマーク比較の診断的再検討が含まれる：

1. 不確かさの定量化: 著者らは、$^4$Heから$^{100}$Snまでの原子核について、統計的、虚時間、有限体積、時間離散化、補助場誘起、および残留符号の不確かさを体系的に定量化している。
2. 非摂動的スピン軌道処理: 以前のアプローチではスピン軌道力が摂動的な挿入として扱われていたが、LAT-OPT1は、スピン軌道密度をユークリッド時間投影ハミルトニアンの中に直接埋め込んでいる。
3. レギュレータの一貫性チェック: 論文では、以下の異なる計算形式を厳密に比較している：
   • 転送行列 vs. ハミルトニアン: 有限の時間格子間隔（$a_t$）を持つ転送行列計算と、$a_t \to 0$のハミルトニアン極限との差異の分析。
   • 境界条件: 核物質における有限体積のシェル効果に対する、周期境界条件（PBC）と平均ツイスト境界条件（ATBC）の影響の調査。
   • 繰り込み: 連続体極限と短距離相互作用の繰り込みの検討。

主要な貢献と結果

• 系統的誤差の定量的制御:
  本研究は、二重魔法原子核（$^{100}$Snまで）の基底状態エネルギーにおける累積的な多体計算不確かさが、パーセントレベルを十分に下回っていることを示している。

  • 虚時間および有限体積: 無限投影時間（$\tau \to \infty$）および無限体積（$L \to \infty$）への外挿は安定している。検討された原子核に対して、ボックスサイズ $L=11$（$a=1.32$ fm）が十分に収束していることが示されている。
  • 時間離散化: 有限の $a_t$ を持つ転送行列形式の使用は、標的となるハミルトニアン $H$ とは異なる有効ハミルトニアン $H_{eff}$ を誘起する。著者らは、$H_{eff}$ が $a_t$ に依存するものの、エネルギー誤差は $O(a_t^2)$ でスケールするため、精密な制御が可能であることを示している。
  • 補助場誘起項: スピン軌道項に対して符号問題のないサンプリングを維持するために必要な平方完成は、誘起された二体（$\rho_s^2$）および三体項を生成する。これらは定量化されており、数MeV程度の小さな補正であることが確認され、構成の妥当性を裏付けている。
  • 残留符号問題: 奇数質量数（odd-$A$）および奇ー奇原子核では、軽微な残留符号問題が存在する。しかし、平均位相は1に近い値（0.94–0.98）を維持しており、統計的不確かさは偶ー偶系と同程度である。

• スピン軌道結合の非摂動的役割:
  スピン軌道相互作用は、ユークリッド時間投影において非摂動的に含まれる必要があるという点が重要な知見である。これをSU(4)対称な波動関数に対する摂動的な挿入としてのみ扱うと、中質量原子核における数から十MeVのシェル構造への寄与を見逃してしまう。投影過程に組み込むことで、正しいシェル構造が動的に生成される。

• ベンチマークの不一致の解決:
  本論文は、異なる格子理論間の概念的および技術的な誤りを特定することで、最近の批判（文献[1]）に対処している：

  • 不等価な比較: $a_t \to 0$ のハミルトニアンからのハートリー＝フォック（HF）境界と、有限の $a_t$ を持つ転送行列からのAFQMCの結果を比較することは無効である。有限の $a_t$ 転送行列は、誘起相互作用を持つ、異なるレギュレートされた理論を定義する。同一のレギュレートされた理論を一貫して使用する場合、変分原理は成立する（HFエネルギーは投影エネルギーを下に抑える）。
  • 相関エネルギー: 「Essential Elements (EE)」相互作用における中性子物質の小さな相関エネルギーは、自明な波動関数を意味するものではない。既約密度の詳細な分析により、エネルギー差は小さいものの、顕著なマルチモーダル超流動性（s波、p波、およびクアテット相関）が存在することが明らかになった。
  • 核飽和: 高密度における相互作用エネルギーの平坦化は、格子充填のアーティファクトではない。ガリレイ不変なレギュレータを持つ連続空間における解析的導出により、相互作用エネルギーの平坦化は、運動量空間のレギュレーションの自然な帰結であることが示されている。NLEFT相互作用の飽和特性は、格子サイト占有率ではなく、スミアリングによるレギュレータ・スケールに従う。
  • 連続体極限: 格子相互作用が連続体極限において系の崩壊を招くという主張は、否定されている。批判で使用されたスケーリング引数は、誤ったボルンレベルの繰り込みに基づいている。物理的な散乱観測量への適切な繰り込みは、崩壊を防ぐ異なるスケーリング挙動（$C_c \propto a$）をもたらす。

意義と主張
本論文は、完全なカイラル力を用いたNLEFTシミュレーションにおける完全な誤差分析のための強固な基礎を確立したと主張している。符号問題のない手法が重い原子核に対してもサブパーセントの精度を達成できることを厳密に定量化することで、本研究はNLEFTが信頼できる多体手法のベンチマークであることを検証している。

著者らは、一貫したパラメータ、体積、および境界条件を持つ等価なレギュレートされた理論の間で比較を行うことで、近年の文献における明白な不一致は解決されると断言している。LAT-OPT1作用は、単なる現象論的なツールではなく、本質的なシェル構造のダイナミクスを捉える、厳密に解ける非摂動的なベンチマークモデルとして提示されている。この枠組みにより、真の核力の欠陥と多体ソルバーの誤差を分離することが可能となり、中質量および重い原子核におけるシェルの進化、分光、および高密度物質の予測的研究への道が開かれる。