[論文レビュー] Simulations of the Glasma in 3+1D

David I. Müller|arXiv (Cornell University)|Apr 1, 2019

High-Energy Particle Collisions Research参考文献 127被引用数 2

ひとこと要約

本学位論文は、相対論的重イオン衝突におけるグレースマの3+1次元数値フレームワークを、ブースト不変性近似を超えて構築した。有限核モデルと、新規のゲージ共変性を有する粒子-格子法を用い、主要な結果として、圧力の異方性は依然として大きいが、明示的なブースト対称性の破れによりエネルギー密度がラピディティに依存するようになり、単純なパラメータで実験データと一致する。

ABSTRACT

Das Glasma ist ein gluonischer Zustand, welcher in relativistischen Schwerionenkollisionen erzeugt werden kann und nur für sehr kurze Zeit existiert, bevor er in das Quark-Gluon-Plasma zerfällt. Die Existenz des Glasmas ist eine Vorhersage des Farbglaskondensats (engl. „color glass condensate“ (CGC)). Das CGC ist eine klassische effektive Theorie, welche direkt aus der fundamentaleren Theorie der Quantenchromodynamik abgeleitet werden kann. In vielen analytischen und numerischen Rechnungen im Rahmen des CGCs kommt die boost-invariante Näherung zur Anwendung. In dieser Näherung nimmt man an, dass die dünne longitudinale Ausdehnung von Atomkernen (also der Lorentz-kontrahierte Durchmesser entlang der Strahlachse bzw. Bewegungsrichtung) infinitesimal ist. Eine Konsequenz dieser Näherung ist, dass das erzeugte Glasma invariant unter Lorentz-Boosts wird und somit effektiv in 2+1 Dimensionen beschrieben werden kann. Observablen, also im Prinzip beobachtbare Größen wie die Energiedichte, die Druckkomponenten und die Gluonenbesetzungszahl, sind dadurch per constructionem unabhängig von der Rapiditätskoordinate. Das Thema dieser Dissertation ist eine neue Methode zu entwickeln, mit der man die Annahme der Boost-Invarianz lockern und umgehen kann. Ich beginne mit einer Diskussion über die physikalischen Eigenschaften des Glasmas und des CGCs im boost-invarianten Fall. Als einfaches Modell für große, schwere Atomkerne kommt das McLerran-Venugopalan-Modell (MV) zum Einsatz. Die Yang-Mills-Gleichungen, welche die Dynamik des Glasmas bestimmen, können im Allgemeinen nicht analytisch gelöst werden. Numerische Lösungsmethoden sind somit oft der einzige verlässliche Weg, um das Glasma zu untersuchen. Daher wird Echtzeit-Gittereichtheorie benötigt, welche die Standardmethode zur numerischen Lösung der Yang-Mills-Gleichungen darstellt. Nach dieser Einführung in die Standardwerkzeuge, welche verwendet werden, um das boost-invariante Glasma zu beschreiben, liegt der Fokus auf der Entwicklung einer numerischen Methode für den Fall, dass Boost-Invarianz nicht mehr gilt, also wenn man für relativistische Atomkerne eine kleine, aber endliche Ausdehnung entlang der Bewegungsrichtung annimmt. Diese kleine Änderung führt dazu, dass viele Annahmen und Vereinfachungen, die noch im boost-invarianten Fall verwendet werden konnten, nicht mehr gültig sind. Inbesondere muss die Kollision im Laborbzw. Schwerpunktsystem in drei räumlichen Dimensionen beschrieben werden, anstelle des sich mit dem Glasma mitbewegenden Koordinatensystems im boost-invarianten Szenario. Dieser Koordinatensystemwechsel erfordert unter anderem, dass die Farbladungen der Atomkerne explizit berücksichtigt werden müssen. In numerischen Simulationen gelingt das mit der Particle-in-Cell-Methode, verallgemeinert auf Farbladungen. Die neue numerische Methode wird getestet, indem Kollisionen von Kernen mit endlicher Dicke simuliert werden. Als Anfangsbedingung für diese Simulationen dient ein erweitertes MV-Modell, welches einen neuen Parameter für longitudinale Ausdehnung besitzt. Es wird gezeigt, dass die neue Methode das boost-invariante Szenario als Grenzfall beschreiben kann. Weiters wird auch die Anisotropie der Druckkomponenten des dreidimensionalen Glasmas untersucht, wobei nur wenige Unterschiede zum boost-invarianten, zweidimensionalen Glasma festgestellt werden können. Anders verhält es sich mit der Brechung der Boost-Invarianz: betrachtet man die Energiedichte des Glasmas im lokalen Ruhesystem, kann eine starke Rapiditätsabhängigkeit festgestellt werden, welche durch die Dicke der kollidierenden Kerne beeinflusst wird. Im Vergleich mit experimentellen Resultaten von echten Kollisionsexperimenten zeigt sich, dass mit diesem sehr einfachen Modell realistische Rapiditätsprofile erzeugt werden können. Die numerische Methode, welche für dreidimensionale Kollisionssimulationen entwickelt wurde, ist auf die Wahl der Simulationsparameter sensibel und kann in gewissen Fällen instabil werden. Die Ursache dieser numerischen Instabilität wurde identifiziert und eine Erweiterung der ursprünglichen Methode entwickelt, welche sich als stabil erweist. Es wird gezeigt, dass diese neue Methode eichkovariant ist und das Gaußsche Gesetz während der Simulation auch für große Zeitschritte erfüllt bleibt.

