Numerical simulations of black hole-neutron star mergers with equal and near-equal mass ratios

기술 요약: 등질량 및 근등질량비를 가진 블랙홀-중성자별 병합의 수치 시뮬레이션

문제 및 동기
질량비 $q \approx 1/2.6$인 블랙홀-중성자별(BHNS) 병합으로 해석되는 중력파 이벤트 GW230529의 검출은, 이전에 예상했던 것보다 더 대칭적인 BHNS 계가 존재할 가능성을 시사한다. 표준 항성 진화 모델은 $q \lesssim 1/3.5$의 질량비를 예측하지만, 밀집 성단에서의 역학적 형성 채널이나 2세대 병합은 $q \sim 1/2$ 또는 심지어 $q \approx 1$인 계를 생성할 수 있다. 이러한 대칭적인 계들은 이론적으로 거대한 강착 원반을 생성하여 킬로노바와 같은 밝은 전자기적 과도 현상을 유도할 것으로 기대된다. 그러나 등질량 및 근등질량비 매개변수 공간은 수치 상대론(NR) 시뮬레이션에 의해 채워진 데이터가 부족한 상태이다. 이러한 데이터의 부재는 중력파 파형 모델과 잔해 특성에 대한 해석적 피팅 공식의 정확도를 제한하며, 이는 실제 사건의 오분류와 전자기적 대조군 식별 기회의 상실로 이어질 수 있다.

방법론
저자들은 BAM 코드를 사용하여 비회전 블랙홀-중성자별 병합에 대한 12개의 NR 시뮬레이션을 수행하였다. 시뮬레이션은 질량비 $q \in \{1, 1/2, 1/3\}$와 중성자별(NS) 질량 $m_{NS} \in \{0.8, 1.2, 1.4, 2.2\} M_\odot$를 다룬다. 두 가지 상태 방정식(EOS)인 연한 SLy와 단단한 DD2가 사용되었다. 시뮬레이션은 파스셜 폴리트로픽(piecewise-polytropic) 표현법을 사용하여 EOS를 나타내며, 공형 Z4 정식화(conformal Z4 formulation)를 통한 시공간과 발렌시아 정식화(Valencia formulation)를 통한 유체역학을 사용하여 약 9번의 궤도 동안 시스템을 진화시킨다. 초기 데이터는 FUKA 코드로 구축되었으며, 이심률은 $10^{-3}$ 미만으로 감소되었다.

본 연구는 다음을 포함한다:

1. 중력파 파형 분석: 여러 반경에서 파형을 추출하고, 2차 외삽을 수행하며, 질량 중심 표류를 보정한다. 이를 IMRPhenomNSBH 및 SEOBNRv4_ROM_NRTidalv2_NSBH 파형 모델과 비교한다.
2. 잔해 특성 검증: 동적 방출물(dynamical ejecta)과 강착 원반의 질량을 계산하고, 병합 후 10ms 시점에서의 잔해 블랙홀(BH) 질량과 스핀을 계산한다. 이들은 [28, 31, 40, 41]의 해석적 피팅 공식과 비교된다.
3. 원반 역학: 유의미한 원반 질량을 가진 구성에 대해 원반의 형태, 회전 법칙 및 전역 밀도 진동을 분석한다.
4. 강착 모델링: 강착 행동을 분류하고 원반 진동에 의한 변조를 식별한다.
5. 킬로노바 모델링: 동적 방출물 프로파일과 추가된 구형 윈드 성분(원반 질량의 30%)을 기반으로 POSSIS 복사 전달 코드를 사용하여 킬로노바 광도 곡선을 모델링한다.

주요 결과

