Prospects of resolving and localising individual supermassive black hole binaries with pulsar timing arrays: the host ranking challenge

기술 요약: 펄서 타이밍 어레이를 이용한 개별 초거대 블랙홀 쌍성계의 해상 및 국지화 전망

문제 정의
펄서 타이밍 어레이(PTA)는 최근 초거대 블랙홀(SMBH) 쌍성계 집단에서 기원하는 것으로 보이는 확률적 중력파(GW) 배경을 검출했다. 다음 주요 이정표는 이 배경 내에서 개별 쌍성계를 해상(resolution)하는 것이며, 이는 다중 신호원 천문학(multi-messenger astronomy)을 가능하게 할 것이다. 그러나 주요 과제는 PTA의 제한된 하늘 국지화(sky localisation) 능력이다. 최초로 해상될 쌍성계의 90% 신뢰도 국지화 영역($A_{90}$)은 수백 평방 도(square degrees)에 달할 것으로 예상되며, 여기에는 수십만 개의 후보 은하가 포함될 수 있다. 또한, 기존의 전자기파(EM) 은하 카탈로그는 별이나 먼지 오염이 심한 영역(예: 은하 평면)에서 불완전하며, 이로 인해 "누락된" 호스트 후보들이 발생한다. 현재는 방대한 수의 후보들을 GW 사후 확률(posteriors) 및 EM 특성에 따라 효율적으로 순위를 매겨 후속 관측의 우선순위를 정할 표준화된 방법이 없다.

방법론
저자들은 holodeck 소프트웨어를 사용하여 NANOGrav 15년 데이터 세트(NG15)에서 측정된 중력파 배경과 일치하는 현실적인 SMBH 쌍성계 집단을 시뮬레이션한다. 그들은 1,000개의 쌍성계 집단 실현(realizations)을 생성하고, 미래의 PTA 구성에서는 검출 가능하지만 NG15에서는 검출 불가능한 가장 강력한 쌍성계를 선택한다.

세 가지 PTA 구성을 미래의 민감도를 나타내기 위해 시뮬레이션한다:

• IPTA_20: 116개의 펄서를 가진 20년 기준선 (다가오는 IPTA 3차 데이터 릴리스와 유사).
• IPTA_25 (표준 모델): 158개의 펄서를 가진 25년 기준선 (IPTA_20에 42개의 새로운 펄서 추가).
• IPTA_30: 200개의 펄서를 가진 30년 기준선.

저자들은 선택된 쌍성계들을 확률적 중력파 배경을 포함하는 시뮬레이션된 타이밍 데이터에 주입하고, QuickCW 파이프라인을 사용하여 쌍성계 파라미터(하늘 위치, 처프 질량, 변형률, 주파수)를 제약한다. 저자들은 "국지화 가능한" 소스를 잘 제약된 하늘 위치를 위해 일반적으로 요구되는 임계값인 신호 대 잡음비(S/N) $\ge 8$을 가진 것으로 정의한다.

호스트를 식별하기 위해, 저자들은 결과물인 $A_{90}$ 국지화 지도를 두 가지 포괄적인 전천(all-sky) EM 카탈로그와 교차 참조한다:

1. ETG 카탈로그: 2MPZ 및 WISExSCOS로부터 유도된 광도 적색편이와 항성 질량을 포함하는 약 1,940만 개의 정적 초기형 은하(ETG)를 포함한다.
2. AGN 카탈로그: W1 밴드 겉보기 등급을 가진 약 450만 개의 활동 은하 핵(AGN) 후보 (Assef et al. 2018의 R90 카탈로그)를 포함한다.

각 후보 은하에 점수($H$)를 부여하는 호스트 순위 시스템이 개발되었다. 점수는 하늘 위치 확률($p_{sky}$)과 은하의 EM 특성이 GW 사후 확률과 일치할 확률($p_{cat}$)의 곱이다.

• ETG의 경우, $p_{cat}$은 은하의 적색편이(광도 거리로 변환됨)와 추정된 SMBH 질량(항성 질량으로부터 스케일링 관계를 통해 유도됨)으로부터 도출된다.
• AGN의 경우, $p_{cat}$은 쌍성계의 총 질량 및 광도 거리의 함수인 겉보기 W1 등급으로부터 도출된다.
  저자들은 EM 카탈로그가 커버하지 못하는 $A_{90}$ 영역의 비율을 추정함으로써 "누락된" 호스트의 수를 명시적으로 계산한다.

주요 결과

