[논문 리뷰] The art of simulating the early Universe -- Part I

Daniel G. Figueroa, Adrien Florio외 2명|arXiv (Cornell University)|2020. 06. 26.

Cosmology and Gravitation Theories참고 문헌 275인용 수 13

한 줄 요약

이 논문은 고정밀 격자 기법을 사용하여 팽창하는 초기 우주의 스칼라 및 게이지 장 역학을 시뮬레이션하기 위한 종합적인 프레임워크를 제시한다. O(δt²)에서 O(δt¹⁰)까지의 심플렉틱 적분기법을 도입하여 기하학적 제약 조건을 기계 정밀도 수준에서 유지함으로써, 비선형 역학의 3+1차원 정확한 시뮬레이션, 특히 자기 일관된 팽창을 포함한 시뮬레이션을 가능하게 한다. 주요 기여는 후기 인플레이션 역학, 예를 들어 SU(2)×U(1) 프리헤이팅을 연구하는 데 사용할 수 있는 공개된 C++ 패키지인 CosmoLattice의 구축이다.

ABSTRACT

We present a comprehensive discussion on lattice techniques for the simulation of scalar and gauge field dynamics in an expanding universe. After reviewing the continuum formulation of scalar and gauge field interactions in Minkowski and FLRW backgrounds, we introduce basic tools for the discretization of field theories, including lattice gauge invariant techniques. Following, we discuss and classify numerical algorithms, ranging from methods of $O(dt^2)$ accuracy like $staggered~leapfrog$ and $Verlet$ integration, to $Runge-Kutta$ methods up to $O(dt^4)$ accuracy, and the $Yoshida$ and $Gauss-Legendre$ higher-order integrators, accurate up to $O(dt^{10})$. We adapt these methods for their use in classical lattice simulations of the non-linear dynamics of scalar and gauge fields in an expanding grid in $3+1$ dimensions, including the case of `self-consistent' expansion sourced by the volume average of the fields' energy and pressure densities. We present lattice formulations of canonical cases of: $i)$ Interacting scalar fields, $ii)$ Abelian $U(1)$ gauge theories, and $iii)$ Non-Abelian $SU(2)$ gauge theories. In all three cases we provide symplectic integrators, with accuracy ranging from $O(dt^2)$ up to $O(dt^{10})$. For each algorithm we provide the form of relevant observables, such as energy density components, field spectra and the Hubble constraint. Remarkably, all our algorithms for gauge theories respect the Gauss constraint to machine precision, including when `self-consistent' expansion is considered. As a numerical example we analyze the post-inflationary dynamics of an oscillating inflaton charged under $SU(2) imes U(1)$. The present manuscript is meant as part of the theoretical basis for $CosmoLattice$, a modern C++ MPI-based package for simulating the non-linear dynamics of scalar-gauge field theories in an expanding universe, publicly available at www.cosmolattice.net

연구 동기 및 목표

• 팽창하는 Friedmann-Lemaître-Robertson-Walker (FLRW) 우주에서 스칼라 및 게이지 장 역학을 시뮬레이션하기 위한 수치적 안정성과 고정밀 격자 통합 기법을 개발하는 것.
• 특히 장의 에너지와 압력에 의해 유도되는 자기 일관된 팽창 상황에서도 게이지 불변성과 제약 조건 유지 문제를 해결하는 것.
• 2차에서 10차 정확도에 이르는 수치 알고리즘의 체계적 분류 및 구현을 통해 고전적 격자 장 이론에서의 응용을 가능하게 하는 것.
• CosmoLattice 패키지의 이론적 및 계산적 기반을 마련하여, 예를 들어 프리헤이팅과 같은 비선형 초기 우주 현상의 고정밀 시뮬레이션을 가능하게 하는 것.

제안 방법

