새로운 물리학으로 본 GRB 221009A 초고에너지 광자 현상

초록

GRB 221009A에서 LHAASO와 Carpet‑3이 각각 10 TeV와 100 TeV 이상의 광자를 관측했다. 표준 전자‑양전자 쌍생성 흡수 모델로는 이러한 고에너지 광자를 설명할 수 없으며, 저자들은 이를 두 가지 새로운 물리 시나리오, 즉 아시온‑유사 입자(ALP)와 2차 서브루미널 라오렌츠 불변성 위반(LIV) 모델에 적용해 공동 플루언스 스펙트럼을 구축하고 통계적 적합을 수행하였다. 결과는 ALP 혼합이 전체 데이터에 대해 유의미한 개선을 보이며, LIV는 제한된 파라미터 영역에서만 미미한 개선을 제공한다는 결론을 제시한다.

상세 분석

본 논문은 GRB 221009A에서 관측된 초고에너지 광자(>10 TeV, >100 TeV)가 우주배경광(EBL, CMB)과의 γγ→e⁺e⁻ 쌍생성에 의해 거의 전부 흡수될 것으로 기대되는 상황에서, 실제로 검출된 사실을 새로운 물리학적 메커니즘으로 해석하려는 시도이다. 저자들은 먼저 두 실험(LHAASO, Carpet‑3)의 플럭스 데이터를 시간‑에너지 구간별로 정리해 전체 플루언스 스펙트럼을 구축하였다. LHAASO는 0–2000 s 구간을 4개의 서브셋(두 검출기·두 시간 구간)으로 나누어 Gaussian 통계로, Carpet‑3은 4536 s 후 단일 이벤트를 Poisson 통계로 처리하였다.

표준 물리 시나리오에서는 Saldana‑Lopez 2021 EBL 모델과 CMB를 이용해 광자 흡수 optical depth τ를 계산했으며, 이 경우 250 TeV 이벤트는 거의 0 확률로 예측된다. 따라서 두 가지 대안 모델을 도입한다.

1. ALP 혼합: 광자는 GRB 호스트 은하와 은하계 자기장 내에서 ALP와 혼합될 수 있다. 저자들은 밀도 행렬 방식을 사용해 전파 방정식을 풀고, 강한 혼합이 예상되는 파라미터 영역(마스 m≈10⁻⁶–10⁻⁷ eV, 결합 g_{aγγ}≈10⁻¹¹–10⁻¹⁰ GeV⁻¹)을 탐색했다. 호스트 은하와 은하계 자기장의 불확실성을 y₀와 φ 두 매개변수에 대해 20번 무작위 샘플링하고 χ²를 평균하였다. 최적 파라미터는 m=5.16×10⁻⁷ eV, g_{aγγ}=6×10⁻¹¹ GeV⁻¹이며, 이는 표준 모델 대비 Δχ²=30.48(24 자유도)이라는 큰 개선을 보인다. ALP 혼합은 5 TeV 근처의 KM2A 데이터에서 관측된 ‘딥’도 재현하고, 10 TeV 이상에서는 강한 혼합 영역에서 광자 플럭스를 약 1/3 수준으로 억제해 Carpet‑3 포인트와도 일치한다.

2. LIV: 라오렌츠 불변성 위반을 2차 서브루미널 형태(E_{LIV,2})로 파라미터화하고, 수정된 γγ→e⁺e⁻ 단면적(정확식 Eq. 2)으로 optical depth를 재계산했다. χ² 곡선은 E_{LIV,2}≈4×10¹² GeV에서 얕은 최소값을 보이며, Δχ²=12.99(25 자유도) 정도의 미미한 개선만을 제공한다. 또한 이 파라미터는 기존 항공·우주 입자 실험이 제시하는 상한(E_{LIV,2}>10¹³–10¹⁴ GeV)과 겹쳐, 독립적인 제약과 충돌한다.

통계적 비교에서는 ALP 모델이 전반적인 플루언스 스펙트럼을 가장 잘 설명하며, 특히 5 TeV ‘딥’과 10 TeV 이상에서의 억제 효과를 동시에 재현한다. LIV는 Carpet‑3 포인트만을 부분적으로 설명하지만, 전체 χ² 개선이 제한적이며 기존 제한과도 모순된다.

저자들은 또한 EBL 모델 교체(Franceschini & Rodighiero 2017)와 자기장 강도 ±30% 변동을 테스트했으며, ALP 최적 파라미터는 68% 신뢰구역 내에 머무르고 LIV 결과는 변하지 않았다. 따라서 결론은 “ALP 혼합이 현재 관측된 초고에너지 GRB 광자를 설명하는 가장 유력한 새로운 물리 후보”라는 것이다.