마르카르시안 501 블라자르 제트의 필라멘트 구조에서 프로톤 가속과 고에너지 중성미자 생성

초록

본 논문은 마르카르시안 501(Mrk 501) 블라자르 제트에 전자기 필라멘트가 존재한다는 가정 하에, 다단계 확산 충격 가속 메커니즘을 적용해 프로톤이 EeV(10¹⁸ eV) 수준까지 가속될 수 있음을 보이고, 이들 고에너지 프로톤이 광자와 상호작용해 생성하는 중성미자 플럭스를 계산한다. 에너지 손실·탈출·홀이 조건 등을 종합적으로 고려한 결과, IceCube와 같은 km³ 규모 탐지기의 감도 내에서 검출 가능성이 있음을 제시한다.

상세 분석

이 연구는 블라자르 제트 내부에 강한 자기장이 얽힌 필라멘트 구조가 형성될 수 있다는 최신 MHD 시뮬레이션 결과를 출발점으로 삼는다. 필라멘트는 전도성 플라즈마가 얇은 고밀도 스트림으로 나뉘어, 각 스트림 사이에 강한 전류와 자기장이 집중되는 형태를 띤다. 이러한 구조는 전통적인 균일 자기장 모델보다 훨씬 높은 자기장 세기와 작은 스케일을 동시에 제공하므로, 입자 가속 효율을 크게 향상시킬 수 있다. 저자들은 다단계 확산 충격 가속(DSA)을 적용한다. 첫 단계는 제트 내부의 미세 충격에서 시작해, 입자가 필라멘트 내부에서 반복적으로 충돌하면서 에너지를 획득한다. 두 번째 단계는 제트 외부의 대규모 충격(예: 재단성 충격)으로, 여기서는 입자들이 필라멘트 간 경계면을 넘어가며 추가 가속을 받는다. 각 단계마다 가속 시간 τ_acc ≈ η r_g c⁻¹ (η는 난류 계수, r_g는 회전반경)으로 표현되며, η는 필라멘트의 비균일성에 따라 1–10 사이로 추정된다.

에너지 제한 요인으로는 (1) 동기복사 손실, (2) 광자-프로톤 상호작용에 의한 광전산(pγ) 손실, (3) 입자 탈출, (4) 홀이 제한(자기장 규모와 강도) 등이 있다. 동기복사 손실은 B≈0.1–1 G, γ_p≈10⁹–10¹⁰ 범위에서 τ_syn ≈ 6π m_p³c⁵/(σ_T m_e² B² γ_p) 로 계산되며, EeV 수준에서는 τ_syn이 가속 시간보다 길어 손실이 지배적이지 않다. 광전산 손실은 제트 내부의 X‑ray/γ‑ray 광자밀도 n_γ≈10⁹ cm⁻³를 가정했을 때, τ_pγ ≈ (κ σ_pγ n_γ c)⁻¹ 로 추정되며, 이 역시 다단계 가속이 충분히 빠르면 무시할 수 있다. 탈출 시간은 필라멘트 반경 R_f≈10¹⁶ cm에 대해 τ_esc≈R_f²/D (D는 확산계수) 로 정의되며, D≈(1/3) c r_g (δB/B)⁻² 로 표현된다. 최종적으로, 모든 제한을 종합하면 최대 에너지 E_max≈e B R_f (η⁻¹) ≈ 10¹⁸ eV 수준이 가능함을 보인다.

고에너지 프로톤이 주변 광자와 충돌해 pγ → Δ⁺ → n π⁺/p π⁰ 과정을 일으키면, π⁺는 μ⁺와 ν_μ를, μ⁺는 e⁺와 ν_e, (\barν_μ)를 생성한다. 저자들은 프로톤 스펙트럼을 dN_p/dE ∝ E⁻²·¹ 로 가정하고, 광자 스펙트럼을 관측된 SSC 모델에 맞춰서 설정한다. 이를 바탕으로 중성미자 플럭스 Φ_ν(E)≈(1/4π d²) (κ f_π) E_p² (dN_p/dE_p) 로 계산했으며, 여기서 f_π는 광전산 효율(≈0.1)이다. 결과는 0.1–10 PeV 범위에서 E²Φ_ν≈10⁻⁹–10⁻⁸ GeV cm⁻² s⁻¹ 정도이며, IceCube의 현재 감도(≈10⁻⁹ GeV cm⁻² s⁻¹)와 비교했을 때 검출 가능성이 있다. 또한, 시간 가변성(플라즈마 블롭의 이동에 따른 광자 밀도 변화)과 방향성(제트가 관측자와 거의 일치) 효과를 고려하면 순간적인 플럭스 상승이 기대된다.

결론적으로, 필라멘트형 제트 구조와 다단계 DSA는 기존 균일 자기장 모델보다 높은 프로톤 최대 에너지를 제공하며, 이에 따른 고에너지 중성미자 생산도 충분히 강력해 현재 및 차세대 km³ 규모 탐지기로 검출 가능성을 열어준다.