The multimessenger view of Pulsar Timing Array black holes with the Horizon-AGN simulation

문제 정의
초거대 블랙홀(MBH)의 기원과 진화는 여전히 불확실하며, 특히 국부 우주(Local Universe)에서의 병합 과정이 그러하다. 최근 펄서 타이밍 어레이(PTA)는 거대 블랙홀 쌍성계(MBHB)의 확률론적 중첩인 중력파 배경(GWB) 신호를 검출했으나, 관측된 진폭과 이론적 예측 사이에는 지속적인 긴장이 존재한다. 이론적 예측은 대개 더 낮은 진폭을 지지하기 때문이다. 또한, 데이터셋의 유한한 기간과 배경 위에 서 있는 개별적으로 분해 가능한 연속파(CW) 소스의 존재 가능성은 GWB 스펙트럼의 해석을 복잡하게 만든다. 이론적 모델을 PTA 데이터와 비교하는 현재의 방법들은 편향을 초래할 수 있는 단순화된 가정을 사용하는 경우가 많다. 본 논문은 Horizon-AGN 우주론적 시뮬레이션을 활용하여 PTA 대역에 기여하는 MBHB의 특성을 연구하고, GWB 및 CW 예측을 위한 현실적인 모집단을 생성하며, 타이밍 잔차로부터 GWB 스펙트럼 특성을 추론하는 데 사용되는 방법론을 비판적으로 평가함으로써 이러한 문제들을 해결하고자 한다.

방법론
본 연구는 표준 $\Lambda$CDM 우주론을 따르는 $142 \, \text{cMpc}^3$ 부피를 모델링하는 Horizon-AGN 유체역학적 우주론적 시뮬레이션을 활용한다. 저자들은 시뮬레이션 내에서 MBH ID를 추적하고 수치적 병합을 식별함으로써 MBH 병합 카탈로그를 구축하였다. 수치적 병합(약 $4 \, \text{kpc}$ 분리)과 실제 합병 사이의 물리적 지연을 고려하기 위해, 저자들은 다음을 포함하는 후처리 파이프라인을 적용하였다:

1. 역학적 마찰(Dynamical Friction): 등온 구체 근사 내에서의 항성 역학적 마찰을 기반으로 계산되었다.
2. 쌍성 경화(Binary Hardening): 기존의 처방(예: Sesana & Khan 2015; Dotti et al. 2015)을 따라 항성 경화, 원형 궤도 주변 디스크의 점성 토크, 그리고 중력파(GW) 방출을 통해 모델링되었다.
3. 이심률 진화(Eccentricity Evolution): 로스-콘(loss-cone) 산란, 점성 항력, 그리고 GW 방출을 통해 추적되었으며, 이전의 수치적 결과에 기반하여 초기 이심률을 할당하였다.

GWB를 추정하기 위해 저자들은 두 가지 별개의 접근 방식을 사용하였다:

• 가우시안 앙상블(Gaussian Ensemble, GE) 접근법: 소스의 연속적인 분포를 가정하는 결정론적 특성 변형률 스펙트럼 계산으로, 문헌에서 자주 사용된다.
• 모집단(Population) 접근법: 병합 카탈로그를 바탕으로 포아송 분포로부터 입사 중인 쌍성계를 추출하여 2,000개의 "우주 실현(Universe realizations)"을 생성하는 확률론적 방법이다. 이는 우주론적 분산(cosmic variance)과 소스의 이산적 특성을 고려한다.

저자들은 두 가지 방법을 사용하여 이러한 모집단에 대한 PTA 타이밍 잔차를 시뮬레이션하였다:

1. 가우시안 앙상블: 평균화된 펄서 응답을 통해 특성 변형률을 타이밍 잔차 전력 스펙트럼 밀도(PSD)로 변환한다.
2. 모집단 접근법: 모든 쌍성의 지구 항(Earth-term) GW 변형률 기여를 합산하며, 유한한 관측 기간($T_{\text{obs}}$)과 sinc 함수를 통한 스펙트럼 누설 효과를 명시적으로 모델링한다.

저자들은 이러한 이론적 예측을 EPTA(25개 펄서, 10.33년) 및 미래의 SKA(80개 펄서, 10.33년) 감도를 나타내는 "토이(toy)" PTA 데이터셋과 비교하였다. 또한 잠재적인 CW 후보(신호 대 잡음비(SNR) $> 3$인 개별 쌍성계)를 식별하고, AGN 방출 및 호스트 은하의 특성을 모델링하여 이들의 전자기적 특성(광학, NIR, 라디오, X-선)을 분석하였다.

주요 기여 및 결과

