[논문 리뷰] Efficient Evaluation of Cosmological Angular Statistics

Valentin Assassi, Marko Simonović외 1명|arXiv (Cornell University)|2017. 05. 14.

Cosmology and Gravitation Theories참고 문헌 28인용 수 12

한 줄 요약

이 논문은 파워법 함수 기반의 FFTlog 방법을 제안하여, 운동량 공간 통계를 거듭제곱 법칙 기저에 분해함으로써 구면 Bessel 함수의 해석적 통합을 가능하게 하여, 천체물리학적 각 지수 스펙트럼과 삼중 스펙트럼을 효율적으로 계산한다. 이 방법은 매우 진동하는 구면 Bessel 함수의 고비용 수치 통합을 피함으로써 계산 시간을 크게 단축시키며, 특히 각 지수 삼중 스펙트럼의 경우 빠른 성능 향상을 보이며, 다중극수 $\sigma$에 관계없이 속도 향상이 가능하다.

ABSTRACT

Angular statistics of cosmological observables are hard to compute. The main difficulty is due to the presence of highly-oscillatory Bessel functions which need to be integrated over. In this paper, we provide a simple and fast method to compute the angular power spectrum and bispectrum of any observable. The method is based on using an FFTlog algorithm to decompose the momentum-space statistics onto a basis of power-law functions. For each power law, the integrals over Bessel functions have a simple analytical solution. This allows us to efficiently evaluate these integrals, independently of the value of the multipole $\ell$. In particular, this method significantly speeds up the evaluation of the angular bispectrum compared to existing methods. To illustrate our algorithm, we compute the galaxy, lensing and CMB temperature angular power spectrum and bispectrum.

연구 동기 및 목표

• 매우 진동성이 큰 구면 Bessel 함수로 인해 천체물리학 관측량의 각 통계를 평가하는 데 발생하는 계산 블로킹 문제를 해결하기 위해.
• 은하, 렌징, CMB 이방이성과 같은 관측량의 빠르고 정확한 각 지수 스펙트럼과 삼중 스펙트럼 평가를 가능하게 하는 방법을 개발하기 위해.
• Bessel 적분이 해석적 해를 가지는 거듭제곱 법칙 기저로 운동량 공간 통계를 변환함으로써 수치 비용을 감소시키기 위해.
• 다양한 천체물리학 관측량과 상관 함수에 적용 가능한 일반적이고 재사용 가능한 프레임워크를 제공하기 위해.
• 기존의 수치 통합 대비 특히 각 지수 삼중 스펙트럼에서 뚜렷한 속도 향상을 보여주기 위해.

제안 방법

• 이 방법은 운동량 공간 통계(예: 스펙트럼 및 삼중 스펙트럼)를 $k$ 공간에서 거듭제곱 법칙 함수 기저에 투영하기 위해 FFTlog 알고리즘을 사용한다.
• 각 거듭제곱 법칙 성분에 대해, 구면 Bessel 함수 $j_\ell(k\chi)$에 대한 적분은 수치 적분을 피하고 해석적으로 평가된다.
• 전체 각 상관 함수는 이러한 해석적으로 풀 수 있는 성분들의 가중 합으로 계산되며, 문제는 부드러운 시간 적분으로 단순화된다.
• 운동량 공간에서 $n$-점 상관 함수의 분리 가능 형태를 활용하며, $f_i(k_i, z_i)$ 의 인과적 종속성을 가정한다.
• 고 $\ell$ 에서의 수치 안정성을 확보하기 위해 특수 함수 $\mathbf{I}_\ell(\nu, t)$ 의 재귀적 평가를 사용하며, 정밀도를 확보하기 위해 재초기화가 수행된다.
• 빠르게 수렴하는 급수를 가진 변환된 일반화된 하이퍼기하함수를 통해 $t^p \mathbf{I}_\ell(\nu, t)$ 를 포함하는 적분을 효율적으로 평가한다.

실험 결과

연구 질문

• RQ1 매우 진동성이 큰 Bessel 함수의 직접 수치 통합을 피함으로써 천체물리학 관측량의 각 지수 스펙트럼과 삼중 스펙트럼을 더 빠르게 계산할 수 있는가?
• RQ2 관측량의 운동량 공간 통계는 어떻게 분해되어야 하며, 이를 통해 구면 Bessel 함수 적분의 해석적 평가가 가능해지는가?
• RQ3 이 방법은 기존의 수치 통합 대비 어떤 성능 향상을 보이며, 특히 각 지수 삼중 스펙트럼의 경우 어떻게 되는가?
• RQ4 이 방법은 은하 수밀도, 렌징 포텐셜, CMB 온도 이방이성과 같은 다양한 관측량에 균일하게 적용될 수 있는가?
• RQ5 고 다중극수 $\ell$ 에서 정확성과 안정성을 확보하기 위한 수치 전략은 무엇이며, 속도를 희생시키지 않고 어떻게 실현되는가?

주요 결과

• 이 방법은 다중극수 $\ell$ 와 무관하게 각 지수 스펙트럼과 삼중 스펙트럼의 빠른 평가를 가능하게 하며, Bessel 함수의 $\ell$-의존성에 따라 계산 시간이 증가하지 않는다.
• 은하 토모그래피의 경우, 랩탑에서 $N_\ell = 200$ 다중극수를 사용해 30초 내로 각 지수 스펙트럼이 계산된다.
• CMB 렌징 스펙트럼은 $N_\ell = 200$ 일 때 5초 내로 계산되며, 저-$\ell$ 관측량에 대해 매우 효율적임을 보여준다.
• CMB 이방이성의 각 삼중 스펙트럼은 $10^5$ 개의 삼각형에 대해 $N_\nu = 70$, $N_\chi = 200$, $N_r = 80$ 조건에서 10분 내로 계산된다.
• 고 $\ell$ 에서도 빠른 수렴을 보장하는 변환된 하이퍼기하함수 형태를 사용함으로써 높은 정밀도를 확보하였으며, 재귀적 평가를 통한 오차 제어가 가능하다.
• 이 방법은 운동량 공간에서 분리 가능한 $n$-점 함수를 가진 모든 관측량에 일반적으로 적용 가능하며, 교차 상관 함수를 포함한 다양한 응용에 유용하다.