Noisy Timing Behavior is a Feature of Central Compact Object Pulsars

문제 정의
중심 콤팩트 천체(Central Compact Objects, CCOs)는 초신성 잔해(supernova remnants, SNRs)에서 발견되는 독특한 부류의 냉각 중성자별이다. 이들은 안정적인 열적 X선 방출, 펄서 풍 성운(pulsar wind nebulae)의 부재, 그리고 느린 스핀다운율($|\dot{f}| < 10^{-15}$ s$^{-2}$)으로부터 추론된 매우 약한 표면 쌍극자 자기장($B_s < 10^{11}$ G)을 특징으로 한다. 이러한 약한 자기장은 관측된 X선 펄세이션을 생성하는 데 필요한 비균일한 표면 온도(핫스팟)를 설명하기에는 불충분하여 이론적 과제를 제기한다. 또한, 이 객체들의 타이밍 불규칙성의 기원은 불분명하다. 일반적인 회전 동력 펄서(rotation-powered pulsars) 집단에서 글리치 활동과 타이밍 노이즈는 스핀다운율($\dot{f}$)과 강하게 상관관계가 있지만, CCO는 $\dot{f}$ 값이 매우 낮아 이러한 현상을 보일 통계적 가능성이 낮다. 본 논문은 알려진 세 개의 CCO 펄서(1E 1207.4−5209, PSR J0821−4300, PSR J1852+0040)가 타이밍 불규칙성을 보이는지, 그리고 만약 그렇다면 이를 이산적인 글리치(discrete glitches) 또는 연속적인 타이밍 노이즈(continuous timing noise)로 모델링할 수 있는지 다룬다.

연구 방법론
저자들은 Chandra, XMM-Newton, NICER 관측소를 통해 20년 동안 수집된 X선 관측 데이터를 사용하여 세 개의 알려진 CCO 펄서에 대한 타이밍 연구를 수행하였다.

• 데이터 처리: 광자 도착 시간은 DE405 에페메리스를 사용하여 바리센트 역학 시간(TDB)으로 변환되었다. 펄스 주기는 피리어도그램 탐색을 통해 측정되었으며, 특정 펄스 프로파일 정의(1E 1207.4−5209의 경우 최소값, PSR J0821−4300의 경우 최대값, PSR J1852+0040의 경우 골(valley) 중심)를 기반으로 위상-제로 도착 시간(ToAs)이 생성되었다.
• 타이밍 솔루션: 연구팀은 TEMPO 소프트웨어를 사용하여 위상 연결형(phase-connected) 타이밍 솔루션을 시작하였다. 이들은 단순 이차 적합(주파수 $f$ 및 미분 $\dot{f}$)부터 삼차 적합( $\ddot{f}$ 포함), 그리고 글리치 모델(주파수 $f$와 미분 $\dot{f}$의 이산적 단계)에 이르는 다양한 모델을 테스트하였다.
• 비교 분석: 타이밍 잔차(residuals)를 분석하여 이산적인 글리치 이벤트와 연속적인 타이밍 노이즈를 구분하였다. 저자들은 타이밍 노이즈 지표($\sigma_{TN,2}$)를 계산하고, 이를 Shannon & Cordes (2010)가 정의한 유사한 스핀 파라미터를 가진 정형적(canonical) 펄서의 기대값과 비교하였다.

주요 결과

