Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

기술 요약: 미라지 소스(Mirage Sources) 및 대규모 TeV 헤일로-펄서 오프셋

문제 정의
최근 LHAASO(Large High Altitude Air Shower Observatory)와 HAWC의 관측을 통해 중년 펄사 주변에서 확장된 TeV 감마선 방출("펄사 헤일로")이 확인되었다. 이러한 소스를 해석하는 데 있어 중요한 과제는 펄사의 위치와 TeV 방출의 중심점 사이에 빈번하게 관찰되는 큰 공간적 오프셋(spatial offset)이다. 전자-양전자($e^\pm$) 쌍이 성간 매질(ISM) 내에서 등방성 확산을 한다고 가정하는 전통적인 모델은 이러한 오프셋을 자연스럽게 설명할 수 없다. 왜냐하면 등방성 확산은 물리적 규모가 여러 자기 결맞음 길이(magnetic coherence lengths)보다 클 때만 유효한 근사치이기 때문이다. 그러나 매우 높은 에너지의 전자(>100 TeV)의 경우 냉각 속도가 매우 빠르기 때문에, 전파 규모가 난류 성간 자기장의 결맞음 길이보다 작을 수 있다. 결과적으로 입자 수송은 균일하게 확산하기보다는 국부적인 자기력선을 따라 매우 비대칭적으로 일어날 수 있다. 본 논문은 이러한 비대칭적 전파가 가속기로부터 떨어진 곳에 위치한 것처럼 보이는 "미라지 소스(mirage sources)"—즉, 겉보기 방출 피크와 큰 공간적 오프셋—를 생성할 수 있는지 조사한다.

연구 방법론
저자들은 100 TeV 전자를 난류 자기장 내로 주입하여 그 궤적을 추적하기 위해 GPU 가속 기반의 제일 원리(first-principles) 테스트 입자 시뮬레이션을 사용한다.

• 자기장 생성: 난류 자기장은 푸리에 공간 방법을 통해 생성된 네 개의 독립적인 콜모고로프 스펙트럼($n_{in} = -11/3$) 자기장의 중첩으로 모델링된다. 이 다중 스케일 접근 방식(그리드 간격 $4 \times 10^{-9}$ ~ $5 \times 10^{-4}$ kpc 사용)은 현실적인 관성 범위(inertial range)와 발산 없는 자기장($\nabla \cdot B = 0$)을 보장한다. 전체 자기장은 규칙적($B_r$) 성분과 난류($B_t$) 성분을 모두 포함한다.
• 입자 전파: 전자는 로렌츠 힘 하에서의 상대론적 운동 방정식을 풀기 위해 보리스 푸셔(Boris pusher) 알고리즘을 사용하여 전파된다. 시뮬레이션은 싱크로트론 복사 및 우주 마이크로파 배경(CMB)과 성간 복사 장(ISRF)에 의한 역 컴프턴(IC) 산란을 통한 복사 에너지 손실을 고려한다.
• 통계적 샘플링: 초기 에너지 스펙트럼($dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$)을 효율적으로 샘플링하기 위해 저자들은 가중치 부여된 매크로 입자(weighted macro-particle) 접근법을 사용한다. 8,192개의 매크로 입자가 주입되며, 각 입자는 해당 에너지 빈(bin)의 물리적 전자 수에 비례하는 가중치를 가진다.
• 관측 합성: 시뮬레이션된 전자 분포는 IC 산란 행렬을 통해 합성 감마선 지도(25–40 TeV 대역)로 변환된다. 이 지도들은 LHAASO-KM2A를 모사한 가우시안 점 퍼짐 함수(PSF, $\sigma_{PSF} = 0.43^\circ$)와 컨벌루션(convolution)된 후, 포아송 노이즈 배경 모델에 더해진다.
• 소스 식별: MCMC(Markov Chain Monte Carlo)를 사용하는 베이지안 모델 비교 프레임워크가 적용된다. 구성 요소를 추가할 때마다 테스트 통계량(TS)이 증가하는지를 확인하여, 증가하는 수의 가우시안 소스를 가진 모델들을 적합(fitting)시킨다. 소스는 구성 요소를 추가했을 때 TS가 25 이상 증가할 경우($\Delta \log L > 12.5$) "식별된" 것으로 간주된다.