• 중력파 파형: 현재의 BHNS 파형 모델(IMRPhenomNSBH 및 SEOBNRv4)은 대칭적 질량비 영역에서 NR 파형과 낮은 일치도를 보인다. 대부분의 경우 병합 시점에 약 1 라디안의 위상차(dephasing)를 보인다. 모델들은 종종 NR 데이터에 비해 더 이른 병합 시간과 확장된 링다운(ringdown) 단계를 예측한다. 특히 아태양 질량 NS 케이스($m_{NS}=0.8 M_\odot$)의 경우, 타이달 변형성(tidal deformability)이 기존 모델의 유효 매개변수 공간을 초과하며, 외삽 시 모델들은 상당한 위상 및 진폭 불일치를 보인다.
• 잔해 특성:
  • 방출물 및 원반 질량: 동적 방출물 질량과 원반 질량에 대한 해석적 피팅 공식은 시뮬레이션 데이터와 잘 일치한다.
  • 잔해 BH 질량: 해석적 공식은 잔해 BH 질량을 매우 정확하게 예측한다(오차 일반적으로 < 1%).
  • 잔해 BH 스핀: 업데이트된 피팅 공식은 잔해 스핀을 예측하는 데 있어 덜 정확하며, 특히 타이달 변형성이 큰 시스템(DD2 EOS)에서 그러하다. 아태양 NS 케이스의 경우, 스핀 예측이 거의 50% 편차를 보인다. 저자들은 이것이 낮은 타이달 변형성과 더 비대칭적인 질량비를 가진 시스템에 과적합(overfitting)되었기 때문일 수 있다고 제안한다.
• 원반 구조 및 역학:
  • 형태: 등질량($q=1$) 원반는 병합 후 거의 즉시($\sim 3$ ms) 원형화되는 반면, $q=1/2$ 원반는 나선형 밀도파와 낙하 물질(fallback matter)에 의해 더 오래 교란되며, $\sim 8$ ms 후에야 축대칭으로 안정화된다.
  • 회전 법칙: 비특이적 각운동량은 $j \propto \Omega^{-n}$의 멱법칙을 따른다. 지수 $n$은 반경 방향으로 변화하며, ISCO 내부의 일정 각운동량($n=0$)에서 외부 원반의 아케플러(sub-Keplerian) 값으로 전이된다. 이는 장기 축대칭 시뮬레이션에서의 고정된 회전 법칙 가정과 모순된다.
  • 진동: 질량이 $10^{-2} M_\odot$보다 큰 원반에서 전역 포획된 $g$-모드 밀도 진동이 식별된다. 그 주파수는 잔해 BH의 최대 전행(prograde) 에피사이클릭 주파수와 일치한다.
  • 강착: 세 가지 강착 유형이 구분된다. Type B ($q=1/2$)는 전역 $g$-모드 진동에 의해 직접 조절되는 강착률 가변성을 보인다. Type A ($q=1/2$)와 Type C (직접 추락)는 다른 거동을 보인다. 일부 시스템의 $g$-모드 주파수는 감마선 폭발에서 관측되는 고주파 준주기 진동(QPO)과 호환 가능하다.
• 킬로노바 방출: 모델링된 킬로노바는 주로 원반 윈드 방출물에 의해 지배된다. 더 단단한 DD2 EOS를 가진 시스템은 더 부드러운 SLy EOS를 가진 시스템보다 더 질량이 큰 원반을 생성하므로 더 밝은 킬로노바를 생성한다. 200 Mpc 거리에서, 대부분의 시스템은 Vera C. Rubin 천문대가 4일 이내에, 그리고 DECam가 $r$ 및 $i$ 밴드에서 2일 이내에 검출할 수 있을 것이다.

의의 및 주장
본 논문은 대칭적 BHNS 병합에 대한 NR 시뮬레이션 매개변수 공간의 중요한 공백을 메운다고 주장한다. 주요 의의는 현재의 중력파 모델이 등질량 및 근등질량비에서 충분히 정확하지 않음을 입증하여, GW230529와 같은 미래의 검출에 대한 매개변수 추정의 편향을 초래할 수 있음을 보여주는 데 있다. 본 연구는 잔해 질량에 대한 해석적 피팅 공식을 검증하는 동시에, 높은 타이달 변형성을 가진 시스템의 잔해 스핀을 예측하는 데 있어서의 한계를 강조한다. 또한, 원반 역학에 대한 상세한 분석은 복잡한 회전 법칙과 진동에 의한 강착 조절을 드러내며, 이는 $g$-모드가 감마선 폭발 광도 곡선에 관측 가능한 흔적을 남길 수 있음을 시사한다. 마지막으로, 킬로노바 모델링은 대칭적 BHNS가 현재의 광시야 탐사 도구들에 의한 전자기적 후속 관측의 유망한 대상임을 나타내며, 이는 기존 관측을 지배해 온 비대칭적 계들에 비해 더 높은 검출 확률을 제공한다.