Physically motivated AGN emissivity profiles and their effects on quasar microlensing signatures. 1. Multi-epoch accretion disc size inference

기술 요약: 물리적으로 동기화된 AGN 방출 프로파일과 이것이 퀘이사 마이크로렌징 시그니처에 미치는 영향

문제 제기
퀘이사 마이크로렌징은 마이크로 초 단위의 스케일에서 강착 흐름(accretion flows)의 물리적 크기를 탐사하는 주요 도구이다. 그러나 마이크로렌징 통계로부터 추론된 강착 원반의 크기와 표준 Shakura-Sunyaev (1973) 원반 이론이 예측하는 크기 사이에 지속적인 불일치가 존재하며, 이는 흔히 "원반 크기 문제(disc-size problem)"라고 불린다. 이러한 긴장은 연속체 리버베레이션 매핑(continuum reverberation mapping) 캠페인에서 발견된 결과와도 유사하다. 현재의 마이크로렌징 분석은 대개 원형 가우시안 프로파일이나 표준 얇은 원반 모델과 같은 단순화된 처방을 사용하여 광원 표면 밝기 분포를 모델링한다. 이러한 접근 방식은 광범위한 연속체 방출, 구체적으로는 광역선 영역(Broad-Line Region, BLR)에서 기원하는 확산 자유-결합(free-bound) 연속체 방출의 존재를 간과하는 경우가 많다. 저자들은 (컴팩트한 원반 + 확장된 BLR)의 복합적인 방출 구조를 단일 컴팩트 원반으로 해석하는 것이 추론된 광원 크기에 체계적인 편향(systematic biases)을 유발한다고 주장한다.

방법론
저자들은 대규모 확산 방출이 마이크로렌징 시그니처에 미치는 영향을 정량화하기 위해 물리적으로 동기화된 프레임워크를 구축한다. 방법론은 세 단계로 진행된다:

1. 고유 원반 모델링: 저자들은 에너지적으로 자기 일관적인(self-consistent) 방출 프로파일을 생성하기 위해 agnsed 프레임워크(Kubota & Done 2018)를 변형하여 적용한다. 이들은 두 가지 전형적인 강착 흐름 구조를 모델링한다:

   • 표준 Shakura-Sunyaev (1973) 열적 원반.
   • 표준 외곽 원반과 내부의 뜨거운 X선 플라즈마 사이에 따뜻한 컴프턴화 영역(warm Comptonised region, warm corona)을 갖는 수정된 원반.
     이 모델들은 Novikov-Thorne (1973) 온도 프로파일을 활용하며, 원반 양면을 가열하는 X선 조사(irradiation)를 고려한다. 이 모델들은 전체 스펙트럼 에너지 분포(SED)와 특정 대역폭에 대한 2D 방출 맵을 생성한다.

2. BLR 확산 연속체 계산: 저자들은 cloudy 광이온화 코드(Ferland et al. 2017)를 사용하여 BLR에서 발생하는 확산 자유-결합 연속체 방출을 계산한다. 저자들은 쌍원추형(bi-conical) 기하 구조와 약 $f_c \approx 0.3$의 피복률(covering fraction) 및 광이온화 평형 상태를 가정한다. BLR 방출은 agnsed 모델로부터 얻은 이온화 SED를 기반으로 계산되며, 조명된 BLR 표면으로부터의 등방성 복사를 가정한다.

3. 마이크로렌징 시뮬레이션 및 크기 회복:

   • 증폭 맵(Magnification Maps): 저자들은 양의 패리티(최솟값)와 음의 패리티(안장점) 렌즈 구성 모두에 대해 Fast Multipole Method–Inverse Polygon Mapping (FMM–IPM) 기법을 사용하여 광원 평면 증폭 맵을 생성한다.
   • 컨볼루션(Convolution): (원반와 BLR 성분을 포함하는) 2D 방출 맵을 증폭 맵과 컨볼루션하여 확장된 광원 증폭 분포를 생성한다.
   • 모의 관측: 저자들은 Cosmograil 데이터셋을 기반으로 한 관측 제약 조건(시간 분해능, 노이즈, 공백)을 시뮬레이션하기 위해, 증폭 맵을 가로지르는 무작위 궤적을 추출하여 100개의 모의 마이크로렌징 광도 곡선을 생성한다.
   • 크기 추론: 저자들은 모의 광도 곡선의 변동 통계(마이크로렌징 증폭의 표준 편차, $\sigma_{\Delta m}$)를 시도된 반광도 반지름(half-light radii) 그리드와 비교하여 광원의 크기를 회복한다. 회복 작업은 세 가지 모델링 가정 하에 수행된다: (1) 표준 Shakura-Sunyaev 원반만 있는 경우, (2) 따뜻한 컴프턴화 원반만 있는 경우, (3) 전체 복합 따뜻한 컴프턴 + BLR 모델.

