Fermi LAT를 위한 펄서 시뮬레이션 도구

초록

본 논문은 Fermi 대형 면적 망원경(LAT)의 펄서 관측 능력을 극대화하기 위해 개발된 시뮬레이션 패키지 PulsarSpectrum의 구조와 기능을 소개한다. 이 도구는 이론적·현상학적 모델에서 얻은 에너지·위상 분포를 입력으로 하여 감마선 광자를 생성하고, 천체의 위치·속도·바이너리 궤도 등을 고려한 도플러·시간 지연 효과를 적용해 위성 좌표계에서의 도착 시간을 계산한다. 시뮬레이션 결과를 통해 LAT의 감도·시간 분해능·위치 정확도를 검증하고, 실제 관측에 적용 가능한 펄서 데이터 세트를 제공한다.

상세 분석

PulsarSpectrum은 Fermi‑LAT 데이터 분석 파이프라인에 직접 연결될 수 있도록 설계된 모듈형 시뮬레이터이다. 첫 단계에서는 사용자가 선택한 이론 모델(예: 외부 전자 가속 모델, 슬롯 갭 모델, 외부 전자-양성자 혼합 모델)이나 관측 기반의 페노멘올로지 모델에서 에너지 스펙트럼 Φ(E)와 위상 프로파일 P(φ)를 정의한다. 이때 스펙트럼은 일반적인 파워‑로우 혹은 커브드 파워‑로우 형태로 파라미터화되며, 위상 프로파일은 다중 가우시안 혹은 푸리에 급수로 표현된다.

다음으로 시뮬레이터는 각 광자를 무작위 추출하여 에너지와 위상을 할당하고, 광자 방출 시점을 시뮬레이션한다. 여기서는 펄서의 회전 주기와 회전 주기 파생(Ṗ)값을 이용해 시간‑위상 변환을 수행하고, 타이밍 노이즈(랜덤 워크 형태)와 장기 스핀‑다운 변동을 포함한다. 바이너리 펄서의 경우, 궤도 원소(a, e, i, ω, Ω, T₀)를 입력받아 로렌츠 변환과 라이트 트래블 타임(LTT) 보정을 적용한다.

위성 궤도와 관측 윈도우를 고려한 최종 단계에서는 Fermi‑LAT의 실제 관측 조건을 반영한다. 위성의 위치와 속도는 공개된 SPICE 커널을 통해 실시간으로 조회되며, 이를 바탕으로 바리센트릭 보정(지구‑태양‑천체 중심 변환)과 상대론적 도플러 시프트가 계산된다. 또한, LAT의 에너지‑방향 응답 함수(IRF)를 적용해 검출 효율과 포인트 스프레드 함수(PSF)를 모사한다. 결과적으로 시뮬레이션된 광자 리스트는 실제 LAT 데이터 포맷(FT1/FT2)과 동일하게 출력되어, 기존 분석 툴(ScienceTools)로 바로 처리할 수 있다.

핵심적인 기술적 기여는 (1) 복합적인 타이밍 보정 체계, (2) 바이너리 궤도 모델링, (3) 실시간 위성 자세·위치 정보와 IRF 통합, (4) 사용자 정의 모델을 손쉽게 삽입할 수 있는 플러그인 구조에 있다. 이러한 기능들은 LAT의 감도 한계, 주기 탐지 효율, 그리고 펄서 스펙트럼 재구성 정확도를 사전에 평가하는 데 필수적이며, 실제 관측에서 발견된 새로운 펄서 후보의 검증에도 활용될 수 있다.