Three-dimensional GRMHD simulations of jet formation and propagation in self-gravitating collapsing stars

기술적 요약: 자체 중력을 가진 붕괴성에서의 제트 형성 및 전파에 대한 3 차원 GRMHD 시뮬레이션

문제 및 동기
장감마선 폭발 (GRB) 은 일반적으로 빠르게 회전하는 거대 별 (콜랩사) 이 항성 질량 블랙홀로 붕괴하는 것과 밀접하게 연관되어 있습니다. 콜랩사 모델은 장감마선 폭발에 대한 표준 설명이지만, 많은 수치 시뮬레이션은 고정된 블랙홀 질량과 스핀을 가정하거나 progenitor 별의 외피 자체 중력을 무시함으로써 물리를 단순화합니다. 실제로는 항성 외피의 질량이 초기 블랙홀 질량을 수 배 초과하는 경우가 많아, 자체 중력이 강착 흐름, 제트 형성 및 전파에 상당한 영향을 미칠 수 있음을 시사합니다. 저자들의 이전 2 차원 연구는 자체 중력이 강착 변동성을 증가시키고 제트 형성을 방해할 수 있음을 나타냈으나, 제트 전파를 포함한 3 차원 분석은 부재했습니다. 본 연구는 동일한 초기 조건 하에서 자체 중력을 포함하는 모델과 포함하지 않는 모델을 직접 비교하여 자체 중력이 제트 개구각, 출력, 최종 로런츠 인자 및 변동성을 포함한 제트 특성에 미치는 영향을 정량화하는 것을 목표로 합니다.

방법론
저자들은 HARM의 업데이트된 버전인 HARM-SELFG 코드를 사용하여 시간 의존적 3 차원 일반 상대성 자기유체역학 (GRMHD) 시뮬레이션 세트를 수행했습니다. 이 시뮬레이션은 블랙홀의 질량 ($M$) 과 각운동량 ($J$) 변화에 따라 커 (Kerr) 계량 성분이 동적으로 업데이트되는 진화하는 시공간에서 시스템을 진화시킵니다.

• 자체 중력 구현: 계산적으로 불가능한 완전한 아인슈타인 방정식을 풀기 대신, 저자들은 테우콜스키 (Teukolsky) 방정식과 CCK 형식주의에 영감을 받아 커 계량에 섭동 항을 추가하여 자체 중력을 구현했습니다. 이 접근법은 항성 외피의 자체 중력을 계량 질량 계수 $M(t,r)$과 스핀 매개변수 $a(t,r)$을 수정하는 섭동으로 취급합니다. 이러한 수정은 사건의 지평선과 주어진 반경 사이의 체적에 걸쳐 에너지 - 운동량 텐서 성분 ($T^t_t$ 및 $T^r_\phi$) 을 적분하여 계산됩니다.
• 초기 조건: 모델은 씨앗 블랙홀 ($M_0 = 3 M_\odot$, 스핀 $a_0 = 0.8$) 로의 준구면 강착 흐름을 가정하며, 이는 항성 핵 ($M_\star = 25 M_\odot$) 으로 둘러싸여 있습니다. 가스는 다방정식 상태 방정식 ($\gamma = 4/3$) 을 따릅니다. 각운동량은 가장 안쪽 안정 원형 궤도 (ISCO) 에 대해 스케일링됩니다. 자기장은 $2 r_g$에서 가스 압력과 자기 압력의 비율인 플라즈마 매개변수 $\beta_0$가 1 또는 10 이 되도록 초기화되었으며, 이는 약 $10^8$ G 에서 약 $10^{15}$ G 까지의 자기장 세기에 해당합니다.
• 시뮬레이션 매개변수: 네 가지 모델이 실행되었습니다. 자체 중력을 포함하는 두 모델 (Model-1-SG, Model-2-SG) 과 포함하지 않는 두 모델 (Model-1-NSG, Model-2-NSG) 로, 주로 초기 자화 ($\beta_0 = 1$ 대 $\beta_0 = 10$) 에서 차이가 있습니다. 시뮬레이션은 $5 \times 10^4$ 중력 시간 단위 ($t_g$) 동안 실행되었습니다.

주요 기여 및 결과
이 연구는 자체 중력을 가진 콜랩사에서의 제트 형성에 대한 최초의 3 차원 GRMHD 시뮬레이션을 제시합니다. 주요 발견 사항은 다음과 같습니다:

