Modeling quasielastic lepton-nucleus interactions with ab initio spectral functions from infinite nuclear matter

技術要約：Ab Initio スペクトル関数を用いた準弾性レプトン・原子核相互作用のモデリング

問題提起
次世代のニュートリノ振動実験（DUNE、T2HKなど）には、系統誤差を制御するために、ニュートリノ・原子核相互作用の精密な理論モデリングが必要である。中心となる課題は、単一の枠組みの中で異なる相互作用チャネルと最終状態相互作用（FSI）を整合的に扱うことである。スペクトル関数（SF）形式は広く用いられているが、FSIはしばしば無視されるか、あるいは光学ポテンシャルなどを通じて近似的に扱われる。さらに、既存の微視的な電弱断面積の計算は、グリーン関数モンテカルロ法や結合クラスター理論を用いた包括的（inclusive）な観測量に限定されており、粒子セクターにおける一貫したFSIの取り扱いを欠いていることが多い。ニュートリノ実験に関連する低・中運動量転移領域において、準弾性（QE）散乱の記述を改善するために、初期状態（正孔）と最終状態（粒子）の両方に対してab initio スペクトル関数を組み込んだ枠組みが必要とされている。

手法
著者らは、包括的なレプトン・原子核散乱をモデル化するために、局所密度近似（LDA）とインパルス近似（IA）に基づいた枠組みを提案している。

1. 理論的基礎:

   • スペクトル関数 (SF): 本研究では、自己一貫グリーン関数（SCGF）理論のGorkovフレームワークを用い、無限対称核物質（NM）から導出されたab initio スペクトル関数を利用する。これらの計算には、カイラル有効場理論に基づく現実的な核ハミルトニアン（具体的には $\Delta$N2LOGO(394)、$\Delta$N2LOGO(450)、および N2LOsat(450)）と、"ADC(3)-D" 切断レベルにおける代数的図式構成（ADC）法を用いている。
   • LDAの実装: 有限核（$^{12}$C, $^{16}$O, $^{40}$Ca）は、無限核物質のSFを核密度プロファイル $\rho(r)$ にわたって積分することによってモデル化される。この枠組みには、初期束縛核子を記述する正孔SF（$S_h$）と、FSIを考慮した粒子SF（$S_p$）の両方が含まれている。
   • 2つのモデル: 以下の2つの運動学的領域を検討している。
     • IA+FSI: 正孔セクターと粒子セクターの両方にab initio 核物質SFを使用する。
     • PWIA (平面波インパルス近似): ab initio 正孔SFを使用するが、粒子SFを相対論的フェルミガスモデル（パウリブロッキングのみを考慮）に置き換える。これは、非相対論的なab initio SFが信頼性を欠く高運動量転移領域に適している。

2. ニューラルネットワークによる補間:

   • SCGF計算はエネルギー、運動量、密度、および離散的な点においてSFを提供するため、著者らはこれらのデータを連続的に補間するための人工ニューラルネットワーク（ANN）を開発した。
   • 3つのハミルトニアンと2つのアイソスピン（陽子/中性子）の組み合わせごとに、計6つの独立したネットワークを訓練した。
   • 入力特徴量は、エネルギー ($E$)、運動量 ($|p|$)、密度 ($\rho$)、および平滑化因子 ($\Gamma$) である。出力はスペクトル関数の値である。
   • 正規化和則（$\int dE S(E) = 1$）を強制するためのペナルティ項を組み込んだ、物理学に基づいた損失関数が採用された。

3. 断面積の計算:

   • ハドロンテンソルは、正孔と粒子のスペクトル関数の積を核体積にわたって積分することにより構築される。
   • 本枠組みは、電磁（EM）応答および、ニュートリノ散乱のフラックス平均全断面積を含む荷電流（CC）応答を計算する。

主な貢献

• 一貫したFSIの取り扱い: 現象論的な光学ポテンシャルに頼るのではなく、ab initio から導出された粒子スペクトル関数を用いて、FSIの一貫した取り扱いを導入した。
• ニューラルネットワークによる補間: 高次元のab initio スペクトル関数データを補間するためのANNの開発により、任意の運動学的変数に対する断面積計算が可能となった。
• ベンチマーキング: モデルは、$^{12}$C, $^{16}$O, $^{40}$Ca の電子散乱データ、および $^{16}$O に対する ab initio LIT-CC（ローレンツ積分変換を用いた結合クラスター法）計算と比較して検証されている。
• ニュートリノ予測: 本枠組みを $^{16}$O の荷電流応答および、$\nu_\mu$-$^{12}$C のフラックス平均全断面積の計算に適用し、LSND の実験データと比較した。

結果

• 電子散乱:
  • $^{12}$C について、IA+FSI モデルは、様々な運動量転移における $\omega \lesssim \omega_{QE}$ での縦波応答（$R_L$）のローゼンブース分離された実験データと極めて良好な一致を示した。
  • FSI の導入は、PWIA モデルと比較してピーク位置と強度を補正する。PWIA はQEピークを高エネルギー側にシフトさせ、強度を過大評価する傾向がある。
  • 横波応答（$R_T$）は系統的に過小評価されており、これは計算における二体電流の除外と一致している。
  • $^{16}$O について、IA+FSI の結果は QE ピーク付近で LIT-CC 計算とよく一致しているが、LDA アプローチは、LIT-CC で見られる詳細な殻構造の閾値効果を欠いている。
  • $^{40}$Ca について、モデルは実験データを良好に再現しており、縦波および横波の両方のチャネルで良好な一致を示している。
• ニュートリノ散乱:
  • IA+FSI モデルで計算された $\nu_\mu$-$^{12}$C のフラックス平均全断面積は、LSND の実験データと一致している。
  • PWIA モデルは、IA+FSI よりも約20%低い断面積を与えており、これは低運動量において QE ピークが積分窓の外側にシフトすることに起因する。
  • 理論的不確定性（ハミルトニアンの選択および平滑化因子 $\Gamma$ によるもの）は、運動量転移 $|q| \approx 0.2\text{--}0.6$ GeV における QE ピークにおいて、約7〜15%と評価されている。

意義と主張
著者らは、本枠組みが「我々の近似の妥当性の範囲内（具体的には低・中運動量転移領域）において実験データと良好な一致を示す」と主張している。その意義は以下の通りである：

1. 改善された一貫性: LDA 内で、初期状態と最終状態の両方に対して微視的な ab initio 記述を提供すること。これは、詳細な殻構造が広範なエネルギーフラックスによって平滑化されるニュートリノ振動実験にとって特に適している。
2. 拡張性: 本枠組みは、追加の動力学的メカニズム（種々のインメディウム・ピオン生成など）や他の核ハミルトニアンへと容易に拡張できるように設計されている。
3. 検証: 電子散乱データおよび LIT-CC 計算に対する成功裏のベンチマーキングは、アイソスピン対称な原子核における包括的観測量に対する本アプローチの妥当性を立証している。

論文では、限界についても控えめに述べている。すなわち、LDA は詳細な殻構造を犠牲にしており（これにより、排他的な $(e, e'p)$ データには適さなくなる）、また、相対論的補正や QE を超えるメカニズム（ピオン生成など）が支配的となる非常に高い運動量転移ではモデルが破綻する。今後の課題として、SCGF 計算を（DUNE に関連する $^{40}$Ar のために）アイソスピン非対称物質へと拡張すること、および本枠組みを核内カスケードモデルと結合させることが挙げられている。