MAXI 조사로 본 지역 AGN X선 광도 함수 재조명

초록

MAXI/GSC 4–10 keV 전천구조에서 37개의 비블레이저 AGN(0.002 < z < 0.2)을 이용해 지역( z≈0 ) X선(2–10 keV) 광도 함수를 새롭게 구축하였다. 흡수된 AGN 비율이 광도에 따라 급격히 감소함을 확인하고, 이중 전력법칙 형태의 광도 함수를 도출했다. 2–10 keV와 14–195 keV(Swift/BAT) 대역의 광도 함수를 비교한 결과, 평균 스펙트럼 인덱스 Γ가 L_X≈10^44 erg s⁻¹를 기준으로 변화한다는 새로운 해석을 제시한다.

상세 분석

본 연구는 MAXI/GSC가 제공하는 4–10 keV 전천구조 데이터를 기반으로, 비블레이저 AGN 37개(식별 완전도 > 97 %)를 선별하여 지역( z ≈ 0 ) X선 광도 함수를 재구성하였다. 샘플 선정 과정에서 Cen A와 복합 소스는 제외하고, 4–10 keV 플럭스 한계 1.5 × 10⁻¹¹ erg cm⁻² s⁻¹에 해당하는 34,000 deg² 영역을 사용하였다. 각 AGN에 대해 기존 문헌(ASCA, XMM‑Newton, Suzaku 등)에서 얻은 흡수 컬럼 N_H와 전력 지수 Γ를 활용해 4–10 keV 카운트율을 2–10 keV 흡수 보정 광도 L_X로 변환하였다. 흡수 보정 과정에서는 Ω = 2π인 냉각 반사 성분을 포함시켜 실제 스펙트럼을 재현하였다.

N_H 함수는 L_X와 z에 대한 확률 분포 f(L_X,z;N_H)로 정의하고, 20 ≤ log N_H ≤ 24 구간을 4개의 동일 폭 구간으로 이산화하였다. 흡수 비율 ψ(L_X)는 log N_H = 22–24 구간에 속하는 AGN 비율로 설정하고, ψ(L_X) ∝ (L_X)^{-β} 형태를 채택해 광도 의존성을 정량화하였다. 최소화된 likelihood 함수를 통해 ψ(L_X)와 광도 함수 Φ(L_X) 두 파라미터를 동시에 추정했으며, Φ(L_X)는 부드럽게 연결된 이중 전력법칙(γ₁ = 0.84 고정, γ₂ = 2.0 ± 0.2, L_* = 10^{43.3 ± 0.4} erg s⁻¹)으로 모델링하였다.

결과적으로, L_X ≈ 10^{42–43.5} erg s⁻¹ 구간에서는 흡수된 AGN 비율이 0.73 ± 0.25이며, L_X ≈ 10^{43.5–45.5} erg s⁻¹에서는 0.12 ± 0.09로 급격히 감소한다는 기존의 “흡수‑광도 반전” 현상을 재확인하였다. 또한, MAXI/GSC 결과는 HEAO‑1(2–10 keV)과 L_X > 10^{43.5} erg s⁻¹ 구간에서 일치하지만, 낮은 광도(L_X < 10^{43.5} erg s⁻¹)에서는 더 높은 AGN 수밀도를 보인다.

Swift/BAT(14–195 keV)와의 비교에서는 동일한 AGN 샘플에 대해 L_X(2–10 keV)와 L_X(14–195 keV) 사이에 비선형 관계가 존재함을 발견했다. 구체적으로, L_X ≲ 10^{44} erg s⁻¹에서는 평균 스펙트럼 인덱스 Γ ≈ 1.7, L_X ≳ 10^{44} erg s⁻¹에서는 Γ ≈ 2.0으로, 광도에 따라 X선 스펙트럼이 부드러워지는 경향을 보인다. 이는 10 keV 이하와 이상에서 각각 도출된 광도 함수가 겉보기에 상이해 보이는 원인을 설명한다.

이러한 결과는 지역 AGN의 X선 특성이 광도에 따라 변한다는 점을 강조하며, 특히 낮은 광도 AGN가 더 평탄한 스펙트럼을, 높은 광도 AGN가 더 가파른 스펙트럼을 보인다는 점이 AGN 진화 모델에 중요한 제약을 제공한다.