Stellar black hole binaries from two common envelope evolution phases in triple stellar systems

기술 요약: 삼중 성계 내 두 차례의 공통 외층 진화 단계를 거친 항성 블랙홀 쌍성계

문제 제기
중력파 검출기(LIGO-Virgo-KAGRA)는 약 300개의 블랙홀 쌍성(BBH) 병합 이벤트를 관측했다. 핵심적인 관측적 특징은 유효 인스파이럴 스핀 파라미터, $\chi_{\text{eff}}$의 분포이다. 데이터는 비대칭적인 분포를 보여준다: 평균값은 약간 양수($\langle \chi_{\text{eff}} \rangle \approx 0.054$)이지만, 음의 $\chi_{\text{eff}}$를 가진 계들의 유의미한 꼬리(tail)가 존재한다.
기존의 형성 채널들은 고립된 상태에서 이러한 비대칭성을 재현하는 데 어려움을 겪는다:

• **고립된 쌍성 진화(Isolated binary evolution)**는 일반적으로 $\chi_{\text{eff}}$가 0 근처에 좁게 집중된 분포를 예측한다.
• 순수 역학적 조립(Pure dynamical assembly) (예: 구상 성단 또는 AGN 디스크 내)은 0을 중심으로 대칭인 분포를 예측하며, 관측된 양수의 평균값을 설명하지 못한다.
• 하이브리드 시나리오들이 존재하지만, 음의 꼬리와 양수의 평균을 동시에 설명하는 구체적인 기원은 여전히 미해결 과제로 남아 있다.

방법론
저자들은 핵붕괴 초신성(CCSNe)을 위한 **지터링 제트 폭발 메커니즘(Jittering-Jets Explosion Mechanism, JJEM)**의 틀 안에서 특정 진화 채널을 제안하고 모델링한다. 이 시나리오는 계층적 삼중 항성계 내에서의 두 가지 뚜렷한 공통 외층 진화(Common Envelope Evolution, CEE) 단계를 포함한다.

1. 진화 시나리오:

   • 초기 상태: 질량이 큰 주성($M_1$), 질량이 큰 부성($M_2$), 그리고 저질량 삼중성($M_3 \lesssim 5 M_\odot$)으로 구성된 계층적 삼중 항성계.
   • 제1차 CEE: 주성이 주계열 단계에서 벗어남에 따라 삼중성을 집어삼킨다. 삼중성은 주성의 핵으로 나선형으로 진입하며, 이를 통해 핵을 조석 파괴하고 스핀을 가속시킨다. 이러한 빠른 회전은 JJEM을 통해 성공적인 CCSN을 거쳐 첫 번째 블랙홀(BH1)을 형성하는 데 필수적이다. BH1의 스핀($\vec{\chi}_1$)은 내부 궤도 각운동량($\vec{L}_{\text{orb,in}}$)과 정렬되는데, 이 내부 궤도는 외부 궤도와 경사져 있다.
   • 제2차 CEE: 부성이 적색 초거성(RSG)으로 진화하며 BH1을 집어삼킨다. BH1은 내부로 나선형 진입하여 부성의 핵을 회전시키며 스핀을 가속시킨다. 이후 부성은 붕괴하여 두 번째 블랙홀(BH2)을 형성한다. BH2의 스핀($\vec{\chi}_2$)은 쌍성 궤도 각운동량($\vec{L}_{\text{orb}}$)과 정렬된다.
   • 결과: 최종 BBH 계는 $\vec{\chi}_2$는 궤도와 정렬되어 있지만, $\vec{\chi}_1$은 초기 내부 궤도와 외부 궤도 사이의 경사로 인해 크게 어긋나 있거나(misaligned) 심지어 반정렬(anti-aligned)될 수 있는 스핀 분포를 갖는다. 두 번째 CEE 동안의 강착(accretion)이 $\vec{\chi}_1$을 부분적으로 정렬시킬 수 있으나, 상당한 불일치가 남아 있어 음의 $\chi_{\text{eff}}$ 꼬리를 생성할 수 있다.

