Population synthesis of active galactic nuclei based on the radiation-regulated unification model

기술 요약: 복사 조절 통일 모델에 기반한 활동은하핵의 인구 합성

문제 제기
활동은하핵 (AGN) 은 초대질량 블랙홀 (SMBH) 에 의해 구동되며, 그 진화는 은하의 특성과 연결되어 있다. X 선 관측은 AGN 의 복사, 차폐, 그리고 강착률에 대한 중요한 제약을 제공하지만, 이전의 인구 합성 모델들은 반사와 흡수를 독립적인 요인으로 취급해 왔으며, 종종 자기 일관적인 물리 모델링 대신 현상론적 템플릿에 의존해 왔다. 또한, 컴프턴 두꺼운 (CTK, $N_H \ge 10^{24} \text{ cm}^{-2}$) AGN 의 고유 비율은 여전히 논쟁의 대상이며, 모델 예측치는 10% 에서 50% 사이로 다양하다. 주요 과제는 차폐와 광도 사이의 반비례 관계를, 방향성만이 AGN 분류를 결정한다고 주장하는 '통일 모델'과 조화시키는 것이다. 최근 제안된 복사 조절 통일 모델은 복사압이 핵주위 환경을 비워냄에 따라 차폐 물질의 덮개 비율이 광도가 아닌 에딩턴 비율 ($\lambda_{\text{Edd}}$) 에 의해 결정된다고 제안한다. 본 연구는 국소 AGN 인구의 진화와 먼지 원환체의 기하학적 구조를 자기 일관적으로 제약하면서, 이러한 복사 조절 모델이 광범위한 X 선 관측 가능량과 일관성을 갖는지 검증하는 것을 목표로 한다.

방법론
저자들은 AGN 인구의 기하학적 및 물리적 특성을 제약하기 위해 시뮬레이션 기반 추론 (SBI) 접근법을 사용한다. 방법론은 다음과 같은 단계를 포함한다:

1. 인공 인구 생성: X 선 광도 함수 (XLF) 에서 표본을 추출하는 대신, 인공 인구는 $z < 0.3$에 대해 Ananna 등 (2022, A22) 이 유도한 고유 국소 활동 블랙홀 질량 함수 (BHMF) 와 에딩턴 비율 분포 함수 (ERDF) 를 표본 추출하여 생성된다. 단순한 진화 처방이 적용되는데, 여기서 활동 SMBH 의 수 밀도는 적색편이 ($e(z)$) 에 따라 진화하는 반면, BHMF 와 ERDF 의 형태는 일정하게 유지된다.
2. 복사 조절 기하학: RefleX 광선 추적 코드를 사용하여 복잡한 기하학적 모델이 구축된다. 이 모델은 다음을 포함한다:
   • $10 r_g$에 위치한 이중 점광원인 '램프-포스트' X 선 원.
   • 평평한 고리 형태의 강착 원반.
   • 고정된 덮개 비율 0.2 를 가진 차가운 먼지가 없는 쐐기 형태로 모델링된 광선 영역 (BLR).
   • 내부 가장자리가 BLR 의 외부 경계 ($R_{\text{BLR,out}}$) 에 설정된 먼지 원환체. 핵심적으로, 원환체의 덮개 비율 ($C_F$) 은 $\lambda_{\text{Edd}}$의 시그모이드 함수로 모델링되어 강착률이 증가함에 따라 감소하며, 이는 복사 조절 모델과 일치한다.
   • 복사 스펙트럼 (고에너지 차단이 있는 전력 법칙) 은 $\lambda_{\text{Edd}}$와의 경험적 관계에 따라 할당된다.
3. 수치 시뮬레이션: RefleX 를 사용하여 저자들은 컴프턴 산란, 광전 흡수, 그리고 형광을 고려하여 핵주위 물질을 통과하는 X 선 광자 (0.1 keV – 1 MeV) 에 대한 몬테카를로 시뮬레이션을 수행한다. 기하학적 및 스펙트럼 매개변수의 격자 전반에 걸쳐 $6 \times 10^6$개 이상의 스펙트럼이 생성된다.
4. 추론: 모델 매개변수는 Sequential Neural Posterior Estimation (SNPE) 방법을 사용하여 시뮬레이션된 인구를 관측 제약 조건과 비교함으로써 추론된다. 제약 조건에는 다음이 포함된다:
   • 우주 X 선 배경 (CXB) 스펙트럼 (2–100 keV).
   • 세 가지 에너지 대역 (14–195 keV, 8–24 keV, 2–10 keV) 에서의 미분 AGN 개수.
   • 관측된 국소 AGN 의 $N_H$ 분포.
   • $\lambda_{\text{Edd}}$에 따른 차폐된 AGN 의 비율.
   • 제한 플럭스에 따른 관측된 CTK AGN 의 비율.

