Global Analysis of Nuclear Parton Distribution Functions using Markov Chain Monte Carlo

あらゆる原子の内部では、陽子や中性子は固形かつ不可分な球体ではなく、クォークやグルオンと呼ばれるより小さな粒子がひしめき合う都市のようなものです。これらの構成要素は総称してパルトンと呼ばれ、強い核力によって結び付けられながら、光速に近い速度で飛び回っています。恒星の中心から大型ハドロン衝突型加速器（LHC）内部の衝突に至るまで、最も極限的な環境において物質がどのように振る舞うかを理解するために、物理学者はこれらのパルトンが親粒子の運動量をどのように共有しているのかを正確に描き出さなければなりません。この地図はパルトン分布関数と呼ばれます。陽子や中性子が重い原子核の中に閉じ込められると、そのルールは変わります。混雑した核環境はパルトンの地図を歪ませ、科学者が核パルトン分布関数と呼ぶ複雑で修正された景観を作り出します。これらの地図を決定することは高エネルギー実験の結果を予測する上で極めて重要ですが、データが乏しく、最適な適合を見つけるために用いられる数学的な景観がしばしば起伏に富んでいるため、非常に困難であることで知られています。

何十年もの間、研究者たちはこの地形をナビゲートするための標準的な数学的ショートカットに頼ってきました。ヘッシアン法として知られるこの手法は、最適解の周囲の景観が、完璧な放物線のような、滑らかで対称的なボウル状であることを前提としています。このボウルの側面の傾斜を測定することで、科学者は自分たちの地図がどの程度不確実であるかを推定できます。しかし、この仮定はデータが豊富であり、変数間の関係が単純である場合にのみ最もよく機能します。重い原子核の場合、実験データが限られており、異なる種類の粒子間のつながりが絡み合っているため、この滑らかなボウルという仮定はしばしば崩壊します。真の景観は、デコボコしていたり、傾いていたり、あるいは解が隠れている可能性のある複数の離れた谷を含んでいたりするかもしれません。もし科学者が、実際にはギザギザであるにもかかわらず景観が滑らかであると仮定してしまうと、彼らの不確実性の推定は誤解を招くものとなり、原子核の真の複雑さを覆い隠してしまう可能性があります。

これに対処するため、ナシム・デラクシャニアン率いる研究チームは、マルコフ連鎖モンテカルロ法と呼ばれる手法を用いて、これらの核の地図を探索する新しい方法を開発しました。この手法は、景観を滑らかなボウルと仮定する代わりに、探索を広大で未知の領域を巡る「誘導された歩行」として扱います。研究者は、数百万回のランダムなステップを踏みながら、各ステップが利用可能な実験データにどれだけ適合するかを常にチェックするコンピュータプログラムを構築しました。時間が経つにつれて、プログラムは景観の真の形状を学習し、データがよく適合する深い谷には多くの時間を費やし、データが適合しない高い峰には少ない時間を費やすようになります。これにより、彼らは確率分布全体を直接サンプリングすることができ、従来のメソッドでは見逃されてしまうであろう非対称性、複数の解、そして複雑な相関関係を捉えることができます。

研究者はこの新しいアプローチを、実験で使用される最も重く複雑な原子の一つである鉛の原子核に適用しました。彼らは、自分たちの知見の堅牢性をテストするために、2つの異なるタイプの分析を行いました。第一に、彼らは鉛核のみのデータのみを調べ、外部からの影響を受けない、適合の固有の構造を見るための「鉛のみ」の地図を作成しました。この分析は、極めて複雑で非自明な構造を明らかにしました。コンピュータによるサンプリングは、価クォーク（陽子のアイデンティティを与える核となる粒子）を記述するパラメータが、単一の滑らかな谷に位置しているわけではないことを明らかにしました。その代わりに、景観には適合が非常に良好となる複数の局所的な極小値、すなわち二次的な谷が存在していました。これらの解の分布は歪んでおり、非対称でした。つまり、不確実性はあらゆる方向に対して同一ではありませんでした。単一の対称な谷を想定する標準的なヘッシアン法は、このマルチピークの現実を捉えることができず、事実上、複雑な構造を滑らかにしてしまい、不確実性の不完全な姿を提供してしまいました。

より多くの情報が追加されたときにこの複雑な図がどのように変化するかを見るために、研究者は次に、鉛だけでなく、幅広い軽い原子核からのデータを含む第二の分析を行いました。異なる原子質量間のデータを標準的な数学的関係を通じて結びつけることで、鉛の振る舞いを他の軽い元素と結びつける新しい制約を導入しました。この「マルチ核（多核）」分析は、景観を大きく変えました。鉛のみの分析で顕著であった二次的な谷は大部分が抑制され、パラメータ空間はより規則的になり、単一の支配的な谷に近くなりました。軽い原子核のデータの包含は強力なアンカーとして機能し、クォーク分布への制約を強め、鉛の地図における不確実性を減少させました。しかし、これらの追加の制約があっても、研究者は景観が完全に左右対称ではないことを見出しました。特に価クォークの分布において、微妙な非ガウス的特徴や非対称な裾（テール）が残存しており、これは標準的なヘッシアン法であれば依然として見逃してしまうものでした。

両アプローチの直接比較は、データの扱いにおける決定的な違いを浮き彫りにしました。鉛のみの分析において、新手法はクォーク分布の不確実性の幅が、旧手法が示唆するものよりもはるかに広く、かつ非対称であることを示しました。滑らかなボウルの仮定に依存する旧手法は、起こりうる可能性の真の範囲を過小評価していました。マルチ核データが追加されると、不確実性は縮小しましたが、新手法は依然として、残存する非対称性を考慮した、より現実的で制約された不確実性を明らかにしました。興味深いことに、強い力を媒介する粒子であるグルオンの分布は、両方の分析を通じて比較的安定していました。グルオンの分布は、LHCでの鉛の衝突による直接的なデータによってすでに十分に制約されているため、軽い原子核からのデータを追加しても、クォーク分布が新しい制約に対して敏感であったのと同様に、グルオンの図を大きく変えることはありませんでした。

この研究はまた、これらの複雑な結果を他の科学者が利用可能にするための実用的な方法も示しました。新手法は、全確率分布を代表する数千の個別のサンプルを生成しますが、これを直接共有するにはデータ量が大きすぎます。研究者は、本質的な統計情報を失うことなく、この膨大なデータセットを管理可能な小さなレプリカのセットに圧縮する方法をテストしました。彼らは、中心値を得るためには一定の間隔で少数のサンプルを選択するだけでも十分に機能することを発見しましたが、サンプルサイズが小さくなりすぎると不確実性の幅を歪めてしまう可能性があることも発見しました。より洗練された圧縮アルゴリズム、すなわち、データセット全体の形状と相関関係を保持するようにサンプルを慎重に選択する手法は、最終的なセットが極端に削減されない限り、分布の複雑な特徴を保持する上でより優れた性能を示しました。

最終的に、この研究は、核パルトン分布関数の統計的景観が、これまで想定されていたよりもはるかに複雑であることを確立しました。解空間を滑らかで対称的なボウルとして近似する標準的な手法は、特に鉛のような重い原子核を扱う場合、全容を捉えるには不十分です。全景観を直接サンプリングする手法を用いることで、研究者は、複数の競合する解や非対称な不確実性といった、以前は見えなかった隠れた構造を明らかにしました。このアプローチは、核に関する未知の領域を定量化するための、より信頼性が高く透明性の高い方法を提供し、将来の高エネルギー物理学実験の予測が、データの複雑さを真に反映した基盤の上に築かれることを保証するものです。