상대론적 MHD 시뮬레이션이 밝힌 사지점성 Sgr A의 비밀

초록

이 논문은 3차원 일반 상대론적 자기유체역학(GRMHD) 시뮬레이션과 두 온도 싱크로트론 방출 모델을 결합해, 우리 은하 중심 초대질량 블랙홀 Sgr A*의 1.3 mm 및 0.87 mm 파장에서의 관측 데이터를 재현한다. 장거리 간섭계(VLBI)와 스펙트럼 지수 측정 결과를 이용해, 면관측(얼굴) 각이 20도 이하인 모델은 3σ 수준에서 배제된다. 최적 모델은 경사각 i≈50°(±15–35°), 위치각 ξ≈−23°(−22° ~ +97°), 전자 온도 Tₑ≈5.4×10¹⁰ K(±3.0×10¹⁰ K), 질량 유입률 Ṁ≈5×10⁻⁹ M☉ yr⁻¹(−2 ~ +15)이며, 90% 신뢰구간을 가진다. 최적 모델에서는 블랙홀 그림자가 가려지지 않아, 칠레‑캘리포니아·애리조나·멕시코 간의 1.3 mm·0.87 mm 베이스라인으로 직접 검출이 가능하다. 변광은 내부 자기 난류에 의해 유발되며, 광학 얇은 주파수에서 변동은 질량 유입률과 강하게 상관된다.

상세 분석

이 연구는 최근 Event Horizon Telescope(EHT)와 같은 초고해상도 mm‑파 관측이 제공하는 정밀 데이터를, 최신 3차원 일반 상대론적 자기유체역학(GRMHD) 시뮬레이션과 결합함으로써, Sgr A*의 물리적 구조와 방사 메커니즘을 정량적으로 규명하려는 시도이다. 저자들은 두 온도 플라즈마 모델을 채택했는데, 이는 전자와 이온이 서로 다른 온도를 유지한다는 가정에 기반한다. 전자 온도는 시뮬레이션 내부에서 직접 계산되지 않으므로, 전자-이온 비율을 파라미터화하여 관측과 일치하도록 조정한다. 이는 특히 1.3 mm(230 GHz)와 0.87 mm(345 GHz) 대역에서의 싱크로트론 방출이 전자 온도에 민감하게 반응한다는 점을 활용한다.

시뮬레이션은 회전하는 흑색체 주변의 마그네틱 레이디에이션 전단( MRI )에 의해 유발되는 난류와 플라즈마 흐름을 재현한다. 저자들은 세 가지 서로 다른 초기 자기장 구성을 사용했으며, 모두 비슷한 전반적 구조를 보였다: 내측( r ≲ 10 GM/c² )에 고밀도, 고온의 전자 플라즈마가 존재하고, 이 영역이 mm‑파 방출의 주된 원천이다. 특히, 이 영역은 거의 등온성을 유지하며, ν/ν_c ≈ 1–20 범위의 싱크로트론 임계 주파수 ν_c와 비교해 광대역에 걸쳐 방출한다. 이는 전자 온도가 Tₑ ≈ 5×10¹⁰ K 수준으로, 열 싱크로트론이 광학 얇은 상태에서 효율적으로 방출됨을 의미한다.

관측과의 비교는 두 단계로 진행된다. 첫째, 시뮬레이션에서 생성된 이미지의 푸리에 변환을 통해 복잡한 가시도(amplitude)와 폐쇄 위상(closure phase)을 계산하고, 이를 실제 VLBI 베이스라인(칠레–캘리포니아, 칠레–애리조나, 칠레–멕시코 등)과 비교한다. 둘째, 스펙트럼 지수(α, S ∝ ν^α)와 전체 플럭스 밀도(≈2.4 Jy at 230 GHz)를 재현하도록 전자 온도와 질량 유입률을 조정한다. 통계적 검증을 위해 마르코프 체인 몬테카를로(MCMC) 샘플링을 수행했으며, 결과적으로 i ≤ 20°인 면관측 모델은 3σ 수준에서 배제되었다. 이는 기존의 “얼굴‑온” 블랙홀 그림자 모델이 실제 데이터와 불일치함을 강력히 시사한다.

최적 파라미터는 i ≈ 50°(±15–35°), ξ ≈ −23°(−22° ~ +97°)이며, 이는 블랙홀 회전축이 우리 시야에 대해 중간 정도 기울어져 있음을 의미한다. 전자 온도 Tₑ ≈ 5.4×10¹⁰ K와 질량 유입률 Ṁ ≈ 5×10⁻⁹ M☉ yr⁻¹는 각각 ±3.0×10¹⁰ K, −2 ~ +15×10⁻⁹ M☉ yr⁻¹의 불확실성을 가진다. 이러한 값들은 이전 연구에서 제시된 범위와 일치하지만, 특히 전자 온도는 기존의 “극히 뜨거운” 모델보다 낮은 편이다. 이는 전자-이온 비열 교환 효율이 상대적으로 높아 전자 온도가 이온 온도에 근접한다는 물리적 해석을 가능하게 한다.

변광 측면에서는, 시뮬레이션이 30–60 분 주기의 플럭스 변동을 재현한다. 변동은 주로 내부 마그네틱 압력의 급격한 변동, 즉 MRI‑유도 난류에 의해 전자 온도와 밀도가 순간적으로 상승하면서 발생한다. 흥미롭게도, 광학 얇은 주파수(230 GHz, 345 GHz)에서의 변동은 질량 유입률 Ṁ와 높은 상관관계를 보이며, 이는 관측된 “mm‑플레어”가 실제로는 내측 흡입 흐름의 변동을 반영한다는 강력한 증거가 된다. 또한, 시뮬레이션은 플레어가 발생할 때 이미지 중심이 약간 이동하고, 비대칭성이 증가함을 보여준다. 이는 향후 고해상도 VLBI가 폐쇄 위상 변화를 통해 플레어 발생 메커니즘을 직접 탐지할 수 있음을 시사한다.

마지막으로, 저자들은 블랙홀 그림자(광학적 실루엣)가 최적 모델에서 전혀 가려지지 않으며, 1.3 mm와 0.87 mm 대역에서 베이스라인 길이가 3000–5000 km 정도인 경우 가시도 최소값이 0.1–0.2 Jy 수준으로 떨어진다. 이는 현재 EHT가 이미 충분히 감지 가능한 수준이며, 특히 폐쇄 위상 측정이 그림자 구조를 명확히 구분하는 데 유리함을 강조한다. 따라서 향후 관측 캠페인에서는 더 많은 베이스라인과 높은 동시성(시간 해상도) 확보가 핵심 과제가 된다.