[논문 리뷰] Marginalizing the likelihood function for modeled gravitational wave searches

Sanjeev Dhurandhar, B. Krishnan외 1명|arXiv (Cornell University)|2017. 07. 25.

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

한 줄 요약

이 논문은 4개의 외재적 중력파 매개변수—진폭, 초기 위상, 편광 각도, 기울기—를 회전군 행렬 원소를 사용하여 복소 진폭으로 재구성함으로써 이들의 우도 함수에 대한 근사적 적분을 수행하는 분석적 방법을 제시한다. 주요 기여는 고신호대비비 비율 영역에서 B-통계량에 대한 폐쇄형 근사식을 도출한 것으로, 이는 일시적 이진별 융합과 중성자별에서 발생하는 연속 중력파에 대한 효율적인 베이지안 탐지 가능하게 한다.

ABSTRACT

Matched filtering is a commonly used technique in gravitational wave searches for signals from compact binary systems and from rapidly rotating neutron stars. A common issue in these searches is dealing with four extrinsic parameters which do not affect the phase evolution of the system: the overall amplitude, initial phase, and two angles determining the overall orientation of the system. The F-statistic maximizes the likelihood function analytically over these parameters, while the B-statistic marginalizes over them. The B-statistic, while potentially more powerful and capable of incorporating astrophysical priors, is not as widely used because of the computational difficulty of performing the marginalization. In this paper we address this difficulty and show how the marginalization can be done analytically by combining the four parameters into a set of complex amplitudes. The results of this paper are applicable to both transient non-precessing binary coalescence events, and to long lived signals from rapidly rotating neutron stars.

연구 동기 및 목표

• 베이지안 탐지에서 4개의 외재적 중력파 매개변수(A, ϕ₀, ψ, ι)에 대한 적분 문제를 해결하기 위한 계산적 과제를 다루기.
• 노이즈 매개변수를 회전군 행렬 원소 Tℓmn를 사용하여 재구성함으로써 B-통계량 적분의 비가역성 문제를 해결하기.
• 고신호대비비 비율 영역에서 B-통계량에 대한 고급 근사식을 개발하여, 일시적 및 연속 중력파 신호에 대한 실용적인 베이지안 탐색을 가능하게 하기.
• 단순화된 이진별 융합과 고속으로 회전하는 중성자별에서 발생하는 장기간 지속되는 신호에 모두 적용 가능한 통합 프레임워크 제공하기.
• 분석적 적분을 통해 계산 가능성을 유지하면서도 천체물리학적 사전 확률을 탐지에 활용할 수 있도록 하기.

제안 방법

• S³ 위의 호프 좌표를 사용하여 4개의 외재적 매개변수(A, ϕ₀, ψ, ι)를 복소 진폭으로 재구성함으로써, SU(2) 회전군을 매개변수화하기.
• 이 복소 진폭을 기반으로 우도 함수를 표현함으로써, 적분을 점점 더 점근적 분석에 적합한 형태로 변환하기.
• 가장 급격한 경로 방법을 적용하여 고신호대비비 비율에서 B-통계량의 근사값을 도출함으로써, 데이터 벡터 Y와 노이즈 가중 내적을 포함하는 폐쇄형 표현식을 도출하기.
• B-통계량에 대해 두 가지 다른 분석적 표현식을 유도함: 하나는 신호가 순수하게 원형 편광이 아닐 경우에 유효한 식(식 104), 다른 하나는 한 복소 진폭 성분이 0이 되는 특수한 경우에 유효한 식(식 111).
• 이로 인해 노이즈 매개변수에 대한 수치적 적분이 필요 없어지며, 계산 비용이 크게 감소함.
• 비특이 케이스에서의 지수 인자와 데이터 벡터 B̂에 대한 의존성은 F-통계량과 일치함을 보여 검증함. 특이 케이스에서는 비물리적 발산을 피하기 위해 다른 분석적 표현식을 사용함.

실험 결과

연구 질문

• RQ1 4개의 외재적 중력파 매개변수에 대한 우도 함수의 적분을 분석적으로 수행할 수 있는가? 이를 통해 효율적인 베이지안 탐지가 가능해지는가?
• RQ2 계산 비용이 큰 수치적 적분을 피하면서도 고신호대비비 비율 영역에서 B-통계량을 어떻게 근사할 수 있는가?
• RQ3 한 복소 진폭 성분이 0이 되는 특이 케이스(예: 순수한 원형 편광 신호)에서 B-통계량의 행동은 어떠한가?
• RQ4 현재 주로 주파수 모드에 국한된 파형에 제한되어 있지만, 이 형식은 고급 모드를 가진 비예약 시스템이나 예약하는 이진별에도 확장 가능한가?
• RQ5 천체물리학적 사전 확률이 포함된 경우, 제안된 방법이 표준 F-통계량과 탐지 능력과 계산 효율성 측면에서 어떻게 비교되는가?

주요 결과

• 논문은 고신호대비비 비율에서 B-통계량에 대한 폐쇄형 근사식을 도출하였으며, 이는 일시적 이진별 융합과 연속 중력파 신호 모두에 유효하다.
• B-통계량은 데이터 벡터 Y와 노이즈 가중 내적 B̂†N B̂를 기반으로 표현되며, 비특이 케이스에서는 F-통계량의 지수 의존성과 동일한 기능적 형태를 가짐.
• 특이 케이스(예: 순수한 원형 편광 신호)에서는 복소 진폭이 0이 되는 경우를 정규화하기 위해 (ζ² − k²)^1/4 요소를 포함하는 다른 분석적 표현식(식 111)을 사용함으로써 비물리적 발산을 피함.
• 결과적으로 도출된 B-통계량은 계산적으로 다룰 수 있으며, 진폭과 방향에 대한 천체물리학적 사전 확률을 활용할 수 있어 F-통계량보다 탐지 능력이 향상됨.
• 사전 확률이 자연의 진짜 분포와 일치할 경우, 이 방법이 네이만-피어슨 최적 탐지 기준과 일관됨을 입증함.
• 베이지안 프레임워크는 파rameter mismatch를 수용할 수 있으므로, 고급 모드나 예약하는 시스템으로의 확장이 자연스럽게 가능함.