[논문 리뷰] Rapid inference of gravitational-wave signals in the time domain using a heterodyned likelihood

Neha Sharma, A. Vijaykumar외 1명|arXiv (Cornell University)|2026. 01. 16.

Pulsars and Gravitational Waves Research인용 수 0

한 줄 요약

시간 영역 이종화된 가능도를 도입하여 파동 모드의 다운샘플링과 미리 계산된 요약 데이터를 사용함으로써 긴 중력파 신호에 대한 베이지안 매개변수 추정을 극적으로 가속화하고, 포스터리는 구별 불가능한 수준으로 유지된다.

ABSTRACT

Parameter estimation of gravitational wave signals is computationally intensive and typically requires millions of likelihood evaluations to construct posterior probability distributions. This computational cost increases significantly in the time domain, which requires non-diagonal covariance matrices to compute the likelihood. Consequently, parameter estimation of long-duration gravitational wave signals, such as binary neutron star mergers, becomes computationally infeasible in time domain. In this work, we detail a framework for the heterodyned likelihood that enables rapid inference in the time domain. Our method is applicable to signals with arbitrary mode content, and leverages the smoothness of the ratio of complex-valued waveform modes, approximating the ratio as a linear function within appropriately chosen time bins. This allows downsampling of the waveform modes and a reformulation of the likelihood, such that it depends only on the bin edges. We demonstrate that this likelihood recovers posteriors that are indistinguishable from those obtained using the standard likelihood in the time domain. We also observe dramatic improvement in speed - for a 128 seconds-long gravitational wave signal, our method is at least $\sim 400$ times faster than the standard time-domain analysis, reducing the wall clock time to just a few hours. We also demonstrate the reliability and unbiasedness of the likelihood using percentile-percentile tests for binary black hole and binary neutron star injections. We use the Gohberg-Semencul representation of the inverse of Toeplitz covariance matrix to accelerate matrix-vector products, which has potential applications even in non heterodyned time-domain inference.

연구 동기 및 목표

• 시간 영역 GW 분석에서 높은 계산 비용으로 인한 신속한 매개변수 추정의 필요성에 대한 동기 부여.
• waveform 모드 비율의 선형 근사를 시간 구간 내에서 사용하는 이종화된 가능도 프레임워크를 개발합니다.
• 가능도 평가를 가속화하되 정확성을 유지하기 위한 요약 데이터와 빙닝 전략을 형식화합니다.
• BBH 및 BNS 시나리오에서 주입 및 p-p 검사들을 통해 정확도와 속도 향상을 입증합니다.

제안 방법

• 가능도를 fiducial waveform h^{3lm}_o에 대한 복소 파동 모드 h^{3lm}의 비율의 관점으로 표현하고 이 비율을 시간 구간 내에서 선형 근사(이종화)합니다.
• 필드 구간 끝점에서만 모드들을 다운샘플링하고 필요할 때 보간으로 전체 파동을 재구성합니다.
• 샘플링 도중 큰 행렬 곱셈을 피하기 위해 fiducial 모드, 검출기 노이즈 및 데이터에 의존하는 요약 데이터를 미리 계산합니다.
• 로그 가능도를 비율 r^{3lm}와 요약 데이터를 사용하여 표현하고, 매 반복마다의 행렬 역행렬 및 곱을 피합니다.
• 진화기에 따라 적응적 시간 빙닝 방식을 도입하여 Inspiral 구간에서 빙을 확장하고 합성 near merger 구간에서 확장 또는 축소하여 정확도와 속도 사이의 균형을 맞춥니다.

실험 결과

연구 질문

• RQ1 이종화된 시간 영역 가능도가 전체 시간 영역 가능도와 동일한 포스터리를 재현할 수 있는가?
• RQ2 시간 영역 이종화로 다양한 지속 시간의 GW 신호에 대해 어떤 측면의 속도 향상을 얻을 수 있는가?
• RQ3 특히 합성 근처에서 정확도를 유지하면서 계산을 최소화하기 위해 시간 빙은 어떻게 선택해야 하는가?
• RQ4 부분지배 모드와 프리세싱 파형으로의 접근 방식이 확장 가능한가?
• RQ5 BBH 및 BNS 주입에서 분위수-주입Tests이 방법의 바이어스 없는 포스터리를 어떻게 검증하는가?

주요 결과

GW 신호 지속 시간 [초]	빈 수	요약 데이터 계산 시간 [초]	이종화된 가능도 [ms]	전체 가능도 [ms]	속도 상승 요인	빈 수 (주파수 영역)	이종화된 가능도 [ms]	전체 가능도 [ms]	(주파수 도메인 속도 상승)
2	191	8.6989	2.290	42.9296	120	119	1.445	4.596
4	193	53.9353	2.418	89.97518	121	121	1.724	8.662
8	302	116.1305	3.860	216.558	121	121	1.812	14.14
16	382	561.3549	6.378	791.557	124.10740	121	1.801	27.35
32	476	967.96419	8.580	884.42	103.07940	121	1.848	51.18
64	483	4103.1448	9.140	4163.54	455.5295	121	1.817	105.3
128	486	20012.4119	9.191	3129.884	340.5379	121	1.849	209.6

• 이종화된 시간 영역 가능도에서의 후진 분포는 테스트된 주입에 대해 전체 가능도 포스터리와 구별될 수 없었다.
• 2초 신호의 경우 이종화가 전체 시간 영역 방법에 비해 가능도 평가를 약 19배 빠르게 만들었다.
• 16초 및 128초 신호의 경우 단일 가능도 평가가 각각 약 6.4 ms 및 약 9.2 ms로, 속도 상승은 각각 약 124배 및 약 340배였다.
• 속도 상승은 신호 지속 시간에 비례하여 확장되며, 주간 단위에서 수주 내로 분석을 완료할 수 있게 한다.
• 적응적 빙 전략은 근처 합성 근처에서 정확도와 파형 평가 수를 급격히 줄이면서도 정확도를 유지한다.
• BBH 및 BNS 주입에 대한 p-p 검정은 편향되지 않은 포스터리를 보여 주의성을 확인시킨다.