Joint Multi-Period Fermi-LAT and LHAASO Constraints on Axion-Like Particles from Mrk 421 Using Profile Likelihood with Gaussian Copula Correlation

기술 요약: Mrk 421에 대한 Fermi-LAT 및 LHA also의 결합 다중 주기 Fermi-LAT 및 LHAASO의 액시온 유사 입자(ALP) 제약

문제 정의
액시온 유사 입자(ALP)는 표준 모델의 잘 동기 부여된 확장으로서, 외부 자기장에서 광자-ALP 진동을 유도하여 먼 거리의 광원으로부터 오는 고에너지 감마선 스펙트럼을 변화시킬 수 있다. 이전 연구들이 ALP 파라미터를 제약해 왔으나, 여전히 제약되지 않은 상당한 영역이 파라미터 공간에 남아 있다. 이러한 탐색에서의 주요 과제는 계통 오차, 특히 감마선을 감쇠시키는 외계 배경 복사(EBL)를 처리하는 것이다. 전통적인 접근 방식은 EBL 흡수와 광자-ALP 진동을 독립적인 곱셈 인자로 취급하여, 이들의 비자명한 상호작용을 간과하곤 한다. 또한, 다중 에포크(multi-epoch) 분석은 관측 기간들이 통계적으로 독립적이라고 가정하는 경우가 많으며, 공유된 소스 가변성, 기기 계통 오차 및 모델 가정에서 기인하는 약한 잔류 상관관계를 무시한다. 이러한 단순화는 편향된 제약이나 인위적으로 좁혀진 신뢰 구간을 초래할 수 있다. 추가로, 블레이자(blazar) Mrk 421은 GeV(Fermi-LAT)와 TeV(LHAASO) 방출 사이의 약한 상관관계를 보이며, 이는 서로 다른 물리적 기원을 시사하여 광대역 데이터 해석을 복잡하게 만든다.

방법론
저자들은 2021년 3월부터 2024년 3월까지의 Mrk 421에 대한 다섯 개의 동시 관측 에포크 데이터를 사용하여, Fermi-LAT과 LHAASO를 결합하여 ALP 파라미터($m_a$ 및 $g_{a\gamma}$)를 제약하기 위한 포괄적인 프레임워크를 제시한다.

1. 전파 형식론(Propagation Formalism): 인수 분해 근사법과 달리, 본 연구는 EBL 흡수가 광자-ALP 혼합 해밀토니안 내에 흡수 항($i/2\lambda_\gamma$)으로 직접 포함되는 자기 일관적(self-consistent) 광자 전파 프레임워크를 채택한다. EBL 강도는 불확실성을 고려하기 위해 가우시안 사전 확률($\mu=1.0, \sigma=0.2$)을 가진 섭잡수 파라미터(nuisance parameter)인 정규화 계수($K_{EBL}$)로 취급된다.
2. 소스 모델링: Mrk 421의 약한 GeV-TeV 상관관계를 해결하기 위해 두 가지 별개의 제트 방출 모델을 테스트한다:
   • 컴팩트한 내부 영역(플레어링 VHE 방출)과 확장된 외부 영역(정적 방출)을 특징으로 하는 이중 존 하이브리드 모델(two-zone hybrid model).
   • 약한 상관관계와 뚜렷한 가변성 진폭을 고려할 때 물리적으로 타당한 단일 존 하드론 모델(single-zone hadronic model).
     두 모델 모두 헬리컬(helical) 및 엉킨(tangled) 자기장 성분을 활용하며, 광자를 소스, 은하간 매질(IGMF), 은하 자기장(GMF)을 통해 전파시킨다.
3. 통계적 프레임워크: 저자들은 결합 프로파일 가능도(joint profile-likelihood) 분석을 개발한다. 에포크들을 독립적으로 취급하는 한계를 해결하기 위해, 보수적인 상관 계수($\rho = 0.03$)를 가진 가우시안 코풀라(Gaussian copula) 형식을 도입한다. 이는 변환된 $\Delta\chi^2$ 과잉 분포를 통해 에포크 간의 약한 상관관계를 모델링한다. 분석은 ALP가 없는 가설 하에서 섭잡수 파라미터(고유 스펙트럼 파라미터 및 $K_{EBL}$)를 전역적으로 최적화한 후 ALP 질량과 결합 상수를 프로파일링한다.
4. 데이터 처리: 광자 생존 확률은 gammaALPs 코드를 사용하여 4차원 그리드(에너지, ALP 파라미터, EBL 정규화) 상에서 사전 계산된 후, 가능도 평가 중에 보간된다. 고유 스펙트럼 모델은 유한한 샘플 크기에 대해 수정된 아카이케 정보 기준(AICc)을 사용하여 에포크별로 선택된다.

