[논문 리뷰] Testing Higher--Order Lagrangian Perturbation Theory Against Numerical Simulations -- 2. Hierarchical Models

Adrian L. Melott, Thomas Buchert외 1명|ArXiv.org|1994. 04. 08.

Cosmology and Gravitation Theories참고 문헌 1인용 수 14

한 줄 요약

이 논문은 음수의 파워스펙트럼 지수를 가진 계층적 우주론 모델에 대해 N-body 시뮬레이션과 함께 두 번째 단계 라그랑주 변형 이론(Lagrangian perturbation theory, LPT)을 시험한다. 두 번째 단계 LPT는 제르드로비치 근사(Zel'dovich approximation, ZA)에 비해 역학적 정확도를 크게 향상시키며, 특히 중력 역학을 탐색할 때 두드러진다. 이는 후기 시점과 소규모 스케일에서도 안정되어 있어, 냉각, 난류, 혼합 어두운 물질 우주론 모델을 고성능으로 모델링하는 데 이상적이다.

ABSTRACT

We present results showing an improvement of the accuracy of perturbation theory as applied to cosmological structure formation for a useful range of scales. The Lagrangian theory of gravitational instability of Friedmann--Lema\^ıtre cosmogonies investigated and solved up to the third order in the series of papers by Buchert (1989, 1992, 1993), Buchert \& Ehlers (1993), Buchert (1994), Ehlers \& Buchert (1994), is compared with numerical simulations. In this paper we study the dynamics of hierarchical models as a second step. In the first step (Buchert, Melott and Weiß 1994) we analyzed the performance of the Lagrangian schemes for pancake models, i.e., models which initially have a truncated power spectrum. We here explore whether the results found for pancake models carry over to hierarchical models which are evolved deeply into the non--linear regime............We find that for spectra with negative power--index the second--order scheme performs considerably better than TZA in terms of statistics which probe the dynamics, and slightly better in terms of low--order statistics like the power--spectrum. In cases with much small--scale power the gain from the higher--order schemes is small, but still measurable. However, in contrast to the results found for pancake models, where the higher--order schemes get worse than TZA at late non--linear stages and on small scales, we here find that the second--order model is as robust as TZA, retaining the improvement at later stages and on smaller scales. In view of these

연구 동기 및 목표

• Zel'dovich 근사 이외의 고차수 라그랑주 변형 이론이 우주론적 구조 형성 모델링에서 얼마나 잘 작동하는지 평가하기 위해.
• 파워 스펙트럼이 잘린 패닉 모델(pancake models)에서 관찰된 개선 사항이 다양한 소규모 파워를 가진 계층적 모델으로까지 확장되는지 확인하기 위해.
• 크로스상관관계를 최대화하기 위해 초기 조건 필터링(필터 유형 및 스케일)을 최적화하기 위해.
• 특히 표준 어두운 물질 우주론 모델에 대해 두 번째 단계 LPT가 깊이 비선형 영역에서도 얼마나 견고한지 평가하기 위해.
• 대규모 구조 시뮬레이션과 은하 환경 모델링에서 두 번째 단계 LPT의 계산적 및 물리적 이점 설정하기 위해.

제안 방법

• Buchert(1989–1994)가 유도한 세 번째 단계 라그랑주 변형 이론 해를 활용하여, 일차 및 이차 단계 해법에 집중한다.
• 가우시안, 스타프-k 등 다양한 필터 유형을 사용하여 초기 밀도 및 속도 장을 가공하여 파워 스펙트럼을 잘라내고 분석 모델 성능을 최적화한다.
• 파워 법칙 지수 -3에서 +1 사이의 계층적 모델에 대해 N-body 시뮬레이션을 수행하며, 다양한 필터 스케일을 가진 초기 조건을 사용한다.
• 크로스상관 통계량(파워 스펙트럼 및 역학 진단 지표 포함)을 사용하여 분석적 LPT 예측과 N-body 결과를 비교한다.
• LPT와 N-body 결과 간의 크로스상관 계수를 최대화함으로써 최적의 필터 스케일을 도출한다.
• 후기 시점과 소규모 스케일에서 이차 해의 행동을 분석하고, ZA 및 고차수 해법과 비교한다.

실험 결과

연구 질문

• RQ1 음수의 파워스펙트럼 지수를 가진 계층적 모델에서 두 번째 단계 라그랑주 변형 이론이 Zel'dovich 근사보다 역학적 특성을 더 잘 포괄하는가?
• RQ2 특히 후기 시점과 소규모 스케일에서 두 번째 단계 LPT의 성능은 깊이 비선형 영역에서 어떻게 변화하는가?
• RQ3 초기 조건 필터링에서 다양한 필터 유형과 스케일에 대해 ZA에 비해 향상된 두 번째 단계 LPT 성능이 견고한가?
• RQ4 강한 비선형성이 나타나는 고소규모 파워를 가진 모델에서도 두 번째 단계 보정의 정확도를 유지할 수 있는가?
• RQ5 두 번째 단계 해법이 CDM, HDM, MDM 등의 표준 우주론 모델에서 얼마나 효과적인가?

주요 결과

• 두 번째 단계 라그랑주 변형 이론은 Zel'dovich 근사에 비해 특히 중력 역학을 탐색하는 통계량에서 역학적 정확도를 크게 향상시킨다.
• 두 번째 단계 보정의 개선 효과는 깊이 비선형 상태와 소규모 스케일에서도 견고하며, 패닉 모델에서 고차수 해법이 악화되는 것과는 정반대이다.
• 음수의 파워 법칙 지수를 가진 스펙트럼에서는 두 번째 단계 LPT가 ZA보다 역학 통계에서 훨씬 우수하며, 파워 스펙트럼 같은 저차수 통계에서는 약간의 개선에 그친다.
• 고소규모 파워를 가진 모델에서도 두 번째 단계 보정의 이점은 측정 가능하나, 크기는 작다.
• 두 번째 단계 잘라낸 라그랑주 모델은 Zel'dovich 근사와 유사한 CPU 비용을 가지며, 계산적으로 효율적이다.
• 이 방법은 냉각, 난류, 혼합 어두운 물질 우주론 모델의 대규모 구조 모델링에 특히 적합하며, 은하 군 질량 규모(약 10^13 M☉)까지 적용 가능하며, 대규모 시뮬레이션에서 노이즈 주파수에 가까운 잘라내기 스케일에 도달해도 여전히 효과적이다.