Study of Flat Spectrum Radio Quasars and BL Lacertae Objects as Sources of Diffusive Ultra High-Energy Cosmic Rays

기술적 요약: 확산 초고에너지 우주선의 원천으로서 평탄 스펙트럼 전파 퀘이사와 BL Lacertae 천체의 연구

문제 제기
초고에너지 우주선 (UHECRs, $E \gtrsim 10^{17}$ eV) 의 기원은 현대 천체물리학에서 여전히 해결되지 않은 중요한 문제입니다. 스펙트럼 특징은 은하 외 기원을 시사하지만, 은하 및 은하 외 자기장에 의한 편향으로 인해 정확한 원천을 식별하기 어렵습니다. 피에르 오제 관측소 (Pierre Auger Observatory) 와 텔레스코프 어레이 (TA) 의 최근 고정밀 데이터는 특정 스펙트럼 특징 (발목 부위에서의 경화, 급격한 감소, GZK 억제) 과 조성 경향 (에너지 증가에 따른 가벼운 원자핵에서 무거운 원자핵으로의 전이) 을 드러냈습니다. 인근 전파 은하가 강력한 후보이지만, 플럭스, 조성, 비등방성에 관한 관측적 제약 조건에 대해 엄격하게 검증해야 할 활성 은하핵 (AGN), 특히 평탄 스펙트럼 전파 퀘이사 (FSRQs) 와 BL Lacertae (BL Lac) 천체의 잠재적 기여도가 필요합니다.

방법론
저자들은 표준 $\Lambda$CDM 우주론 내의 포괄적인 전파 모델을 사용하여 FSRQ 와 BL Lac 를 확산 UHECR 의 가능한 원천으로 조사합니다.

1. 원천 진화 모델:

   • 본 연구는 페르미-LAT 데이터에서 도출된 관측된 감마선 광도 함수 (LF) 에 기반한 현실적인 광도 의존 밀도 진화 (LDDE) 함수를 활용합니다.
   • FSRQ: 원천 밀도의 적색편이 피크가 광도에 의존 ($z_c \propto L^\alpha$) 하는 LDDE 함수를 사용하여 모델링됩니다.
   • BL Lac: 광도 상관관계를 가진 광자 지수의 고유 분포 (가우시안) 와 동일한 LDDE 형식을 결합하여 모델링됩니다.
   • 이러한 모델들은 $(1+z)$ 의 멱함수 등 단순한 현상론적 진화 형태를 관측된 감마선 데이터에 기반한 물리적으로 동기를 부여한 진화로 대체합니다.

2. 전파 물리:

   • UHECR 은 우주 마이크로파 배경 (CMB) 및 은하 외 배경 광 (EBL) 과의 상호작용으로 인한 에너지 손실을 고려하여 은하 외 공간을 통해 전파됩니다.
   • 포함된 과정: 베트 - 하이틀러 쌍생성, 광파이온 생성 (양자의 경우), 광분열 (원자핵의 경우). 본 연구는 광분열에 대해 연속 에너지 손실 (CEL) 근사를 사용하여 2 차 캐스케이드를 명시적으로 풀지 않고 1 차 원자핵 (p, He, N, Si, Fe) 의 생존 확률을 추적합니다.
   • 자기 확산: 전파는 RMS 강도 $B_{rms} = 1$ nG 와 결맞음 길이 $l_c = 0.5$ Mpc 인 균일한 난류 은하 외 자기장 (IGMF) 내에서 발생합니다. 에너지 의존 확산 계수 $D(E)$ 는 콜모고로프 난류에 기반하여 계산되며, 공명 및 비공명 확산 영역이 분리됩니다.

3. 플럭스 계산 및 피팅:

   • 확산 플럭스는 적색편이와 광도에 걸쳐 원천 방출률을 적분하여 계산하며, 균일하게 분포된 합성 원천들의 기여를 합산하기 위해 전파 정리를 적용합니다.
   • 모델은 자유 매개변수로 취급된 분율을 가진 혼합 원자핵 조성 (p, He, N, Si, Fe) 을 포함합니다.
   • 원천 방출률 정규화는 바리온 로딩 인자 $\kappa = L_{CR}/L_\gamma$ 를 통해 감마선 광도와 연결됩니다.
   • 매개변수 (원천 분리 거리 $d_s$, 가속 효율 $\eta$, 원자핵 분율 $\epsilon_Z$) 는 Auger 와 TA 의 관측 데이터 (플럭스 및 평균 샤워 깊이 $\langle X_{max} \rangle$) 에 맞춰 피팅됩니다.
   • 통계적 모델 비교는 축소된 $\chi^2$, 아카이케 정보 기준 (AIC), 베이지안 정보 기준 (BIC) 을 사용하여 수행됩니다.
   • 적분 플럭스 ($10-100$ EeV) 와 쌍극자 비등방성의 전천도 (Full-sky maps) 는 HEALPix 프레임워크를 사용하여 생성됩니다.

