Full Configuration Interaction Quantum Monte Carlo for Accurate $\textit{Ab Initio}$ Nuclear Structure Calculations

技術要約：正確なアブイニシオ核構造計算のためのフル構成相互作用量子モンテカルロ法

問題提起
強く相互作用する原子核に対するアブイニシオ多体手法の開発は、現代の核理論における中心的な課題である。カイラル有効場理論（EFT）は、量子色力学から核相互作用を導出するための系統的な枠組みを提供するが、中質量および重い原子核に対する多体系シュレディンガー方程式の求解は、計算量的に極めて困難である。

• 無心殻モデル (NCSM): 軽い原子核には成功しているものの、NCSMはモデル空間の次元の指数関数的な増大によって制限されており、現在のスーパーコンピュータの計算能力（$\sim 10^{12}$ 次元の次元に達する）を急速に超えてしまう。その結果、NCSMはより大きな系に対しては、主に外挿に依存することになる。
• 打ち切り多体系展開: 多体摂動論（MBPT）、イン・メディウム・シミラリティ・リノーマライゼーション・グループ（IMSRG）、および結合クラスター（CC）理論などの手法は、多項式スケールを実現しており、より重い原子核の計算を可能にしている。しかし、これらは打ち切り誤差を伴う。高い精度を得るには、これらの展開をより高次の項まで進める必要があり、数学的な複雑さと計算コストが増大する。さらに、特に強いクラスタリングや集団的相関を持つ系において、打ち切られたスキームでは、信頼できる多体不確かさの推定を行うことが困難である。

全構成相互作用（FCI）の精度を維持しつつ、ハミルトニアンを明示的に格納することなく、打ち切られた展開手法を検証し導くためのベンチマークとして機能する手法が切実に求められている。

手法
著者らは、核構造の領域にフル構成相互作用量子モンテカルロ（FCIQMC）法を導入する。FCIQMCは、虚時間シュレディンガー方程式を解く確率論的な射影法である：
$$-\frac{d}{d\tau} |\Psi(\tau)\rangle = (\hat{H} - E_0) |\Psi(\tau)\rangle$$
実装における主な特徴は以下の通りである：

• 確率的サンプリング: 波動関数は、全構成空間内のスレーター行列式に分布する「ウォーカー」（符号付き粒子）のアンサンブルによって表される。ある行列式の係数は、その行列式上のウォーカーの正味の符号付き集団によって決定される。
• 伝播ステップ: 虚時間ステップ（$\Delta\tau$）ごとに、ウォーカーのアンサンブルは3つの確率的ステップを経て進化する：
  1. スポーニング (Spawning): 接続された行列式上に新しいウォーカーを作成するために、非対角項のハミルトニアン結合（$H_{ij}$）をサンプリングする。
  2. 死滅/クローニング (Death/Cloning): パラメータ $S$ によってシフトされた対角項（$H_{ii}$）をサンプリングし、現在の行列式上のウォーカーの集団を調整する。
  3. 消滅 (Annihilation): 同じ行列式上に存在する反対の符号を持つウォーカーを打ち消す。これは、固定ノード近似やガイド波動関数を必要とせずにフェルミオン符号問題を克服するために極めて重要である。
• イニシエーター近似: 大きなモデル空間における符号問題を制御するために、著者らは適応型シフト補正を組み合わせたイニシエーター近似を採用している。これにより、ウォーカー数（$N_w$）が増加するにつれて、FCI極限の制御可能かつ改善可能な推定が可能となる。
• システム設定: 計算は、球形調和振動子基底（$\hbar\omega = 20$ MeV）内において、最適化されたカイラル二体相互作用（N2LOopt）を用いて、$^4$He、$^8$Be、$^{12}$C、および $^{16}$O に対して行われた。これらのFCIQMC計算では、構成の打ち切りは適用されていない。

主要な結果
本研究は、選択された原子核について、サブパーセントレベルの多体不確かさを伴う基底状態エネルギー（$E$）および電荷半径（$R_{ch}$）を提示している。

• 小空間における検証: 小さなモデル空間（$e_{max}=2$）において、FCIQMCの結果は直接対角化によって得られた厳密なFCI結果と一致し、確率的サンプリングの正確性が検証された。
• 大きなモデル空間での性能: 大きなモデル空間（$e_{max}=10$）において、$^{16}$O のヒルベルト空間次元が $\sim 10^{33}$ に達し、直接対角化が不可能な場合でも、FCIQMCは正確な結果を提供した。
  • $^4$He について、FCIQMCの結果は、事実上厳密なNCSMの外挿と完全に一致した。
  • NCSMが収束に失敗するより重い原子核（$^8$Be、$^{12}$C、$^{16}$O）において、FCIQMCは精密な結果をもたらした。
• 打ち切り手法のベンチマーク:
  • MBPT: FCIQMCの結果からの偏差が最大であった。
  • IMSRG: 一般には妥当であるが、$^8$Be において顕著な偏差を示した。ここでは高次補正（IMSRG(3f2)+T）が厳密な値からさらに逸脱しており、欠落している集団的相関の存在を示唆している。
  • 結合クラスター (CC): CCSD から CCSDT-3 にかけて系統的な改善を示した。しかし、$^{16}$O の電荷半径については、主要な反復トリプル（iterative triples）の導入がFCIQMCと比較して大きな偏差をもたらしており、特定の観測量における収束の難しさを浮き彫りにした。
• クラスタリング効果: 顕著な分子状／クラスタリング構造を持つ $^8$Be の結果は、参照系に基づく展開（単一参照IMSRGなど）が、高次の集団的相関を省略しているために収束に苦戦することを実証した。あらゆる励起ランクの相関を基底内に含めることにより、FCIQMCはこれらのクラスタリング状態に対して精度を維持した。

意義と主張
本論文は、FCIQMCがアブイニシオ核構造計算における新しい標準を確立すると主張している：

1. 確率的フル空間ソルバー: 従来のNCSMの到達範囲を超える系を、全構成空間における厳密な波動関数を直接サンプリングすることによって扱うことができる。
2. 定量化された不確かさ: 不確かさの推定が困難な打ち切り展開手法とは異なり、FCIQMCは定量化されたサブパーセントの多体不確かさを伴う結果を提供する。
3. 理論のベンチマーク: その結果は、カイラル相互作用の精度および多項式スケールの多体系手法（MBPT、IMSRG、CC）の収束を評価するための厳密なベンチマークとして機能する。本研究は、打ち切られたスキームでしばしば欠落している高次の相関の決定的な役割を明示的に示している。
4. 将来への影響: 著者らは、この手法がより高精度な次世代の多体系手法の開発を促進し、カルシウム鎖における電荷半径の進化や、超選択的フェルミ崩壊におけるアイソスピン破壊補正など、高次の多体効果が中心となる核理論における長年の謎を解決する助けになると断言している。

本研究は、強い相互作用下にある有限核のアブイニシオ計算にFCIQMCを適用した最初の事例であり、前例のない精度で核力と基本対称性の理解へと向かう重要な一歩である。