Proton and neutron electromagnetic form factors from lattice QCD in the continuum limit

技術要約：連続極限における陽子および中性子の電磁形因子

問題提起
陽子および中性子の電磁形因子は、ハドロンの内部構造を記述する基本的な量である。陽子の形因子は実験的に非常によく測定されている一方で、中性子の形因子は自由中性子標的の欠如により、電子・重水素または電子・ヘリウム散乱からの間接的な抽出に頼らざるを得ないため、精度が低い。格子量子色力学（Lattice QCD）は、これらの量を計算するための第一原理的なアプローチを提供する。しかし、従来の格子計算には、連続極限への外挿が欠けていたり、非物理的なクォーク質量からのカイラル外挿に依存していたり、あるいは（アイソスカラー項を含む）非連結寄与（disconnected contributions）を十分な精度で含めていなかったりするという限界があった。本研究は、物理的なクォーク質量を持つアンサンブルを用いて、物理的なパイ中間子質量において直接的に形因子を計算し、連続極限を取ることで、これらの課題に対処する。

手法
著者らは、$N_f=2+1+1$ のツイスト質量クローバー改善フェルミオンを用いた格子QCDを用いている。シミュレーションでは、$a=0.080$ fm、$0.068$ fm、および $0.057$ fm の格子間隔を持つ3つのゲージアンサンブルを使用しており、これらはすべて「物理点（physical point）」に調整されている。

• 相関関数: 核子行列要素は、2点相関関数および3点相関関数から抽出される。励起状態の混入を抑制するために、著者らは 0.5 fm から 1.5 fm の範囲の複数のシンク・ソース時間間隔を利用している。
• 励起状態解析: 厳密なマルチステート・フィット解析が行われる。標準的なプラトー・フィットとは異なり、この手法では、電磁電流が補間場（interpolating field）とは異なる多粒子状態（例：パイ中間子・核子状態）に結合する可能性を考慮し、2点関数と3点関数で異なる励起状態エネルギーを許容する。
• 非連結寄与: アイソスカラー寄与（クォークループ、すなわち非連結図を含む）は、統計的ノイズを軽減するために高統計手法を用いて計算される。これには、低モード・デフレーション（low-mode deflation）、カラー・スピン希釈（color-spin dilution）、および階層的プロービング（hierarchical probing）が含まれる。非連結寄与は局所ベクトル電流を用いて計算され、これには繰り込みが必要であるが、連結寄与は保存ベクトル電流を用いる。
• 繰り込み: フレーバー非シングレットおよびシングレットのベクトル演算子の繰り込み因子は、RI'-MOM スキームを用いて決定される。解析により、シングレットと非シングレットの因子が一致することが確認されており、これによりアイソスカラー電流の繰り込みが簡略化される。
• 観測量の抽出: 形因子は、ダイポールおよび $z$-展開のアンサッツ（Ansatz）の両方を用いて抽出される。著者らは、「ワンステップ」のアプローチを採用し、$Q^2$ 依存性と連続極限への外挿（$a^2$ に関して）を同時にフィッティングする。異なるフィット範囲やパラメータ化の結果を組み合わせるためにモデル平均化が用いられ、最終的な統計誤差および系統誤差を導出するために赤池情報量基準（AIC）に基づいた重み付けが行われる。

主要な貢献

1. 物理点における連続極限: 本研究は、物理的なパイ中間子質量を持つアンサンブルを用いて、陽子の電磁形因子の連続極限を達成した最初期の計算の一つであり、これによりカイラル外挿の必要性を排除している。
2. 包括的な励起状態制御: 異なる時間間隔を用いたマルチステート・フィットの使用、および2点関数と3点関数で異なる励起状態エネルギーを許容することにより、従来の研究よりも堅牢な基底状態行列要素の抽出を実現している。
3. 非連結寄与の包含: アイソスカラー形因子に対する高統計な非連結寄与の解析を含んでおり、階層的プロービングやデフレーションなどの高度なノイズ低減技術を用いることで、連結寄与に匹敵する精度を達成している。
4. 系統誤差の推定: フィット範囲、モデル依存性（ダイポール vs $z$-展開）、および連続極限の外挿から生じる系統的な不確かさの詳細な内訳を提供し、モデル平均化を用いてこれらの影響を定量化している。

結果
著者らは、連続極限における陽子および中性子の電気半径、磁気半径、および磁気モーメントを提示している。結果は、統計誤差および系統誤差（後者はフィットおよび連続極限の外挿による不確かさを含む）と共に記載されている：

• 陽子の電気半径: $\sqrt{\langle r^2_E \rangle_p} = 0.860(38)(23)$ fm
• 中性子の電気半径: $\langle r^2_E \rangle_n = -0.147(48)$ fm$^2$
• 陽子の磁気半径: $\sqrt{\langle r^2_M \rangle_p} = 0.870(53)(15)$ fm
• 中性子の磁気半径: $\sqrt{\langle r^2_M \rangle_n} = 0.913(67)(19)$ fm
• 陽子の磁気モーメント: $\mu_p = 2.849(92)(52)$
• 中性子の磁気モーメント: $\mu_n = -1.819(76)(29)$

さらに、ハイパーファイン分裂の計算に関連するゼマッハ半径およびフライヤー半径も計算している：

• 陽子のゼマッハ半径: $r_Z^p = 1.091(34)$ fm
• 中性子のゼマッハ半径: $r_Z^n = -0.0405(81)$ fm

半径および磁気モーメントに関する結果は、統計的および系統的な誤差の範囲内で、粒子データグループ（PDG）の値と一致している。計算された陽子の電気半径の不確かさは、歴史的な電子・陽子散乱による大きな値と、ミューオン的水素分光法による小さな値の両方と重なっている。

意義と主張
本論文は、これらの結果が格子QCDの定式化のベンチマークとなり、非摂動計算を実験データに対して検証するものであると主張している。物理的なパイ中間子質量において直接的に、かつ連続極限で形因子を計算することで、本研究は、特に中性子のように実験的に知られていないハドロン構造の量を評価する格子QCDの能力を実証している。

著者らは、自らの結果が実験データと一致しているものの、陽子の電気半径の精度は現在、統計的および系統的な誤差によって制限されていることを控えめに述べている。彼らは、個々のアンサンブルの精度を高め、より小さな格子間隔およびより大きな物理的体積（$Q^2=0$ に近い運動量にアクセスするため）を持つアンサンブルを含めることが、これらの誤差をさらに減少させ、陽子半径に関する異なる実験的決定の間の不一致を解決する可能性があると示唆している。また、本研究は、格子QCDにおける非連結寄与の包含および厳密な励起状態解析を行うための手法の検証としても機能している。