주요 결과

1. 플럭스 및 조성:

   • FSRQ 와 BL Lac 시나리오 모두 혼합 조성을 사용하여 관측된 UHECR 에너지 스펙트럼과 $\langle X_{max} \rangle$ 데이터를 재현할 수 있습니다.
   • 최적 피팅 가속 효율은 $\eta \approx 0.1$ 로, 상대론적 제트 가속의 수치 시뮬레이션과 일치합니다.
   • 필요한 바리온 로딩 인자는 $\kappa \approx 5.65$ 로, 우주선 (CR) 파워가 에너지적으로 실현 가능함을 나타냅니다 (감마선 광도의 몇 배 수준).
   • FSRQ 시나리오: 플럭스를 피팅하기 위해 원천 분리 거리 $d_s \approx 8 \pm 1$ Mpc 가 필요합니다. 이는 페르미-LAT 광도 함수에서 도출된 관측적 추정치보다 훨씬 높은 원천 밀도를 의미합니다.
   • BL Lac 시나리오: $d_s \approx 160 \pm 21$ Mpc 가 필요하며, 이는 BL Lac 밀도에 대한 관측적 추정치와 일치합니다.

2. 통계적 비교:

   • 결합된 Auger 와 TA 데이터에 대한 축소된 $\chi^2$ 값은 FSRQ 의 경우 2.37, BL Lac 의 경우 2.55 입니다.
   • 모델 선택 기준 (AIC 및 BIC) 은 BL Lac 모델을 선호합니다 ($\Delta$AIC = 4.17, $\Delta$BIC = 4.83). 이는 BL Lac 가 FSRQ 보다 데이터를 더 잘 설명한다는 중간 정도의 증거를 제공합니다.

3. 비등방성 제약:

   • BL Lac: 플럭스 피팅 매개변수 ($z_{max}=0.5$) 를 사용하여 예측된 쌍극자 비등방성은 Auger 관측과 잘 일치합니다.
   • FSRQ: 동일한 매개변수를 사용할 때 예측된 쌍극자 진폭은 관측치보다 현저히 높습니다. 관측된 비등방성과 일치시키기 위해서는 원천 집단을 인위적으로 저적색편이 영역 ($z_{max} \approx 0.01$) 으로 제한해야 하는데, 이는 플럭스 피팅 요구 사항과 FSRQ 의 잘 알려진 강한 우주론적 진화와 모순됩니다.

4. 광분열:

   • 광학 두께와 흡수된 분율의 계산은 무거운 원자핵 (예: Fe) 이 가벼운 원자핵 (예: He) 에 비해 더 낮은 에너지에서 상당한 광분열을 겪음을 보여주며, 이로 인해 최고 에너지에서 무거운 성분이 점진적으로 억제됩니다.

의의 및 결론
본 논문은 FSRQ 와 BL Lac 모두 관측된 UHECR 플럭스를 재현할 수 있지만, 관측적 제약 조건 전체를 고려할 때 BL Lac 시나리오가 강력하게 선호됨을 결론지었습니다.

• FSRQ 는 확산 UHECR 의 지배적 원천으로서 불리합니다. 플럭스를 피팅하기 위해 비현실적으로 높은 원천 밀도 (작은 $d_s$) 가 필요하며, 원천 집단을 저적색편이에서 인위적으로 단절하지 않는 한 관측 한계를 초과하는 쌍극자 비등방성을 예측하기 때문입니다.
• BL Lac 는 관측과 일관성을 유지합니다. 플럭스, 조성, 쌍극자 비등방성에 대한 자기 일관된 설명을 제공하며, 이는 페르미-LAT 데이터와 일치하는 원천 밀도 및 적색편이 진화를 가집니다.

저자들은 FSRQ 가 다중 메신저 채널 (중성미자 및 감마선) 에 여전히 기여할 수 있지만, 확산 UHECR 플럭스에는 BL Lac 천체에 비해 하위 기여자일 가능성이 높다고 제안합니다. 본 연구는 UHECR 기원을 제약할 때 단순화된 현상론적 형태 대신 현실적이고 광도 의존적인 원천 진화 모델을 사용하는 것의 중요성을 강조합니다.