A Stellar Role Reversal: Multiple Features in the Mass and Mass Ratio Distributions of Merging Binary Black Holes from Stable Mass Transfer

기술 요약: 별의 역할 역전: 안정적 질량 이동으로부터 기인한 병합 이중 블랙홀의 질량 및 질량비 분포 내 다중 특징

문제 제기
중력파(GW) 관측 사례가 늘어나고 있음에도 불구하고, 병합되는 이중 블랙홀(BBH) 집단의 기원은 여전히 불확실하다. 관측을 통해 BBH 질량 함수(예: $10 M_\odot$ 및 $35 M_\odot$ 근처의 피크)와 질량비 분포에서 특징들이 드러났으나, 측정 가능한 성질들의 퇴행성과 항성 진화 및 이중성 상호작용의 불확실성으로 인해 서로 다른 형성 채널의 기여를 분리해내는 것은 매우 어렵다. 본 연구는 모든 로슈 로브 초과(RLOF) 단계가 안정적으로 유지되는 '안정적 질량 이동(Stable Mass Transfer, SMT-only)' 채널을 구체적으로 조사하며, 공통 외층(Common Envelope, CE) 채널은 제외한다. 저자들은 질량 이동 흡수 효율($\beta$)에 대한 가정이 병합되는 BBH의 일차 질량(primary mass) 및 질량비 분포에 어떻게 뚜렷한 특징을 각인시키는지 규명하는 것을 목표로 한다.

방법론
본 연구는 이진 집단 합성(BPS)과 해석적 프레임워크를 결합한 이중 접근 방식을 사용한다:

1. 집단 합성: 저자들은 cosmic 코드(v4.1.0)를 사용하여 병합되는 BBH의 우주론적 집단을 시뮬레이션한다.

   • 초기 조건: 시뮬레이션은 30개의 금속 함량($10^{-4}$ ~ $0.03$),$0.08$~$100 M_\odot$ 범위의 초기 질량 함수(Kroupa 2001), 그리고 0.7의 이중성 분율을 포함하는 격자를 다룬다.
   • 물리 법칙: 모델에는 금속 함량 의존적 항성풍(Vink et al. 처방), 조석 원형화, 그리고 공여체의 진화 상태(Case A, B, C, BC)에 따른 특정 질량 이동 안정성 기준이 포함된다.
   • 컴팩트 천체 형성: 블랙홀(BH) 질량은 "지연된(delayed)" 초신성 처방(Fryer et al. 2012)을 사용하여 할당되며, 최소 BH 질량은 $3 M_\odot$이다. 쌍 불안정성 초신성(PISN)은 최대 BH 질량을 $44 M_\odot$로 제한하도록 모델링되었다. 탄생 킥(Natal kicks)은 로그 정규 분포에서 추출되며, 낙하(fallback)에 의해 억제된다.
   • 질량 이동: 기준 모델은 항성 흡수체에 대해서는 보존적 질량 이동($\beta=1$)을, 컴팩트 천체 흡수체에 대해서는 에딩턴 한계 흡정($\beta \approx 0$)을 가정한다.
   • 우주론: cosmerge를 통해 IllustrisTNG 시뮬레이션에 기반한 별 형성사(SFH)를 적용하여 $10^6$개의 병합 시스템 카탈로그를 생성하고, 허블 시간 내에 병합되는 SMT-only 시스템을 필터링하였다.

2. 해석적 프레임워크: 저자들은 흡수 효율($\beta_1$), 임계 질량비($q_{crit}$), 그리고 초기 질량비($q_{ZAMS}$)를 결과적인 최소 BH 질량과 연결하기 위해 기존의 해석적 모델(van Son et al. 2022)을 수정 및 적용한다. 수정 사항에는 더 나은 cosmic 시뮬레이션 일치를 위한 가변적 핵 질량 분율($f_{core}$)과 특정 초신성 질량 손실 함수가 포함된다.

주요 기여