2. 수치적 접근법:

   • 저자들은 해당 채널의 생존 가능성을 추정하기 위해 전체 인구 합성(population synthesis) 대신 **반-해석적 접근법(semi-analytical approach)**을 사용한다.
   • 추출 확률 ($P_{\text{draw}}$): 질량, 주기, 질량비, 이심률을 포함하여 경험적으로 보정된 질량 분포(Moe & Di Stefano 2017)를 사용하여 계산한다. 시스템이 진화 요구 사항(예: $M_1 \geq 25 M_\odot$, 특정 궤도 주기 범위)을 충족하도록 제약 조건이 적용된다.
   • 성공 분율 ($\beta$): ZAMS(주계열 전 단계) 조건을 만족하는 시스템이 전체 시퀀스(두 번의 CEE, 두 번의 CCSN, 그리고 최종적인 좁은 궤도)를 성공적으로 거쳐 허블 시간 내에 병합될 조건부 확률로 정의된다.
   • 율(Rate) 정규화: 계산된 병합율은 관측된 국지적 CCSN율($R_{\text{CCSN}} \approx 6.9 \times 10^4 \text{ Gpc}^{-3} \text{ yr}^{-1}$) 및 관측된 BBH 병합율($R_{\text{BBH}} \approx 17.4 - 31.3 \text{ Gpc}^{-3} \text{ yr}^{-1}$)에 대해 정규화된다.
   • 안정성 검증: 역학적 안정성은 Aarseth & Mardling (2001) 기준을 사용하여 강제된다. 또한, 저자들은 제안된 진화 시퀀스가 역학적으로 실현 가능한지 확인하고 허용된 내부 궤도 주기 범위를 정밀화하기 위해 Multiple Stellar Evolution (MSE) 코드를 사용한다.

주요 결과

• 파라미터 제약: 생존 가능한 파라미터 공간은 $M_1 \in [25, 150] M_\odot$, $M_2 \in [18, M_1] M_\odot$, $M_3 \in [1.25, 4.5] M_\odot$, $\log(P_{\text{in}}/\text{day}) \in [0.7, 1.5]$, $\log(P_{\text{out}}/\text{day}) \in [2.3, 3.5]$를 요구한다.
• 추출 확률: 의례적인(fiducial) 삼중성 파라미터에 대해, 이러한 초기 특성을 가진 시스템을 추출할 확률은 $P_{\text{draw}}^{\text{triple}} \approx 1.5 \times 10^{-3}$이다.
• 성공 분율 ($\beta$):
  • 의례적 제약($M_1 \geq 25 M_\odot$) 하에서, 관측된 BBH 율을 맞추기 위해 필요한 성공 분율은 $\beta_{\text{triple}} \lesssim 0.67$이다.
  • 주성 질량 제약을 이론적 BH 형성 하한선($M_1 \geq 20 M_\odot$)으로 완화하면, 필요한 성공 분율은 $\beta_{\text{triple}} \approx 0.47$로 떨어진다.
• 쌍성 채널과의 비교: 삼중성 채널의 성공 분율은 그와 유사한 고립 쌍성 채널($\beta_{\text{binary}} \approx 0.023$)보다 현저히 높지만, 삼중성 채널은 "다중성 페널티(multiplicity penalty)"(특정한 파라미터를 가진 안정적인 삼중성을 형성할 낮은 확률)로 인해 더 희귀하다.

의의 및 주장
본 논문은 이 삼중성 진화 채널이 병합하는 BBH 인구의 생존 가능하며 무시할 수 없는 기여자라고 주장한다.

• 질적 적합성: 이 메커니즘은 관측된 비대칭적 $\chi_{\text{eff}}$ 분포를 자연스럽게 설명한다. 삼중계의 기하학적 구조로 인한 첫 번째 탄생 BH의 스핀 불일치는 음의 꼬리를 생성하며, 두 번째 탄생 BH의 스핀 정렬은 양의 평균에 기여한다.
• 율의 생존 가능성: 계산된 성공 분율($\beta < 1$)은 이 채널이 율 제약에 의해 배제되지 않음을 나타낸다. 저자들은 이 채널이 이론적으로 전체 BBH 인구를 만들어낼 수 있음에도 불구하고(만약 $\beta$가 1에 가깝다면), 다른 채널들과 함께 중요한 기여를 할 가능성이 높다고 명시한다.
• 한계점: 저자들은 이 채널이 BBH 인구를 지배한다고 주장하지 않는다. 현재의 연구가 정량적인 $\chi_{\text{eff}}$ 분포를 도출하거나 이 채널이 공급하는 인구의 정확한 비율을 결정하지는 않는다는 점을 언급한다. 높은 요구 성공 분율($\sim 0.5-0.7$)은 이 채널이 타당할 수는 있지만, 유일한 형성 메커니즘은 아닐 수 있음을 시사한다.

요약하자면, 본 연구는 관측된 BBH 스핀의 양수 평균과 음의 꼬리 사이의 긴장을 해결하는 물리적으로 동기 부여된 삼중성 경로를 제안하며, 이러한 계들이 현재의 중력파 관측과 일치하는 비율로 형성될 수 있음을 입증한다.