방사선 구동 분수 모델을 통한 AGN 토러스 피복 인자와 에디션 비율 관계 연구

초록

본 연구는 2차원 방사선 유체역학 시뮬레이션 결과를 바탕으로 XSTAR를 이용한 1차원 광전이온 평형 계산을 수행하여, 중성 수소(H I)와 이온화 수소(H II)의 각도별 분포와 그에 따른 Compton‑thin(C₍₂₂₎) 및 Compton‑thick(C₍₂₄₎) 피복 인자를 추정한다. 결과는 C₍₂₂₎(전체 가스)는 에디션 비율에 무관하게 약 70 %이며, 중성 수소만 고려한 C₍₂₂₎와 C₍₂₄₎는 약 30 %로 일정함을 보여준다. 관측과 비교했을 때 고에디션 비율에서는 이론값이 관측값보다 낮으며, 이를 해소하려면 현재 시뮬레이션 영역 바깥 10 pc 이상까지의 Compton‑thin 가스가 필요함을 제시한다.

상세 분석

이 논문은 AGN 토러스의 피복 인자를 물리적으로 해석하기 위해 두 단계의 계산을 결합한다. 첫 단계는 Y. Kudoh et al. (2024)에서 수행한 2차원 축대칭 방사선 유체역학(RHD) 시뮬레이션으로, SMBH 질량 10⁷ M☉, 먼지‑가스 비율 1 %를 가정하고 로그 Eddington 비율(log R_Edd) = –3, –2, –1, 0 네 경우에 대해 r = 10⁻⁴–2 pc, z = ±2 pc 범위의 가스 밀도와 속도장을 도출한다. 시뮬레이션 결과는 적도면에 log n_H ≈ 8 (cm⁻³) 수준의 고밀도 원반과, 극지방으로 갈수록 log n_H ≈ 4 (cm⁻³) 수준의 흐름을 보여준다. 평균 수소 컬럼 밀도 N_H(θ)는 θ ≈ 40–50°에서 10²² cm⁻²에 도달하고, 65–75° 사이에서는 10²⁵ cm⁻²를 초과해 거의 완전 중성 상태가 된다.

두 번째 단계는 XSTAR를 이용한 1차원 광전이온 평형 계산이다. 각 RHD 셀을 독립적인 광원으로 간주하고, 내부 셀에서 얻은 스펙트럼을 외부 셀의 입력으로 순차 전파한다. 입력 파라미터는 셀의 내·외부 반경, n_H, 1–1000 Ryd 구간의 X‑ray 광도, 그리고 SED(디스크+코로나)이다. 디스크 성분은 등방성으로 가정했으며, 코로나의 광자 지수 Γ는 Γ = 0.32 log R_Edd + 2.27 로 설정하였다. 이렇게 얻은 ionization parameter ξ = L/(n_H r²)와 온도, 전자 밀도 등을 바탕으로 H I와 H II의 비율을 계산하였다. 결과적으로 H I가 차지하는 영역은 전체 가스의 약 30 %에 해당하며, 이 영역이 10²² cm⁻²와 10²⁴ cm⁻² 임계값을 모두 초과하는 경우가 동일하게 나타나 C₍₂₂₎(H I) ≈ C₍₂₄₎(H I) ≈ 0.3 로 수렴한다.

관측과의 비교에서는 Ricci et al. (2017a,b)의 X‑ray 기반 피복 인자와 일치하는 부분과 차이가 있는 부분이 명확히 드러난다. 특히 고에디션 비율(log R_Edd = 0)에서는 관측된 C₍₂₂₎가 30–40 % 수준인 반면, 이론값은 30 %에 머물러 차이가 발생한다. 저에디션 비율에서는 이론값이 관측값과 거의 일치한다. 저자들은 이 차이를 메우기 위해 현재 시뮬레이션이 2 pc까지밖에 확장되지 않은 점을 지적하고, 최소 10 pc 이상의 Compton‑thin(10²² < N_H < 10²⁴ cm⁻²) 가스가 추가로 존재해야 한다고 주장한다. 이는 실제 AGN 주변에 광범위한 저밀도 가스가 존재한다는 최근 적외선·분자 가스 관측과도 일맥상통한다.

이 연구는 기존의 방사선‑구동 분수 모델에 광전이온 평형을 결합함으로써, 토러스의 물리적 구조와 피복 인자를 직접적으로 연결하는 중요한 시도를 보여준다. 다만, 2차원 축대칭 가정, 등방성 디스크 방사, 그리고 시뮬레이션 영역 제한 등 몇몇 근사점이 결과 해석에 영향을 미칠 수 있기에, 향후 3차원 전방위 시뮬레이션과 더 넓은 공간 스케일을 포함한 모델링이 필요하다.