주요 기여

• 통합 전파 모델: 본 논문은 광자-ALP 혼합과 EBL 흡수를 전파 방정식 내에서 자기 일관적으로 처리하는 모델을 구현하여, 곱셈적 인수 분해에 따른 편향을 피한다.
• 상관 다중 에포크 분석: 관측 기간 사이의 약한 상관관계를 처리하기 위해 가우시안 코풀라를 도입함으로써, 다중 에포크 데이터셋을 결합하는 통계적으로 일관된 방법을 제공하여 인위적으로 좁혀진 신뢰 구간을 방지한다.
• 광대역 제약: 본 연구는 Fermi-LAT(GeV)와 LHAASO(multi-TeV to PeV)의 독특한 광대역 커버리지를 활용하여 넓은 질량 범위(0.1–500 neV)에 걸쳐 ALP 파라미터를 조사한다.
• 강건한 모델링: 두 가지 모델(이중 존 및 단일 존 하드론)을 테스트함으로써, 본 연구는 복잡한 방출 물리(emission physics)를 고려하여 제약이 단일 방출 가정에 과도하게 의존하지 않도록 보장한다.

결과
결합 프로파일 가능도 프레임워크를 사용하여, 저자들은 광자-ALP 결합 상수($g_{a\gamma}$)에 대한 95% 신뢰 수준(CL) 상한선을 도출한다:

• 저질량 영역 ($m_a \lesssim 1$ neV): 제약은 매우 엄격하며, 하한값은 $5.0 \times 10^{-14}$ GeV$^{-1}$, 상한값은 $8.33 \times 10^{-14}$ GeV$^{-1}$이다. 이 영역은 CAST 제약에 민감하다.
• 중간 및 고질량 영역 ($0.1 - 500$ neV):
  • 이중 존 하이브리드 모델의 경우, 보수적인 상한값은 $g_{a\gamma} < 6.50 \times 10^{-12}$ GeV$^{-1}$이다.
  • 단일 존 하드론 모델의 경우, 보수적인 상한값은 $g_{a\gamma} < 7.30 \times 10^{-12}$ GeV$^{-1}$이다.
• 거동: 제외 윤곽선(exclusion contours)은 특징적인 "점프-플래토(jump-plateau)" 형태를 보인다. $m_a \approx 1$ neV 부근에서 상한값이 약 10배 증가하는 급격한 전이가 발생하며, 이후 높은 질량에서는 점진적으로 완화된다. 하한값은 전체 질량 범위에 걸쳐 안정적으로 유지된다.
• 강건성: 대안적인 상관 계수($\rho = 0, 0.15, 0.3$)를 사용한 테스트 결과, 제약의 변화가 10배 미만 내에 머물러 결과의 안정성을 확인하였다.

의의
본 논문은 이러한 결과가 LHAASO의 광대역 초고에너지 커버리지와 Mrk 421에 대한 장기 모니터링으로부터 큰 혜택을 받았다고 주장한다. 개발된 프레임워크는 차세대 감마선 관측소들을 위한 미래의 다중 메신저 ALP 탐색에 적합한 강건하고 확장 가능한 파이프라인을 제공한다. 약한 상관관계를 명시적으로 모델링하고 EBL 흡수를 자기 일관적으로 처리함으로써, 본 연구는 기존의 인수 분해 접근 방식에 존재하는 계통적 편향을 해결하고, ALP 파라미터 공간에 대해 더 신뢰할 수 있고 보수적인 제약을 제공한다. 본 연구는 Mrk 421의 약한 GeV-TeV 상관관계가 정확한 ALP 제약을 위해 복잡한 방출 모델(이중 존 또는 하드론)을 필요로 함을 입증하며, 제안된 통계적 방법론이 관련 불확실성을 효과적으로 처리함을 보여준다.