[논문 리뷰] Searching for gravitational wave burst in PTA data with piecewise linear functions

Heling Deng, B. Bécsy외 3명|arXiv (Cornell University)|2023. 06. 29.

Pulsars and Gravitational Waves Research인용 수 4

한 줄 요약

이 논문은 펄서 타이밍 어레이(PTA) 데이터에서 알려지지 않은 버스트 웨이브폼을 조각별 선형 함수로 모델링함으로써, 전체 베이지안 프레임워크를 제안한다. 이는 웨이브폼 파rameter에 대한 분석적 마진화를 가능하게 하며, 탐지 문제를 천체 위치와 신호 강도의 세 가지 핵심 파rameter로 축소시킨다. 신호가 존재할 경우 효율적인 추론과 정확한 웨이브폼 재구성 가능하며, 신호 대 잡음비(SNR) ≈ 14.7 및 ≈ 6.5인 시뮬레이션 데이터에서 이를 입증하였다.

ABSTRACT

Transient gravitational waves (aka gravitational wave bursts) within the nanohertz frequency band could be generated by a variety of astrophysical phenomena such as the encounter of supermassive black holes, the kinks or cusps in cosmic strings, or other as-yet-unknown physical processes. Radio-pulses emitted from millisecond pulsars could be perturbed by passing gravitational waves, hence the correlation of the perturbations in a pulsar timing array can be used to detect and characterize burst signals with a duration of $\mathcal{O}(1 ext{-}10)$ years. We propose a fully Bayesian framework for the analysis of the pulsar timing array data, where the burst waveform is generically modeled by piecewise straight lines, and the waveform parameters in the likelihood can be integrated out analytically. As a result, with merely three parameters (in addition to those describing the pulsars' intrinsic and background noise), one is able to efficiently search for the existence and the sky location of {a burst signal}. If a signal is present, the posterior of the waveform can be found without further Bayesian inference. We demonstrate this model by analyzing simulated data sets containing a stochastic gravitational wave background {and a burst signal generated by the parabolic encounter of two supermassive black holes.

연구 동기 및 목표

• 알려지지 않은 웨이브폼을 가진 펄서 타이밍 어레이(PTA) 데이터에서 일시적인 중력파 버스트를 효율적으로 탐지하기 위한 베이지안 방법을 개발한다.
• 시간 도메인에서 웨이브폼을 두 개의 조각별 선형 함수로 모델링함으로써, 웨이브폼 파rameter에 대한 분석적 통합을 가능하게 한다.
• 웨이브폼에 대해 마진화함으로써 자유도 파rameter 수를 천체 위치와 신호 강도의 단지 세 개로 줄여 계산 효율성을 향상시킨다.
• 신호가 탐지된 경우 추가적인 베이지안 추론 없이도 직접적으로 버스트 웨이브폼을 재구성할 수 있도록 한다.
• 스토하스틱 중력파 배경과 초거대 블랙홀의 타원형 접근에서 발생하는 버스트를 포함한 시뮬레이션된 PTA 데이터를 대상으로 방법을 검증한다.

제안 방법

• 버스트 웨이브폼은 시간 도메인에서 각각 두 개의 조각별 선형 함수로 모델링되며, 각각 + 및 × 편광에 대응한다.
• 웨이브폼 파ram터는 공분산이 단일 하이퍼파rameter q로 제어되는 가우시안 사전분포를 부여받으며, 이는 가능하게 하여 우도에 대한 분석적 마진화를 실현한다.
• 마진화된 우도는 q(신호 강도), θ 및 φ(천체 위치)를 포함하여 총 세 개의 파rameter에만 의존하며, 내재적 및 공통 잡음 파rameter도 포함된다.
• 마르코프 체인 몬테카를로(MCMC) 샘플링을 통해 q, θ, φ의 사후분포를 탐색하며, 전체 웨이브폼은 사후 샘플들로부터 재구성된다.
• 실제 PTA 설정(스토하스틱 중력파 배경 포함)을 반영한 세 개의 시뮬레이션 데이터 세트를 대상으로 방법을 테스트하였다: 두 개는 버스트가 포함되어 있으며, 신호 대 잡음비(SNR)는 각각 ≈ 14.7 및 ≈ 6.5이고, 하나는 버스트가 없는 경우이다.
• 베이즈 요인을 통한 모델 비교를 통해 버스트 신호 존재 여부에 대한 증거를 평가한다.

실험 결과

연구 질문

• RQ1 알려지지 않은 웨이브폼을 가진 PTA 데이터에서 조각별 선형 모델링을 사용하여 베이지안 프레임워크가 중력파 버스트를 효율적으로 탐지할 수 있는가?
• RQ2 웨이브폼 파arameter를 분석적으로 마진화할 수 있으며, 이로 인해 추론 문제의 파arameter 수가 단지 세 개로 줄어드는가?
• RQ3 중간에서 높은 신호 대 잡음비를 가진 버스트 신호의 천체 위치와 웨이브폼을 얼마나 잘 복원할 수 있는가?
• RQ4 신호가 약할 경우 배경 잡음과의 구별에서 방법의 성능은 어떠한가?
• RQ5 더 나은 웨이브폼 해상도를 확보하기 위해 그리드 포인트의 수와 간격을 적응적으로 선택할 수 있도록 모델을 확장할 수 있는가?

주요 결과

• 신호 대 잡음비(SNR) ≈ 14.7인 강한 버스트 신호의 경우, 사후분포에서 천체 위치 파arameter θ 및 φ가 진짜 값 근처에 피크를 이룬다. 이는 소음 전용 모델 대비 신호 존재에 대한 강력한 지지를 나타낸다.
• MCMC 샘플들로부터 재구성된 웨이브폼은 진짜 주입된 웨이브폼과 매우 유사하여 정확한 신호 복원 능력을 입증한다.
• 신호 대 잡음비(SNR) ≈ 6.5인 더 약한 신호의 경우에도 천체 위치 사후분포는 여전히 진짜 값 근처에 피크를 이룬다. 이는 웨이브폼의 불확실성에도 불구하고 탐지 가능성은 유지됨을 시사한다.
• 버스트가 없는 경우, 버스트 모델과 소음 전용 모델 간의 베이즈 요인은 ≈ 1에 가까워, 모델이 소음 전용 데이터와 일관됨을 확인한다.
• 스토하스틱 중력파 배경의 파aram터는 모든 시뮬레이션 케이스에서 진짜 값과 호환되는 것을 성공적으로 복원하였다.
• 분석적 마진화 덕분에 계산이 효율적이며, 반복적인 웨이브폼 피팅이 필요 없이 빠른 추론이 가능하다.