研究の動機と目的

• ブースト不変性近似を緩和した3+1次元におけるグレースマの数値シミュレーション手法を開発すること。
• 衝突する核の有限な縦方向厚さをカラーグラスコンデンエート（CGC）フレームワークに組み込み、ブースト不変性を破ること。
• 大きな時間ステップを用いた時間発展演算における数値的安定性とゲージ共変性を保証すること。
• 有限核の厚さがグレースマの観測量（エネルギー密度や圧力の異方性など）に与える影響を調査すること。
• 重イオン衝突におけるエネルギー密度のラピディティプロファイルを実験データと照合してモデルを検証すること。

提案手法

• 共動座標系から実験系または質量中心系に移行することで、時間発展演算の3+1次元への適応を実現した。
• 有限厚さを持つ核の色電荷を明示的に取り入れるために、カラーパarticle-in-cell（CPIC）法を実装した。
• パラメータ化されたMVモデルを用いて、有限核の厚さを修正作用に組み込むことで、標準的なヤン・ミルズ方程式を拡張した。
• 時間ステップが大きくてもガウス制約を保存する、新しい半陽的でゲージ共変性を持つ数値スキームを開発した。
• 平面波アンザッツと分散解析を用いて、CFL条件下での新スキームの安定性を証明した。
• 核の厚さを変化させた3+1次元のシミュレーションを実施し、エネルギー密度のラピディティ依存性と圧力の分析を行った。

実験結果

リサーチクエスチョン

• RQ1 有限核の厚さによるブースト不変性の破れによって、グレースマの構造はどのように変化するか？
• RQ2 3+1次元グレースマにおける圧力成分は、ブースト不変性を仮定した2+1次元の場合と比べてどの程度ずれるか？
• RQ3 局所静止系におけるエネルギー密度プロファイルは、有限核の厚さによってラピディティ依存性を示すか？
• RQ4 3+1次元時間発展演算における数値的安定性に、シミュレーションパラメータがどのように影響を与えるか？
• RQ5 本モデルから得られるラピディティプロファイルは、RHICおよびLHCの実験データを定量的に再現できるか？

主な発見

• 新しい3+1次元数値スキームは、核の厚さがゼロに近づく極限でブースト不変性グレースマを正確に再現した。
• 3+1次元グレースマにおける圧力の異方性は依然として大きく、2+1次元ブースト不変系と顕著な差異を示さなかった。
• 局所静止系におけるエネルギー密度は、明示的なブースト不変性の破れにより強くラピディティ依存性を示した。
• エネルギー密度のラピディティプロファイルの幅は、核の縦方向厚さによって制御され、実験データと良い定性的一致を示した。
• 数値スキームに以前に特定されていなかった不安定性が見つかったが、これは時間微分の取り扱いに起因し、新しいゲージ共変性で制約を保存するアルゴリズムにより解決された。
• 新しいスキームはCFL条件下で安定であり、大きな時間ステップでもガウス制約を保存するため、長期的な信頼性が保証された。