스핀 효과가 이끄는 중성자별 병합의 새로운 모습

초록

본 연구는 SFHo 방정식 상태와 중성자별의 초기 스핀을 다양하게 설정한 3+1 일반 상대론적 유체역학 시뮬레이션을 통해, 동등·비동등 질량 이진 중성자별 병합에서 스핀이 동역학, 잔류물 구조, 물질 방출, 중력파 및 중성미자 방출에 미치는 영향을 체계적으로 조사한다. 특히, 고스핀(χ=0.67) 경우에 동등 질량 시스템에서 최대 0.06 M⊙의 물질이 방출되고, 고질량 시스템에서는 차원 없는 블랙홀 스핀이 0.92에 달함을 보고한다.

상세 분석

이 논문은 최신의 유한 온도, 조성 의존적 SFHo 방정식 상태(EoS)를 사용하고, 중성미자 방출·흡수를 포함한 M0+Leakage 모델을 적용한 3+1 차원 일반 상대론적 유체역학(GRHD) 시뮬레이션을 수행하였다. 초기 데이터는 Kadath 라이브러리의 Fuka 브랜치를 이용해 다양한 스핀 구성(정렬, 반정렬, 혼합)과 질량비(q≤2.05)를 갖는 12개의 모델을 생성했으며, 전체 질량은 2.55 M⊙, 3.05 M⊙, 4.10 M⊙ 세 가지 경우로 나누었다. 시뮬레이션은 WhiskyTHC 코드와 CTGamma의 Z4c 제약 억제 포뮬레이션을 결합해 높은 정확도의 중력파 위상과 진폭을 확보했으며, 7단계 AMR을 통해 최세밀 격자 해상도는 222 m(LR은 308 m)까지 도달하였다.

동등 질량(1.27 M⊙)의 무스핀 모델은 장시간 지속되는 초중성자별(remnant NS)을 형성했지만, 질량을 1.52 M⊙와 2.05 M⊙로 증가시키면 즉시 블랙홀로 붕괴하였다. 스핀을 정렬(χ>0)으로 증가시키면 궤도 이탈(hang‑up) 효과가 나타나 병합 전 단계에서 인스피라션 시간이 연장되고, 결과적으로 블랙홀 형성 시점이 지연된다. 반정렬(χ<0) 경우에는 반대로 병합이 가속화되어 조기 붕괴가 관찰된다. 특히 χ=0.67인 고스핀 동등 질량 모델에서는 동역학적 물질 방출량이 약 0.06 M⊙에 달했으며, 이는 기존 연구에서 보고된 최대값보다 두 배 이상 큰 수치이다. 방출 물질의 전자분율(Y_e)은 스핀 정렬에 따라 약간 상승하여 r‑process 핵합성에 유리한 조건을 제공한다.

블랙홀 형성 모델에서는 최종 블랙홀의 차원 없는 스핀 파라미터가 χ≈0.92에 이르렀으며, 이는 현재까지 보고된 BNS 병합 블랙홀 스핀 중 최고값이다. 이는 초기 스핀이 높은 경우 각운동량이 효율적으로 보존되어 블랙홀에 전달되기 때문이다. 또한, 디스크 질량은 스핀 정렬 모델에서 평균 0.12 M⊙, 반정렬 모델에서는 0.07 M⊙ 수준으로 차이를 보였으며, 디스크 바람(ejecta) 역시 스핀에 민감하게 반응한다. 중성미자 방출 측면에서는 정렬 스핀 모델이 평균 에너지와 광도를 약 15 % 상승시켰으며, 이는 디스크와 바람의 온도 상승에 기인한다.

중력파 스펙트럼 분석 결과, 기본 (l,m)=(2,2) 모드의 피크 주파수 f₂는 스핀 정렬 시 낮아지고 반정렬 시 높아지는 경향을 보였다. f₁ 모드(1,2) 역시 스핀에 따라 약 100 Hz 정도 이동했으며, 이는 관측 가능한 후방 파형 변조에 기여한다. 전체적으로 스핀은 질량 비율, 방출 물질 양·조성, 디스크·블랙홀 스핀, 중력파 및 중성미자 신호에 복합적인 영향을 미치며, 특히 고스핀 정렬 경우는 관측 가능한 전자기·중성미자 신호를 크게 강화한다는 점이 강조된다.