전이성 백색왜성 중성자별 병합 전일반상대론 시뮬레이션을 통한 회전 TZlO 형성 연구

초록

본 논문은 전일반상대론(GR) 시뮬레이션을 이용해 백색왜성‑중성자별( WD‑NS ) 이진계의 인스파이럴 및 병합 과정을 조사한다. 실제 백색왜성 대신 크기를 축소한 ‘의사백색왜성(pWD)’ 모델을 사용해 0.98 M⊙ WD와 1.4 M⊙ NS의 병합을 모사하였다. 결과는 회전하는 Thorne‑Zytkow‑like Object(TZlO)와 거대한 디스크가 형성되며, 총 질량이 차가운 EOS가 지탱할 수 있는 최대 질량을 초과함에도 불구하고 즉각적인 블랙홀 붕괴는 일어나지 않는다. 열냉각을 적용해도 회전 지지에 의해 붕괴가 억제되며, 장기적으로는 점성·자기장에 의한 각운동량 재분배가 지연 붕괴를 일으킬 것으로 예상한다.

상세 분석

이 연구는 WD‑NS 이진계의 물리적 규모 차이가 10⁴ ~ 10⁶에 달해 전일반상대론 수치 시뮬레이션이 실질적으로 불가능한 점을 극복하기 위해 ‘의사백색왜성(pWD)’이라는 스케일링 기법을 도입한 것이 핵심이다. pWD는 실제 WD와 동일한 질량을 유지하면서 반경을 인위적으로 축소시켜, 압축비(R_NS : R_pWD)를 5:1 ~ 20:1 수준으로 조정한다. 이렇게 하면 충돌 속도 v_c와 음속 c_s의 비인 마하 수가 보존되며, 충돌 시 발생하는 충격 가열과 회전·중력 에너지 비율(T/|W|)도 동일하게 유지된다. 따라서 pWD 기반 시뮬레이션 결과를 실제 WD에 직접 스케일링할 수 있다.

EOS는 6‑parameter 조각다항식(polytropic) 형태로 설계되어, 실제 핵물리 EOS가 보여주는 질량‑중심밀도 관계와 불안정 구간을 재현한다. 특히 최대 NS 질량을 1.8 M⊙, 최대 WD 질량을 1.43 M⊙(챈드라세카르 한계)로 맞추어, 두 별의 구조적 특성을 동시에 만족한다. 이 EOS는 압축비 조절을 위해 κ₁, κ₂, κ₃와 지수 n₁, n₂, n₃를 자유롭게 설정할 수 있어, pWD의 반경을 원하는 비율로 줄이면서도 물리적 일관성을 유지한다.

시뮬레이션 코드는 Illinois AMR GRHD 코드를 기반으로 하며, BSSN 포멀리즘으로 중력장을, 고해상도 적응격자망으로 물질 흐름을 해결한다. 초기 데이터는 로체 한계(Roche limit)에서 원형 궤도를 갖는 WD‑NS 이진을 구성하고, 질량비 q ≈ 0.7 > q_crit(≈0.5)인 경우를 선택해 불안정 질량 전달(UMT) 상황을 재현한다.

핵심 결과는 두 단계로 나뉜다. 첫 번째는 기존 연구