A Boundary-Consistent Two-Zone Electron Kernel for Distant Pulsar Contributions to Positron Flux and Anisotropy

기술 요약: 원거리 펄서의 양전자 플럭스 및 이방성 기여를 위한 경계 일관적 2-존 전자 커널

문제 제기
국부 우주선 전자 및 양전자 스펙트럼은 인근 고에너지 가속기, 특히 펄사를 조사하는 정밀한 탐침 역할을 한다. 제밍가(Geminga)나 모노젬(Monogem)과 같이 성숙한 펄서들은 TeV 규모의 렙톤 플럭스의 자연스러운 후보이지만, 핵심적인 질문이 남아 있다: 국부 렙톤 플럭스는 반드시 가장 가까운 펄서들에 의해 지배되어야 하는가, 아니면 더 먼 펄서들이 집단적으로 기여할 수 있는가? 이 질문은 두 가지 요인에 의해 복잡해진다: (1) 급격한 복사 냉각(싱크로트론 및 역 컴프턴 산란)은 고에너지 입자의 전파 지평을 제한하며, (2) Fermi-LAT와 AMS-02의 이방성 한계는 단일 근접 소스가 너무 넓게 지배하는 모델을 제약한다. 기존 연구들은 먼 소스들이 저에너지 전자를 재현할 수 있다고 시사했으나, 근접 소스의 느린 확산(slow diffusion)이 존재하는 상황에서 기하학적 소스 개수 계수(geometric source-count factor)를 정량화할 수 있는 안정적인 반-해석적 프레임워크는 결여되어 있었다. 더욱이 은하계의 가장 안쪽 수백 파섹(Local Bubble)은 낮은 분자 가스 밀도와 펄서의 출생 시 킥(natal kick)으로 인해 연속적인 펄서 밀도를 가정하기 어려운 물리적 불연속성을 나타낸다.

방법론
저자들은 구형 2-존 확산 모델 내에서 전자 및 양전자의 전파를 위한 반-해석적 급수 해(series solution)를 개발한다.

• 2-존 커널: 모델은 소스를 둘러싼 반지름 $R_1$의 느린 확산 영역(확산 계수 $D_1(E)$)과 표준 성간 영역($R_1 < r < R_2$, 확산 계수 $D_2(E)$)을 정의한다. 두 계수는 동일한 에너지 의존성($D \propto f_D(E)$)을 공유하여, 해석적인 고유함수 해(eigenfunction solutions)를 가능하게 한다.
• 경계 조건: 해는 $R_1$에서의 밀도 및 확산 플럭스의 연속성을 보장하며, $R_2$(20 kpc로 설정)에서의 흡수 경계를 적용한다. 이 정식화는 2-존 모델의 직접 적분 평가에서 발생하는 진동적 상쇄(oscillatory cancellations)를 피하여, 날카로운 복사 냉각 경계를 보존한다.
• 냉각: 복사 손실(CMB, 적외선, 광학 장에 대한 싱크로트론 및 클라인-니시나 역 컴프턴 산란)은 에너지 특성(energy characteristics)을 통해 처리되며, 이는 주입 에너지와 소스의 연령을 연결한다.
• 집단 모델링: 커널은 연속적인 디스크 형태의 펄서 집단에 적용된다. 저자들은 Local Bubble과 펄서의 고유 운동을 고려하여 $R < 0.1$ kpc에 빈 공동(empty cavity)을 설정한다.
• 통계 분석: AMS-02 양전자 플럭스 데이터, CALET 전전자 쌍극자 이방성 측정값, 그리고 Fermi-LAT 이방성 상한치를 피팅하기 위해 "솔버-인-더-루프(solver-in-the-loop)" MCMC 분석이 수행된다. 피팅에는 연속적인 스핀다운 주입 이력, 깨진 거듭제곱 법칙(broken-power-law) 쌍 주입 스펙트럼, 그리고 이차적 성격의 양전자 기저와 매끄러운 거듭제곱 법칙 배경에 대한 비음수 선형 진폭에 대한 프로파일이 포함된다.

핵심 기여

