[논문 리뷰] Spectral calculations of 3D RMHD simulations of super-Eddington accretion onto a stellar-mass black hole

Brianna S. Mills, Shane W. Davis외 2명|arXiv (Cornell University)|2023. 04. 17.

Astrophysical Phenomena and Observations인용 수 4

한 줄 요약

이 연구는 6.62 M⊙ 블랙홀에 대한 초-에드링턴 질량 합성 시뮬레이션의 몬테카를로 방사선 전달 후처리를 통해 합성 X선 스펙트럼을 계산한다. 다중군 콜만 스캐터링 처리 방식이 회색 근사보다 더 딱딱한 스펙트럼과 높은 온도를 산출하며, 이는 NGC 1313 X-1과 같은 ULX에서 관측된 딱딱한 X선 尾를 매우 잘 재현한다.

ABSTRACT

We use the Athena++ Monte Carlo (MC) radiation transfer module to post-process simulation snapshots from non-relativistic Athena++ radiation magnetohydrodynamic (RMHD) simulations. These simulations were run using a gray (frequency integrated) approach but were also restarted and ran with a multi-group approach that accounts for Compton scattering with a Kompaneets operator. These simulations produced moderately super-Eddington accretion rates onto a 6.62 $M_\odot$ black hole. Since we only achieve inflow equilibrium out to 20-25 gravitational radii, we focus on the hard X-ray emission. We provide a comparison between the MC and RMHD simulations showing that the treatment of Compton scattering in the gray RMHD simulations underestimates the gas temperature in the regions above and below the accretion disk. In contrast, the restarted multi-group snapshots provides a treatment for the radiation field that is more consistent with the MC calculations, and result in post-processed spectra with harder X-ray emission compared to their gray snapshot counterparts. We characterize these MC post-processed spectra using commonly employed phenomenological models used for spectral fitting. We also attempt to fit our MC spectra directly to observations of the ultraluminous X-ray source (ULX) NGC 1313 X-1, finding best fit values that are competitive to phenomenological model fits, indicating that first principle models of super-Eddington accretion may adequately explain the observed hard X-ray spectra in some ULX sources.

연구 동기 및 목표

• 초-에드링턴 질량 합성에 의한 별질량 블랙홀에 대한 3D 비상대론적 RMHD 시뮬레이션에서 합성 X선 스펙트럼을 계산하기 위해.
• 콜만 스캐터링 처리 방식(회색 vs. 다중군)이 시뮬레이션된 스펙트럼과 기체 온도에 미치는 영향을 평가하기 위해.
• 후처리된 몬테카를로 스펙트럼을 ULX인 NGC 1313 X-1의 관측 데이터와 비교하기 위해.
• 현상학적 모델에 의존하지 않고도 원리 기반 시뮬레이션이 ULX에서 관측된 딱딱한 X선 성분을 재현할 수 있는지 평가하기 위해.

제안 방법

• Athena++ 비상대론적 RMHD 시뮬레이션의 스냅샷을 몬테카를로 방사선 전달 모듈을 사용해 후처리한다.
• 초기 시뮬레이션에서 주파수 통합 투과도와 방사선 힘을 사용한 회색 방사선 전달 방법을 적용한다.
• 콜만 스캐터링을 더 정확하게 모의하기 위해 콘스탄틴 연산자( Kompaneets operator)를 사용한 다중군 접근법을 사용해 시뮬레이션을 재시작한다.
• 흑체 가정을 사용해 초기 콜만 냉각 추정치를 얻고, 몬테카를로 광자 전달을 통해 스펙트럼 에너지 분포를 계산한다.
• 결과로 얻은 합성 스펙트럼을 표준 현상학적 모델(예: SIMPL, POW)로 피팅하고, XMM-Newton 및 NuSTAR 관측치와 직접 비교한다.
• 몬테카를로 스펙트럼에서 XSPEC 테이블 모델을 생성하여 관측 데이터와 직접 피팅할 수 있도록 한다.

Figure 1: The mass accretion rate $\dot{M}$ as a function of gravitational radius $r_{\rm g}$ from the black hole for each Athena++ RMHD simulation snapshot (see Table 1 ). ULX4a is the dot-dash line, ULX4b is the dashed line, ULX2.5 is the dotted line, and ULX1.3 is the solid line.

실험 결과

연구 질문

• RQ1 회색 근사 대비 다중군 RMHD 시뮬레이션에서 콜만 스캐터링 처리 방식이 예측된 X선 스펙트럼과 기체 온도에 어떤 영향을 미치는가?
• RQ2 3D RMHD 시뮬레이션의 후처리 몬테카를로 스펙트럼이 NGC 1313 X-1과 같은 ULX에서 관측된 딱딱한 X선 꼬리와 얼마나 잘 일치하는가?
• RQ3 초-에드링턴 질량 합성의 원리 기반 시뮬레이션이 스펙트럼 피팅에 사용되는 현상학적 모델과 경쟁 가능한 스펙트럼을 생성할 수 있는가?
• RQ4 현재 시뮬레이션 설정의 한계—특히 유입 평형과 방사선장 처리 측면에서—관측된 ULX 스펙트럼을 재현하는 데 어떤 영향을 미치는가?

주요 결과

• 흑체 가정을 사용한 회색 RMHD 시뮬레이션은 흡착 디스크 상하의 터널 영역에서 기체 온도를 과소평가한다. 이는 콜만 냉각 추정치의 정확도가 떨어지기 때문이다.
• Kompaneets 유사 콜만 스캐터링 처리를 적용한 다중군 시뮬레이션은 그 회색 대비 더 높은 온도와 상당히 딱딱한 X선 스펙트럼을 산출한다.
• 다중군 스냅샷에서 유도된 스펙트럼은 딱딱한 X선 대역에서 Γ ≈ 2–3의 파wr-법칙 기울기를 보이며, NGC 1313 X-1 및 홀름버그 IX X-1의 관측치와 일치한다.
• 회색 RMHD 스냅샷의 스펙트럼을 피팅할 때 SIMPL 모델이 파워-법칙(POW) 모델보다 딱딱한 X선 성분에 더 잘 맞는다. 이 경우 Γ ≈ 2–4이다.
• 몬테카를로 스펙트럼을 NGC 1313 X-1 관측치에 직접 피팅한 결과, 10 keV 이하에서는 좋은 피팅이 이루어졌지만, 소스 성분에 따라 딱딱한 X선 꼬리 부분에서 약간 너무 기울어지거나 평탄한 경향이 있었다.
• 최적 피팅 정규화는 관측된 것보다 낮은 빛의 강도를 암시하며, 이는 일반 상대성 이론 효과가 누락되었거나, 유입 평형이 완전하지 않아 발생한 오차일 가능성이 있다.

Figure 2: Azimuthally averaged gas density $\rho$ (g cm -2 ) of the Athena++ simulation showing the inner 25 $r_{\rm g}$ with a small inset plot that shows the gas density out to 500 $r_{\rm g}$ .