On gamma rays as predictors of UHECR flux in AGNs

기술 요약: AGN의 $\gamma$선에 대한 UHECR 플럭스 예측 연구

문제 제기
활동 은하핵(AGN)은 초고에너지 우주선(UHECR)의 주요 후보입니다. 그러나 관측된 매우 높은 에너지(VHE) $\gamma$선 플럭스를 UHECR 광도($L_{CR}$)의 직접적인 대리 지표(proxy)로 사용하는 최근의 현상론적 모델들은 피에르 오제 관측소(Pierre Auger Observatory)의 실험 데이터와 일치하는 데 실패하고 있습니다. 구체적으로, 이러한 모델들은 에너지 스펙트럼, 질량 조성, 그리고 도착 방향 지도를 동시에 맞추는 데 어려움을 겪습니다. 결정적인 실패 지점은 블레이자(blazar)인 Mkn 421의 기여도를 과대평가하여, 관측 결과와 일치하지 않는 핫스팟과 쌍극자(dipole) 방향을 생성한다는 점입니다. 또한, 이러한 모델들은 서로 다른 AGN 클래스(예: 라디오 은하와 블레이자)를 동일한 분석에 혼합할 때 "가짜(spurious)" 쌍극자 방향과 좋지 못한 스펙트럼 형태 적합도를 나타내는 경우가 많습니다.

방법론
저자들은 이러한 불일치가 이전 연구들의 암묵적인 가정, 즉 UHECR 방출이 상대론적 제트 내의 $\gamma$선 방출과 동일하게 비밍(beamed)된다는 가정에서 비롯된다고 제안합니다. 논문은 $\gamma$선이 제트 내의 "덩어리(blobs)"에서의 상대론적 운동으로 인해 강한 도플러 부스팅($D^4$)의 영향을 받는 반면, UHECR은 소스 환경을 탈출하기 전에 등방화(isotropization, decollimation) 과정을 거칠 가능성이 높다고 주장합니다.

본 연구는 다음 단계로 진행됩니다:

1. 이론적 검토: 저자들은 AGN 제트에서의 가속 메커니즘(충격파, 재결합, 전단 가속)과 제트 로브(lobe)를 통한 입자의 전파를 검토합니다. 그들은 제트 로브와 소스 코쿤(cocoon) 내의 자기장이 UHECR을 산란시켜, 비록 이들이 상대론적 흐름에서 기원할지라도 방출을 효과적으로 등방화시킨다고 주장합니다.
2. 대리 지표 교정: 저자들은 관측된 $\gamma$선 광도($L^{obs}_{\gamma}$) 대신 고유 $\gamma$선 광도($L^{int}_{\gamma}$) 또는 복사 제트 출력(radiative jet power)을 $L_{CR}$의 대리 지표로 사용할 것을 제안합니다. 이는 관측된 플럭스를 도플러 인자($D$)의 $q$ 승으로 보정하는 과정을 포함합니다 ($L^{int}_{\gamma} = D^{-q}L^{obs}_{\gamma}$). 두 가지 시나리오가 테스트됩니다: $q=4$ (assuming $L_{CR} \propto L^{int}_{\gamma}$)와 $q=2$ (assuming $L_{CR} \propto$ radiative jet power).
3. 시뮬레이션 및 피팅: 저자들은 CRPropa3 프레임워크를 사용하여 다섯 가지 주요 핵종(H, He, N, Si, Fe)에 대해 에너지 및 거리 격자에 걸쳐 1차원 전파 시뮬레이션을 수행합니다. 그들은 시뮬레이션된 플럭스를 피에르 오제 데이터(스펙트럼 및 $X_{max}$ 모멘트) 및 도착 방향 지도와 비교합니다.
4. 비교 분석: 세 가지 케이스가 비교됩니다:
   • Case (a): $L^{obs}_{\gamma}$를 대리 지표로 사용 (표준 접근법).
   • Case (b): $D^{-2}L^{obs}_{\gamma}$를 대리 지표로 사용.
   • Case (c): $D^{-4}L^{obs}_{\gamma}$를 대리 지표로 사용.