1. 안정적인 급수 해: 본 논문은 복사 냉각을 포함하는 2-존 확산 문제에 대한 수치적으로 안정적인 반-해석적 급수 해를 제공한다. 이 방법은 계산 효율적이며, 특히 중간 거리 영역에서 브루트 포스 유한 부피 방법(finite-volume methods)이나 진동 적분이 갖는 수치적 불안정성을 피한다.
2. 기하학적 정량화: 본 연구는 "쉘 계수(shell factor, 소스 개수 대 거리)"가 전파 커널과 어떻게 상호작용하는지 정량화한다. 3D 균일 분포는 10–100 GeV에서 킬로파섹 규모의 쉘에 의해 지배되는 반면, 디스크 기하 구조($h=0.2$ kpc)는 기여를 안쪽으로 이동시키지만, 1 kpc 너머의 소스들이 여전히 10–100 GeV 플럭스의 37–47%를 공급한다는 것을 입증한다.
3. 느린 헤일로(Slow Halo)의 역할: 저자들은 느린 확산 헤일로의 포함이 1-존 모델에 비해 먼 소스의 플럭스를 강화하는 것이 아니라, 오히려 모든 거리로부터 오는 절대적 플럭스를 감소시키는 에너지 의존적 냉각 페널티를 도입하며, 고에너지에서 더 강한 억제를 일으킨다는 점을 명확히 한다.
4. 이방성 제약: 본 연구는 CALET 쌍극자 진폭과 Fermi-LAT 상한치를 사용하여 공통 헤일로 파라미터의 타당성을 테스트함으로써 이방성 제약을 모델에 직접 통합한다.

결과

• 플럭스 분해: (R < 0.1 kpc를 제외한) 피팅된 연속 펄서 집단에서, 0.3 kpc 너머의 소스들의 기여는 모든 에너지 범위(10 GeV ~ TeV 스케일)에서 국부 0.1–0.3 kpc 고리(annulus)의 기여를 압도한다. 10 GeV에서 $R > 0.3$ kpc의 기여는 국부 고리보다 약 10배 더 크다.
• 헤일로 파라미터: MCMC 분석 결과, 느린 헤일로 반지름 $R_1$의 사후 분포 중심은 $0.0709^{+0.0510}_{-0.0407}$ kpc ($\approx 71$ pc)이며, 확산 계수는 $D_1(100 \text{ TeV}) \approx 2.08 \times 10^{28}$ cm$^2$ s$^{-1}$이다. 이 결과는 제밍가와 모노젬의 TeV 헤일로 형태학 연구와 일치한다.
• 이방성 호환성: 도출된 헤일로 파라미터를 가진 피팅 모델은 현재의 이방성 제한(CALET 및 Fermi-LAT)과 호환성을 유지한다. 쌍극자 진폭 예측은 95% 상한치 범위 내에 있다.
• 에너지 의존성: 먼 소스들이 서브-TeV 플럭스를 지배하지만, 복사 냉각으로 인해 기여가 TeV 스케일 근처로 국지화된다. 1 TeV에서 디스크 모델의 플럭스 중 약 90%는 1 kpc 이내에서 기인한다.

의의 및 주장
본 논문은 먼 펄서(즉각적인 Local Bubble 근처를 벗어난)의 집단적 기여가 수백 GeV 미만의 양전자 및 전자 플럭스에 대해 타당하고 지배적인 기여자가 될 수 있음을 결과로 보여준다. 피팅된 펄서 성분은 가장 가까운 소스들이 아니라 0.3 kpc 너머로 확장되는 집단에 의해 지배된다.

그러나 저자들은 자신들의 발견이 유일하다는 점에 대해서는 신중한 입장을 유지한다:

• 플럭스 및 이방성 데이터만으로는 헤일로 크기($R_1$)나 확산 계수($D_1$)를 유일하게 결정할 수 없다. 주입 형태, 배경 모델링, 방출 이력, 카탈로그 불완전성과의 퇴행성(degeneracy)이 존재한다.
• 도출된 "제밍가 규모"의 헤일로 크기는 특정 집단 모델과 데이터 세트에 기반한 조건부 제약이지, 물리적 아확산(sub-diffusion) 헤일로에 대한 확정적 측정값이 아니다.
• 단순한 일관성 테스트에서 느린 확산 영역의 결정으로 나아가기 위해서는, 개별 헤일로의 TeV 역 컴프턴 표면 밝기 프로파일, 상세한 펄서 집단 데이터(고유 운동, 탄생 위치), 그리고 전체 벡터 이방성 가능도(likelihood)와 같은 추가적인 관측치가 필요하다.

본 연구는 현재의 데이터가 공통 헤일로 파라미터 공간을 테스트할 수 있음을 보여주는 전파 제어적 입증 사례로서, 미래의 집단 추론을 위한 안정적인 커널을 제공하는 동시에 완전한 해결책을 위해 다중 메신저 입력이 필수적임을 강조한다.