An 18-year Fermi-LAT stacking limit on GeV $γ$-ray emission from particle-accelerating colliding-wind binaries

기술 요약: 입자 가속 충돌-바람 쌍성계의 GeV $\gamma$-선 방출에 대한 18년 차 Fermi-LAT 스태킹 한계치

문제 제기
질량이 큰 충돌-바람 쌍성계(CWBs)는 초음속 항성풍이 충돌하여 발생하는 충돌-바풍 영역(WCR)을 보유하고 있으며, 이곳에서는 강한 충격파가 생성되어 입자를 상대론적 에너지까지 가속할 수 있다. 비열적 라디오 싱크로트론 방출은 이러한 계에서 입자 가속이 일어나고 있음을 확인시켜 주며, $\eta$ Carinae는 GeV $\gamma$-선원으로 명확하게 검출된다. 그러나 집단 수준에서의 GeV 방출은 알려진 바가 없다. 제한된 노출량으로 수행된 이전의 전용 조사들은 개별 시스템으로부터 유의미한 방출을 발견하지 못했다 (Vela 영역의 혼잡함으로 인해 $\gamma^2$ Velorum의 검출이 $\eta$ Car보다 덜 확실하다는 점을 제외하면). 본 연구가 다루는 핵심 질문은 입자 가속 충돌-바람 쌍성계(PACWBs) 전체 집단이 집단적으로, 개별적으로는 임계값 미만인 GeV 신호를 생성하는지 여부이다.

방법론
저자들은 1 GeV 이상의 17.8년치(2008년 8월부터 2026년 5월까지) Fermi-LAT 데이터를 분석하였다. 에너지 임계값은 1 GeV 미만에서의 넓은 응답과 달리, 점 퍼짐 함수(PSF)가 크게 좁아져($68\%$ 포함 범위에서 $\lesssim 0.8^\circ$) 쌍성계를 은하 평면의 인접 천체들과 분리할 수 있다는 점을 고려하여 선택되었다.

1. 샘플 선택: 연구에는 De Becker & Raucq 카탈로그에 등재된 61개의 확정된("list A") PACWB를 활용하였다.

   • 제외 대상: 개별적으로 검출된 두 시스템($\eta$ Car 및 $\gamma^2$ Vel)은 검출되지 않은 집단의 집단적 신호를 탐색하기 위해 별도로 분류하였다.
   • 클리닝: 밝은 소스 주변에 엄격한 $1.0^\circ$ 마스크를 적용하였다(FL16Y, Signif$_{Avg} \ge 5$). 이 마스크 내에 있는 시스템들은 PSF가 관련 없는 밝은 소스(예: 항성 형성 지역, 은하 중심 소스)와 섞이기 때문에 제외되었으며, 이는 개별 소스에 대한 플럭스 측정을 신뢰할 수 없게 만들기 때문이다.
   • 병합: 상관관계가 있는 스태킹을 피하기 위해 $1.5^\circ$ 미만으로 떨어진 시스템들은 병합하였으며, 가장 높은 테스트 통계량(TS)을 가진 멤버를 유지하였다.
   • 최종 샘플: 이 과정을 통해 상호 고립된 고위도 시스템인 $N=6$개의 "깨끗한" 샘플을 얻었다.

2. 스태킹 분석:

   • 광자 지수($\Gamma \in [1.0, 3.5]$)와 적분 플럭스($\log_{10} F \in [-13, -6]$)에 대해 2차원 우도 스캔(likelihood scan)을 수행하였다.
   • 결합 테스트 통계량($TS_{stack}$)은 $N$개 소스의 개별 우도 비율의 합으로 계산되었다.
   • 제어 필드(Control Fields): 영 분포(null distribution)를 설정하기 위해 200개의 제어 필드를 생성하였다. 이들은 은하 위도($|b|$), 국부 소스 밀도(FL16Y의 $2^\circ$ 내 합산된 Signif$_{Avg}$), 그리고 은하 확산 강도($> 1$ GeV) 측면에서 타겟 샘플과 일치하도록 매칭되었다. 결정적으로, 제어 분포가 밝은 소스와의 우연한 일치를 고려할 수 있도록 동일한 $1.0^\circ$ 카탈로그 마스크를 적용하였다.

주요 결과

