R-process nucleosynthesis from magnetar giant flares in neutron star--white dwarf mergers: A unified picture for peculiar long gamma-ray bursts

기술 요약: 중성자별-백색왜성 병합에서의 마그네타 거대 플레어로부터의 r-과정 핵합성

문제 제기
GRB 211211A 및 230307A로 대표되는 일련의 "특이한" 긴 감마선 폭발(LGRB)은 표준적인 GRB 분류의 이분법에 도전 과제를 던진다. 이 사건들은 일반적인 코어 붕괴형 LGRB의 특징인 긴 지속 시간($T_{90} > 2$ s)을 보이지만, 동시에 컴팩트 쌍성 병합에서 기원하는 단기 GRB(SGRB)의 특징인 다중 성분 프롬프트 방출(전구체, 주 폭발, 연장 방출)과 X선 플래토(plateau), 그리고 킬로노바 신호를 보유하고 있다. 기존 모델들은 다음과 같은 상당한 난관에 봉착해 있다:

1. 블랙홀 잔해 모델: 블랙홀 엔진을 남기는 NS–BH 또는 NS–NS 병합 모델은 세 가지 뚜의 방출 에피소드(주 폭발[MB], 연장 방출[EE], X선 플래토)와 특정한 MB–EE 사이의 트러프(trough)를 단일 중심 엔진으로 설명하는 데 어려움을 겪는다.
2. 표준 NS–WD 병합 모델: 중성자별-백색왜성(NS–WD) 병합은 초가속(hyperaccretion) 및 프로펠러 단계를 통해 프롬프트 방출 구조를 재현할 수 있으나, 기존의 유체역학 시뮬레이션은 조석 파쇄된 WD 파편으로부터 형성된 디스크가 양성자가 풍부한(proton-rich) 상태임을 시사한다. 이러한 환경은 강력한 r-과정 중원소를 생성하기에 부적합하며, 이는 이 사건들에서 관측된 킬로노바 신호(특히 GRB 230307A에서 확인된 r-과정 원소의 분광학적 증거)를 설명하는 데 결핍을 초래한다.

방법론
저자들은 기존의 마그네타($B_s \gtrsim 10^{14}$ G)와 질량이 큰 백색왜성($M_{WD} \gtrsim 1 M_\odot$)의 병합에 기반한 통합된 그림을 제안한다. 방법론은 다음을 통합한다:

• 가속 역학: 일정한 엔트로피를 가진 강착 디스크(CENW/CEEW 규정)에 대한 유체역학 시뮬레이션(Kaltenborn et al. 2023)을 활용하여 초가속 단계($\dot{M} \gtrsim 10^{-4} M_\odot \text{s}^{-1}$)와 그에 따른 자기 프로펠러 단계로의 전이를 모델링한다.
• 마그네타 물리학: 병합 과정 중 토로이달 자기장($B_\phi$)을 $10^{15} - 5 \times 10^{16}$ G까지 증폭시키는 다이너모 메커니즘을 적용한다. 이는 반복적인 성크웨이크(starquakes)와 마그네타 거대 플레어(GF)로 이어진다.
• 핵합성 메커니즘: Cehula et al. (2024)의 프레임워크를 채택하여, 거대 플레어의 초기 하드 스파이크(hard spikes)가 팽창하는 $e^\pm$-$\gamma$ 파이어볼을 구동한다고 가정한다. 이 파이어볼은 중성자별 지각(crust)을 충격하여 바리온 질량을 방출하며, 높은 전자 분율($Y_e \gtrsim 0.4$)에도 불구하고 $\alpha$-풍부 동결( $\alpha$-rich freeze-out) 메커니즘을 통해 r-과정 원소를 합성한다.
• 관측 모델링: GRB 211211A 및 230307A의 다파장 데이터(X선, 가시광선, 적외선, 라디오)에 대한 결합 MCMC 피팅을 수행한다. 이 모델은 표준 외부 충격 후광 모델과 $^{56}$Ni 대신 r-과정 가열을 사용하는 준해석적 마그네타 강화 킬로노바 정식화(formalism)를 결합한다.

주요 기여