1. 제트 억제 및 정적 간격: 자체 중력을 포함하는 모델에서 저자들은 제트 분사가 일시적으로 억제되거나 완전히 중단되는 "정적 간격 (Quiescent interval)"을 관찰했습니다. Model-1-SG 에서 제트는 발사되지만 질량 축적과 급격한 강착으로 인해 이후 억제되어 휴면 상태의 중심 엔진을 초래합니다. Model-2-SG(더 약한 자화) 에서는 자체 중력이 제트 형성을 완전히 방해하여 "실패한 GRB"를 초래합니다. 이러한 억제는 동일한 조건 하에서 자체 중력을 포함하지 않는 대응 모델에서는 관찰되지 않습니다.
2. 제트 에너지 및 변동성: 자체 중력은 제트가 방출하는 총 에너지를 크게 감소시킵니다. 자체 중력 모델과 비자체 중력 모델 사이의 방출 에너지 비율은 약 0.17~0.21 입니다. 자체 중력 모델의 최종 로런츠 인자 ($\Gamma_\infty$) 는 일반적으로 더 낮고 변동성이 크며, 급격한 초기 스파이크 followed by 감소를 보입니다. 반면, 비자체 중력 모델은 더 안정적이고 고출력의 제트를 유지합니다.
3. 제트 개구각: 자체 중력은 더 좁은 제트 개구각을 초래합니다. Model-1-SG 에서 각도는 제트가 정적 단계에 접근함에 따라 약 $10^\circ$에서 약 $4^\circ$로 감소합니다. 이는 자체 중력 유체가 제트 깔때기에 가하는 추가 압력 때문입니다. 비자체 중력 모델은 더 넓고 팽창하는 제트 (최대 약 $22^\circ$) 를 보입니다.
4. 블랙홀 진화: 자체 중력 모델은 비자체 중력 모델에 비해 훨씬 빠른 질량 강착률 (일시적으로 $10^3 M_\odot \text{s}^{-1}$ 초과) 을 보입니다. 결과적으로 블랙홀 질량은 급격히 증가하여 SG 모델에서는 약 $28 M_\odot$에 도달하는 반면, NSG 모델에서는 약 $5\text{--}16 M_\odot$에 도달합니다. 블랙홀 스핀 진화 또한 다르며, 비자체 중력 모델은 블랜드포드 - 즈나젝 (Blandford-Znajek) 과정을 통해 더 많은 회전 에너지를 추출하여 더 낮은 최종 스핀을 보이는 반면, 자체 중력 모델은 효율적인 에너지 추출의 방해로 인해 더 높은 스핀을 유지합니다.
5. MAD 상태 안정성: 효율적인 제트 발사에 중요한 자기적으로 정지된 원반 (MAD) 상태는 자체 중력 모델에서 유지하기 어렵습니다. 자체 중력에 의해 주도되는 급격한 질량 유입은 강착 흐름을 막을 만큼 자기 플럭스가 충분히 축적되는 것을 방지하여 차원 없는 자기 플럭스 ($\phi_{MAD}$) 가 낮아지고 제트 효율 ($\eta$) 이 감소합니다. 비자체 중력 모델에서는 MAD 상태가 더 오랜 기간 안정적으로 유지됩니다.

의의 및 주장
저자들은 계산 결과가 자체 중력이 제트 전파와 특성을 크게 수정하는 결정적 요소임을 보여준다고 주장합니다. 구체적으로 그들은 다음과 같이 논합니다:

• 자체 중력은 질량 축적으로 인한 일시적인 제트 억제 메커니즘을 통해 GRB 즉각 방출의 특정 특징, 특히 정적 간격과 변동성을 설명할 수 있습니다.
• 자체 중력을 포함하면 성공적인 상대론적 제트를 생성할 수 있는 초기 조건 (자기장 세기, 회전, 스핀) 의 범위가 좁아져, 비자체 중력 모델이 성공할 시나리오에서도 실패한 GRB로 이어질 수 있습니다.
• 제트의 시간 척도, 변동성 및 개구각은 자체 중력이 포함되는지에 근본적으로 의존하므로, 이 효과를 무시한 이전 모델들은 제트의 안정성과 효율성을 과대평가했을 가능성이 있습니다.

이 연구는 자체 중력이 초기 제트 발사 시간을 반드시 지연시키지는 않지만 (제트가 feasible 한 경우), 그 이후의 진화를 극적으로 변화시켜 종종 간헐적이거나 실패한 제트 활동을 초래한다고 결론지으며, 이는 장감마선 폭발에서 관찰되는 복잡한 광도곡선 구조에 대한 보완적 설명을 제공합니다.