1. 비검출: 6개 시스템의 클린 스택은 누적 TS 3.91을 나타냈다. 이는 환경이 매칭된 제어 필드 영 분포(중앙값 TS = 10.0)와 일치하며, p-값은 0.83($-0.9\sigma$)이다. PACWB 집단으로부터의 집단적인 GeV 방출에 대한 증거는 없다.
2. 상한치: 광자 지수 $\Gamma = 2$를 가정할 때, 평균 소스당 플럭스의 95% 상한치는 다음과 같다:
   $$F(> 1 \text{ GeV}) \lesssim 1.1 \times 10^{-11} \text{ ph cm}^{-2} \text{ s}^{-1}$$
   이는 광도 한계 $L_\gamma \lesssim 1.3 \times 10^{31} (d/\text{kpc})^2 \text{ erg s}^{-1}$에 해당한다.
3. 효율 제약: 대표적인 바람 운동 에너지($L_{wind} \sim 3.3 \times 10^{36} \text{ erg s}^{-1}$)에 대해, $\gamma$-선 생성 효율은 다음과 같이 제한된다:
   $$\eta = L_\gamma / L_{wind} \lesssim 4 \times 10^{-6} (d/\text{kpc})^2$$
   샘플의 중앙 거리($\sim 1.7$ kpc)에서, 이는 $\eta$ Car의 효율($\eta \sim 10^{-3}$)보다 약 두 자릿수 낮다.
4. 아티팩트 식별: 카탈로그 마스크를 제거하고 전체 집단(59개 시스템)을 스택했을 때, $\gtrsim 3.5\sigma$의 가짜 과잉(spurious excess)이 나타났다. 저자들은 이것이 은하 평면의 밝은 카탈로그 소스와 PACWB 사이의 우연한 일치에 의한 아티팩트임을 입증하였으며, 이는 혼잡한 영역에서의 스태킹에 대한 중요한 방법론적 주의 사항을 강조한다.

입자 가속에 대한 제약
단일 존 모델(single-zone model)을 사용하여, 저자들은 플럭스 한계치를 입자 가속에 대한 제약으로 변환하였다:

• 전자: 이 한계치는 라디오 데이터로부터 추론된 자기장($B$, 수 mG에서 $\sim 1$ G)에 대해 전자 가속 효율($\eta_e$)을 $\sim 10^{-3}$로 제한한다. 표준적인 효율($\eta_e \sim 0.01–0.1$)을 달아가려면, 자기장이 일반적으로 CWBs에서 추론되는 값보다 높은 자기-광자 등가값($B \gtrsim 3–7$ G)과 비슷하거나 그보다 높아야 한다.
• 양성자: 이 한계는 양성자가 $\gamma$-선을 방출하기 위해 $pp \to \pi^0$ 붕괴를 거치기 전에 WCR에서 탈출할 가능성이 높기 때문에 양성자 가속을 훨씬 더 약하게($\eta_p \lesssim 0.11$) 제약한다.

의의 및 주장
본 논문은 $\eta$ Car가 일반적인 GeV 방출 CWB 집단의 가장 밝은 구성원이 아니라, 극단적인 바람 출력과 이심률을 가진 특이한 객체라고 결론짓는다. 비검출 결과는 대다수의 PACWB가 바람 에너지를 GeV $\gamma$-선으로 전환하는 데 있어 현저히 효율이 낮거나, LAT 대역에서의 방출을 억제하는 스펙트럼 특성(예: 컷오프)을 가지고 있음을 시사한다.

이 연구는 스태킹을 수행할 때 은하 위도뿐만 아니라 국부 소스 밀도와 확산 강도까지 매칭된 제어 필드를 사용하는 것이 필수적임을 강조한다. 타겟과 제어 필드에 동일한 마스킹을 적용하지 않을 경우, 스태킹 분석은 카탈로그의 우연한 일치에 의해 가짜 양성(false positive)을 만들어낼 수 있음을 보여준다. 결과로 얻은 제약 조건은 고위도에서 오염되지 않은 하위 집단에 대해서는 견고하지만, 마스킹으로 인해 탐사가 불가능해진 가장 유망한 후보들(밀집된 항성 형성 지역에 위치한 것들)은 여전히 조사되지 않은 상태로 남아 있다.