1. 1E 1207.4−5209: 저자들은 2015년 에포크 이후로 타이밍 솔루션을 확장하였다. 데이터는 글리치 전 에페메리스로부터 유의미하게 벗어난다. 글리치 모델(주파수와 미분의 이산적 도약)은 $\Delta f/f \approx 1.6 \times 10^{-9}$의 글리치 크기와 함께 데이터를 잘 설명한다. 그러나 글리치가 없는 대안 모델(이차 또는 삼차 적합) 또한 데이터를 잘 설명하지만, $\chi^2$ 값은 더 높다. 주목할 점은, 글리치 후의 주파수 미분($\dot{f}$)이 글리치 전 값보다 약 9% 덜 음수(less negative)라는 것인데, 이는 일반적인 펄서 글리치와 비교했을 때 부호와 크기 면에서 이례적인 변화이다.
2. PSR J0821−4300: 이 펄서는 유사하게 노이즈가 많은 타이밍 거동을 보인다. 저자들은 2015년경에 타이밍 불규칙성을 확인하였다. 글리치 모델은 $\Delta f/f \approx 7 \times 10^{-10}$의 주파수 도약을 시사하는 반면, 더 이른 시기(약 2011년경)에 발생한 "안티글리치(antiglitch)" 모델(주파수 감소)이 데이터를 거의 비슷하게 잘 설명한다. 1E 1207.4−5209와 마찬가지로, 이벤트 발생 후 주파수 미분은 덜 음수가 되었다(약 7%).
3. PSR J1852+0040: 12년간의 모니터링 공백으로 인해, 저자들은 전체 기간에 걸쳐 위상 연결형 솔루션을 유지할 수 없었다. 2021–2023년의 데이터는 불확실성 범위 내에서 2004–2009년의 에페메리스를 외삽한 결과와 일치한다. 결과적으로, 본 연구는 다른 두 CCO에서 나타난 정도의 글리치나 타이밍 노이즈를 탐지할 만큼의 민감도를 갖지 못했다.
4. 타이밍 노이즈 크기: 이차 적합을 사용하여 타이밍 노이즈로 분석했을 때, 1E 1207.4−5209와 PSR J0821−4300의 잔차는 유사한 $\dot{f}$를 가진 정형적 펄서의 평균값을 각각 약 157배와 57배 초과하였다. 이는 노이즈 지표의 자연 로그 값에서 각각 5 표준 편차와 4 표준 편차의 편차를 나타낸다.

의의 및 주장
본 논문은 노이즈가 많은 타이밍 거동이 CCO 펄서의 특징임을 확립하며, 특히 PSR J0821−4300이 1E 1207.4−5209와 양적으로 유사한 불규칙성을 보이는 두 번째 알려진 CCO임을 확인하였다.

저자들은 희소한 샘플링과 라디오 펄서에 비해 낮은 ToA 정밀도로 인해 이산적인 글리치와 연속적인 타이밍 노이즈를 구분하는 것이 어렵지만, 불규칙성의 크기가 유사한 $\dot{f}$를 가진 일반 펄서 집단에 비해 극단적이라고 주장한다.

결정적으로, 본 논문은 관측된 $\dot{f}$의 변화(덜 음수가 됨)와 극단적인 타이밍 노이즈가 높은 내부 온도, 강력한 매몰된 자기장, 초유체 역학 등 다른 젊고 에너지가 넘치는 펄서들과 공유되는 내부적 특성에서 기인할 가능성이 높다고 상정한다. 저자들은 CCO의 안정적인 열적 X선 방출과 펄서 풍의 부재를 고려할 때, 이러한 내부 메커니즘이 자기권 효과(magnetospheric effects)보다 더 유력한 동력원이라고 명시적으로 밝히면서도, 1E 1207.4−5209에서 약한 라디오 펄싱이 최근 발견된 점은 자기권의 기여를 재고할 필요가 있음을 인정하였다. 대안적으로, 잔류 낙하 원반(fallback disk)으로부터의 낮은 수준의 강착(accretion) 또한 타이밍 노이즈의 원인으로서 배제되지 않았다.

본 연구는 CCO가 타이밍 불규칙성을 이해하기 위한 독특한 실험실임을 결론지으며, 이들의 "노이즈가 많은" 거동이 내부 초유체와 자기장 사이의 상호작용과 관련하여 일반적인 펄서 집단의 글리치 및 타이밍 노이즈 이론에 정보를 제공할 수 있음을 시사한다.