[논문 리뷰] Jet Formation in Black Hole Accretion Systems II: Numerical Models

Jonathan C. McKinney|arXiv (Cornell University)|2005. 06. 16.

Gamma-ray bursts and supernovae인용 수 3

한 줄 요약

이 논문은 페어 생성을 고려한 일반 상대성 이론적 자기유체역학(그림: GRMHD) 시뮬레이션을 통해 감마선 폭발(GRB), 활성 은하핵(AGN), X선 이중성계에서의 블랙홀 순환의 원리를 모델링한다. Poynting-주도 제트는 자기 에너지 추출로 인해 콜라프서 모델에서 Γ ∼ 100–1000의 로렌츠 인자에 도달하는 것으로 나타났으며, 보존된 양성자로 오염된 Poynting-보존 제트는 Γ ∼ 1.5–3에 도달한다. 제트 가속은 큰 반경에서 발생하며, 시간적 변동성은 토로이드형 자기장 불안정성에 의해 주도된다.

ABSTRACT

In a companion theory paper, we presented a unified model of jet formation. We suggested that primarily two types of relativistic jets form near accreting black holes: a potentially ultrarelativistic Poynting-dominated jet and a Poynting-baryon jet. We showed that, for the collapsar model, the neutrino-driven enthalpy flux (classic fireball model) is probably dominated by the Blandford-Znajek energy flux, which predicts a jet Lorentz factor of $Γ\sim 100-1000$. We showed that radiatively inefficient AGN, such as M87, are synchrotron-cooling limited to $Γ\sim 2-10$. Radiatively efficient x-ray binaries, such as GRS1915+105, are Compton-drag limited to $Γ\lesssim 2$, but the jet may be destroyed by Compton drag. However, the Poynting-baryon jet is a collimated outflow with $Γ\sim 1-3$. Here we present general relativistic hydromagnetic simulations of black hole accretion with pair creation used to simulate jet formation in GRBs, AGN, and x-ray binaries. Our collapsar model shows the development of a patchy ``magnetic fireball'' with typically $Γ\sim 100-1000$ and a Gaussian structure. Temporal variability of the jet is dominated by toroidal field instabilities for $\gtrsim 10^2$ gravitational radii. A broader Poynting-baryon jet with $Γ\sim 1.5$ could contribute to a supernova.

연구 동기 및 목표

• 다양한 블랙홀 순환 시스템, 즉 GRB, AGN, X선 이중성계에서의 상대론적 제트 형성의 물리적 메커니즘을 통합하는 것.
• 최소한의 물리 원칙을 사용하여 제트 에너지 기원, 조성, 빔화 정도, 최종 로렌츠 인자에 대한 원인을 규명하는 것.
• McKinney(2005b)의 이론 프레임워크를 페어 생성과 뉴트론 확산을 포함한 수치 시뮬레이션을 통해 검증하는 것.
• 복사 효율이 높은 시스템에서 콜럼브 드래그와 싱크로트론 냉각이 제트의 로렌츠 인자를 어떻게 제한하는지 평가하는 것.
• 미래의 모델링에 사용하기 위해, 제트 구조의 보편적 특성, 예를 들어 초음속 유동과 자기 유사성 행동을 식별하는 것.

제안 방법

• 수치 시뮬레이션은 페어 생성에 대한 효과적 모델과 뉴트론의 푸아르 확산을 고려한 일반 상대성 이론적 자기유체역학(GRMHD) 방정식을 해결한다.
• 시뮬레이션은 케러 메트릭 스페이스타임을 사용하며, 전자기장의 스트레스-에너지 텐서가 $ T^{ m EM}_{ ueta} = \frac{b^2}{2}(u_\nu u_\beta + P_{\nu\beta}) - b_\nu b_\beta $ 로 주어지는 이상 자기유체역학 방정식을 해결한다. 여기서 $ b^\mu $ 는 자기장 4벡터이다.
• 모델은 GRB에서 '자기 화염공' 시나리오를 포함한다. 여기서 페어 생성은 터널 영역에서 발생하며, 뉴트론 확산에 의해 제한되어 최종 로렌츠 인자가 결정된다.
• 제트 역학은 알프베인, 패스트, 정지 표면과 같은 특성 표면을 사용하여 분석되며, 알프베인 점이 에르고스피어 내부에 있을 경우에만 에너지 추출이 가능하다.
• 시뮬레이션은 내부 순환 디스크에서 큰 반경까지 제트의 진화를 추적하며, 외부 패스트 표면 이외의 영역에서 가속 영역을 식별한다.
• 곡선 피팅을 통해 제트 구조의 조각별 자기 유사성을 추출하여 향후 제트 역학 모델링에 활용한다.

실험 결과

연구 질문

• RQ1 블랙홀 순환 시스템에서의 상대론적 제트의 최종 로렌츠 인자는 어떤 물리적 메커니즘에 의해 결정되는가?
• RQ2 페어 생성과 뉴트론 확산은 GRB에서 Poynting-주도 제트의 조성과 최대 로렌츠 인자에 어떤 영향을 미치는가?
• RQ3 반경 100 중력 반경 이상에서 토로이드형 자기장 불안정성은 상대론적 제트의 시간적 변동성에 어떤 역할을 하는가?
• RQ4 콜럼브 드래그와 싱크로트론 냉각은 X선 이중성계와 AGN와 같은 복사 효율이 높은 시스템에서 로렌츠 인자를 어떻게 제한하는가?
• RQ5 GRB, AGN, X선 이중성계의 제트 구조는 어느 정도 자기 유사성 또는 보편적 특성을 보이는가?

주요 결과

• 콜라프서 모델에서 Poynting-주도 제트는 관측된 단단한 감마선 폭발의 특성과 일치하는 최종 로렌츠 인자 Γ ∼ 100–1000에 도달한다.
• 제트 가속은 블랙홀 근처가 아니라 넓은 반경 범위에서 발생하며, 외부 패스트 표면이 주요 가속 영역을 나타낸다.
• 반경 ∼100 중력 반경 이상에서 제트의 시간적 변동성은 토로이드형 자기장 불안정성에 의해 주로 지배된다.
• Poynting-보존 제트는 콜라프서 제트로의 빔화된 유동을 형성하며, 최종 로렌츠 인자는 Γ ∼ 1.5에 도달한다. 이는 초신성에 의해 유도된 유동과 일치한다.
• 제트 구조는 조각별 자기 유사성을 보이며, 향후 제트 역학 모델링에 활용하기 위해 곡선 피팅 결과가 제공된다.
• 모델은 복사 효율이 높은 시스템에서 콜럼브 드래그가 Poynting-주도 제트를 파괴할 수 있으며, 이로 인해 Γ ≲ 2로 제한된다는 것을 확인한다. 또한 싱크로트론 냉각은 AGN 제트의 로렌츠 인자를 Γ ∼ 2–10으로 제한한다.