전이형 밀리초 펄서 PSR J1023 0038의 궤도와 펄스 지연 최신 관측

초록

본 연구는 2021∼2023년 Aqueye+ 광학 데이터와 2023년 NICER X‑ray 데이터를 이용해 전이형 밀리초 펄서 PSR J1023 0038의 상승노드 통과 시각(Tasc)의 장기 변화를 측정하고, 동시 관측을 통해 광학 펄스와 X‑ray 펄스 사이의 위상 지연을 0.067 ± 0.018(112 ± 31 µs)로 새롭게 제시한다. Tasc는 연간 약 20 초씩 증가하며, 2017년 이후 포물선 형태의 증가 추세를 보인다. 이는 궤도 주기가 늘어나고 이진 분리도가 확대되고 있음을 의미하며, 비보존적 로체-리브 오버플로우와 질량 손실이 가속화되고 있음을 시사한다.

상세 분석

본 논문은 전이형 밀리초 펄서(Transitional Millisecond Pulsar, tMSP)인 PSR J1023 0038의 이진 궤도와 펄스 메커니즘을 정밀 타이밍 분석을 통해 재조명한다. 먼저, 2021년 1월, 2022년 1∼2월, 2023년 1월 및 12월에 Aqueye+ 광학 포톤 카운터를 이용해 얻은 5개의 관측 세트를 전처리하였다. 각 광자 도착 시각은 TEMPO2와 JPL DE405를 사용해 태양계 중심 좌표계(TDB)로 변환하고, Archibald et al. (2013)의 궤도 파라미터(a sin i = 0.343356 s, P_orb = 0.1980963155 d)를 적용해 이진 운동 보정을 수행하였다.

Tasc(ascending node passage time)의 변화를 추적하기 위해 각 관측 세트마다 에포크 폴딩(epoch folding) 탐색을 수행했으며, Tasc만을 자유 파라미터로 변동시켜 χ² 최대값을 찾았다. 탐색 범위는 예상값을 중심으로 ±20 s, 스텝 0.5 s였으며, 최적값은 가우시안 피팅으로 결정하고 공분산 행렬을 이용해 불확도를 산출하였다. 결과는 Table 1에 제시된 바와 같이 2021년 44.35 ± 0.17 s, 2023년 1월 81.05 ± 1.23 s, 2023년 12월 97.33 ± 1.39 s의 오프셋을 보이며, 연간 약 20 s씩 증가한다.

이러한 Tasc 변화를 시간에 대한 2차 다항식 ΔTasc(T)=A+BT+½CT² 로 모델링했을 때, C = 1.5 ± 0.1 × 10⁻⁵ s MJD⁻² (또는 C′ = 3.52 ± 0.23 × 10⁻¹¹ s s⁻¹) 로, 궤도 주기의 도함수(dP_orb/dt)와 일치한다. 이는 2017년 이후 관측된 Tasc가 포물선 형태로 증가함을 의미하며, 이전에 보고된 불규칙적인 변동과는 대조적이다.

광학 펄스와 X‑ray 펄스의 위상 차이를 측정하기 위해 2023년 1월 21일에 동시 관측된 Aqueye+와 NICER 데이터를 사용하였다. NICER 데이터는 0.5–12 keV 범위로 필터링하고, 동일한 이진 보정 파라미터와 Tasc = 80.25 ± 0.89 s (X‑ray 기준)으로 정렬하였다. 두 데이터 세트를 동일한 기준 시각(t₀ = 59965.0 MJD)과 회전 주기(P ≈ 1.687987447 ms)로 폴딩한 뒤, 2차 조화 성분의 위상을 비교하였다. 결과, 광학 펄스가 X‑ray 펄스보다 112.3 ± 30.7 µs(위상 0.067 ± 0.018) 늦게 도착함을 확인하였다. 이 위상 차이는 0–0.15 사이의 작은 값으로, 광학과 X‑ray 펄스가 동일한 물리적 메커니즘, 즉 펄서 풍과 디스크 사이의 충돌(shock)에서 발생하는 동기화된 싱크로트론 방출에 기인한다는 기존 모델을 강력히 지지한다.

마지막으로, 관측된 Tasc 증가와 궤도 분리도 확대는 비보존적 로체-리브 오버플로우(non‑conservative Roche‑lobe overflow) 시나리오와 일치한다. 질량 전달률보다 훨씬 높은 질량 손실률을 보이는 동반성은 펄서 풍에 의해 물질이 탈출하면서 시스템의 총 각운동량이 감소하고, 이에 따라 궤도 반경이 증가한다. 이러한 과정은 장기적인 궤도 주기 증가와 Tasc 포물선 증가를 동시에 설명한다.

요약하면, 본 연구는 (1) 최신 광학·X‑ray 타이밍을 통해 Tasc가 연간 ~20 s씩 증가하고 포물선 형태를 보임을 확인, (2) 광학‑X‑ray 위상 차이를 112 µs로 정밀 측정해 동일한 충격 가속 싱크로트론 방출 모델을 뒷받침, (3) 비보존적 질량 손실에 의한 궤도 팽창을 제시함으로써 전이형 밀리초 펄서의 이진 진화 메커니즘을 한층 구체화하였다.