상대론적 MHD 시뮬레이션을 통한 사지타리우스 A와 M87의 복사 모델

초록

이 논문은 일반 상대론적 MHD(GRMHD) 시뮬레이션에서 얻은 물리량에 전자 분포와 복사 모델을 결합해, Sagittarius A와 M87의 mm‑VLBI 관측과 비교한다. 시뮬레이션 기반 이미지와 스펙트럼을 이용해 Sgr A의 시야각·전자 온도·질량 흡수율을 제약하고, M87에서는 저각도에서의 제트와 디스크 복합 모델을 통해 제트 발사 영역의 크기와 블랙홀 그림자 검출 가능성을 예측한다.

상세 분석

본 연구는 세 가지 핵심 구성요소—GRMHD 시뮬레이션, 전자 분포 모델, 그리고 전자 복사 전이(레이트레이싱)—를 체계적으로 연결한다. 먼저, 두 개의 코드(Cosmos++와 3D‑HARM)로 생성된 네 개의 3차원 GRMHD 시뮬레이션을 사용해, 강자성장 불안정(MRI)에 의해 유도된 난류와 자기장 구조를 재현한다. 시뮬레이션은 반지름 0.25 M까지의 준정상 상태 흐름을 제공하지만, 방사와 전자-이온 온도 차이를 직접 포함하지 않는다. 따라서 전자 온도는 압력·밀도에서 유도된 이온‑전자 온도 비(Ti/Te) 혹은 자기 에너지에 비례하는 비열 비율(nnth∝b²)로 지정한다. Sgr A*는 전자들이 거의 열적이라고 가정하고, 일정한 Ti/Te 비율을 통해 전자 온도를 계산한다. 반면 M87는 디스크와 제트가 동시에 기여할 가능성을 열어두고, 제트 영역을 b²/ρc²>1 로 정의해 비열 전자 분포(전력법칙)와 자기 에너지 비례 스케일링을 적용한다. 복사 전이는 synchrotron 방출만을 고려하고, inverse‑Compton은 무시했지만, 이는 M87 고에너지 스펙트럼에 영향을 미칠 수 있다. 레이트레이싱은 grtrans와 geokerr를 이용해 광선이 null geodesic을 따라 이동하면서 복사 전이 방정식을 적분한다. 시뮬레이션 데이터를 2000–4000 M 시간 구간에 걸쳐 시간 평균하거나, 개별 프레임을 이용해 ‘영상 영화’를 생성한다.

Sgr A*에 대해서는 1.3 mm VLBI 가시도(visibility)와 0.4 mm·1.3 mm 총 플럭스를 동시에 피팅한다. 자유 파라미터는 질량 흡수율(ṁ), Ti/Te, 시야각(i), 그리고 스카이 방향(ξ)이며, 베이지안 혹은 χ² 최소화 방법을 적용한다. 결과는 i≈60°±15°, ξ≈−70°(+86°/−15°), 전자 온도 Te≈6±2×10¹⁰ K, 질량 흡수율 ṁ≈3⁺⁷₋₁×10⁻⁹ M⊙ yr⁻¹ 로, 기존 RIAF 모델과 일치한다. 현재 VLBI 배열(SMTO, CARMA, APEX 등)만으로는 물리적 흐름 모델을 구분하기 어렵지만, 향후 ALMA·LMT 등 북‑남 방향 긴 기저선이 추가되면 블랙홀 그림자(≈3 Gλ에서 최소값) 검출이 가능해진다.

M87에 대해서는 디스크/제트(DJ1)와 제트 전용(J2) 두 모델을 모두 적용했으며, 대규모 제트 관측에서 얻은 시야각(i<40°)과 스카이 방향(−115°<ξ<−75°)을 고정한다. 제트 모델에서는 반대 제트(counter‑jet)의 synchrotron 복사가 주된 mm 방출원이며, 제트 발사 영역이 반경 ≈10 M 이내에 국한될 경우 이미지 크기는 33–44 µas(≈4–5 Schwarzschild 반경)로 예측된다. 이는 현재 Sgr A* 이미지와 거의 동일한 규모이며, 향후 북‑남 기저선이 포함된 VLBI 관측에서 블랙홀 그림자(광자 궤도 근처 밝은 링)가 3000 Gλ 부근에서 최소값을 보이며 검출 가능성이 제시된다.

전반적으로 이 논문은 3차원 GRMHD 시뮬레이션을 기반으로 전자 물리와 복사를 일관되게 연결함으로써, mm‑VLBI 관측이 블랙홀 주변 흐름의 시야각·전자 온도·제트 구조 등을 직접 제약할 수 있음을 실증한다. 또한, 현재 배열로는 그림자 검출이 어려우나, ALMA·LMT와 같은 고감도·북‑남 기저선이 추가될 경우 최초의 블랙홀 사건지평선(그림자) 관측이 실현될 전망이다.