[論文レビュー] Bayesian estimation of the low-energy constants up to fourth order in the nucleon-nucleon sector of chiral effective field theory

Isak Svensson, A. Ekström+1 more|arXiv (Cornell University)|Jun 16, 2022

Nuclear physics research studies参考文献 62被引用数 11

ひとこと要約

本稿では、核子-核子相互作用における次々に高次の有効場理論（χEFT）の低エネルギー定数（LECs）をN3LOまで推定するためのベイズ推論フレームワークを提示する。npおよびpp散乱データに条件付け、共役事前分布を用いてχEFTの切り捨て誤差を組み込むことで、LECsの31次元の後方確率密度関数（PDF）をサンプリングし、実験的nn散乱長データを用いてそれを拡張し、アイソスピン破れ効果を評価する。主な結果として、1S0チャンネルにおけるアイソスピン破れ効果は、次に次に高次の項（NNLO）でのみ明確に検出可能であることが示され、χEFTにおける系列的誤差処理の重要性が強調される。

ABSTRACT

We use Bayesian methods and Hamiltonian Monte Carlo (HMC) sampling to infer the posterior probability density function (PDF) for the low-energy constants (LECs) up to next-to-next-to-next- to-leading order (N3LO) in a chiral effective field theory ($\chi$EFT) description of the nucleon-nucleon interaction. In a first step, we condition the inference on neutron-proton and proton-proton scattering data and account for uncorrelated $\chi$EFT truncation errors. We demonstrate how to successfully sample the 31-dimensional space of LECs at N3LO using a revised HMC inference protocol. In a second step we extend the analysis by means of importance sampling and an empirical determination of the neutron-neutron scattering length to infer the posterior PDF for the leading charge-dependent contact LEC in the $^{1}S_0$ neutron-neutron interaction channel. While doing so we account for the $\chi$EFT truncation error via a conjugate prior. We use the resulting posterior PDF to sample the posterior predictive distributions for the effective range parameters in the $^{1}S_0$ wave as well as the strengths of charge-symmetry breaking and charge-independence breaking. We conclude that empirical point-estimate results of isospin breaking in the $^{1}S_0$ channel are consistent with the PDFs obtained in our Bayesian analysis and that, when accounting for $\chi$EFT truncation errors, one must go to next-to-next-to-leading order to confidently detect isospin breaking effects.

研究の動機と目的

• 核子-核子系における次々に高次の有効場理論（χEFT）の低エネルギー定数（LECs）をN3LOまで推定するための堅牢なベイズ推論プロトコルの開発。
• LECsの31次元の後方確率密度関数（PDF）をハミルトニアン・モンテカルロ（HMC）を用いてサンプリングする課題に、見直されたサンプリングプロトコルを適用することによる対処。
• npおよびpp散乱データから得られたLECsの後方PDFを、中性子-中性子（nn）散乱長の不確実性を組み込むことで拡張し、アイソスピン破れ効果の推論を可能にする。
• χEFTの切り捨て誤差が1S0部分波におけるアイソスピン破れ効果の検出可能性に与える影響を、切り捨て誤差に共役事前分布を用いることで評価する。
• 標準的な切り捨て誤差モデルの妥当性を、散乱長および有効範囲の実験データと比較することで検証する。

提案手法

• LECsの31次元の後方PDFをN3LOで効率的に探索できるよう、パrameter空間の幾何的性質を活用した見直し済みプロトコルを用いたハミルトニアン・モンテカルロ（HMC）を採用。
• エネルギー範囲Tlab ∈ [0, 290] MeVにおけるnpおよびpp散乱データに条件付け、相関のないχEFTの切り捨て誤差を考慮した尤度モデルを用いる。
• EFT展開の順次収束パターンに基づいて情報化された共役事前分布を用いて、χEFTの切り捨て誤差をモデリングし、一貫した誤差伝搬を実現する。
• 重要度サンプリングを用いて、npおよびpp LECの後方PDFを実験的nn散乱長ann^expの情報に拡張し、1S0チャンネルにおける主要な電荷依存接触項LECsの推論を可能にする。
• 拡張された後方分布からのサンプルを用いて、有効範囲パラメータおよびアイソスピン破れ観測量（例：電荷対称性破れおよび電荷不変性破れ）の後方予測分布（PPDs）を構築する。
• 実用的用途を想定し、LECsの後方分布の多変量正規近似（平均ベクトルおよび分散-共分散行列）を提示し、補足資料に格納する。

実験結果

リサーチクエスチョン

• RQ1 ハミルトニアン・モンテカルロ（HMC）は、χEFTのN3LOにおける31次元の後方確率密度関数（PDF）を効果的にサンプリングするために適応可能か？
• RQ2 実験的中性子-中性子（nn）散乱長データの組み込みが、1S0チャンネルにおける主要な電荷依存LECsの後方PDFにどのように影響を与えるか？
• RQ3 χEFTの切り捨て誤差は、1S0部分波におけるアイソスピン破れ効果の検出可能性にどの程度影響を与えるか？
• RQ4 1S0チャンネルにおける実験的点推定値のアイソスピン破れ効果は、ベイズ推論により得られた完全な後方PDFと整合的か？
• RQ5 χEFTで一般的に用いられる切り捨て誤差モデルは、N3LOに拡張され、実験的nnデータと組み合わせた場合にも依然として頑健か？

主な発見

• 見直されたHMCプロトコルは、N3LO LECの31次元の後方PDFを効果的にサンプリングでき、高次元パrameter空間における計算スケーリングと収束性が管理可能であることが示された。
• 1S0チャンネルにおける主要なnn接触項LECsの後方平均は、N3LOでeCnn_1S0 = −0.1283+0.0007−0.0007 × 10^4 GeV⁻²と推定され、95％信用区間は−0.1306から−0.1290 fm⁻¹にわたる。
• 1S0チャンネルにおけるアイソスピン破れ効果の実験的点推定値は、ベイズ推論により得られた完全な後方PDFと整合的であり、不確実性下での点推定値の使用が正当化される。
• χEFTの切り捨て誤差を考慮した場合、1S0チャンネルにおけるアイソスピン破れ効果は、NLOでは検出可能ではなく、NNLOでのみ明確に検出可能であることが示され、高次の補正項の必要性が強調される。
• NLO、NNLO、N3LOにおける散乱長および有効範囲の後方予測分布は、実験データと良好に一致しており、理論的および実験的不確実性がNNLOでほぼバランスが取れている。
• 切り捨て誤差に共役事前分布を用いることで、推論が正則化され、特にNNLOとN3LOの間で収束パターンが不規則な場合に予測の頑健性が向上した。