Search for neutrino emission from blazar $γ$-ray flares accounting for possible neutrino time delays

기술 요약: 가능한 중성미자 시간 지연을 고려한 블레이저 $\gamma$-선 플레어에서의 중성미자 방출 탐색

문제 정의
IceCube 중성미자 관측소가 고에너지 중성미자($E_\nu \gtrsim 1$ TeV)의 존재를 입증했음에도 불구하고, 확산 플럭스(diffuse flux)에 기여하는 구체적인 소스는 특정 이벤트가 블레이저(예: TXS 0506+056)나 활동 은하(예: NGC 1068)와 연관되어 있다는 증거가 있음에도 여전히 명확히 식별되지 않은 상태로 남아 있다. 블레이저와 IceCube 중성미자 사이의 상관관계를 찾는 기존의 탐색들은 일반적으로 중성미자 도착과 주요 전자기적(특히 $\gamma$-선) 플레어 사이의 시간적 일치를 가정했다. 그러나 이론적 고려 사항에 따르면 이러한 일치가 성립하지 않을 수 있다. 만약 중성미자가 $p\gamma$ 상호작용을 통해 생성된다면, 제트 프레임에서의 양성자 가속 시간($t'_{acc}$)과 양성자 에너지 손실 시간($t'_{p\gamma}$)은 전자기 플레어를 일으키는 전자들의 시간보다 수십 배에서 수백 배 더 길 수 있다. 결과적으로, 고에너지 중성미수($\gtrsim 100$ TeV)는 관련 $\gamma$-선 플레어의 정점으로부터 수개월에서 수년 후에 도착할 수 있다. 본 논문은 이러한 잠재적인 시간 지연을 고려한 탐색이 부족하다는 점을 다룬다.

