[논문 리뷰] Simulations of the Glasma in 3+1D

David I. Müller|arXiv (Cornell University)|2019. 04. 01.

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차원으로 적응시켰다.
• 유한 두께 핵의 색 전하를 명시적으로 포함하기 위해 색 전하가 있는 입자-격자(CPIC) 방법을 구현하였다.
• 표준 양-밀스 방정식을 수정된 작용으로 확장하여 매개변수화된 MV 모델을 통해 유한한 핵 두께를 포함시켰다.
• 큰 시간 단계에서도 여전히 가우스 제약 조건을 유지하는 새로운 반-암시적, 게이지 고유성 수치 기법을 개발하였다.
• 플레인 웨이브 안사와 분산 분석을 통해 신규 기법이 CFL 조건 하에서 안정함을 증명하였다.
• 다양한 핵 두께로 3+1차원에서 시뮬레이션을 수행하여 에너지 밀도와 압력의 빠르기 의존성을 연구하였다.

실험 결과

연구 질문

• RQ1 유한한 핵 두께에 의해 보정 불변성이 깨지면 글라스마의 구조는 어떻게 변화하는가?
• RQ2 3+1차원 글라스마에서 압력 성분은 보정 불변성 2+1차원 경우와 얼마나 다를까?
• RQ3 로컬 정지 프레임에서의 에너지 밀도 프로파일은 유한한 핵 두께로 인해 빠르기 의존성을 보일 수 있는가?
• RQ4 3+1차원 실시간 격자 양자 게이지 이론 시뮬레이션에서 시뮬레이션 매개변수가 수치적 안정성에 어떻게 영향을 미치는가?
• RQ5 이 모델에서 유도된 빠르기 프로파일은 라이크(RHIC)와 대형 강입자 충돌기(LHC)의 실험 데이터를 정량적으로 일치시킬 수 있는가?

주요 결과

• 새로운 3+1차원 수치 기법은 핵 두께가 0에 수렴할 때 보정 불변성 글라스마를 제한적인 경우로 정확히 재현한다.
• 3+1차원 글라스마에서의 압력 이방성은 여전히 크며, 2+1차원 보정 불변성 경우와 유의미한 차이를 보이지 않는다.
• 로컬 정지 프레임에서의 에너지 밀도는 명시적 보정 불변성 파괴로 인해 강한 빠르기 의존성을 보인다.
• 에너지 밀도의 빠르기 프로파일 너비는 핵의 종방향 두께에 의해 제어되며, 실험 데이터와 양적 일치를 보인다.
• 수치 기법에서 이전에 알려지지 않은 불안정성이 발견되었으며, 이는 시간 도함수 처리 방식에 기인했고, 새로운 게이지 고유성, 제약 조건 유지 알고리즘으로 해결되었다.
• 새로운 기법은 CFL 조건 하에서 안정함이 증명되었으며, 큰 시간 단계에서도 가우스 제약 조건을 유지하여 장기적인 신뢰성을 확보한다.