X‑선 이진성, 우주선 ‘무릎’ 스펙트럼을 지배하는 새로운 파워플랜트

초록

본 논문은 은하계에 존재하는 수백 개의 X‑ray binary( XRB )가 PeV 수준의 입자를 가속해 ‘무릎(knee)’ 에너지대의 우주선 스펙트럼에 주요 기여를 할 수 있음을 제시한다. 저자들은 최신 XRB 카탈로그를 기반으로 은하 중심에서의 방사형 분포 함수를 도출하고, 이를 DRAGON 전파 모델에 삽입해 에너지 스펙트럼, 이방성, 그리고 확산 감마선 방출을 계산한다. 결과는 XRB가 PeVatrons 로서 SNR보다 더 큰 비중을 차지할 가능성을 보여주며, 향후 LHAASO·HAWC 관측을 통해 검증될 수 있음을 강조한다.

상세 분석

본 연구는 X‑ray binary( XRB )를 은하계 PeVatrons 로 설정하고, 그 집합적 기여가 ‘무릎(knee)’ 영역의 우주선 스펙트럼을 지배할 수 있음을 정량적으로 검증하려는 시도이다. 먼저 저자들은 Neumann et al.(2023)와 Avakyan et al.(2023)에서 제공된 172개의 HMXB와 360개의 LMXB 데이터를 사용해 은하 평면에 투영하였다. 거리‑의존 선택 효과를 보정하기 위해 Sun‑GC 중심으로 동일 간격의 원형 격자를 구축하고, 각 격자 셀에서 관측된 XRB 표면 밀도를 구해 가중 평균을 취하였다. 이 과정에서 Poisson 오차를 고려한 σρij = √Nij/Nij 식을 적용했으며, 최종 방사형 분포는 ρ(R)=A·exp(−aR/R⊙)+B·exp(−bR²/R⊙²) 형태로 피팅되었다. 피팅 파라미터 a=2.67±0.21, b=39.33±3.88, A=1.54±0.29, B=9.74±1.99는 XRB가 은하 중심에 집중되면서도 디스크 외곽까지 확산된 구조를 나타낸다.

이후 DRAGON 전파 코드에 위 분포 함수를 source term으로 삽입하고, 표준 확산 파라미터(Dxx∝β·(R/4 GV)δ, δ≈0.33)와 대류·재가열 항을 포함한 3D 전파 모델을 구축하였다. XRB가 방출하는 입자 스펙트럼은 파워‑law with exponential cutoff 형태(E−γ exp(−E/Ecut))를 가정하고, γ≈2.2, Ecut≈3 PeV 로 설정하였다. 결과적으로 얻어진 전체 우주선 스펙트럼은 1 PeV 이하에서는 전통적인 SNR 기여와 거의 겹치지만, 3 PeV 근처에서 XRB 기여가 50 % 이상을 차지한다는 점이 핵심이다. 또한, 모델은 관측된 우주선 이방성(특히 10‑30 TeV 대역에서의 작은 위도 편차)과 일치하도록 조정되었으며, 확산 감마선 배경 역시 LHAASO가 보고한 확산 감마선 플럭스와 좋은 일치를 보였다.

하지만 몇 가지 한계점도 존재한다. 첫째, XRB의 가속 효율 ηCR을 10 % 정도로 고정했는데, 실제 제트 물리와 충돌 가스 밀도에 따라 크게 변동할 수 있다. 둘째, HMXB와 LMXB를 동일한 스펙트럼 파라미터로 취급했으며, 이는 각각의 제트 구조와 전자‑양성자 비율 차이를 무시하는 것이다. 셋째, DRAGON 내부의 은하 자기장 모델은 단순화된 2‑scale 파라미터화에 머물러 있어, 실제 복잡한 나선 팔 구조와의 상호작용을 충분히 반영하지 못한다. 마지막으로, 관측된 XRB 수가 실제보다 낮을 가능성이 있는데, 은하 중심의 은폐와 감도 제한으로 인해 수백 개에서 수천 개까지 추정치가 크게 차이한다. 이러한 불확실성을 감안하면, XRB가 ‘무릎’ 영역을 지배한다는 결론은 아직 가설 단계에 머물며, 향후 LHAASO·HAWC의 장기 모니터링과 고해상도 제트 시뮬레이션이 필요하다.