낮은 각운동량 흐름으로 본 우리 은하 중심 블랙홀 Sgr A 의 변동성

초록

이 연구는 Sgr A* 주변 별들의 풍선에서 추정된 특정 각운동량과 전체 에너지를 기반으로, 2차원 수소역학 시뮬레이션을 수행한다. 낮은 각운동량을 가진 물질이 비정상적으로 흡입되면서 수십백여 개 슈바르츠시트 반경(Rs) 안에서 불규칙하게 진동하는 충격파가 형성되고, 이로 인해 광도와 질량 유출이 각각 수 퍼센트5배, 2~7배 정도 변동한다. 흐름은 고도로 advective이며 복사 효율은 10⁻⁵–10⁻³ 수준에 불과한다. 두 온도 모델과 동기트론 복사를 포함하면, 입력 질량 흡입률 4 × 10⁻⁶ M⊙ yr⁻¹가 관측된 10³⁶ erg s⁻¹의 광도를 재현한다.

상세 분석

본 논문은 Sgr A* 주변의 젊은 별들(특히 IRS 13E 등)에서 방출되는 강풍이 블랙홀에 공급하는 물질의 특정 각운동량(λ)과 특정 에너지(ε)를 직접 측정한 최신 관측 결과를 토대로, 저각운동량(λ≈1.5–2.0 R_g c) 흐름이 어떻게 흡수·방출 과정을 겪는지를 2차원 비점성, 비자기수소역학 모델로 재현한다. 시뮬레이션은 축대칭 좌표계에서 라디얼·극 방향으로 512 × 256 격자를 사용했으며, 초기 조건은 원형 원판 형태의 가스 구름을 가정하고, 경계조건은 외부에서 일정한 질량 유입률(4 × 10⁻⁶ M⊙ yr⁻¹)과 자유 방출(아웃플로우)으로 설정하였다.

핵심 결과는 다음과 같다. 첫째, 입력된 λ와 ε가 특정 범위(λ≈1.6–1.9 R_g c, ε≈0.001–0.003 c²)에 있을 때, 흐름은 정상적인 원형 원판이 되지 못하고, 중심부(30–160 R_s)에서 강한 충격파가 주기적으로 형성·소멸한다. 이 충격파는 전형적인 Rankine‑Hugoniot 조건을 만족하지만, 비정상적인 압력·밀도 구배와 고속 인플로우가 결합되어 비선형 진동을 일으킨다. 충격 전후의 온도 차이는 약 10배 이상이며, 포스트‑쇼크 영역에서는 전자와 이온이 서로 다른 온도로 유지되는 두 온도 구조가 형성된다.

둘째, 충격 진동에 따라 질량 유출률(𝑀̇_out)은 0.02–0.99 𝑀̇_in 사이에서 크게 변동한다. λ가 증가할수록 포스트‑쇼크 가스가 원심력에 의해 더 많이 탈출하게 되며, 최악의 경우 입력 질량의 거의 전부가 바깥으로 방출된다. 이는 관측된 Sgr A*의 저밀도, 저광도 상태와 일맥상통한다.

셋째, 전체 흐름은 고도로 advective이며, 방사 효율 η = L/Ṁc²가 10⁻⁵–10⁻³ 수준에 머문다. 이는 전통적인 표준 얇은 원판 모델(η≈0.1)과는 크게 차이가 있다. 저효율의 원인은 두 온도 플라즈마에서 전자와 이온 간 열 교환이 비효율적이며, 동기트론 복사가 주된 냉각 메커니즘이지만, 자기장 강도가 약해(β≈10) 복사 손실이 제한적이기 때문이다.

넷째, 시뮬레이션에서 관측된 변동 시간 스케일은 1 시간에서 10 일에 이른다. 이는 충격 파동이 중심부를 왕복하는 자유 낙하 시간(t_ff≈R/c_s)과 일치한다. 이러한 시간 변동은 Sgr A*에서 보고된 X‑ray·NIR 플레어(수시간~수일 주기)와 직접적인 연관성을 제시한다.

마지막으로, 저각운동량 흐름이 정상적인 원판 형태를 이루지 못하고 충격‑진동 구조를 보이는 이유는, λ와 ε가 회전하는 무점성·단열 흐름에서 정상적인 충격이 존재하는 파라미터 영역에 위치하거나 그 근처에 있기 때문이다. 따라서, 작은 파라미터 변화만으로도 흐름이 완전히 다른 동역학(정상 충격 vs. 비정상 진동)으로 전이될 수 있음을 보여준다.

이러한 결과는 Sgr A*와 같은 저활성 은하핵의 물리적 상태를 이해하는 데 중요한 힌트를 제공한다. 특히, 저각운동량, 저에너지 유입이 고전적인 ADAF(Advection‑Dominated Accretion Flow) 모델과는 다른, 충격‑유도 변동성을 내포한다는 점은 향후 고해상도 관측 및 MHD 시뮬레이션에서 검증될 필요가 있다.