Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses

문제 정의
LIGO–Virgo–KAGRA의 네 번째 중력파 카탈로그(GWTC-4.0)가 발표됨에 따라, 확신 있게 검출된 이진 블랙홀(BBH) 병합의 수가 150개를 넘어섰다. 초기 기대와 달리, 소수의 검출만으로 단일 주요 형성 채널을 식별하기에는 충분하지 않을 것이라는 예측이 있었으나, 현재 데이터는 단일 메커니즘으로는 쉽게 설명할 수 없는 복잡한 인구 집단을 보여준다. 스핀 기울기(spin tilts)를 이용해 하위 집단을 식별하려는 이전의 시도들은 측정의 어려움과 이론적 불확실성으로 인해 난관에 봉착했다. 또한, 매개변수 간의 전역적 상관관계가 존재함에도 불구하고 그 천체물리학적 해석은 여전히 쉽지 않은 문제이다. 본 연구에서 다루는 핵심 문제는, 일차 질량($m_1$)에서 명확한 전이 경계를 가지며 서로 다른 질량비($q$) 및 스핀 크기($\chi$) 분포를 특징으로 하는 구별된 BBH 병합 하위 집단이 존재하는지 여부이다.

방법론
저자들은 누적된 GWTC-4.0 데이터를 사용하여(이차적 모호성으로 인한 GW190814 및 파형 계통 오차와 극단적 스핀에 대한 우려로 인한 GW231123 제외) BBH 집단을 모델링하기 위해 현상론적 접근 방식을 채택한다. 이들은 다양한 복잡도 수준을 가진 모델들을 비교하기 위해 계층적 베이지안 추론을 사용한다.

주요 방법론적 특징은 다음과 같다:

• 매개변수화: 일차 질량 분포($p(m_1)$)는 두 개의 가우시안 피크를 가진 깨진 거듭제곱 법칙(broken power law)으로 모델링된다. 스핀 크기($p(\chi)$)와 기울기($p(\cos \theta)$)는 각각 절단된 정규 분포(truncated normals)와 등방성/정렬 성분의 혼합 모델을 사용한다. 적색편이($p(z)$)는 거듭제곱 법칙을 따른다.
• 전이 질량: 핵심적인 혁신은 $q$, $\chi$, $\cos \theta$의 분포가 최대 세 개의 전이 질량에서 형태가 변하도록 허용함으로써, 결과적으로 최대 네 개의 하위 집단을 생성하는 것이다. 분포의 함수 형태는 전이를 가로질러 동일하게 유지될 수 있으나, 형태를 제어하는 하이퍼파라미터는 서로 다르게 설정된다.
• 모델 선택: 저자들은 과적합을 피하기 위해 옥캄의 면도날(Occam penalty)을 적용하여 베이즈 인자(Bayes factors)를 통한 모델 선택을 수행한다. 이들은 "다중 집단(multipop)" 모델을 기준이 되는 강하게 매개변수화된 모델(gwtc4-fiducial) 및 특정 특징이 제거된 변형 모델들과 비교한다.

주요 기여 및 결과
본 분석은 두 개의 전이 질량 $m_t^{AB} = 27.7^{+4.1}_{-3.4} M_\odot$ 및 $m_t^{BC} = 40.2^{+4.7}_{-3.2} M_\odot$에 의해 구분되는, 최소 세 개의 구별된 병합 이진 블랙홀 하위 집단에 대한 강력한 증거를 제공한다. "다중 집단" 모델은 자연 로그 베이즈 인자 12.3로 기준 모델보다 선호된다.

세 하위 집단은 다음과 같이 특징지어진다:

1. 하위 집단 A ($m_1 \lesssim 27.7 M_\odot$):

   • 질량비: 거의 평탄한 분포($\beta_0 \approx -0.8$)를 보인다.
   • 스핀: 작은 스핀 크기를 특징으로 하며, 주로 $\chi \leq 0.5^{+0.1}_{-0.1}$ 범위 내에서 지지된다.
   • 해석: 효율적인 각운동량 수송을 포함한 고립 쌍성계 등 여러 형성 채널의 기여를 포함할 가능성이 높다.

2. 하위 집단 B ($27.7 M_\odot \lesssim m_1 \lesssim 40.2 M_\odot$):

   • 질량비: 등질량 쌍 형성($q \approx 1$)에 대한 훨씬 뚜렷한 선호를 보인다 ($\beta_1 \approx 3.6$). 이 독특한 쌍 형성 분포는 일차 질량 분포에서 관찰되었던 $\sim 35 M_\odot$ 부근의 피크를 해소하며, 이는 해당 피크가 모델 오지정(misspecification)에 의한 인공물이었음을 시사한다.
   • 스핀: 스핀 크기 및 기울기 분포는 하위 집단 A와 통계적으로 구별되지 않는다.
   • 해석: 저자들은 화학적 균질 진화(CHE) 또는 인구 III 항성을 잠재적 기원으로 제안하지만, 표준 CHE 예측과 스핀 차이가 없는 점을 조화시키는 데 따르는 어려움을 언급한다. 스핀 기울기 분포의 부재를 고려할 때 역학적 형성은 가능성이 낮은 것으로 간의된다.

3. 하위 집단 C ($m_1 \gtrsim 40.2 M_\odot$):

   • 스핀: 상당한 수준의 $\chi > 0.5$ 지지(약 $\chi \sim 0.7$에서 피크)를 포함하는 훨씬 넓은 스핀 크기 분포를 특징으로 한다.
   • 질량비: $q \approx 0.5$를 선호하는 경향에 대한 잠재적 지지를 보인다.
   • 해석: 속성은 넓은 스핀 분포와 $q \approx 0.5$ 선호를 고려할 때 계층적 병합(이전 세대 블랙홀의 병합)과 대체로 일치한다. 그러나 이 그룹 내에 작은 스핀이 존재한다는 점과 스핀 기울기의 전이가 나타나지 않는다는 점은 밀집 성단 내의 계층적 병합이 이 하위 집단을 완전히 설명하지 못할 수도 있음을 시사한다.

의의 및 주장
본 논문은 데이터가 이 세 하위 집단이 단일 연속적인 집단이 아니라 공존하는 상태로 존재한다는 해석을 강력히 선호한다고 주장한다. 질량비와 스핀 크기 분포에서 발견된 날카로운 전이 경계의 식별은 $\sim 25 M_\odot$ 이상에서 BBH 매개변수의 집단 분포에 대한 분리 가능성(separability) 가정이 무너짐을 의미한다.

저자들은 자신들의 결과가 강하게 매개변수화된 모델에 의존하며 모델 오지정의 가능성이 있으나, 신호의 통계적 강도와 다른 문헌 분석에서의 독립적인 전이 질량 확인은 이러한 하위 집단의 존재가 실재함을 시사한다고 강조한다. 결론적으로 서로 다른 형성 경로가 서로 다른 질량 범위에서 지배적일 가능성이 높으며, 하위 집단 B는 CHE 또는 인구 III 항성과 연결될 수 있고, 하위 집단 C는 계층적 병합과 연결될 수 있으나, AGN 디스크 내의 강착과 같은 대안적 시나리오도 여전히 가능하다. 본 연구는 $10 M_\odot$ 피크 주변의 전이에 대한 증거를 찾지 못했음을 명시하며, 해당 하위 구조는 향후 조사를 위해 남겨둔다.