1. 특이한 LGRB에 대한 통합적 기원: 본 논문은 단일 NS–WD 병합 시나리오가 프롬프트 방출 구조(MB, EE, 트러프), X선 플래토, 그리고 킬로노바 신호를 모두 설명할 수 있는 자기 일관적인 프레임워크를 제시한다.
2. r-과정 결핍 문제 해결: 저자들은 NS–WD 병합 내에서 r-과정 핵합성이 일어나는 새로운 채널을 식별하였다. WD 파편 디스크(양성자가 풍려함)에 의존하는 대신, 마그네타 거대 플레어에 의해 구동되는 지각 방출물(crustal ejecta) 내에서 r-과정 원소가 합성된다고 제안한다. 이는 NS–WD 병합 모델과 GRB 230307A의 r-과정 원소 분광학적 증거 사이의 모순을 해결한다.
3. "플레어의 숲(Forest of Flares)": 주 폭발(Main Burst)을 단일 연속적인 제트가 아니라 약 800개에서 2,000개의 연속적인 GF 초기 스파이크들의 집합체로 재해석한다. 이 해석은 관측된 최소 변동성 시간 규모($\tau_{min} \sim 9-16$ ms)를 은하계 및 인근 은하에서 관측된 개별 GF 스파이크의 지속 시간과 일치시킨다.
4. 정량적 제약 조건: 관측된 킬로노바를 구동하는 데 필요한 방출물 질량($M_{ej} \sim 10^{-5} - 10^{-3} M_\odot$), 팽창 속도 및 불투명도에 대한 구체적인 제약 조건을 제공하며, 이는 유체역학 및 핵 반응 네트워크 계산과 일치함을 보여준다.

결과

• 프롬프트 방출: 모델은 주 폭발(MB)을 GF 스파이크들의 중첩으로, 연장 방출(EE)을 프로펠러 단계의 결과로 성공적으로 재현한다. 강착에서 프로펠러로의 전이는 MB와 EE 사이에서 관측되는 트러프를 자연스럽게 설명한다.
• 킬로노바 특성: 결합 모델링 결과, 빔 효과를 고려했을 때 $M_{ej} \sim 10^{-2} M_\odot$의 방출물 질량을 산출하였으며, 이는 시뮬레이션과 일치하는 속도를 보인다. GRB 230307A의 경우, 높은 불투명도($\kappa \sim 4.28 \text{ cm}^2 \text{ g}^{-1}$)는 두 번째 및 세 번째 피크 r-과정 핵종의 합성에 의해 물리적으로 정당화된다.
• 마그네타 매개변수: 추론된 병합 후 마그네타 표면 자기장($B_{s,post} \sim 4 - 7 \times 10^{14}$ G)과 초기 회전 주기는 병합 중 토로이달 자기장의 증폭에 따른 예상치와 일치한다.
• 에너지 예산: GF 스파이크의 누적 에너지는 알려진 은하계 마그네타 플레어보다 2~4 자릿수 더 높으며, 이는 진화된 마그네타보다 훨씬 강력한 자기장($B_\phi \sim \text{few} \times 10^{16}$ G)을 암시한다.

의의 및 주장
본 논문은 GRB 211211A 및 230307A와 같은 특이한 LGRB의 기원에 관한 오랜 논쟁을 해결하는 "통합된 그림"을 제공한다고 주장한다. 사전 병합 마그네타를 포함하는 단일 NS–WD 병합 사건에 프롬프트 방출, 후광, 킬로노바를 연결함으로써, 모델은 다음 사항들에 대한 일관된 설명을 제공한다:

• 프롬프트 방출의 다중 성분 시간 구조.
• 마그네타 스핀다운에 의해 구동되는 X선 플래토의 존재.
• r-과정 중원소 생성 및 그에 따른 킬로노바 생성(기존 NS–WD 병합 모델의 한계를 극복).

저자들은 모델이 견고하지만, 활발한 강착에도 불구하고 GF 메커니즘을 지원할 만큼 중성자별 지각이 충분히 차갑게 유지되어야 한다는 가정에 의존하고 있으며, 이는 향후 연구에서 추가적인 조사가 필요한 한계점임을 인정한다.