[论文解读] 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维数值框架，用于在非洛伦兹不变近似下模拟相对论重离子碰撞中的Glasma，采用有限核模型和一种新型规范协变粒子-网格方法。关键结果表明，尽管压强各向异性依然显著，但由于显式破坏了洛伦兹对称性，能量密度表现出快度依赖性，且通过简单参数即可与实验数据匹配。

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维中模拟Glasma的数值方法，以摆脱标准的洛伦兹不变性近似。
• 将碰撞核的有限纵向厚度纳入颜色玻璃凝聚体（CGC）框架中，从而破坏洛伦兹不变性。
• 在使用大时间步长的真实时间格点规范理论模拟中，确保数值稳定性和规范协变性。
• 研究有限核厚度对Glasma可观测量（如能量密度和压强各向异性）的影响。
• 将模型与重离子碰撞中能量密度的实验快度分布进行对比测试。

提出的方法

• 通过从共动参考系转换到实验室参考系或质心参考系，将真实时间格点规范理论推广至3+1维。
• 实现了彩色粒子-网格（CPIC）方法，以在模拟中显式包含有限厚度核的颜色电荷。
• 通过参数化MV模型，将有限核厚度引入标准杨-米尔斯方程的修正作用量中。
• 开发了一种新型半隐式、规范协变的数值格式，即使在大时间步长下也能保持高斯约束的守恒。
• 采用平面波假设和色散分析，证明了该新格式在CFL条件下的稳定性。
• 在3+1维中进行模拟，通过改变核厚度，研究能量密度的快度依赖性及压强特性。

实验结果

研究问题

• RQ1 当有限核厚度破坏洛伦兹不变性时，Glasma结构如何变化？
• RQ2 在3+1D Glasma中，压强分量与2+1D洛伦兹不变情况相比偏离程度如何？
• RQ3 由于有限核厚度的存在，局部静止参考系中的能量密度分布是否表现出快度依赖性？
• RQ4 模拟参数如何影响3+1D真实时间格点规范理论中的数值稳定性？
• RQ5 该模型生成的快度分布能否在定量上与RHIC和LHC的实验数据相匹配？

主要发现

• 当核厚度趋近于零时，新开发的3+1D数值格式能成功复现洛伦兹不变Glasma作为极限情况。
• 3+1D Glasma中的压强各向异性依然显著，与2+1D洛伦兹不变情况相比无明显偏离。
• 在局部静止参考系中，能量密度表现出强烈的快度依赖性，这是由于显式破坏了洛伦兹不变性所致。
• 能量密度快度分布的宽度由核的纵向厚度控制，与实验数据具有良好的定性一致性。
• 发现该数值格式存在一种此前未识别的不稳定性；该问题被追溯至时间导数的处理方式，并通过一种新型规范协变、约束保持的算法得以解决。
• 该新格式在CFL条件下被证明是稳定的，即使在大时间步长下也能保持高斯约束守恒，确保了长期模拟的可靠性。