거미형 이진계의 테라전자볼트 방출 예측: 충돌충격 기하학의 영향

초록

본 논문은 스파이더(거미) 이진계에서 전자·양전자 쌍이 충돌충격(intrabinary shock)에서 가속되어 발생하는 역컴프턴(TeV) 방출을 예측한다. 최신 싱크로트론 커널, 다양한 입자 주입 스펙트럼 형태, 그리고 구형·윌킨·칸토 등 세 가지 충격 기하학을 코드에 구현해 위상‑의존 스펙트럼과 에너지‑의존 광곡선을 산출하였다. 결과는 충격 형태와 관측 각도에 따라 TeV 플럭스가 크게 달라짐을 보여, 차세대 체레노프 망원경으로 검증 가능함을 시사한다.

상세 분석

본 연구는 기존 UMBRELA 코드의 한계를 보완하기 위해 네 가지 주요 개선점을 도입하였다. 첫째, 싱크로트론 방출 계산에 MacLeod(2000)의 체비쉐프 다항식 기반 커널을 적용함으로써 수정된 Bessel 함수 K5/3의 적분을 테이블 보간 대신 정확히 평가하였다. 이는 스펙트럼의 고에너지 절단부에서 부드러운 전이와 미세한 정밀도를 제공한다. 둘째, 입자 주입 스펙트럼을 단순 파워‑로우에서 세 가지 형태(단순 파워‑로우, 부드러운 절단을 포함한 브레이크 파워‑로우, 지수·초지수 절단 파워‑로우)로 확장하였다. 각 형태는 가속 효율 ηₚₐᵣₜ, 쌍 다중성 Mₚₐᵢᵣ, 최소·최대 에너지(E_min, E_max)와 연계되어 전류·광도 보존 조건을 만족하도록 정규화된다. 셋째, 광자 에너지 그리드의 로그 스케일을 확장하고, 도플러 변환 후 라보 프레임에 매핑할 때 단순 인덱싱이 아닌 선형 보간을 적용해 에너지 분포의 누락을 최소화하였다. 넷째, 충격 기하학을 구형(‘Sphere’), 구형+원뿔(‘Sphere+Cone’), 윌킨(‘Wilkin’) 및 칸토(‘Canto’) 네 가지 모델로 구현하였다. 각 모델은 입자 흐름의 아크 길이(s)와 속도 베타(β) 분포를 다르게 설정해, 입자 밀도와 도플러 부스트가 위상에 따라 어떻게 변하는지를 정량화한다. 특히 칸토 모델은 두 풍의 질량·운동량 비율 β′에 따라 충격이 퍼스트 혹은 컴패니언을 둘러싸는 형태를 전환시키며, 뉴턴‑라프슨 방법으로 θ₁을 계산한다. 이러한 기하학적 차이는 PSR J2339‑0533을 대상으로 한 시뮬레이션에서, 180° 위상에서의 최대 TeV 플럭스 위치와 광도 변동 폭을 크게 바꾸었다. 예를 들어, ‘Sphere’ 모델은 고속 베타γ_max≈7인 영역에서 강한 도플러 부스트를 보이지만, ‘Wilkin’ 모델은 충격 곡률이 커져 베타가 감소함에 따라 플럭스가 약 30% 감소한다. 또한, 관측 각도(i) 변화에 따라 광도 변조가 비선형적으로 반응함을 확인했으며, 이는 실제 관측 시 시스템의 기울기와 충격 형태를 역추정하는 데 중요한 제약조건이 된다. 전반적으로, 입자 주입 스펙트럼의 형태와 충격 기하학이 TeV 역컴프턴 방출의 위상‑에너지 특성을 결정짓는 핵심 파라미터임을 명확히 제시한다.