Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions

宇宙で最も激しい衝突の中心、原子核が光速に近い速度で衝突する場所では、物質は通常の条件下ではあり得ないような振る舞いを見せます。大型ハドロン衝突型加速器のような施設で陽子や鉛の原子核が衝突すると、クォーク・グルーオン・プラズマとして知られる、一瞬で消え去る超高温の亜原子粒子のスープが生成されます。ほんのわずかな間、このプラズマは流体のように流れ、膨張し、冷却された後、最終的に検出器に記録される粒子へと凍結します。このプロセスを理解するために、物理学者は流体の動きや変化を記述する「流体力学」と呼ばれる一連の数学的規則に依拠しています。しかし、これらの衝突の極限的な相対論的速度にこれらの規則を適用することは、長年の課題でした。この流体をモデル化するために使用される標準的な手法では、数学的な安定性を保ち、方程式が破綻するのを防ぐために、追加の目に見えない変数が必要になることがよくあります。これらの追加変数はセーフティネットのような役割を果たしますが、それは複雑さを増大させ、どの特性が流体の実在する性質であり、どれが単なる数学的モデルの産物（アーティファクト）であるのかを判別することを困難にします。

研究チームは、これらのシミュレーションを実行するための新しい方法を開発し、不要なセーフティネットを取り除くことで、流体の挙動をよりクリーンで直接的なモデルとして作成しました。最近の研究において、ジャン＝フランソワ・パケ、マヤンク・シン、デレク・ティーニーは、「密度フレーム（Density Frame）」と呼ぶ手法を紹介しました。流体の流れとその内部応力を別々の進化する対象として追跡するのではなく、この新しいアプローチでは、流体のエネルギーと運動量のみが時間の経過とともに変化するものとして扱います。流れに抵抗する内部摩擦を含む他のすべての要素は、その瞬間の流体の動きに基づいて即座に計算されます。これにより、以前のモデルを煩雑にしていた、非物理的な追加変数の必要性が排除されました。研究者たちは、真空中で粒子がどのように振る舞うかという既知の解と比較したり、この分野で最も広く使われているMUSICと呼ばれるシミュレーション・ソフトウェアと並行して実行したりすることで、この新しいエンジンをテストしました。その結果、新しい手法は古い手法よりも流体の微細で鋭い波立ちをより速く平滑化するものの、衝突から出現する粒子の最終的な結果は驚くほど似ていることが分かりました。このことは、古いモデルにおける複雑な追加ステップが、クォーク・グルーオン・プラズマがどのように流れるかという全体像を理解する上で、必ずしも必要ではない可能性を示唆しています。

この研究の核心は、研究者が流体内のエントロピー、すなわち無秩序をどのように扱うかという点にあります。実際のあらゆる流体において、摩擦は熱を生み出し、無秩序を増大させますが、これは継続的に起こらなければならないプロセスです。新しい手法は、「散逸関数（dissipative function）」と呼ばれる数学的ツールを中心に構築されており、これが無秩序のマスター台帳として機能します。コンピューターコードを設計する際、シミュレーションの各ステップがこの無秩序を常に一貫した方法で増加させるようにすることで、研究者たちは、流体が非常に薄くなったり、あるいは乱流状態になったりしても、モデルが安定し続けることを保証しました。この安定性により、数学的なクラッシュを防ぐための人工的な「レギュレーター（調整器）」を必要とせずに、衝突のシミュレーションを行うことが可能になりました。テストにおいて、彼らは衝突直前の原子核の混沌とした状態を模した「IP-Glasma」と呼ばれる特定の初期条件から開始しました。そして、温度と流速がどのように進化するかを観察しながら、シミュレーションを時間軸に沿って進めました。

新しい密度フレームの結果を標準的なMUSICコードと比較した際、両者が流体の内部摩擦をどのように扱うかにおいて明確な違いがあることに気づきました。標準的なコードでは、流体の流れに対する抵抗（剪断応力として知られるもの）は、重い物体が動き出したり止まったりするのに時間がかかるのと同様に、変化に対して適応するまでに少しの時間を要します。新しい密度フレームでは、この抵抗は即座に適応します。この違いは、体積粘性（流体が押しつぶされたり引き伸ばされたりすることへの抵抗）においてさらに顕著でした。標準的なモデルでは、特に流体が熱いプラズマからより冷たいガスへと転移する点付近において、この抵抗は時間の経過とともにゆっくりと緩和されます。一方、新しいモデルでは、それは即座に反応します。プロセスの途中におけるこれらの違いにもかかわらず、最終的な結果は驚くほど一致していました。研究者が衝突の終わりに生成される粒子（生成された粒子の総数、平均速度、および異なる方向への広がりなど）を調査したとき、両方の手法はほぼ同一の結果を示しました。

また、研究では、流体が押しつぶされることへの抵抗がないシナリオ、つまり体積粘性がゼロに設定された場合についても調査を行いました。この場合、両手法はほぼ完璧に一致しました。体積粘性が温度に依存する場合（これは、転移点付近での流体の複雑な振る舞いを模したものですが）、衝突領域の外縁部を中心に小さな差異が現れました。しかし、これらの差異は一般に、体積粘性自体の影響よりも小さいものでした。研究者たちは、古いモデルが安定性を維持するために依存している複雑な非流体力学的モードは、衝突の最終的な観測可能な特性を決定する上で主要な役割を果たしていないと考えています。この発見は重要です。なぜなら、現在の流体の特性の測定における不確実性は、流体自体の物理学によるものではなく、むしろシミュレーション・ツールが抱えている余分な数学的荷物によるものである可能性を示唆しているからです。

これらの余分な変数を取り除くことで、密度フレームはクォーク・グルーオン・プラズマの根本的な性質を理解するための、より直接的な経路を提供します。研究者たちは、自身のコードを公開し、文書化しており、他の科学者が既存のシミュレーション・パッケージの流体力学エンジンとしてこれを使用できるようにしました。これにより、数学的モデルに組み込まれた恣意的な時間遅延ではなく、粘性や状態方程式といった流体の実際の物理的特性のみを調整パラメータとする、新しい種類の解析への扉が開かれます。著者らは、このよりクリーンなアプローチが、流体が十分に発達していない、より小さな衝突系を研究する際に特に有用であると示唆しています。そこでは、これらの余分な数学的モードの影響がより顕著になる可能性があるためです。最終的に、この研究は、高エネルギー衝突型加速器で収集されたデータから、クォーク・グルーオン・プラズマの真の輸送特性を抽出するための、より堅牢で透明性の高いツールを提供し、宇宙で最も極端な流体の完全な理解へと私たちを近づけるものです。