• 자연 단위와 평탄한 FLRW 계량을 채택하여, 팽창 우주에서 스칼라 및 게이지 장의 연속장 방정식을 유도하며, U(1) 및 SU(2) 상호작용을 포함한다.
• 장 방정식을 공간 격자에 이산화하고, 이산 푸리에 변환을 사용하여 격자 미분 연산자와 운동량 표현을 구현함으로써, 스태그러드 및 윌슨 유형의 형식을 통해 게이지 불변성을 확보한다.
• 이계에서 10차까지의 수치 적분기법 계층을 구현: 이차 스태그러드 리프플로그 및 버렛 방법에서 시작하여, 4차 명시적 룬게-쿠타 방법과 10차까지의 암시적 가우스-레지앙드 및 요시다 방법에 이르기까지.
• 게이지 이론의 경우, 기하학적 제약 조건을 수치적으로 강제하기 위해 기하학적 제약 조건을 기계 정밀도 수준에서 유지하는 심플렉틱 적분기법을 구성한다. 이는 자기 일관된 시뮬레이션에서도 성립한다.
• 에너지 밀도 성분, 장 스펙트럼, 허블 제약 조건 등의 관측 가능량을 격자 정규화 기대값을 통해 계산하는 프레임워크를 포함한다.
• 실제 적용 사례를 통해 방법을 검증: 전하를 가진 인플라톤이 있는 SU(2)×U(1) 프리헤이팅을 시뮬레이션하여, 다양한 적분기종에서 정확도와 제약 조건 유지가 일관되게 확보됨을 입증한다.

실험 결과

연구 질문

• RQ1 팽창 우주에서 비선형 스칼라 및 게이지 장 역학을 최소한의 수치 오차로 시뮬레이션하기 위해 고차 심플렉틱 적분기법을 어떻게 설계하고 구현할 수 있는가?
• RQ2 특히 장 자체의 에너지에 의해 유도되는 팽창 상황에서, 비아벨 게이지 이론의 격자 시뮬레이션에서 기하학적 제약 조건 위반을 어느 정도 억제할 수 있는가?
• RQ3 3+1차원 격자 장 이론 시뮬레이션에서 명시적 적분기법(예: 룬게-쿠타)과 암시적 적분기법(예: 가우스-레지앙드) 간의 성능 및 정확도 간 상호 교환 관계는 어떠한가?
• RQ4 에너지 보존 및 게이지 불변성과 같은 핵심 물리적 제약 조건을 유지하는 방식으로 상호작용하는 스칼라 및 게이지 장의 격자 형식을 체계적으로 구성할 수 있는가?
• RQ5 2차에서 10차까지의 다양한 적분기종(O(δt²)에서 O(δt¹⁰))이 후기 인플레이션 역학에서 장 스펙트럼과 에너지 밀도의 진화에 어떤 영향을 미치는가?

주요 결과

• 저자들은 3+1차원 격자 시뮬레이션에서 스칼라 및 게이지 장에 대해 O(δt¹⁰) 정확도의 심플렉틱 적분기법을 성공적으로 구현하였으며, 시간 단계 가중치에 대한 명시적 공식을 표 1~3에 제시하였다.
• 모든 게이지 이론 적분기법은 자기 일관된 팽창 시뮬레이션에서도 기하학적 제약 조건을 기계 정밀도 수준에서 유지한다. 이는 물리적 일관성을 유지하는 데 있어 핵심적인 성취이다.
• 가우스-레지앙드 및 요시다 방법은 기하학적 제약 조건을 기계 정밀도 수준에서 유지하는 안정적이고 대칭적인 시간 스텝 스킴을 제공하며, s = 2에서 5 단계까지의 부처 표기법을 통해 검증되었다.
• 전하를 가진 인플라톤이 있는 SU(2)×U(1) 프리헤이팅 시뮬레이션은 다양한 적분기종에서 에너지 및 스펙트럼 진화가 일관되게 유지됨을 보여주며, 수치적 안정성과 정확도를 확인한다.
• 이 프레임워크를 기반으로 한 CosmoLattice 패키지의 구현은 공개된 고정밀도 초기 우주 역학 시뮬레이션을 가능하게 하며, 비아벨 게이지 장과 자기 일관된 팽창을 포함한다.
• 장 스펙트럼과 에너지 밀도 성분을 계산하는 데 높은 정확도를 달성하였으며, 다양한 적분기종에서의 수렴 테스트를 통해 오차 제어가 검증되었다.