방법론
저자들은 IceCat-1 공개 카탈로그의 고에너지 중성미자 경보와 두 가지 서로 다른 블레이저 소스 샘플 간의 교차 매칭 분석을 수행하였으며, 광범위한 제트 프레임 시간 지연($t'_{delay}$) 범위에 대해 상관관계를 스캔하였다.

1. 데이터 세트:

   • 중성미자: IceTop 신호와 중첩되지 않는 IceCat-1 카탈로그의 348개 경보 ($E_\nu \gtrsim 100$ TeV).
   • 소스 목록:
     • H21+: 알려진 도플러 인자($D$), 적색편이($z$), 그리고 Fermi-LAT 광도 곡선을 가진 MOJAVE 데이터베이스(Homan et al. 2021)의 AGN 294개 샘플.
     • R24+: 스펙트럼 에너지 분포 모델링을 통해 도출된 도플러 인자를 가지며, Fermi-LAT 광도 곡선과 교차 매칭된 CGRaBS 카탈로그(Rodrigues et al. 2024b)의 블레이저 295개 샘플.
   • 광도 곡선: Fermi-LAT 광도 곡선 저장소(LCR)의 월 단위 빈(bin) 처리된 광자 플럭스 광도 곡선.

2. 통계적 프레임워크:

   • 공간적 연관성: 중성미자가 확장된 각도 불확실성 영역(계통적 불확실성 $\Delta = 0.78^\circ$ 포함) 내에 위치할 경우 해당 블레이저와 연관된 것으로 간주하였다.
   • 가중치 체계: 각 경보에 대해 신호 강도, 큰 각도 불확실성에 대한 페널티 인자, 그리고 시간적 인자 $G_{\nu a}$를 기반으로 가중치 $w_\nu(t'_{delay})$를 계산하였다.
   • 시간적 인자: 시간적 인자는 중성미자 도착 시간($t_\nu$)을 도착 전 블레이저의 $\gamma$-선 광도 곡선의 전역 최댓값($t^{\gamma \max}_a$)과 연결한다. 이 인자는 시간 차이가 예상되는 관측자 프레임 지연($t_{delay} = t'_{delay}(1+z)/D$)과 빈 크기 내에서 일치할 경우 최대 가중치를 부여하며, 편차에 대해서는 가우시안 페널티를 적용한다.
   • 검정 통계량 (TS): 전역 TS는 모든 경보에 대한 개별 가중치의 합으로 정의된다.
   • 유의성 검정: 실제 데이터의 TS를 중성미자의 적경과 도착 시간을 무작위화한 49,999개의 모의 데이터 세트에서 생성된 TS 분포와 비교하였다. 사전 시험(pre-trial) $p$-값은 계산되었으며, 이후 $t'_{delay}$를 스캔함에 따라 시행 교정(trials correction)을 거쳤다.

주요 기여

• 새로운 탐색 전략: 본 연구는 느린 양성자 가속 또는 에너지 손실로 인해 발생하는 보편적인 제트 프레임 시간 지연($t'_{delay}$)을 구체적으로 고려하여, IceCube 중성미자와 블레이저 $\gamma$-선 플레어 사이의 연관성을 탐색한 최초의 시도이다.
• 이중 소스 샘플: 제트 운동학의 불확실성에 대한 상관관계의 견고성을 테스트하기 위해, 도플러 인자를 결정하는 방법이 서로 다른(밝기 온도 논거 vs. SED 모델링) 두 개의 독립적인 블레이저 카탈로그를 활용하였다.
• 체계적 지연 스캔: 분석은 $\gamma$-선 방출보다 중성미자 생성이 뒤처질 수 있는 이론적 범위를 포괄하도록 1개월에서 약 $3.5 \times 10^4$일까지 $t'_{delay}$를 스캔한다.

결과

• 통계적 유의성: 본 탐색은 IceCat-1 중성미자와 H21+ 또는 R24+ 소스 목록 중 어느 쪽에서도 상관관계에 대한 통계적으로 유의미한 증거를 발견하지 못했다.
  • H21+ 목록의 경우, 최소 사전 시험 $p$-값은 0.037 ($t'_{delay} \approx 7.7 \times 10^3$ d에서)이었다.
  • R24+ 목록의 경우, 최소 사전 시험 $p$-값은 0.026 ($t'_{delay} \approx 1.9 \times 10^3$ d에서)이었다.
  • 시행 교정 후, 이 최솟값들은 각각 0.14 (1.06$\sigma$) 및 0.11 (1.21$\sigma$)에 해당하는 사후 시험 $p$-값이 되었으며, 이는 귀무 가설과 일치한다.
• 최상위 후보: 가장 유의미한 연관성은 3C 454.3의 2010년 플레어와 연결된 IC120523B 및 3C 279의 플레어와 연결된 IC150926A에 의해 주도되었다. 이 특정 이벤트들은 예상된 지연과 시간적 정렬을 보여주었으나, 전체 데이터셋에 걸친 전역적 상관관계를 확립하기에는 불충분했다.
• 도플러 인자의 불일치: 최적의 지연 시간은 H21 및 R24 방법론 간의 3C 454.3 및 3C 279에 대한 도플러 인자 추정치 차이에 의해 크게 영향을 받아, 두 소스 목록 간에 4배의 차이를 보였다.

의의 및 주장
저자들은 예상되는 시간 지연을 고려하더라도, 자신들의 결과가 연구된 특정 블레이저 집단과 유의미한 상관관계를 뒷받침하지 않는다는 결론을 조심스럽게 내렸다. 저자들은 이러한 영 결과(null results)에 대해 다음과 같은 몇 가지 비상호 배타적인 이유를 제시한다:

1. 비효율성: 블레이저는 $\sim 100$ TeV 중성미자의 비효적인 소스일 수 있으며, IceCube 플럭스의 극히 일부($\sim 0.5\%$)만을 기여할 수 있다.
2. 에너지 규모: 만약 블레이저가 효율적인 중성미자 소스라면, 그 방출의 정점은 현재 IceCube 경보 샘플의 주요 민감도 범위를 넘어서는 훨씬 높은 에너지($\gtrsim 10$ PeV)에 있을 수 있다.
3. 모델 가정: 보편적인 제트 프레임 시간 지연($t'_{delay}$) 가정이 가속 메커니즘 및 자기장의 소스별 차이로 인해 유효하지 않을 수 있다.
4. 디커플링(Decoupling): 양성자 가속과 에너지 소산이 전자 과정과 공간적 또는 시간적으로 분리되어 있어, 중성미자 생성이 관측된 $\gamma$-선 플레어와 인과적으로 연결되지 않을 수 있다.

본 논문은 향후 감도 개선을 위해 더 큰 데이터 세트(예: 차기 IceCat-2 및 CAZ 카탈로그)와 $\gamma$-선 플레어보다 지연되는 것으로 알려져 중성미자 생성 시간 척도와 더 잘 상관될 수 있는 라디오 광도 곡량의 포함이 필요할 것이라고 상정한다.