중성자별 병합 디스크 유출에 대한 r‑프로세스 가열 피드백 연구

초록

**
본 논문은 중성자별 병합 후 형성되는 디스크에서 발생하는 r‑프로세스 핵반응 가열을 수치유체역학에 직접 결합하는 새로운 방법을 제시한다. 온도와 전자분율 (Y_e) 이력에 기반한 가열율을 트레이서 입자와 Cloud‑in‑Cell 기법으로 격자에 전달한다. 장기 점성 수치 시뮬레이션에 적용한 결과, r‑프로세스 가열이 무결합 디스크 유출 질량을 약 10 % 증가시키고, (Y_e<0.25)인 물질의 방사속도를 최대 두 배까지 끌어올리며, 경계선에 가까운 대류 유출을 억제한다는 것을 확인하였다.

**

상세 분석

**
이 연구는 기존의 후병합 시뮬레이션이 알파 입자 재결합에 의한 핵에너지 방출만을 고려하고, 비평형 r‑프로세스에 의한 추가 가열을 무시한다는 한계를 극복하고자 한다. 저자들은 Lippuner & Roberts(2015)의 핵반응망 결과를 이용해 온도와 초기 전자분율 (Y_{e,,T6}) (6 GK에서의 (Y_e))에 따라 달라지는 가열율 테이블을 구축하였다. 핵에너지의 50 %는 중성미자로 탈출하고 나머지는 유체 내부에 열로 전환된다고 가정한다.

핵가열 정보를 격자에 전달하기 위해 두 종류의 트레이서를 사용한다. ‘일반 트레이서’는 유동을 추적하고, ‘메모리 트레이서’는 각 격자 셀을 충분히 샘플링하도록 9.9 × 10⁵개를 배치해 (Y_{e,,T6}) 이력을 기록한다. 트레이서는 매 시간 단계마다 현재 온도와 (Y_e) 값을 받아 (T_{\rm heat}=4) GK 이하가 되면 가열을 적용한다(입구 흐름에는 음의 부호를 주어 냉각 효과를 부여). 가열율은 (1) (Y_{e,,T6}) 가 가장 가까운 두 테이블 값 사이를 선형 보간하고, (2) 각 (Y_{e,,T6}) 에 대해 로그 스케일로 온도 보간을 수행한 뒤, 최종적으로 두 보간값을 다시 선형 보간한다.

수치 실험은 2.65 M⊙, 스핀 (a=0.8) 인 블랙홀과 0.1 M⊙ 질량의 등각동량 토러스(초기 (Y_e=0.1), 엔트로피 8 k_B / 바리온)를 초기조건으로 설정하고, α‑디스크 모델(α=0.03, 0.06)과 Helmholtz EOS, 3‑종류 누출 중성미자 스키마를 적용하였다. 결과는 다음과 같다.

1. 질량 증가: r‑프로세스 가열을 포함한 경우 무결합 디스크 유출 질량이 기준(알파 재결합만) 대비 약 10 % 증가한다. 이는 가열이 디스크 외부 압력을 상승시켜 더 많은 물질이 탈출하도록 만든다.
2. 속도 변화: (Y_e<0.25)인 저‑(Y_e) 물질은 방사속도가 최대 두 배까지 증가한다. 고‑(Y_e) 물질은 상대적으로 변화가 적으며, 이는 가열율이 저‑(Y_e)에서 더 오래 지속되기 때문이다.
3. 대류 억제: 경계선에 위치한 약간 결합된 대류 유출은 가열에 의해 내부 압력이 상승하면서 재결합되며, 결과적으로 전체 유출의 비대류 성분 비율이 높아진다.
4. 임계 온도 민감도: (T_{\rm heat}=6) GK로 설정하면 가열 적용 시점이 앞당겨져 질량 증가와 속도 상승 효과가 약 30 % 정도 감소한다. 이는 r‑프로세스가 실제로 시작되는 온도 구간에 대한 불확실성을 반영한다.

저자들은 또한 기존 연구와 비교한다. Just et al.(2015)와 Wu et al.(2016)에서 온도 의존성만을 고려한 가열은 결과가 크게 달랐으며, Klion et al.(2022)와 Kawaguchi et al.(2021)은 시간‑의존적 파라메트릭 가열을 사용했지만, 구성 의존성을 무시했다. 본 연구는 (Y_e) 이력을 직접 추적함으로써 이러한 한계를 보완하고, 가열이 물질의 궤적에 미치는 비선형 효과를 정량화한다.

한계점으로는 (i) 핵가열 테이블이 고정된 엔트로피와 팽창시간(30 k_B, 10 ms)으로만 생성되었으며, 실제 디스크 흐름의 다양성을 완전히 포착하지 못한다는 점, (ii) 중성미자 탈출 비율을 50 %로 고정한 단순화, (iii) 2‑차원 축대칭 시뮬레이션이 3‑차원 마그네틱 터뷸런스와 비대칭 효과를 놓친다. 향후 연구에서는 전자분율‑의존적 가열 테이블을 다중 파라메터(엔트로피, 팽창시간)로 확장하고, 3‑D MHD와 연계한 실시간 핵반응망 결합을 시도할 필요가 있다.

**