On the effective spin-mass ratio relation of binary black hole mergers that evolved in isolation

기술 요약: 고립 진화된 쌍성 블랙홀 병합의 유효 스핀-질량 비 관계에 대하여

문제 제기
LIGO-Virgo-Kagra(LVK) 네트워크에 의한 지상 중력파(GW) 검출 결과, 관측된 집단에서 두 가지 흥미로운 특이점이 드러났다: 유효 스핀 파라미터($\chi_{\text{eff}}$) 분포의 양의 편향(positive bias)과 병합 질량비($q$)와 $\chi_{\text{eff}}$ 사이의 명백한 반비례 관계(anti-correlation)이다. 역학적 형성 채널은 등방성 스핀 분포(0을 중심으로 대칭)를 예측하는 반면, 관찰된 양의 편향은 고립 쌍성 진화의 기여를 시사한다. 그러나 $\chi_{\text{eff}}-q$ 반비례 관계의 기원에 대해서는 논쟁이 계속되고 있으며, 최근의 일부 데이터는 오히려 양의 상관관계(positive correlation)를 시사하기도 한다. 본 연구는 고립된 질량 큰 쌍성 진화, 구체적으로 안정적 질량 전이(SMT) 및 공통 외층(CE) 채널이 AGN과 같은 이례적인 환경을 도입하지 않고도 어떻게 자연스럽게 관측된 스핀 및 질량 특성을 재현할 수 있는지 조사한다.

방법론
저자들은 고립된 질량 큰 쌍성의 인구 통계적 합성을 수행하기 위해 업데이트된 버전의 BSE(Hurley et al. 2002) 신속 이진 진화 코드를 활용한다. 주요 방법론적 업데이트는 다음과 같다:

• BH 스핀 구현: 코드는 Wolf-Rayet(WR) 전성기(progenitor)의 조석 스핀업(tidal spin-up)에 기반한 두 번째 탄생 블랙홀($a_{\text{BH,sec}}$)의 출생 스핀 처방을 포함한다. 두 가지 해석적 모델이 테스트된다: 조밀한 WR-BH 쌍성에서 강체 회전과 완전한 조석 동기화를 가정하는 Be20 모델(Belczynski et al. 2020), 그리고 차동 회전(differential rotation)을 고려하며 WR 질량과 궤도 주기에 모두 의존하는 Ba21 모델(Bavera et al. 2021)이다.
• 질량 전이 물리: SMT 채널을 탐구하기 위해, 공통 외층(CE) 발달을 위한 임계 질량비($q_{\text{crit}}$)를 인위적으로 높은 값($q_{\text{crit}} = 8.0$)으로 설정하여 CE 진화를 억제하고 SMT 하위 채널을 분리한다. 이는 CE 진화를 선호하는 기본 BSE 설정과 대조적이다.
• 물리적 가정: 연구는 금속함량($Z = 0.0001$ ~ $0.02$), 잔존 질량 모델(Fryer et al. 2012의 "rapid" vs "delayed" 엔진), 그리고 CE 진화에서의 HG-별 기증자(donor) 처리를 변화시킨다.
• 인구 통계적 합성: 각 금속함량 및 모델 구성에 대해 $2 \times 10^6$ 개의 쌍성을 진화시킨다. 생성된 병합 집단은 우주 별 형성사 및 금속함량 진화에 의해 가중치가 부여되어 적색편이 $z \approx 0$에서의 현재 관측 가능한 분포를 도출한다.
• 섭동(Perturbations): 또한 역학적으로 조립된 BBH(등방성 스핀)의 혼합 효과와 "BH 스핀 토싱(BH spin tossing, 출생 시 킥으로 인한 전성기 별의 스핀과 BH 스핀 사이의 미정렬)"의 영향을 모델링한다.

주요 기여 및 결과