주요 기여 및 결과
본 연구는 펄사 거리($d$), 전자 주입 스펙트럼 지수($\alpha$), 총 자기장 강도($B$), 난류 결맞음 길이($L_c$), 그리고 규칙적 자기장과 난류 자기장의 비율($B_r/B_t$)을 변화시키며 파라미터 공간을 체계적으로 탐색한다.

1. 미라지 소스 형성: 시뮬레이션은 비대칭적 전파가 단일 가속기로부터 여러 개의 뚜렷한 감마선 소스("미라지 소스")를 생성할 수 있음을 확인했다. 이는 전자가 자기력선을 따라 이동할 때, 투영 효과(projection effects)로 인해 주입 지점으로부터 멀리 떨어진 곳에서 기원하는 것처럼 보일 때 발생한다.
2. 소스 다중성 ($N$):
   • 식별된 소스의 개수는 자기장 강도($B$)에 의해 가장 강력하게 영향을 받는다. $B$가 높을수록 냉각이 빨라져 전자가 더 작은 규모에 갇히게 되며, 이는 여러 개의 뚜렷한 소스를 형성할 가능성을 낮춘다(낮은 $N$).
   • 소스 개수는 거리($d$) 및 결맞음 길이($L_c$)와 음의 상관관계를 가지며, $B^2$와는 양의 상관관계를 가진다.
3. 소스 분리 및 플럭스 비율:
   • 분리(Separation): 식별된 소스 간의 최대 분리 거리는 $L_c$와 양의 상관관계를 가지며, $B^2$와는 음의 상관관계를 가진다.
   • 플럭스 비율 ($\chi$): 주 소스 플럭스와 미라지 소스 플럭스의 비율은 $B^2$가 증가함에 따라 감소한다. 강한 자기장 영역에서는 미라지 소스가 때때로 주 소스보다 더 밝을 수 있다.
   • 규칙적 vs 난류 자기장: 분리 거리는 $B_r/B_t$에 대해 비단조적(non-monotonic) 의존성을 보인다. 최대 분리는 $B_r/B_t \approx 1$일 때 발생하며, 이때 규칙적 자기장이 방향성을 제공하는 동시에 난류가 충분한 산란을 일으켜 오프셋을 생성한다. 순수 난류이거나 규칙적 자기장이 지배적인 경우에는 분리 거리가 작아진다.
4. 공간적 오프셋:
   • 펄사와 가장 가까운 식별된 소스 사이의 오프셋은 $B_r/B_t \approx 1$이고 $L_c$가 클 때 최대가 된다.
   • 결맞음 길이가 큰 영역(예: 은하 팔 사이 영역 또는 은하 헤일로)에서는 오프셋이 최대 50 pc에 달할 수 있다.

의의 및 주장
본 논문은 LHAASO가 관측한 소스들 중 TeV 헤일로가 연관된 펄사와 공간적으로 분리되어 보이는 현상에 대해 정량적인 설명을 제공한다고 주장한다. 구체적으로, 저자들은 1LHAASO J0216+4232u, J0212+4254u, J0206+4302u와 같은 소스들이 별개의 가속기가 아니라, 단일 밀리초 펄서(PSR J0218+4232)에 의해 구동되는 미라지 헤일로일 수 있다고 제안한다.

본 연구는 특히 고에너지 렙톤의 경우 냉각 시간이 자기 결맞음 시간보다 짧기 때문에, 성간 매질 내 가속기 근처의 입자 수송을 모델링하는 데 있어 등방성 확산 근사가 불충분함을 강조한다. 투영 효과가 난류 자기장에서 비대칭적 방출을 일으켜 오프셋과 다중 소스 형태를 자연스럽게 만들어낼 수 있음을 입증함으로써, 본 논문은 관측된 TeV 하늘을 표준 펄서 윈드 네뷸라(PWN) 모델과 조화시킬 수 있는 메커니즘을 제시한다. 저자들은 자신들의 합성 난류가 (MHD의 비등방적 알펜파 난류와 달리) 등방성을 가정하고 있으나, 비대칭적 방출을 일으키는 ISM 난류의 대규모 필라멘트 구조를 질적으로 잘 포착하고 있다고 언급하였다.