Combined KM3NeT/ARCA and ANTARES search for compact neutrino sources

기술 요약: KM3NeT/ARCA와 ANTARES의 결합을 통한 컴팩트 중성미자 소스 탐색

문제 및 목적
본 연구는 점원(point-like) 천체 물리학적 소스에서 기원하는 고에너지 중성미자를 탐색하는 문제를 다룬다. 주요 목적은 두 개의 서로 다른 중성미자 검출기인 ANTARES 망원경(2008~2022년 프랑스 툴롱 인근에서 운영)과 KM3NeT/ARCA 검출기(현재 시칠리아 인근에서 건설 중)의 데이터를 결합할 수 있는 잠재력을 입증하는 것이다. 두 장비 모두 남반구 하늘과 은하 중심부를 관측하기에 최적의 가시성을 제공한다. 본 분석은 단독 분석에 비해 민감도와 발견 잠재력을 향상시키기 위해 빈 기반 가능도(binned likelihood) 접근 방식을 사용하는 통합 탐색 프레임워크를 구축하는 것을 목표로 한다.

방법론
본 분석은 네 가지 별개의 구성 요소로 이루어진 결합 데이터셋을 활용한다:

1. ANTARES 트랙 데이터셋: 4541일의 생존 시간(livetime)을 가지며, 겹치지 않는 트랙 및 샤워 데이터셋으로 분할됨.
2. ANTARES 샤워 데이터셋: 4541일의 분할 데이터에 포함됨.
3. ARCA19 트랙: 19개의 검출 라인을 가진 KM3NeT/ARCA의 트랙 유사 이벤트(48.4일).
4. ARCA21 트랙: 21개의 검출 라인을 가진 KM3NeT/ARCA의 트랙 유사 이벤트(287.4일).

정의: "트랙(Tracks)"은 뮤온 생성을 포함하는 이벤트를 의미하며, "샤워(Showers)"는 모든 중성미자 맛깔(flavor)의 중성 전류 상호작용 및 전자 중성미자의 전하 전류 상호작용에 의해 유도된다.

통계적 프레임워크
본 분석은 빈 기반 가능도 접근 방식을 채택한다. 관측된 이벤트 히스토그램($N$)은 참조 플럭스($\Phi_0 = 10^{-8} \text{ GeV cm}^{-2} \text{ s}^{-1}$)에 기반한 확장 가능한 신호 추정치($S$) 및 배경($B$)과 비교된다. 가능도 함수는 다음과 같이 정의된다:
$$\ln L = \sum_{i}^{N_{bins}} N_i \ln(B_i + \zeta S_i) - (B_i + \zeta S_i)$$
여기서 $\zeta$는 신호 강도 파라미터이다. 검정 통계량 $\lambda$는 최적 적합 신호 강도($\hat{\zeta}$)와 배경-전용 가설($\zeta=0$) 사이의 로그 가능도 비이다.

• 발견 잠재력(Discovery Potential): 50%의 의사 실험(pseudo-experiments)에서 $5\sigma$ 초과를 달야내는 데 필요한 신호 강도($\zeta_{5\sigma}$)로 정의된다.
• 민감도(Sensitivity): 네이만(Neyman) 접근법을 사용하여 계산된 $\zeta$에 대한 중앙값 90% 신뢰 수준(CL) 상한값으로 정의된다.

배경 모델링 및 계통 불확실성
배경 모델링과 관련된 계통 불확실성을 다루기 위해, 적위 띠(declination bands)에 걸쳐 네 가지 독립적인 분석 케이스가 수행되었다:

• 케이스 1: ARCA 배경은 완전히 데이터 기반이며, ANTARES 생성은 데이터 기반, 평가는 몬테카를로(MC) 방식이다.
• 케이스 2: 케이스 1과 유사하지만, ANTARES 생성 배경 데이터가 무작위화(균일한 방위각 분포)되었다.
• 케이스 3: 두 생성 및 평가 단계 모두에서 모든 데이터셋이 DEFT 알고리즘을 사용한다.
• 케이스 4: ARCA는 데이터 기반이며, ANTARES는 전적으로 MC 기반이다.

네 가지 케이스에 걸친 결과의 확산(spread)은 계통 불확실성 밴드를 정의한다. 히스토그램은 각도 거리($\log_{10} \alpha$), 재구성된 에너지($\log_{10} E$), 그리고 ANTARES 트랙의 경우 각도 오차 추정치($\log_{10} \beta$)에 대해 작성된다.

주요 결과
본 분석은 $E^{-2}$ 에너지 스펙트럼과 1:1:1 맛깔 비율을 가정한다.

• 성능 향상: ANTARES 데이터가 결합 분석 결과를 크게 주도하지만, ARCA 데이터와의 결합은 ANTARES 단독 분석에 비해 성능을 약 10% 향提升시킨다.
• 후보 소스 MG3 J225517+2409: 이 소스는 $(1.2 \pm 0.4) \cdot 10^{-4}$의 가장 낮은 사전 시험 p-값(pre-trial p-value)을 기록했으며, 이는 $(3.69 \pm 0.09)\sigma$의 유의성에 해당한다.
  • 최적 적합 플럭스: $E^2\phi = (1.39 \pm 0.02) \cdot 10^{-8} \text{ GeV cm}^{-2} \text{ s}^{-1}$.
  • 사후 시험 확률(Post-trial probability): 이항 공식을 통해 계산된 값은 $0.013 \pm 0.004$ ($(2.26 \pm 0.13)\sigma$)이다. 시뮬레이션 기반 접근법은 $0.018 \pm 0.006$ ($(2.12 \pm 0.14)\sigma$)의 사후 시험 p-값을 산출했다.
  • 1255개의 객체로부터의 초기 선택 및 이후 카탈로그 축소를 고려한 단순화된 사후 시험 추정치는 이러한 후보를 관측할 확률이 $\approx 0.1$임을 시사한다.
• 기타 후보들: 본 분석은 3C403 ($p \approx 1.5 \cdot 10^{-3}$) 및 TXS 0506+056을 포함한 알려진 다른 소스들을 평가하였으며, 106개의 카탈로그 소스에 대해 90% CL 상한값을 제공하였다.
• 스카이 맵(Sky Map): MG3 J225517+2409 근처의 스카이 맵을 보면, 하나의 ARCA21 이벤트가 소스의 위치를 포함하는 오차 박스를 가지고 있으며, 신호 가능도가 배경-전용 시나리오보다 6.3배 높음을 보여준다.

의의 및 주장
본 논문은 ANTARES와 KM3NeT/ARCA 간의 결합 분석의 타당성과 잠재력을 입증한다고 주장한다. 저자들은 현재 결과가 성숙한 ANTARES 데이터셋에 의해 지배되고 있지만, ARCA 데이터의 통합이 민감도에 측정 가능한 개선을 제공한다는 점을 확립하였다. 저자들은 KM3NeT/ARCA의 크기와 노출량이 계속 증가함에 따라 몇 년 내에 ANTARES의 성능 수준에 도달할 것으로 예상되며, 이에 따라 결합 데이터 분석이 미래의 발견에 점점 더 중요해질 것이라고 언급하였다. 본 연구는 여러 배경 모델링 전략을 통해 계통 불확실성을 처리하는 프레임워크를 제공하며, 광범위한 천체 물리학적 후보들에 대한 상한값을 설정한다.