주요 기여 및 결과

• 블레이자 지배력의 억제: 교정을 통해 (Cen A, M87, Fornax A와 같은) 라디오 은하에 비해 (Mkn 421과 같이) 고도로 비밍된 블레이자의 가중치를 크게 줄였습니다. 표준 모델에서는 Mkn 421이 플럭스를 지배하지만, 교정된 모델에서는 그 기여도가 2~5 차수(orders of magnitude) 감소하여 무시할 수 있는 수준이 됩니다.
• 개선된 스펙트럼 및 조성 적합도: 고유 광도 대리 지표를 사용할 때 에너지 스펙트럼과 조성 데이터의 결합 적합도가 크게 향상됩니다. 환산된 $\chi^2/NDF$ 값은 $L^{obs}_{\gamma}$를 사용할 때의 4.6에서 $D^{-2}L^{obs}_{\gamma}$를 사용할 때의 3.1로 감소합니다. 이러한 개선은 멀리 떨어진, 고도로 비밍된 소스의 과잉 기여로 인한 스펙트럼 왜곡을 피함으로써, 근거리 소스와 원거리 소스의 상대적 기여를 더 현실적으로 반영했기 때문입니다.
• 쌍극자 긴장 해소: 모델링된 쌍극자 방향과 관측된 방향 사이의 긴장이 완화되었습니다. $E > 8$ EeV 에너지 영역에서 오제 관측 영역의 쌍극자 불일치는 5.9$\sigma$에서 3.5$\sigma$로 감소했습니다. $E > 32$ EeV 영역에서는 2.1$\sigma$에서 1.1\sigma}로 감소했습니다.
• 핫스팟 재정렬: 도착 방향 지도는 Mkn 421에 의해 지배되던 상태에서 Cen A, Fornax A, M87 및 특정 원거리 블레이자(Mkn 180, 1ES2344+514) 주변의 핫스팟을 보여주는 형태로 이동했습니다. 이 새로운 핫스팟들은 피에르 오제 관측소가 보고한 두 개의 핫스팟 및 텔레스코프 어레이(Telescope Array)가 보고한 하나의 핫스팟과 더 잘 일치합니다.

의의 및 주장
본 논문은 $\gamma$선과 UHECR의 비밍 가정을 분리함으로써 기존의 AGN 기반 UHECR 소스 모델에서 발견된 모순들에 대한 해결책을 제시한다고 주장합니다. 저자들은 다음과 같이 단언합니다:

1. $\gamma$선과 UHECR의 연관성은 약하지만 존재합니다. $\gamma$선은 입자 가속 영역을 추적하지만, UHECR은 탈출하기 전 소스의 자기장 로브에 의해 등방화될 가능성이 높습니다.
2. 관측된 $\gamma$선 플럭스를 대리 지표로 사용하는 것은 UHECR이 비밍되어 있다는 가정을 내포하며, 저자들은 이것이 물리적으로 일어날 가능성이 낮다고 주장합니다.
3. 도플러 부스팅을 교정하는 것($L^{int}_{\gamma}$ 또는 복사 제트 출력)은 스펙트럼, 조성, 그리고 도착 방향 데이터를 일치시키는 더 견고한 대리 지표를 제공합니다.

저자들은 자신들의 모델이 단순화된 결합 피팅(전체 $X_{max}$ 분포나 상세한 자기장 편향을 고려하지 않음)임에도 불구하고, 결과적으로 AGN이 UHECR 소스라는 가설을 강화하며, 다양한 AGN 클래스에 대한 다중 신호(multi-messenger) 데이터를 해석하는 데 있어 더 일관된 틀을 제공한다고 결론짓습니다. 그들은 향후 연구에서 에너지 의존적 탈집속(decollimation)과 정규 은하 및 외계 은하 자기장의 구체적인 효과를 다루어야 한다고 언급했습니다.