거울 같은 고에너지 광원: 테라 전자 미라지와 펄서 오프셋의 비대칭 전파

초록

본 연구는 펄서 주변의 난류 자기장 속에서 100 TeV 전자가 비대칭적으로 전파되는 과정을 GPU 기반 시뮬레이션으로 재현하고, 이로 인해 관측되는 ‘미라지’ 소스가 펄서 위치와 크게 어긋나는 현상을 규명한다. 전자 주입 스펙트럼, 거리, 자기장 세기·코히어런스 길이, 규칙·난류 비율 등 파라미터 변화를 체계적으로 탐색하여 LHAASO가 보고한 이상적인 오프셋을 설명한다.

상세 분석

이 논문은 100 TeV 전자(주로 20–30 TeV γ선으로 방출)의 비대칭 전파를 정밀히 모델링하기 위해, 3차원 난류 자기장을 Fourier 방법으로 생성하고, 다중 스케일(δ₁~δ₄) 구성을 통해 10⁵ pc에 이르는 관측 스케일을 포괄한다. 전자 궤적은 보리스 푸셔(Boris pusher)를 이용해 시간 단계당 gyration 주기의 1/100으로 적분하며, synchrotron 및 inverse Compton 손실을 실시간으로 적용한다. 초기 전자 스펙트럼은 dN/dE ∝ E⁻ᵅ exp(−E/E_cut) 형태로 설정하고, 로그 에너지 구간 8192개에 각각 하나의 매크로 입자를 배치해 가중치를 부여함으로써 통계 효율을 극대화한다.

시뮬레이션 결과는 전자들이 한 코히어런스 길이(L_c) 내에서 평균 자기장 방향을 따라 나선형으로 이동하고, L_c를 초과하면 방향이 급격히 바뀌는 특성을 보인다. 관측선과 일치하는 방향으로 자기장이 정렬될 경우, 전자 밀도가 시각적으로 집중되어 ‘미라지’ 소스가 형성되며, 이는 실제 펄서 위치와 수십 파섹에서 수백 파섹까지 오프셋될 수 있다. 파라미터 스터디에서는 (1) 전자 주입 지수 α가 작을수록(즉, 하드 스펙트럼) 고에너지 전자가 더 오래 살아 남아 멀리 전파하므로 미라지 비율이 증가하고, (2) 자기장 세기 B_t가 클수록 라디에이션 손실이 강화되어 전자 전파 거리가 감소하지만, 동시에 코히어런스 길이가 동일하면 미라지 형성 확률이 약간 상승한다는 복합 효과를 발견했다. (3) 코히어런스 길이 L_c가 5 pc에서 40 pc로 증가할수록 미라지 소스 수 ξ가 비선형적으로 증가하며, 특히 L_c ≥ 20 pc에서 ξ가 2배 이상 급증한다. (4) 규칙 자기장 대비 난류 비율(B_r/B_t)이 높을수록 전자 흐름이 보다 무작위화되어 미라지 소스가 다수 나타나지만, 동시에 전체 γ선 밝기가 감소한다는 트레이드오프가 존재한다.

관측 측면에서는 LHAASO‑KM2A의 PSF(σ≈0.43°)와 25–40 TeV 에너지 밴드를 적용해 합성 γ선 맵을 생성하고, 베이지안 MCMC와 TS(테스트 통계량) 기반 모델 선택을 통해 최소 2개의 가우시안 소스를 식별할 경우 미라지 존재를 확정한다. 시뮬레이션에 기반한 미라지 비율(Mirage Ratio)은 거리(d)와 B_t, L_c, α에 따라 0%~30% 범위로 변동하며, 특히 d≈2 kpc, B_t≈3 μG, L_c≈20 pc, α≈2.0인 경우 15% 이상의 미라지 소스가 기대된다. 이러한 결과는 LHAASO가 보고한 펄서와 γ선 중심 사이의 수십 파섹 오프셋을 자연스럽게 설명한다.

결론적으로, 비대칭 전파와 시각적 투영 효과가 고에너지 펄서 할로의 형태를 크게 왜곡시킬 수 있음을 입증했으며, 이는 기존의 등방성 확산 모델이 과소평가한 현상이다. 향후 관측에서는 미라지 소스의 다중성 및 스펙트럼 차이를 통해 난류 코히어런스 길이와 자기장 구조를 역추정할 수 있는 새로운 진단 도구가 될 전망이다.