• GWB 스펙트럼 추론 및 편향: 본 연구는 이론적 GWB 모델을 PTA 데이터와 비교하는 데 있어 중요한 방법론적 주의 사항을 강조한다. 무한한 관측 기간을 가정하는 표준 가우시안 앙상블 접근법은 유한한 $T_{\text{obs}}$로 인한 스펙트럼 누설을 설명하지 못한다. 현실적인 PTA 분석(타이밍 모델 마진화 포함)과 함께 모집단 접근법을 사용할 경우, 추론된 GWB 스펙트럼은 실제 기저 스펙트럼에 비해 더 완만한 기울기와 더 높은 진폭을 갖는 방향으로 편향된다. 저주파 소스로부터 발생하는 스펙트럼 누설에 의해 유발되는 이 편향은 이론적 예측과 관측된 PTA 진폭 사이의 긴장을 부분적으로 설명할 수 있다.
• 확률론적 성질 및 소스 개수: 저주파 영역($\sim 3 \, \text{nHz}$)에서 GWB는 수백 개에서 수천 개의 쌍성에 의해 지배되며, 이는 확률론적 배경 근사의 타당성을 확인시켜 준다. 그러나 고주파 영역($\sim 10\text{--}50 \, \text{nHz}$)에서는 기여하는 소스의 수가 급격히 감소하여(수십 개 수준), 확률론적 성질의 상실과 증가된 우주론적 분산을 초래한다. 이심률이 있는 모집단은 높은 궤도 고조파에서의 방출 덕분에 원형 궤도 모집단보다 더 높은 주파수까지 확률론적 성질을 유지한다.
• 소스 특성 (배경 vs. CW):
  • 배경(Background): GWB는 적색편이 $z \in 0.05 - 1$ 범위에서 질량 척도(chirp mass) $M_c \sim 10^{8.5} - 10^{9.5} \, M_\odot$를 가진 쌍성들에 의해 지배된다. 이들은 주로 질량 $M_\star > 10^{11} \, M_\odot$인 거대 은하와 질량 $\sim 10^{13} \, M_\odot$인 헤일로 내에 위치한다.
  • CW 후보: 개별적으로 분해 가능한 소스(CW 후보)는 일반적인 배경에 비해 더 무겁고 더 낮은 적색편이($z < 0.4$)에 위치하는 경향이 있다. 이들은 주로 거대 그룹 및 클러스터($M_{\text{halo}} > 10^{13} \, M_\odot$)의 중심 은하에 위치한다.
• 전자기적 대응 천체:
  • 성격: 대부분의 CW 후보는 퀘이사라기보다는 활동 은하 핵(AGN)의 형태를 띠며, 낮은 에딩턴 비율($f_{\text{Edd}} < 10^{-2}$)을 특징으로 한다.
  • 탐지 가능성: 광학 및 근적외선(NIR) 대역에서 대부분의 CW 후보는 호스트 은하에 의해 빛이 가려지므로, 광도 변화 탐색는 어려우며 희귀하고 밝은 시스템에 편향될 수 있다.
  • 최적의 파장: 라디오와 X-선 대역이 식별을 위한 가장 좋은 전망을 제공한다. $f_{\text{Edd}} > 10^{-3}$인 후보의 경우, AGN 방출이 호스트 은하를 압도할 수 있다. 낮은 강착 시스템이라 하더라도 낙관적인 추정치에 따르면 최대 80%(라디오) 및 30%(X-선)가 호스트보다 탐지 가능할 수 있다.
• 탐지 확률: 10년의 기준선을 가정할 때, 적어도 하나의 CW 후보를 $S/N > 3$으로 탐지할 확률은 EPTA급 감도에서 4%, SKA급 감도에서 20%로 추정된다.

의의
본 논문은 최첨단 우주론적 시뮬레이션을 사용하여 PTA 신호를 이해하기 위한 통합된 다중 신호(multi-messenger) 프레임워크를 제공한다는 점에서 주요한 의의를 갖는다. 결정적으로, 본 연구는 이론적 GWB 모델을 관측 데이터와 비교하는 표준 방법론에 이의를 제기한다. 유한한 관측 기간과 스펙트럼 누설이 추론된 GWB 진폭을 인위적으로 팽창시키고 스펙트럼 기울기를 완만하게 만들 수 있음을 보여줌으로써, 저자들은 이론적 모델과 PTA 관측 사이의 명백한 긴장이 기초적인 천체물리학의 실패라기보다는 분석 기법의 산물일 수 있음을 시사한다. 또한, 본 연구는 전자기적 후속 관측을 위한 구체적인 지침을 제공하며, 퀘이사 같은 시스템의 광학 변광성에만 의존하기보다 라디오 및 X-선 대역을 우선시하고 거대 은하 그룹 및 클러스터에 집중할 것을 제안한다.