주요 기여 및 결과

• 방출 프로파일 민감도: 본 연구는 컴팩트한 원반의 상세한 방출 프로파일 형태(열적 vs. 따뜻한 컴프턴화)가 마이크로렌징 프로파일에 미치는 영향은 고차 효과(higher-order effect)에 불과하다는 것을 발견했다. 마이크로렌징 신호는 주로 광원의 유효 반광도 반지름에 민감하다.
• 확산 BLR 방출의 영향: 확산 BLR 연속체의 포함은 마이크로렌징 시그니처를 크게 변화시킨다. BLR 방출은 물리적으로 더 큰 규모에서 기원하기 때문에, 소규모 카우스틱(caustic) 네트워크를 평균화하는 평활화(smoothing) 성분 역할을 한다. 이는 컴팩트 전용 모델에 비해 더 좁은 증폭 분포와 감소된 변동 진폭을 초래한다.
• 파장 의존성: 편향의 크기는 관측된 대역폭 내에서 확산 BLR 방출이 전체 SED에서 차지하는 분율에 따라 결정되므로 파장에 따라 강력하게 달라진다. BLR 기여도가 상당한 저에너지 대역(예: $\lambda \approx 7500$ Å, 약 30%)에서는 평활화 효과가 뚜렷하게 나타난다. 반면 BLR 기여도가 최소인 고에너지 대역(예: $\lambda \approx 1650$ Å, 약 1.4%)에서는 그 효과가 미미하다.
• 체계적 크기 편향: 복합(원반 + BLR) 모델로부터 생성된 모의 광도 곡선을 컴팩트 전용 광원 처방을 사용하여 분석할 때, 추론된 반광도 반지름은 실제 컴팩트한 원반의 크기를 체계적으로 과대평가한다. 회복된 크기는 복합 광원의 유효 반광도 반지름을 향하게 된다.
• 편향 메커니즘: 저자들은 이 편향이 단순히 원반와 BLR 반지름의 광도 가중 평균(flux-weighted average)이 아님을 명확히 한다. 대신, BLR 플럭스 분율이 50% 미만인 광학/자외선 영역에서, 추론된 크기는 전체(원반 + BLR) 플럭스의 절반을 포함하는 데 필요한 반지름에 의해 결정된다. BLR은 내부의 컴팩트한 방출에는 크게 기여하지 않으면서 플럭스를 추가하기 때문에, 전체 플럭스의 50%를 포함하는 데 필요한 반지름은 바깥쪽으로 이동하여 겉보기 크기 증가를 유발한다.

의의
본 논문은 "원반 크기 문제"—즉, 마이크로렌징으로 추론된 강착 원반 크기의 체계적인 과잉—의 일부가 복합적인 원반+BLR 방출을 단일 컴팩트 강착 원반의 방출으로 해석하는 데서 기인한다는 결론을 내린다. 확장된 확산 BLR 연속체를 고려하지 않음으로써, 표준 마이크로렌징 분석은 실제 컴팩트한 원반보다 더 큰 "유효" 크기를 측정하게 된다. 이러한 결과는 마이크로렌징과 연속체 리버베레이션 매핑 모두에서 관찰되는 크기 긴장에 대한 물리적 설명을 제공하며, 이 불일치가 반드시 표준 원반 이론의 실패가 아니라 확장된 재처리 영역(reprocessing regions)을 누락한 불완전한 광원 모델링의 결과일 수 있음을 시사한다. 저자들은 자신들의 결과가 물리적으로 동기화된 모델을 통해 편향을 입증하는 수준에서 다소 제한적인 범위임을 강조하지만, 확산 방출의 분율이 크기 추론에 있어 결정적인 매개변수임을 확립하였다.