PSR J2032+4127의 52년 궤도와 새로운 글리치 발견

초록

Fermi‑LAT 16년 관측을 이용해 γ‑레이 바이너리 PSR J2032+4127의 궤도 주기 ≈ 52.3 yr, 이심률 ≈ 0.98, 반대축 a sin i ≈ 25.3 AU, 궤도 경사 ≈ 47°–55°를 정밀하게 측정하였다. 또한 2021년(MJD ≈ 59500)에서 작은 글리치가 추가로 발견되었다.

상세 분석

본 연구는 Fermi‑LAT이 2008년 8월부터 2025년 10월까지 수집한 16년 연속 γ‑레이 데이터를 50일 구간으로 나누어 펄스 도착 시간(Time‑of‑Arrival, TOA)을 추출하고, TEMPO2와 자체 파이썬 코드를 이용해 전 단계 위상 일관 타이밍(FPCT)과 부분 위상 일관 타이밍(PPCT)을 수행하였다. 초기에는 ν와 ˙ν의 2차원 χ² 탐색을 통해 관측된 스핀 주파수 변동이 이진 궤도에 의한 도플러 효과임을 확인하고, 케플러 방정식을 이용해 5개의 케플러 궤도 매개변수를 추정하였다. 이후 도플러 보정 후 남은 순수 스핀 변화를 이용해 TOA를 재구성하고, DD 모델로 최종 궤도와 스핀 파라미터를 동시에 피팅하였다.

핵심 결과는 다음과 같다. 궤도 주기 P_orb = 19110.5 ± 2.1 day (≈ 52.3 yr), 이심률 e = 0.979889 ± 3×10⁻⁶, a sin i = 12613.2 ± 1.0 lt‑s (≈ 25.3 AU)이며, 질량 함수 f_m을 이용해 비정상적인 경사각 i를 47°–55° 범위로 제한하였다. 이는 Be 별 MT91 213의 질량 13.1–17.5 M⊙와 NS 질량 1.4 M⊙를 가정했을 때 일관된다.

글리치 분석에서는 기존에 보고된 2011년(MJD 55810.76) 글리치 외에 2021년(MJD ≈ 59500)에서 Δν = 2.52×10⁻⁷ Hz, Δ˙ν = −1.8×10⁻¹⁶ Hz s⁻¹의 작은 글리치를 새롭게 확인하였다. 두 글리치의 상대적 크기 Δν/ν는 각각 2.74×10⁻⁷, 3.61×10⁻⁸ 수준으로, 전체 펄서 글리치 샘플에 비해 매우 미미하지만 장기 타이밍 모델링에 중요한 교정 요소가 된다.

이러한 정밀 궤도 파라미터는 γ‑레이 바이너리의 복사 메커니즘을 모델링하는 데 필수적이다. 높은 이심률과 긴 주기로 인해 퍼시스트론이 비행 중에 비정상적인 풍선형 충돌 영역을 통과하지 않으며, 따라서 전자-양성자 풍선 모델보다는 펄스 풍선-비상대론적 입자 가속 모델이 더 적합할 것으로 보인다. 또한 PSR B1259‑63와 비교했을 때, PSR J2032+4127은 비슷한 회전 구동 파워(˙E ≈ 2.7×10³⁵ erg s⁻¹)를 유지하면서도, 궤도 거리와 경사각 차이로 인해 TeV‑광도 변동이 거의 없고, GeV‑광도는 안정적인 특성을 보인다.

결론적으로, 16년 연속 γ‑레이 타이밍을 통한 전 단계 위상 일관 분석은 장기 이진 펄서 시스템의 궤도와 스핀 변동을 동시에 정밀하게 측정할 수 있음을 입증했으며, 향후 다파장 관측과 입자 가속 시뮬레이션에 중요한 기준 데이터를 제공한다.