• 아집단 식별: 본 논문은 SMT-only 채널이 자연스럽게 두 가지 관측적으로 구별되는 아집단인 질량비 역전(Mass Ratio Reversed, MRR) 집단과 비-MRR(non-MRR) 집단을 생성함을 입증한다.
• 분리 메커니즘: 저자들은 이러한 아집단의 경계를 다양한 진화 단계(Case A, B, C)에서의 질량 이동 안정성 기준으로 직접 매핑한다. 이들은 서로 다른 공여자 유형에 대한 안정성 기준이 일차 질량-질량비 평면에서 상호 보완적인 영역을 형성함을 보여준다.
• 매개변수 민감도: 본 연구는 흡수 효율($\beta_1$)이 병합율과 질량 분포에 미치는 영향을 체계적으로 탐구함으로써, 이 불확실한 매개변수를 변화시키는 것이 결과적인 BBH 집단을 어떻게 변화시키는지 이해하기 위한 프레임워크를 제공한다.

결과

1. 질량 분포의 뚜렷한 특징:

   • MRR 아집단: 높은 일차 질량($20-40 M_\odot$)과 높은 질량비($q_{BBH} > 0.75$, 0.9 근처에서 피크)를 특징으로 한다. 이 시스템들은 주로 첫 번째 질량 이동(MT1) 동안 Case A 질량 이동을 겪는다.
   • 비-MRR 아집단: 낮은 일차 질량(약 $10 M_\odot$에서 피크)과 더 넓은 질량비 분포($q_{BBH} \sim 0.3 - 0.7$)를 특징으로 한다. 이 시스템들은 주로 MT1 동안 Case B 또는 C 질량 이동을 통해 형성된다.
   • 분리: 기준 모델에서 약 $20 M_\odot$의 일차 질량 컷오프가 이 두 아집단을 효과적으로 분리한다.

2. 흡수 효율($\beta_1$)의 역할:

   • $\beta_1$을 증가시키면(더 보존적인 질량 이동) 안정적 진화를 위한 매개 변수 공간이 축소되어 총 병합율이 감소한다.
   • 높은 $\beta_1$은 MRR 시스템의 분율을 높이고, 두 아집단의 일차 질량 분포를 더 높은 질량 쪽으로 이동시킨다.
   • $\beta_1$이 증가함에 따라 MRR 아집단은 더 두드러지고 더 높은 질량까지 확장되는 반면, 비-MRR 피크는 $10 M_\odot$ 근처에 머물지만 폭이 좁아진다.

3. 경계로서의 안정성 기준:

   • MRR과 비-MRR 집단 사이의 전이는 초기 질량비($q_{ZAMS}$)와 특정 질량 이동 케이스에 대한 안정성 한계($q_{crit}$) 사이의 상호작용에 의해 결정된다.
   • Case A 질량 이동은 안정성을 유지하기 위해 높은 $q_{ZAMS}$와 높은 초기 질량을 요구하는 가장 엄격한 안정성 기준을 가지며, 따라서 고질량 MRR 피크를 형성한다.
   • Case C 질량 이동은 더 낮은 $q_{ZAMS}$와 낮은 질량을 허용하여, 저질량 비-MRR 피크에 기여한다.

의의 및 주장
저자들은 안정적 질량 이동의 물리학이 뚜렷하고 구별 가능한 특징을 중력파 집단에 생성한다고 주장한다. 구체적으로, SMT-only 채널은 일차 질량 분포의 이봉성(bimodality)과 질량비의 이분법을 예측하며, 이는 현재 및 가까운 미래의 GW 검출기 성능과 일치한다. 본 연구는 GWTC-5 카탈로그에서 관측된 특징들(예: 저질량 피크 및 고질량-고질량비 시스템)이 이러한 뚜렷한 SMT 아집단의 징후일 수 있음을 시사한다.

본 논문은 현재의 모델이 퇴행성으로 인해 관측된 사건을 특정 채널에 확정적으로 할당할 수는 없지만, 식별된 특징들(약 $20 M_\odot$에서의 분리와 고-$q$ MRR 집단)이 테스트 가능한 예측으로서 기능한다는 점을 겸허히 결론짓는다. 향-후 관측은 흡수 효율 및 안정성 기준을 포함한 SMT의 물리학을 제약함으로써 거대 이중성 진화에 대한 이해를 높이는 데 도움이 될 수 있다. 저자들은 항성 반경 처방, 풍 손실, 고정된 안정성 기준의 단순화와 관련된 한계점을 인정하며, 이러한 요인들이 예측된 특징의 정확한 위치를 이동시킬 수 있음을 명시하였다.