자기상호작용 암흑물질 스파이크와 최종 파섹 문제: NANOGrav 15년 배경에 대한 베이지안 제약

초록

이 논문은 초대질량 블랙홀(SMBH) 병합 과정에서 자기상호작용 암흑물질(SIDM) 스파이크가 제공하는 동역학 마찰이 1 pc에서 10⁻² pc까지 이진을 수축시켜 최종 파섹 문제를 해결할 수 있음을 제시한다. 코스모로지적 SMBH 군집 모델에 SIDM 마찰과 중력파 방출을 동시에 적용해 나노헤르츠 중력파 배경을 계산하고, NANOGrav 15년 펄서 타이밍 데이터와 베이지안 MCMC(가우시안 프로세스 가속) 분석을 수행해 σ/m(단위 질량당 단면적)과 최대 원형 속도(vₜ) 등에 대한 사후 확률을 도출한다. 결과는 은하 회전 곡선 및 은하단 렌즈링 제한과 일치하면서, SIDM 스파이크가 파괴되지 않고 충분한 마찰을 제공함을 보여준다.

상세 분석

본 연구는 최종 파섹 문제를 SIDM 스파이크가 해결할 수 있다는 가설을 정량적으로 검증한다. 먼저, NFW 프로파일을 기반으로 하여 중심부에 자기상호작용에 의해 형성된 등온 코어와 그 안에 급격히 상승하는 스파이크를 포함하는 3구역 밀도 모델을 구축한다. 코어 반경 r₁과 속도 분산 v₀는 σ/m·t_age 조건과 질량 보존 조건을 동시에 만족하도록 수치적으로 풀며, 스파이크 반경 r_sp=GM_BH/v₀²와 지수 γ는 상호작용 유형(a=0 또는 a=4) 및 전이 속도 v_t에 따라 결정된다. 이때 전이 반경 r_t는 v(r)=v_t을 만족하는 지점으로, v₀<v_t이면 스파이크가 얕은 γ=3/4 구간과 급격한 γ=7/4 구간을 순차적으로 갖고, v₀>v_t이면 전체가 γ=7/4 형태가 된다.

동역학적으로는 이진의 원형 궤도에서 중력파 방출에 의한 에너지 손실 dE_GW/da와 SIDM 마찰에 의한 에너지 손실 dE_SIDM/da를 합산한다. SIDM 마찰력은 ρ(r)·σ/m·v_rel·v_rel 형태로, 코어와 스파이크 각각의 밀도와 상대 속도를 사용해 계산한다. 이 두 손실을 이용해 이진 반경 a(t)의 감쇠율 da/dt를 구하고, 이를 적분해 각 이진이 방출하는 특성 변형 h_c(f) 스펙트럼을 얻는다.

전체 SMBH 군집에 대해 은하-블랙홀 질량 관계와 은하-광역질량 관계를 이용해 M_BH→M_200→NFW 파라미터(ρ_s, r_s)를 역산하고, 각 이진에 대해 위의 과정을 적용해 GWB 스펙트럼 라이브러리를 만든다. 이 라이브러리를 가우시안 프로세스로 보간해 MCMC 샘플링을 가속화하고, NANOGrav 15년 데이터의 파워 스펙트럼과 직접 비교한다. 사후 분포는 σ/m≈0.1–1 cm² g⁻¹와 v_t≈300–1500 km s⁻¹ 범위에 집중되며, 이는 기존 은하 회전곡선·은하단 렌즈링 제한과 일치한다. 또한, 이 파라미터 영역에서는 스파이크가 파괴되지 않고 충분한 마찰을 제공해 이진이 10⁻² pc 이하까지 효율적으로 수축한다는 것이 확인된다.

핵심 인사이트는 (1) SIDM 스파이크가 코어와 연계돼 동역학 마찰을 지속적으로 공급함으로써 별 손실 장벽을 뛰어넘을 수 있다, (2) GWB 관측이 SIDM 파라미터 공간을 독립적으로 제한할 수 있는 새로운 우주론적 탐색 도구가 된다, (3) 가우시안 프로세스 가속 MCMC가 복잡한 천체물리 모델의 베이지안 추론에 실용적임을 보여준다.