BeAGLE++: A modern modular framework for eA physics at the EIC

あらゆる原子の核心深くには、物質の最小の構成要素を結合させるほど強力な力によって支配された世界が存在しています。この領域は量子色力学と呼ばれる理論によって記述されており、クォークやグルーオンとして知られる粒子が、どのようにして陽子や中性子を形成するために結びついているのかを説明しています。科学者たちは、これらの粒子が自由に浮遊しているときの振る舞いについては古くから理解してきましたが、原子核の中に密集しているときの振る舞いは、現代物理学における大きな謎の一つであり続けています。高速の電子が重い原子核に衝突すると、それは単に単一の粒子に跳ね返るだけではありません。それは原子核構造全体に波及する複雑な連鎖反応を引き起こします。エネルギーがどのように広がり、どのように新しい粒子が誕生し、そして原子核がどのように砕け散ったり、あるいは形を変えたりするのかを正確に理解することは、私たちの宇宙の根本的な規則を解明する上で極めて重要です。これらの激しい衝突を研究するために、研究者たちはコリジョン（衝突型加速器）と呼ばれる強力な装置に頼っていますが、収集されたデータを理解するためには、これらのイベントを極めて高い精度でシミュレートできる高度なコンピュータプログラムが必要となります。

長年にわたり、BeAGLEとして知られるコンピュータプログラムが、これら電子・原子核衝突をシミュレートするための主要なツールとして機能してきました。これにより、科学者たちは電子ビームが重い原子からなるターゲットに衝突したときに何が起こるかを予測することができました。しかし、オリジナルのソフトウェアは古い技術に基づいて構築されていたため、粒子が原子核内部でどのように相互作用するかについての新しいアイデアを導入したり、テストしたりすることが困難でした。これに対処するため、物理学者のチームは、次世代の継承モデルであるBeAGLE++を開発しました。この新しいフレームワークは、単なるマイナーアップデートではありません。これは、衝突の直後の混沌とした余波をシミュレーションがどのように扱うかという点における、完全な再構築です。研究者たちは、現代的なコーディング標準を用いてプログラムをゼロから作り直し、原子核の密な環境内での粒子の動きと衝突の追跡方法について、より現実的な手法を導入することを可能にしました。

この新しいアプローチの核心は、電子が原子核内の核子を直撃した直後の、一瞬の出来事をプログラムがどのように処理するかという点にあります。旧モデルでは、シミュレーションは原子核を、粒子が固定された順序で動く「凍結された写真」のような静的なスナップショットとして扱っていました。新しいBeAGLE++は、原子核を「生きている、呼吸しているシステム」として扱います。この更新された視点では、原子核内部の粒子は常に動き、相互作用するダイナミックな環境にあります。衝突が発生すると、プログラムはそれらを世代ごとにグループ化するのではなく、発生した正確な時間順序に従って、あらゆる相互作用を追跡します。これは、もし最初の衝突で生成された粒子が、二番目の衝突が完了する前に別の粒子に衝突するほど高速であった場合、シミュレーションがその事象を考慮することを意味します。この時間順序に基づいたアプローチは、より詳細な粒子相互作用のリストと組み合わさることで、核の崩壊に関するより豊かで正確な描写を生み出します。

最も重要な改善点の一つは、プログラムが新しい粒子の生成をどのように扱うかに関わるものです。衝突によって新しい粒子が生成されるとき、それは即座に他のものに跳ね返ることができる「固形物」になるわけではなく、形成されるまでに極めて短い時間を要します。新しいモデルは、この「形成時間」を注意深く追跡し、粒子が完全に実体化するまで相互作用を許可しないように制御しています。この微細ながらも決定的な詳細は、シミュレーションの結果、特にエネルギーが原子核内にどのように堆積するかという点に変化をもたらします。研究者たちは、FermilabのE665実験で測定されたミューオン・鉛（µ–Pb）衝突の前方中性子多重度データを用いて、この形成時間のパラメータを校正しました。その結果、新しいモデルは、旧バージョンと比較して、衝突後に原子核から放出される中性子の数が多いことを予測することがわかりました。この増加はエラーではなく、むしろ新しいモデルが衝突中のエネルギー転送をより正確に捉えており、その結果、励起状態となった原子核が冷却過程でより多くの粒子を放出していることを示唆しています。

新しいモデルが他の観測量に対しても正しく機能しているかを検証するため、チームは同じE665実験のデータを用い、今回はミューオン・キセノン（µXe）衝突を用いて比較を行いました。彼らは、生成された荷電粒子の数とその空間分布を調査しました。結果として、BeAGLE++は実験データ、特に原子核自体が崩壊する領域において非常によく一致していることが示されました。旧モデルでは、シミュレーションは砕け散った原子核から放出される粒子の数を過小評価することがよくありました。新しいバージョンはこの問題を修正し、現実により近い結果を提供しました。しかし、研究者たちは、最初に衝突を受けた粒子が飛び去っていく領域においては、新モデルでも依然として生成される粒子の数がわずかに過小評価されていることも指摘しました。これは、新しいフレームワークが核の破片に関する多くの問題を解決した一方で、初期の衝突に関連する複雑な相互作用については、まださらなる研究が必要であることを示唆しています。

また、チームは新しいモデルが軽い原子核クラスターの生成を、以前のモデルとは異なる方法で処理していることも発見しました。旧式のシミュレーションでは、原子核から叩き出された陽子は通常、個々の粒子として飛び去っていました。しかし、新しいシミュレーションでは、これらの陽子の一部がヘリウムや重水素のような、小さく安定した粒子のグループへと結合します。この変化は、検出される自由陽子の総数に影響を与え、粒子の分布をE665実験で実際に観察されたものにより適合するようにシフトさせます。この発見は、新しいフレームワークの主要な強みを浮き彫りにしています。すなわち、個々の粒子と小さなクラスターを区別することができ、高エネルギー衝撃の後に物質がどのように再編成されるかについて、より微細な視点を提供できるようになったのです。

将来を見据えて、著者らはBeAGLE++が将来の発見のための柔軟な基盤として設計されていることを強調しています。これは、新しい物理理論が開発されるにつれて容易に更新できるように構築されており、現在建設中の巨大な新施設である電子イオン衝突型加速器（EIC）にとって不可欠なツールとなります。この新しい装置は、前例のない明晰さで原子核を探索することを可能にするものであり、原子核内部の粒子の複雑なダンスを正確にモデル化できるシミュレーションツールを持つことは、データを解釈する上で極めて重要です。現在のバージョンは最も一般的なタイプの衝突に焦点を当てていますが、このフレームワークは、グルーオンの飽和や原子核の変形の影響を含むような、よりエキゾチックなシナリオを取り込む準備ができています。モジュール化され透明性の高いテスト方法を提供することで、BeAGLE++は、新しい衝突型加速器が稼働を開始したとき、科学者が原子核という奇妙で強力な世界を航海するための信頼できる地図となることを保証しています。