주요 기여

• 자기 일관적 모델링: 본 연구는 고유 BHMF 와 ERDF 를 직접 표본 추출하고, 차폐와 반사가 공간 분포를 통해 본질적으로 연결된 통일된 프레임워크 내에서 에딩턴 비율 의존 기하학 및 복사 스펙트럼을 구현한 최초의 인구 합성 연구이다.
• 복사 조절 검증: 본 연구는 복사 조절 통일 모델을 엄격하게 검증하여, 원환체 덮개 비율이 $\lambda_{\text{Edd}}$에 의존하는 기하학이 CXB, 개수, 그리고 흡수 특성을 동시에 재현할 수 있음을 보여준다.
• 고유 CTK 비율: 저자들은 $40 \pm 3\%$의 견고한 고유 CTK 비율을 도출했는데, 이는 이전 추정치보다 현저히 높으며, 이는 BLR 차폐를 포함함으로써 설정된 약 $20\%$의 기본 CTK 비율에 기인한다.
• 매개변수 제약: 본 연구는 원환체 열량 분포 ($\mu, \sigma$), $\lambda_{\text{Edd}}$에 따른 원환체 덮개 비율의 의존성 ($k$), 그리고 활동 SMBH 비율의 진화 ($p_1, z_c$) 를 포함한 주요 물리 매개변수에 대한 사후 분포를 제공한다.

결과

• CXB 및 개수: 모델은 60 keV 이상 플럭스를 약간 과대평가하지만, CXB 스펙트럼의 전체 강도와 형태를 성공적으로 재현한다. 미분 개수는 불확실성 범위 내에서 재현되며, 14–195 keV 대역의 밝은 끝에서 약간의 과대평가와 더 부드러운 대역의 어두운 끝에서 약간의 과소평가가 관찰된다.
• 흡수 특성: 모델은 이모달 $N_H$ 분포와 차폐 비율과 $\lambda_{\text{Edd}}$ 사이의 반비례 관계를 재현한다. 그러나 관측된 차폐 및 비차폐 AGN 의 절대 수는 약 15–25% 과대평가되는데, 이는 인공 인구에서 완전히 포착되지 않은 선택 효과 때문일 가능성이 높다.
• 원환체 기하학: 추론된 원환체의 내부 가장자리는 $-3 < \log \lambda_{\text{Edd}} < 0.5$에 대해 약 0.5 - 8 pc 범위를 가지며, 이는 중적외선 간섭계 측정치와 복사 조절 시나리오와 일치한다. 덮개 비율은 증가하는 $\lambda_{\text{Edd}}$에 따라 감소하여 국소 AGN 에서 관찰된 경향을 회복한다.
• 반사 - 차폐 상관관계: 모델은 차폐된 원천에서 $N_H$가 증가함에 따라 반사가 증가한다는 관측된 경향을 질적으로 재현하며, 비차폐 AGN 에 대한 반사 수준을 이전 모델들보다 관측과 더 잘 일치하도록 예측한다.
• 진화: 추론된 활동 SMBH 비율의 진화는 $z \approx 1$까지 급격히 증가한 후 정체되는 모습을 보여주며, 이는 XLF 연구에서 사용되는 순수 밀도 진화 모델과 일치한다.

의의 및 주장
본 논문은 복사 조절 통일 모델이 매우 포괄적인 X 선 관측 가능량과 일관성이 있음을 주장한다. 현상론적 모델링에서 물리 모델링으로 전환함으로써, 저자들은 활동 SMBH 비율을 조절하는 단순한 진화 처방이 국소적으로 표본 추출된 BHMF 와 ERDF 와 결합될 때 CXB 를 재현하기에 충분함을 보여준다. 도출된 고유 CTK 비율인 $40 \pm 3\%$는 최근의 국소 우주에 대한 NuSTAR 관측 및 이전 인구 합성 예측과 잘 일치한다. 연구는 차폐 및 비차폐 AGN 간의 BHMF 와 ERDF 의 차이가 단일 부모 분포를 표본 추출함으로써 재현될 수 있으며, 그 차이는 오직 차폐를 지배하는 물리적 및 기하학적 특성에서 비롯됨을 결론짓는다. 저자들은 균질한 구성 요소에 대한 가정과 다른 잠재적 차폐체 (예: 극성 먼지, 유출물) 의 부재와 같은 한계를 지적하며, 이를 향후 개선 영역으로 제안한다.