• 검출 가능성 및 국지화:

  • 향후 0~5년 내(IPTA_20)에, 쌍성계를 해상할 확률($S/N \ge 4$)은 약 21%이지만, 잘 제약된 국지화를 달성할 확률($S/N \ge 8$)은 약 0.3%에 불과하다.
  • 향후 5~10년 내(IPTA_25)에, 쌍성계를 해상할 확률은 약 38%로 상승하며, 국지화할 확률은 약 3.8%이다.
  • 향후 10년 내(IPTA_30)에, 해상 확률은 약 51%에 도달하며, 국지화 확률은 약 14.1%에 달한다.
  • 국지화 가능한 소스들은 상대적으로 가까이 있으며(중앙값 $d_{lum} \approx 400$ Mpc), 매우 무거운 질량($M_{tot} \gtrsim 10^{9.8} M_{\odot}$)을 가질 것으로 예상된다.

• 국지화 영역 및 호스트 수:

  • 표본 구성(IPTA_25)에 대해, 잠재적으로 국지화 가능한 38개의 쌍성계 중 30개는 $A_{90} < 2,000$ deg$^2$이다. 중앙값 $A_{90}$은 $408$ deg$^2$이다.
  • 이 영역들은 중앙값으로 약 190,000개의 ETG와 약 40,000개의 AGN을 포함한다.
  • IPTA_30에서는 중앙값 $A_{90}$이 $212$ deg$^2$로 줄어들며, 중앙값 호스트 수가 약 40% 감소한다.

• 불완전한 커버리지:

  • 불완전한 EM 커버리지는 중요한 제한 요소이다. 표본 구성에서 누락된 ETG 호스트의 중앙값은 약 20,000개이며, 누락된 AGN 호스트는 약 6,000개이다.
  • 일부 소스의 경우, $A_{90}$이 카탈로그의 발자국(footprint) 밖에 부분적으로 또는 완전히 걸쳐 있어(예: 은하 평면과 겹침), 새로운 조사 없이는 실제 호스트를 식별할 수 없다.

• 순위 효율성:

  • 순위 방법은 GW 사후 확률과 일치하지 않는 특성을 가진 은하들을 성공적으로 제외한다.
  • ETG: 약 47% (IPTA_25)에서 61% (IPطا_30)의 ETG가 부적절한 적색편이 또는 SMBH 질량으로 인해 제외된다. 순위 효율성($E_{95-50}$, 95 백분위수 점수와 중앙값의 비율)은 IPTA_25에서 약 12.4, IPTA_30에서 약 17.2이다.
  • AGN: 약 0.3% (IPTA_25)에서 3.8% (IPTA_30)의 AGN만이 제외된다. 순위 효율성은 현저히 낮으며 (IPTA_25에서 약 3.7), 이는 겉보기 등급이 질량과 거리의 퇴화된 조합이며, PTA에 의해 처프 질량이 잘 제약되지 않기 때문이다.

의의 및 주장
본 논문은 최초의 PTA 해상 SMBH 쌍성계에 대한 후보 호스트의 수와 불완전한 하늘 커버리지의 영향을 현실적으로 정량화한 첫 사례라고 주장한다. 저자들은 비록 후보의 수가 많지만, GW 사후 확률 및 EM 특성에 기반한 순위 시스템이 특히 적색편이와 항성 질량을 사용할 수 있는 ETG의 경우 탐색 범위를 효과적으로 줄일 수 있음을 입증한다.

연구의 결론은 다음과 같다:

1. 다중 신호원 발견은 가능하지만 도전적이다: 10년 이내에 소스를 국지화할 확률은 무시할 수 없는 수준($\approx 14\%$)이지만, 국지화 영역은 여전히 크다.
2. 카탈로그 완결성이 결정적이다: 조사 공백(예: 은하 평면)으로 인한 수천 개의 누락된 호스트의 존재는 호스트 식별 시 위음성(false negative)의 위험을 초래한다.
3. 데이터 품질이 순위 성공을 결정한다: 순위 방법의 효율성은 독립적인 적색편이 및 SMBH 질량 추정치의 가용성에 크게 의존한다. 현재 AGN을 효율적으로 순위 매기지 못하는 이유는 AGN 카탈로그에 적색편이/질량 데이터가 부족하고, PTA가 처프 질량을 잘 제약하지 못하기 때문이다.
4. 미래의 개선 사항: 저자들은 펄서 거리 측정의 개선(GW 신호에서 펄서 항을 활용하기 위해)과 적색편이 및 질량 추정치를 포함한 AGN 카탈로그의 개발이 국지화 및 순위 능력을 크게 향상시킬 것이라고 제안한다.

본 논문은 최초의 다중 신호원 탐지가 충분히 도달 가능한 범위 내에 있지만, 그 성공은 PTA의 민감도뿐만 {만} EM 카탈로그의 완결성과 균일성, 그리고 후속 관측의 우선순위를 정하기 위한 효율적인 순위 파이프라인의 개발에 달려 있다고 결론짓는다.