[논문 리뷰] Perturbation-theory informed integrators for cosmological simulations

Florian List, Oliver Hahn|arXiv (Cornell University)|2023. 01. 23.

Numerical methods for differential equations인용 수 6

한 줄 요약

이 논문은 라그랑주 편미분 이론(LPT)에 기반한 시간단계 적분기법을 도입하여, 입자 궤적을 LPT로 예측한 것과 일치시키는 우주론적 시뮬레이션을 제안한다. 이는 더 적은 시간단계로도 정확도를 크게 향상시킨다. 제안된 PowerFrog 적분기법은 FastPM에 비해 2LPT 방향으로는 2.5배 이상, 3LPT 방향으로는 5배 이상의 이격도 필드 오차를 감소시키며, O(10–100)개의 시간단계로 2D 및 3D 시뮬레이션에서 표준 방법보다 뛰어난 성능을 보인다.

ABSTRACT

Large-scale cosmological simulations are an indispensable tool for modern cosmology. To enable model-space exploration, fast and accurate predictions are critical. In this paper, we show that the performance of such simulations can be further improved with time-stepping schemes that use input from cosmological perturbation theory. Specifically, we introduce a class of time-stepping schemes derived by matching the particle trajectories in a single leapfrog/Verlet drift-kick-drift step to those predicted by Lagrangian perturbation theory (LPT). As a corollary, these schemes exactly yield the analytic Zel'dovich solution in 1D in the pre-shell-crossing regime (i.e. before particle trajectories cross). One representative of this class is the popular FastPM scheme by Feng et al. 2016, which we take as our baseline. We then construct more powerful LPT-inspired integrators and show that they outperform FastPM and standard integrators in fast simulations in two and three dimensions with $\mathcal{O}(1 - 100)$ timesteps, requiring less steps to accurately reproduce the power spectrum and bispectrum of the density field. Furthermore, we demonstrate analytically and numerically that, for any integrator, convergence is limited in the post-shell-crossing regime (to order 3/2 for planar wave collapse), owing to the lacking regularity of the acceleration field, which makes the use of high-order integrators in this regime futile. Also, we study the impact of the timestep spacing and of a decaying mode present in the initial conditions. Importantly, we find that symplecticity of the integrator plays a minor role for fast approximate simulations with a small number of timesteps.

연구 동기 및 목표

• 라그랑주 편미분 이론(LPT)에서 도출된 해석적 해를 시간단계 기반 시뮬레이션에 통합하여, 우주론적 N체 시뮬레이션의 정확도와 효율성을 향상시키는 것.
• 쉘-크로싱 이전 영역에서 높은 정확도를 유지하면서도 시간단계 수를 최소화하여 계산 비용을 절감하는 것.
• 적은 시간단계로 빠른 근사 시뮬레이션을 수행할 때, 심플렉틱성(보존성)이 필수적인지 여부를 조사하는 것.
• 단지 제1차 라그랑주 편미분 이론(1LPT)이 아닌, 고차수 LPT 해(예: 2LPT, 3LPT)로 수렴하는 적분기법을 개발하는 것.
• 초기 조건에서 시간단계 간격과 감쇠 모드의 영향이 시뮬레이션 정확도에 미치는 영향을 평가하는 것.

제안 방법

• 시간단계 기반 적분기법을 유도하기 위해, 단일 리프플로그/베르레트 드리프트-킥-드리프트 단계가 라그랑주 편미분 이론(LPT)에 의해 예측된 궤적과 일치하도록 설정한다.
• 1차원에서 셸-크로싱 이전에 제1차 라그랑주 편미분 이론(1LPT)에 일치하는 해를 갖는 적분기법의 집합인 Π-적분기법을 구성한다.
• FastPM와 같은 기존 기법보다 고차수 LPT 항(≥2LPT)을 더 잘 포착할 수 있도록 설계된 대표적 기법으로서, PowerFrog 적분기법을 도입한다.
• 초기 시뮬레이션 단계에서 궤적 정확도를 향상시키기 위해 LPT 기반 속도 업데이트를 통합한 수정된 킥 단계를 사용한다.
• ΛCDM 우주론 모델을 기반으로 한 1D, 2D, 3D 시뮬레이션을 수행하여, 기준 해와의 파워 스펙트럼 및 비스펙트럼 비교를 통해 적분기법의 성능을 평가한다.
• 특히 셸-크로싱 이후 영역에서 가속도 필드의 규칙성 상실로 인해 발생하는 에너지 보존 및 수렴 특성에 대해 분석한다.

실험 결과

연구 질문

• RQ1 LPT에 기반한 시간단계 기반 적분기법이 동일한 시간단계 수로도 기존 표준 적분기법보다 더 높은 정확도를 달성할 수 있는가?
• RQ2 LPT 기반 적분기법은 2D 및 3D 시뮬레이션에서 파워 스펙트럼과 비스펙트럼을 FastPM 및 표준 리프플로그 기법과 비교해 얼마나 잘 재현하는가?
• RQ3 O(10–100)개의 시간단계로 빠른 근사 시뮬레이션을 수행할 때, 심플렉틱성은 에너지 보존에 크게 기여하는가?
• RQ4 PowerFrog와 같은 기법이 제1차 라그랑주 편미분 이론(1LPT)이 아닌 고차수 LPT 해(예: 2LPT, 3LPT)를 얼마나 잘 포착할 수 있는가?
• RQ5 셸-크로싱 이후 영역에서 적분기법의 수렴 차수는 무엇에 의해 제한되며, 고차수 기법이 이곳에서 효과가 없는 이유는 무엇인가?

주요 결과

• 2D 시뮬레이션에서 PowerFrog 적분기법은 FastPM의 DKD 기법 대비 2LPT 해로 햖थ한 평균 이격도 필드 오차를 2.5배 이상 감소시키고, 3LPT 해로 향한 오차는 5배 이상 감소시켰다.
• Quijote 및 Camels 우주론적 박스를 사용한 3D 시뮬레이션에서, Π-적분기법(LPTFrog, TsafPM, PowerFrog)은 FastPM 및 표준 리프플로그 기법보다 O(10–100)개의 시간단계에서 교차 파워 스펙트럼과 등변 비스펙트럼을 더 잘 재현하였다.
• 제1차 라그랑주 편미분 이론(1LPT)에 일치하는 적분기법은 셸-크로싱 이전 1차원에서 정확한 해를 달성하여 이론적 기대와 일치함을 확인하였다.
• 셸-크로싱 이후 영역에서 수렴 차수는 가속도 필드의 비규칙성으로 인해 최대 3/2차까지 제한되며, 이로 인해 고차수 적분기법은 효과를 발휘하지 못한다.
• 빠른 근사 시뮬레이션에서 에너지 보존에 대해 심플렉틱성은 미미한 영향을 미치며, 오히려 시간단계의 선택 및 샘플링 방식이 에너지 오차 제어에 더 중요한 역할을 한다.
• PowerFrog 적분기법의 해는 3LPT 해와 매우 유사하여, a=0에서 단일 PowerFrog 단계를 수행함으로써 LPT 기반 초기 조건을 직접 설정하는 것이 가능할 수 있으며, 이로써 LPT 기반 초기 조건의 필요성을 제거할 수 있다.