[논문 리뷰] Hellings and Downs correlation of an arbitrary set of pulsars

B. Allen, Joseph D. Romano|arXiv (Cornell University)|2022. 08. 15.

Pulsars and Gravitational Waves Research참고 문헌 41인용 수 1

한 줄 요약

이 논문은 펄서 팀 링 어레이(PTA)에서 헬링스와 다운스(Hellings and Downs, HD) 상관관계에 대한 최적 추정기(_estimator)를 개발하며, 각도 구간화, 우주적 변동성, 펄서 쌍 간 상관관계를 고려한다. 구간화된 HD 상관관계 추정치의 분산과 공분산을 유도하여, 많은 펄서가 존재할 경우 이들이 우주적 변동성으로 수렴함을 보이고, 관측된 HD 곡선에서의 이격을 검증하기 위해 χ² 적합도 통계량을 제안함으로써, 노이즈가 있는 데이터에서 중력파 배경을 정밀하게 탐지할 수 있게 한다.

ABSTRACT

Pulsar timing arrays (PTAs) detect gravitational waves (GWs) via the correlations they induce in the arrival times of pulses from different pulsars. We assume that the GWs are described by a Gaussian ensemble, which models the confusion noise produced by expected PTA sources. The mean correlation h2μu(γ) as a function of the angle γ between the directions to two pulsars was predicted by Hellings and Downs in 1983. The variance σtot2(γ) in this correlation was recently calculated [B. Allen, Variance of the Hellings-Downs correlation, Phys. Rev. D 107, 043018 (2023)PRVDAQ2470-001010.1103/PhysRevD.107.043018] for a single noise-free pulsar pair at angle γ, which shows that after averaging over many pairs, the variance reduces to an intrinsic cosmic variance σcos2(γ). Here, we extend this to an arbitrary set of pulsars at specific sky locations, with pulsar pairs binned by γ. We derive the linear combination of pulsar-pair correlations which is the optimal estimator of the Hellings and Downs correlation for each bin, illustrating our methods with plots of the expected range of variation away from the Hellings and Downs curve, for the sets of pulsars monitored by three active PTA collaborations. We compute the variance of and the covariance between these binned estimates, and show that these reduce to the cosmic variance and covariance s(γ,γ′) respectively, in the many-pulsar limit. The likely fluctuations away from the Hellings and Downs curve μu(γ) are strongly correlated/anticorrelated in the three angular regions where μu(γ) is successively positive, negative, and positive. We also construct the optimal estimator of the squared strain h2 from pulsar-pair correlation data. Remarkably, when there are very many pulsar pairs, this determines h2 with arbitrary precision because (in contrast to LIGO-like GW detectors) PTAs probe an infinite set of GW modes. To assess if observed deviations away from the Hellings and Downs curve are consistent with predictions, we propose and characterize several χ2 goodness-of-fit statistics. While our main focus is ideal noise-free data, we also show how pulsar noise and measurement noise can be included. Our methods can also be applied to future PTAs, where the improved telescopes will provide larger pulsar populations and higher-precision timing.

연구 동기 및 목표

• 지정된 천구 위치를 가진 임의의 펄서 집합에 대해 HD 상관관계를 최적의 선형 추정기로 추정하는 것.
• 구간화된 HD 상관관계 추정치의 분산과 공분산을 정량화하여, 펄서 분산과 우주적 변동성 간의 차이를 명확히 하는 것.
• 특히 HD 곡선이 양, 음, 양의 순서로 변화하는 영역에서, 우주적 변동성으로 인한 예상되는 변동성의 특성을 기술하는 것.
• 다수의 펄서 쌍이 존재할 경우 임의의 정밀도를 달성할 수 있는, 펄서 쌍 상관관계로부터 제곱 변형도 h²의 강력한 추정기 구축.
• 관측된 HD 곡선에서의 이격이 통계적으로 예측과 일치하는지 평가하기 위해 χ² 적합도 통계량을 제안하고 검증하는 것.

제안 방법

• 각각의 각도 구간에서 HD 상관관계를 추정하기 위해 펄서 쌍 상관관계의 최적 선형 조합을 유도하여 분산을 최소화하는 것.
• 총 분산 σ²_tot(γ)를 펄서 분산과 우주적 변동성 σ²_cos(γ)의 합으로 계산하며, 후자는 스토하스틱 중력파 배경의 단계 변동성에서 기인함.
• 구면 조화 함수 분해와 클렙슈-고르단 계수를 사용하여 서로 다른 각도에서의 구간화된 HD 추정치 간의 공분산을 계산하는 것.
• 많은 수의 펄서가 존재할 경우, 추정기의 분산이 본질적인 우주적 변동성 s(γ, γ′)로 수렴함을 보이는 것.
• 역 공분산 행렬 기반의 χ² 통계량을 도입하여 관측된 상관관계가 HD 곡선과 얼마나 잘 맞는지 적합도를 테스트하는 것.
• 펄서 팀 링 노이즈와 측정 노이즈를 포함한 프레임워크로 확장하여, 현재 및 향후 PTA에 대한 실제 적용 가능성을 보여주는 것.

실험 결과

연구 질문

• RQ1 펄서 쌍이 각도 간격에 따라 구간화될 경우, HD 상관관계에 대한 최적 선형 추정기는 무엇인가?
• RQ2 구간화된 HD 상관관계 추정치의 분산과 공분산은 각도 분포와 펄서 수에 따라 어떻게 달라지는가?
• RQ3 우주적 변동성이 HD 상관관계 측정 정밀도를 제한하는 데 어떤 역할을 하는가? 그리고 펄서 분산과 비교해 보면 어떻게 되는가?
• RQ4 전체 공분산 구조를 반영한 χ² 통계량은 실제 PTA 데이터에서 HD 곡선에서의 이격을 신뢰성 있게 탐지할 수 있는가?
• RQ5 점점 증가하는 펄서 쌍의 수에 따라 제곱 변형도 h² 추정기의 정밀도는 어떻게 향상되는가? 그리고 많은 펄서 수의 극한에서 임의의 정확도에 도달하는가?

주요 결과

• 각각의 각도 구간에서 HD 상관관계에 대한 최적 추정기는 분산을 최소화하며, 많은 펄서가 존재할 경우 우주적 변동성 s(γ, γ′)로 수렴함.
• 모든 펄서 쌍의 방향에 대해 평균을 내면, 구간화된 HD 상관관계 추정기의 분산은 본질적인 우주적 변동성 σ²_cos(γ)로 줄어듦.
• HD 곡선이 양, 음, 양의 순서로 변화하는 세 영역에서, 곡선 주변의 변동성은 강하게 상관되거나 반대상관됨.
• 다수의 펄서 쌍이 존재할 경우, 중력파 모드의 수가 무한히 많아지므로, 제곱 변형도 h² 추정기는 임의의 정밀도에 도달함.
• 제안된 χ² 통계량은 펄서 쌍 간 상관관계를 고려하여, 노이즈가 존재하는 상황에서도 HD 곡선과의 일치성을 신뢰성 있게 평가할 수 있는 강력한 방법을 제공함.
• 이 프레임워크는 펄서 팀 링 노이즈와 측정 오차를 포함한 실제 데이터로 확장 가능하며, 더 큰 펄서 집단을 가진 현재 및 향후 PTA에 적용 가능함.