거대한 블랙홀, 가스는 어떻게 모이는가

초록

다중 규모 SPH 시뮬레이션을 이용해 은하 규모(킬로파섹)부터 0.1 pc 이하까지 가스 흐름을 추적하였다. 가스가 풍부하고 디스크가 지배적인 은하에서는 중력 불안정이 연속적으로 발생해 초대질량 블랙홀에 1–10 M⊙ yr⁻¹의 높은 유입률을 만든다. 10 pc 이상에서는 ‘바 안의 바’와 유사한 구조가 나타나지만, 나선, 고리, 클럼프 등 다양한 형태가 교대로 나타나며, 1–10 pc에서는 m=1 비대칭 모드가 지배해 가스를 서브 파섹 규모로 운반한다. 별 형성과 피드백, 가스·별의 자체 중력이 핵심 역할을 하며, 블랙홀 성장률과 핵별 형성률 사이에 스케일에 따라 다른 상관관계가 존재한다.

상세 분석

본 연구는 은하 전체 규모에서 초소형 원반까지 연속적인 물리 과정을 포착하기 위해 Smoothed Particle Hydrodynamics(SPH) 방식을 100여 차례에 걸쳐 다양한 초기 조건과 서브그리드 모델로 수행하였다. 초기 조건은 대규모 은하 합병 시뮬레이션과 바 불안정 디스크 시뮬레이션을 통해 300 pc 이하의 핵 영역을 재현했으며, 여기에는 가스, 별, 블랙홀, 자체 중력, 별 형성, 그리고 피드백 메커니즘이 모두 포함되었다. 가스가 풍부하고 디스크가 지배적인 경우, 킬로파섹 규모에서부터 수십 파섹까지 연속적인 중력 불안정(바, 나선, 클럼프 등)이 발생한다. 이 불안정은 토러스 형태의 가스 흐름을 형성하고, 효율적인 각운동량 운반을 통해 블랙홀 근처까지 물질을 끌어들인다. 특히 10 pc 이상에서는 전통적인 ‘bars‑within‑bars’ 모델과 유사하게 다중 바 구조가 나타나지만, 시뮬레이션은 이 구조가 일시적이며 전체 시간의 10–30 % 정도만 존재함을 보여준다. 이는 관측적으로 핵 바와 AGN 활동을 직접 연결하기 어려운 이유를 설명한다.

핵심적인 전이 구간은 1–10 pc 스케일이다. 이 영역에서는 블랙홀의 포텐셜이 지배적이 되어 기존의 바‑모드가 사라지고, 대신 m = 1의 비대칭(편심 원반 또는 단일 나선) 모드가 성장한다. 이 모드는 자체 중력과 별 형성 피드백에 민감하게 반응하며, 가스와 별이 동시에 비대칭 구조에 참여함으로써 각운동량을 급격히 감소시켜 서브 파섹 규모로 물질을 운반한다. 따라서 별 형성률과 피드백을 무시하면 m = 1 모드의 성장과 가스 흐름이 과소평가된다.

시뮬레이션 결과는 블랙홀 질량 증가율(BHAR)과 핵별 형성률(SFR) 사이에 스케일 의존적인 상관관계를 제시한다. 100 pc 이하 핵 영역에서는 SFR과 BHAR이 거의 동시 변동하며, 큰 스케일(킬로파섹)에서는 평균적인 가스 공급이 일정해도 순간적인 BHAR 변동이 크게 나타난다. 이는 관측적으로 ‘AGN‑SF 연결’이 스케일에 따라 다르게 나타나는 현상을 이론적으로 뒷받침한다.

전반적으로 이 연구는 가스 흐름을 다중 스케일에서 일관되게 모델링함으로써, 전통적인 바‑안‑바 메커니즘을 확장하고, m = 1 비대칭 모드가 초소형 원반 형성에 핵심적임을 밝힌다. 또한 별 형성·피드백과 자체 중력이 블랙홀 성장에 미치는 정량적 영향을 최초로 제시함으로써, 관측적 AGN‑SF 상관관계 해석에 중요한 이론적 토대를 제공한다.