근접 마이크로퀘이사와 우주선 무릎 현상의 새로운 해석

초록

LHAASO가 100 TeV 이상 γ선을 방출하는 다섯 개 마이크로퀘이사를 발견했지만, 이들만으로는 우주선 스펙트럼의 ‘무릎’(PeV) 영역을 설명할 수 없다. 저자들은 기존 마이크로퀘이사들의 기여를 정량화하고, 반대 은하 중심 방향으로 약 2.6 kpc 떨어진 미확인 마이크로퀘이사(또는 유사 시스템)가 단일 근접원으로 작용하면 관측된 양성자·헬륨·전체 입자 스펙트럼과 이방성(진폭·위상)을 동시에 재현할 수 있음을 제시한다.

상세 분석

본 논문은 LHAASO가 보고한 다섯 개 마이크로퀘이사(SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070, Cygnus X‑1) 중 네 개가 100 TeV 이상 γ선을 방출한다는 사실을 출발점으로 삼는다. 저자들은 먼저 전통적인 초신성 잔해(SNR) 기반의 배경 우주선 모델을 구축한다. 확산 계수 D(R)=D₀ β (R/R₀)^δ와 파워‑브레이크 형태의 주입 스펙트럼(Q∝R^{-ν₁} for R<R_br, ∝R^{-ν₂} e^{-R/R_cut} for R≥R_br)을 사용해 B/C 비율과 AMS‑02, DAMPE, Voyager 등 관측 데이터를 재현한다. 이 과정에서 δ≈0.33, D₀≈2×10²⁸ cm² s⁻¹ 등 표 1의 파라미터가 채택되었다.

다음 단계에서는 LHAASO γ선 측정치를 이용해 네 개 마이크로퀘이사의 개별 입자 주입 스펙트럼을 역추정한다. γ‑ray 데이터는 주로 양성자‑핵 충돌에 의한 중성파이온 붕괴에서 유래한다는 가정 하에, 각각의 총 입자 출력 L₀와 스펙트럼 지수 β를 조정해 모델 γ선과 관측 γ선을 일치시켰다(표 2). 결과적으로 각 마이크로퀘이사의 기여는 전체 우주선 플럭스의 1 % 이하에 불과함을 확인한다.

핵심 가설은 ‘무릎’ 영역(∼4 PeV)에서의 우주선이 단일 근접원, 즉 미확인 마이크로퀘이사에 의해 지배된다는 것이다. 이를 위해 저자들은 순간 주입 스펙트럼 q(R)=q₀ (R/R₀)^{-β} e^{-R/R_cut,loc}를 도입하고, 확산 반경 σ(R,t)=√(2 D(R) t)를 이용해 지구에 도달한 플럭스를 계산한다. 파라미터 탐색은 거리 d, 나이 t, 절단 강성 R_cut,loc, 스펙트럼 지수 β를 포함한다. 관측된 양성자·헬륨·전체 입자 스펙트럼과 동시에 이방성 진폭·위상(반 은하 중심 방향을 가리키는 10⁻³ 수준)을 재현하도록 제약한다.

베이지안 MCMC 결과는 최적 거리 d≈2.6 kpc, 나이 t≈4 × 10⁵ yr, R_cut,loc≈5 PV, β≈2.2 정도가 가장 높은 사후 확률을 보였으며, 이 위치는 반 은하 중심 방향에 해당한다. 이러한 파라미터 집합은 무릎 근처 스펙트럼의 급격한 경사 변화와 질량 구성 감소 현상을 동시에 설명한다. 또한, 이방성 위상이 관측된 에너지 구간(∼10 TeV–1 PeV)과 일치한다는 점에서 모델의 일관성을 뒷받침한다.

논문의 한계점으로는 (1) 마이크로퀘이사 주변 환경(가스 밀도, 자기장 구조)의 불확실성, (2) 단일 근접원 가정이 실제로는 복수의 미확인 소스가 겹쳐 있을 가능성, (3) 확산 계수가 공간·시간에 따라 변할 수 있다는 점을 들 수 있다. 그럼에도 불구하고, 관측 데이터와 이론 모델을 다중 메신저(γ‑ray, 우주선 스펙트럼, 이방성)로 결합한 접근법은 무릎 현상의 새로운 해석을 제시한다는 점에서 의의가 크다.