1. 주 질량 피크: SMT 채널은 일관되게 약 $\approx 10 M_\odot$ (일부 구성에서는 $\approx 9 M_\odot$)에서 주 블랙홀 질량 분포의 두드한 피크를 생성한다. 이는 LVK 데이터에서 관찰된 피크와 일치하며, SMT가 검출된 사건들의 중요한 형성 채널임을 시사한다.
2. $\chi_{\text{eff}}-q$ 반비례 관계: 본 연구는 특히 질량비 역전(mass-ratio-reversed) 병합(두 번째 탄생 BH가 첫 번째 탄생 BH보다 질량이 큰 경우, $M_{\text{second}} > M_{\text{first}}$) 내에서 $\chi_{\text{eff}}$와 $q$ 사이의 뚜렷한 반비례 관계를 식별하였다.
   • 메커니즘: 이 반비례 관계는 질량 비대칭성이 높은(낮은 $q$) 계가 더 광범적인 질량 전이를 겪어 궤도 수축이 더 크게 일어나기 때문에 발생한다. 이는 더 조밀한 WR-BH 쌍성을 생성하여 궤도 주기를 단축시키고, 결과적으로 WR 별의 조석 스핀업을 강화하여 더 높은 스핀을 가진 더 무거운 두 번째 탄생 BH를 생성한다. 결과적으로 이 하위 집단 내의 낮은 $q$ 병합은 더 높은 $\chi_{\text{eff}}$를 나타낸다.
   • 채널 의존성: 반비례 관계는 "rapid" 잔존 질량 모델과 낮은 금속함량을 사용할 때, 특히 SMT 채널 병합에서 가장 두드러진다. 반면, 보완적인 하위 집단($M_{\text{second}} \le M_{\text{first}}$)은 종종 양의 상관관계나 명확한 경향성을 보이지 않는다.
3. 금속함량 의존성: 강한 항성풍이 쌍성 궤도를 넓히고 조석 동기화를 감소시키기 때문에, 질량비 역전 시스템의 비율($f_{\text{rev}}$)과 두 번째 탄생 BH의 평균 스핀은 금속함량이 증가함에 따라 감소한다. 결과적으로 높은 금속함량에서는 $\chi_{\text{eff}}-q$ 반비례 관계가 사라진다.
4. 강건성 및 한계:
   • 역학적 혼합: 등방성 스핀을 가진 역학적 조립 BBH의 비율을 도입하면 반정렬 스핀($\chi_{\text{eff}} < 0$)의 집단이 증가하지만, 뚜렷한 $\chi_{\text{eff}}-q$ 반비례 관계를 지우는 경향이 있다. 전체 집단에서 반비례 관계는 역학적 비율이 $\lesssim 50\%$일 때만 유지된다.
   • 스핀 토싱: "스핀 토싱"(BH 스핀과 전성기 별 스핀 사이의 미정렬)을 허용하면 반비례 관계의 강도가 약해진다. LVK 카탈로그의 $\chi_{\text{eff}} < 0$ 사건 비율을 맞추기 위해서는 상당한 수준의 스핀 토싱($\omega_{b,\text{max}} > 0.5\pi$)이 필요하며, 이는 반비례 신호를 더욱 약화시킨다.
   • 잔존 모델: "delayed" 잔존 질량 모델은 "rapid" 모델에 비해 더 많은 저질량 BH와 더 적은 질량비 역전 시스템을 생성하여, 전체 집단에서 전반적인 반비례 관계를 약화시킨다.

의의 및 주장
본 논문은 SMT 하위 채널을 통한 고립 쌍성 진화가 관측된 BBH 병합 집단의 특징인 양의 $\chi_{\text{eff}}$ 편향과 특징적인 $\chi_{\text{eff}}-q$ 반비례 관계를 재현할 수 있는 실행 가능한 메커니즘이라고 주장한다. 저자들은 이 반비례 관계가 질량비에 따른 궤도 수축 및 로슈 로브(Roche lobe) 크기의 의존성에 의해 구동되는, 질량비 역전 쌍성의 물리에 기인한 자연스러운 결과라고 논한다.

연구는 관찰된 특징들($\sim 10 M_\odot$ 질량 피크 및 $\chi_{\text{eff}}-q$ 반비례 관계)이 질량비 역전 시스템을 포함하는 고립 쌍성 진화의 징표가 될 수 있다고 결론짓는다. 그러나 저자들은 관측 데이터에서 반비례 관계의 존재 여부가 여전히 논쟁 중이며(GWTC-4는 양의 상관관계를 암시함), 자신들의 결과가 특정 물리적 가정(질량 전이 안정성 및 잔존 질량 모델)에 의존한다는 점을 언급하며 신중한 태도를 유지한다. 그들은 모델이 추세를 재현하기는 하지만, 스핀 토싱의 영향과 정확한 질량 전이 안정성 기준을 고려할 때 SMT 대 CE 또는 역학적 채널의 정확한 기여도를 규명하기 위해서는 추가적인 조사가 필요함을 강조한다.