[논문 리뷰] A direct method to compute the galaxy count angular correlation function including redshift-space distortions

J.E. Campagne, S. Plaszczynski외 1명|arXiv (Cornell University)|2017. 03. 08.

Galaxies: Formation, Evolution, Phenomena참고 문헌 28인용 수 4

한 줄 요약

이 논문은 새로운 구면 조화 함수 덧셈 정리(Addition Theorem)를 통해 중간 단계의 파wer 스펙트럼 계산을 피하고, 적색편이 공간에서 은하 수의 각도 상관 함수 ξ(θ,z₁,z₂)를 직접적으로 수치적으로 효율적으로 계산하는 방법을 제시한다. 이 방법은 백분율 이내의 정확도를 달성하고 기존 코드인 CLASSgal 대비 10-20배 빠른 속도를 보이며, 기준 모델에 대한 가정 없이 대규모 구조 데이터를 빠르게 분석할 수 있다.

ABSTRACT

In the near future, cosmology will enter the wide and deep galaxy survey area allowing high-precision studies of the large scale structure of the universe in three dimensions. To test cosmological models and determine their parameters accurately, it is natural to confront data with exact theoretical expectations expressed in the observational parameter space (angles and redshift). The data-driven galaxy number count fluctuations on redshift shells, can be used to build correlation functions $C(θ; z_1, z_2)$ on and between shells which can probe the baryonic acoustic oscillations, the distance-redshift distortions as well as gravitational lensing and other relativistic effects. Transforming the model to the data space usually requires the computation of the angular power spectrum $C_\ell(z_1, z_2)$ but this appears as an artificial and inefficient step plagued by apodization issues. In this article we show that it is not necessary and present a compact expression for $C(θ; z_1, z_2)$ that includes directly the leading density and redshift space distortions terms from the full linear theory. It can be evaluated using a fast integration method based on Clenshaw-Curtis quadrature and Chebyshev polynomial series. This new method to compute the correlation functions without any Limber approximation, allows us to produce and discuss maps of the correlation function directly in the observable space and is a significant step towards disentangling the data from the tested models.

연구 동기 및 목표

• 중간 단계의 파워 스펙트럼 계산에 의존하지 않고 적색편이 공간에서 ξ(θ,z₁,z₂)를 수치적으로 효율적으로 계산하는 직접적 방법을 개발한다.
• 구면 조화 함수 전개와 덧셈 정리를 사용하여 ξ(θ,z₁,z₂)의 폐쇄형 해석적 표현을 제공함으로써 유도 및 구현을 단순화한다.
• Limber 근사 없이도 적색편이 공간 변형(RSD), 밀도, 도플러, 렌즈 효과를 포함한 상관 함수의 정확한 계산을 가능하게 한다.
• CLASSgal과 같은 기존 코드보다 빠르고 메모리 효율성이 뛰어난 천체물리학적 파arameter 추정을 위한 공개 소프트웨어 도구(Angpow)를 개발한다.
• 관측 가능 공간(각도와 적색편이)에서 직접 작동함으로써 ΛCDM-GR 철학에 대한 데이터 주도적이고 천체물리학 모델에 종속되지 않는 검증을 가능하게 한다.

제안 방법

• 다중극자에 대한 직접 합산을 가능하게 하기 위해 구면 조화 함수 전개와 덧셈 정리를 사용하여 ξ(θ,z₁,z₂)의 폐쇄형 표현을 유도한다.
• 유한한 ℓ_max 절단에 의해 발생하는 진동을 피하기 위해 레지드 폴리노미얼의 무한합으로 상관 함수를 표현한다.
• Cℓ(z₁,z₂)를 직접 계산하지 않고 3C-알고리즘을 활용하여 각도 파워 스펙트럼 Cℓ(z₁,z₂)로부터 상관 함수를 효율적으로 재구성한다.
• 조화 공간 표현에서 추가 도함수와 연산자를 통해 RSD, 도플러 이동, 렌즈 효과와 같은 상대론적 효과를 체계적으로 통합한다.
• apodization 없이도 ℓ_max = 3500인 기준 천체물리학 모델을 사용해 CLASSgal 코드와의 정합성을 검증한다.
• 다중 코어 시스템에서 저메모리 및 고속 계산을 최적화한 공개 소프트웨어 Angpow에 이 방법을 구현한다.

실험 결과

연구 질문

• RQ1 중간 단계의 Cℓ(z₁,z₂) 파워 스펙트럼을 계산하지 않고도 ξ(θ,z₁,z₂)에 대한 직접적이고 폐쇄형 표현을 도출할 수 있는가?
• RQ2 적색편이 공간 변형, 도플러 이동, 렌즈 효과와 같은 상대론적 효과를 상관 함수 형식에 어떻게 체계적으로 통합할 수 있는가?
• RQ3 새로운 방법이 CLASSgal과 같은 기존 접근 방식보다 더 높은 수치적 효율성과 정확도를 달성할 수 있는가?
• RQ4 관측 가능 공간에서 직접 작동함으로써 데이터 분석에서 기준 천체물리학 모델이 필요 없어지는 정도는 어느 정도인가?
• RQ5 적색편이 선택 함수의 선택(디랙 δ 함수 vs. 가우시안)이 관측된 상관 함수의 구조, 특히 BAO 와이브스에 어떤 영향을 미치는가?

주요 결과

• 새로운 방법은 CLASSgal 코드와 비교해 백분율 이내의 일치도를 보이며, 다양한 적색편이와 각도 스케일에서 정확도를 입증한다.
• 구현된 코드는 CLASSgal 대비 10-20배 더 빠르게 동작하여, 16개 스레드에서 18초 만에 실행되며 메모리 사용량은 1GB에 불과하지만, CLASSgal은 8개 스레드에서 3분이 걸리고 100GB의 메모리를 소비한다.
• Cℓ가 고다중수에서 급격히 감쇠함에 따라 고ℓ 효과가 자연스럽게 억제되므로, apodization이나 윈도우 함수가 필요 없어진다.
• RSD와 렌즈 효과를 포함함으로써 중심 상관 피크의 수축과 θ ≈ 0, Δz ≈ ±0.05에서의 음의 상관관계가 정확히 재현된다.
• 가우시안 적색편이 선택 함수는 θ ≈ 0에서 BAO 와이브스를 완화시키지만, 더 높은 각도와 작은 Δz에서는 유지되며 이는 이론적 기대와 일치한다.
• Angpow 소프트웨어는 https://gitlab.in2p3.fr/campagne/AngPow 에 공개되어 있으며, 대규모 구조 연구에서 빠르고 천체물리학 모델에 종속되지 않는 매개변수 추정을 가능하게 한다.