From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

기술 요약: 우주론적 시뮬레이션에서 이진 블랙홀 병합까지

문제 제기
이진 블랙홀(BBH) 병합에 대한 관측은 거대 항성의 진화 과정을 이해하는 데 중요한 통찰을 제공한다. 그러나 관측된 병합 특성을 그 전구체 항성과 연결하기 위해서는 정밀한 우주 별 형성사(SFRH) 및 화학적 풍부도 모델이 필요하다. 현재 이 분야에서는 복잡한 기저의 역사(우주론적 시뮬레이션이나 관측으로부터 도출된)를 대변하기 위해 금속 함량 의존적 우주 별 형성율 밀도, $S(Z, z)$에 대한 해석적 적합(analytical fits)을 빈번하게 사용한다. 이러한 해석적 처방들은 종종 제약 조건이 불분명하며, 추론된 BBH 집단에 미치는 영향 또한 충분히 파악되지 않았다. 또한, 우주론적 시뮬레이션의 해상도와 박스 크기가 도출된 $S(Z, z)$ 및 그에 따른 BBH 병합율에 미치는 영향이 체계적으로 정량화되지 않았다.

방법론
저자들은 해상도와 박스 크기가 서로 다른 IllustrisTNG 시뮬레이션 제품군(구체적으로 TNG50-1, TNG100-1, TNG300-1)을 활용하여, 가스 셀 데이터로부터 직접 시뮬레이션 기반의 $S(Z, z)$를 도출한다. 이들은 직접 도출된 시뮬레이션 결과와 널리 사용되는 해석적 적합 모델(L. A. C. van Son 등 2023에 기반함)을 비교한다. 해당 모델은 총 별 형성 밀도 $S(z)$에 대한 Madau & Dickinson (2014)의 함수 형태와 왜곡된 로그-노멀(skewed log-normal) 금속 분포를 결합한 것이다.

$S(Z, z)$ 입력값이 BBH 집단에 미치는 영향을 평가하기 위해, 저자들은 시뮬레이션 유래 $S(Z, z)$ 모델과 해석적 $S(Z, z)$ 모델 모두를 COMPAS(Compact Object Mergers: Population Astrophysics and Statistics) 급속 이진 인구 통계 합성 코드에 결 결합한다. 이들은 적색편이(redshift)에 따른 병합율과 질량 분포를 계산하기 위해 "우주 통합(cosmic integration)" 프레임워크를 채택한다. 본 연구는 구체적으로 다음 사항을 조사한다:

1. 해석적 적합이 다양한 TNG 해상도에서의 $S(Z, z)$의 형태와 정규화(normalization)를 재현하는 충실도.
2. 해석적 적합과 전체 시뮬레이션 데이터를 사용할 때 발생하는 BBH 집단의 병합율 및 일차 질량 분포($M_{BH,1}$)의 편차.
3. 시뮬레이션 해상도(고해상도 TNG50-1부터 저해상도 TNG300-1까지)가 이러한 결과에 미치는 영향.

주요 결과
본 연구는 $S(Z, z)$에 대한 해석적 처방이 BBH 집단을 모델링할 때 상당한 계통 오차를 유발한다는 것을 발견했다:

• $S(Z, z)$의 불일치: 해석적 적합은 시뮬레이션 기반의 금속 분포에 나타나는 핵심적인 구조적 특징을 포착하지 못한다. 구체적으로, 시뮬레이션에 존재하는 고금속 영역($Z \gtrsim 3 Z_\odot$)의 이차적인 "범프(bump)"와 평탄한 저금속 꼬리 부분을 놓친다. 이 적합 모델은 고적색편이($z \gtrsim 6$) 및 고금속성에서의 별 형성율을 과대평가하는 경향이 있는 반면, 저금속성 및 매우 높은 적색편이에서의 기여도는 시뮬레이션 해상도에 따라 과소평가한다.
• 병합율에 미치는 영향: 고적색편이($z \gtrsim 6$)에서, 해석적 적합은 시뮬레이션 기반 병합율에 비해 10에서 $10^4$배까지 BBH 병합율을 과대평가하며, 시뮬레이션 해상도가 낮아질수록 오차의 크기가 증가한다. 저적색편이($z \approx 0$)에서, 적합 모델로부터 도출된 국지적 병합율은 일반적으로 시뮬레이션 기반 병합율과 약 1.5배 이내의 오차로 일치하지만, TNG50-1의 경우 $\sim 15$ Gpc$^{-3}$ yr$^{-1}$의 과대평가를 보인다.
• 질량 분포에 미치는 영향: 해석적 적합을 사용하면 BBH 일차 질량 분포에 가짜 특징(spurious features)이 도입된다. 가장 눈에 띄는 점은 $z \lesssim 2$에서, 적합 모델은 전체 시뮬레이션 기반 $S(Z, z)$를 사용할 때는 나타나지 않는 일차 질량 분포 근처 $8 M_\odot$ 지점의 인위적인 과잉(bump)을 생성한다. 이 특징은 안정적 질량 전이(stable mass transfer) 채널에 기인하며, 고금속성 별 형성 피크를 올바르게 모델링하지 못하는 적합 모델의 한계로 인해 발생한다. 반대로, 적합 모델은 $M_{BH,1} \gtrsim 35 M_\odot$ 영역에서 고질량 특징(근처 $33 M_\odot$ 및 $45 M_\odot$)을 3~5배 정도 억제한다.
• 해상도 의존성: 시뮬레이션 해상도는 $S(z)$의 정규화와 피크 위치에 상당한 영향을 미친다. 고해상도 시뮬레이션(TNG50-1)은 저해상도 실행($z \approx 2.5$)에 비해 더 높은 별 형성율을 보이며, 약간 더 낮은 적색편이($z \approx 2$)에서 피크를 형성한다. 해석적 적합의 정확도는 시뮬레이션 해상도가 낮아질수록, 특히 고적색편이 기울기를 재현하는 데 있어 저하된다.

의의 및 주장
저자들은 해석적 적합이 계산상의 편의를 제공하기는 하지만, 특히 중력파 관측을 해석할 때 시뮬레이션으로부터 도출된 복잡한 기저의 $S(Z, z)$를 대체하기에는 신뢰할 수 없다고 결론짓는다. 본 연구는 단순한 해석적 모델이 시뮬레이션에서 발견되는 고금속성 범프 및 특정 저금속성 꼬리와 같은 결정적인 물리적 특징을 포착하는 데 실패함을 입증했다.

이러한 불일치는 현재의 해석적 적합의 함수 형태가 시뮬레이션에서 발견되는 다봉형(multi-modal) 금속 분포를 수용할 수 없기 때문에 발생한다고 저자들은 설명한다. 결과적으로, 이러한 적합 모델에 의존하는 것은 BBH의 형성 채널 및 진화사에 대한 잘못된 추론으로 이어질 수 있다. 저자들은 중력파 데이터의 해석 과정에서 가짜 특징을 도입하는 것을 방지하기 위해, 우주론적 시뮬레이션, 해석적 적합, 그리고 인구 통계 합성을 신중하게 통합해야 한다고 주장한다. 저자들은 해석적 적합이 완전히 무용하다는 주장을 하는 것이 아니라, 특히 차세대 검출기가 더 높은 적색편이 영역을 탐사함에 따라 이들의 한계를 강조하고 고충실도 시뮬레이션 데이터에 의한 검증의 필요성을 역설하고 있다.