[論文レビュー] Measuring timing properties of PSR B0540-69

Minjun Kim, Hongjun An|arXiv (Cornell University)|Mar 4, 2019

Pulsars and Gravitational Waves Research参考文献 11被引用数 5

ひとこと要約

本研究は、PSR B0540-69のSwift X線観測データ1,126日を再分析し、そのタイミング解と減速指数を精緻化した。以前報告された異常に低い減速指数 $ n \approx 0.031 $ が、タイミングノイズによってバイアスされていたことが判明し、修正された解析により $ n = 0.163 \pm 0.001 $ が得られた。これは顕著ではあるが、極端な値ではなく、以前の乖離は計測不能なタイミングノイズに起因していた。この結果は、ノイズ環境下におけるパルサータイミング測定の信頼性を向上させた。

ABSTRACT

We report on the timing properties of the `Crab twin' pulsar PSR~B0540$-$69 measured with X-ray data taken with the Swift telescope over a period of 1100\,days. The braking index of the pulsar was estimated to be $n = 0.03\pm0.013$ in a previous study performed in 2015 with 500-day Swift data. This small value of $n$ is unusual for pulsars, and a comparison to old measurements of $n\approx2.1$ for the same target determined $\sim$10 years ago suggests a dramatic change in the braking index. To confirm the small value and large change of $n$, we used 1100-day Swift observations including the data used for measuring $n=0.03$. In this study we found that the braking index of PSR~B0540$-$69 is $n=0.163\pm0.001$, somewhat larger than $0.03$. Since the measured value of $n$ is still much smaller than $2.1$, we can confirm the dramatic change in the braking index for this pulsar.

研究の動機と目的

• 以前の研究で用いられた500日ベースラインを上回る、延長されたSwift X線観測データを用いて、PSR B0540-69のタイミング解と減速指数を再評価すること。
• Marshallら（2016年）が報告した異常に低い減速指数 $ n \approx 0.031 $ が、モデル化されていないタイミングノイズによってバイアスを受けていたかどうかを調査すること。
• 残差タイミング残差におけるタイミングノイズの影響をモデル化・補正することで、パルサータイミング解の精度を向上させること。
• 再評価された減速指数が、中性子星の放射メカニズムおよび内部ダイナミクスに与える影響を検討すること。特に、若年齢かつパルサー風ネビュラと関連するこのパルサーの文脈において。

提案手法

• 著者らは、MJD 57070から58068までの期間にわたる、1,126日分のSwift X線望遠鏡（XRT）データを用いた。これは以前の500日研究の2倍の期間をカバーする。
• 標準的なX線イベントスクリーニングおよび光子抽出を、HEASoftソフトウェアパッケージを用いて実施し、良好な時間区間およびグレードのための標準的なフィルタリング基準を適用した。
• 脈動位相を2次微分までの多項式モデルにフィットさせることでタイミング解を導出し、位相は $ \phi(t) = \phi_0 + \nu t + \frac{1}{2}\dot{\nu} t^2 + \frac{1}{6}\ddot{\nu} t^3 $ としてモデル化した。
• 非定常的挙動を検出するためにタイミング残差を分析し、タイミングノイズの影響を分離するために、立方多項式 $ \phi(t) = at^3 + bt^2 + ct + d $ をフィットした。
• タイミングノイズによる $ \dot{\nu} $ および $ \ddot{\nu} $ のバイアスは、フィットされた立方多項式係数と真の微分値を比較することで定量化され、$ \dot{\nu}_{\text{bias}} = 2b $ および $ \ddot{\nu}_{\text{bias}} = 6a $ が得られた。
• 修正されたタイミング解を用いて、減速指数 $ n = \ddot{\nu} \nu / \dot{\nu}^2 $ を計算し、タイミング解の誤差から不確実性を誤差伝搬法で評価した。

実験結果

リサーチクエスチョン

• RQ1 PSR B0540-69の以前報告された減速指数 $ n \approx 0.031 $ は統計的に信頼できるものであったか、それともモデル化されていないタイミングノイズによってバイアスを受けていたか？
• RQ2 500日から1,126日へと観測ベースラインを延長することで、タイミング解および減速指数測定の精度と正確性にどのような影響があったか？
• RQ3 若年齢で明るいパルサー、特にPSR B0540-69において、タイミングノイズの影響（立方多項式位相トレンドとしてモデル化）が $ \dot{\nu} $ および $ \ddot{\nu} $ の推定にどの程度歪みをもたらすか？
• RQ4 再評価された減速指数 $ n \approx 0.163 $ は、磁場の進化や降着を含む理論的モデルと整合するか、あるいは矛盾するか？
• RQ5 混雑した Large Magellanic Cloud 地域におけるパルサー風ネビュラや近傍源からの空間的混入が、フラックス測定に影響を与え、結果としてタイミング解にバイアスをもたらす可能性はあるか？

主な発見

• 1,126日分のSwiftデータを用いた修正されたタイミング解により、減速指数は $ n = 0.163 \pm 0.001 $ となり、以前報告された $ n = 0.031 \pm 0.013 $ より顕著に高い値が得られた。
• 2つの結果の乖離は、定量的にタイミングノイズによって説明できる。最初の500日間にわたる約10 msの立方位相トレンドが、$ \ddot{\nu} $ に $ -2.3 \times 10^{-21} \, \text{s}^{-3} $ のバイアスをもたらしていた。
• タイミングノイズは $ \phi(t) = at^3 + bt^2 + ct + d $ としてモデル化され、係数は $ a = -3.82 \times 10^{-22} \, \text{s}^{-3} $、$ b = 1.73 \times 10^{-15} \, \text{s}^{-2} $、$ c = 2.22 \times 10^{-8} \, \text{s}^{-1} $ であった。このモデルは、初期の解で観測されたバイアスを説明した。
• 修正された $ \dot{\nu} $ および $ \ddot{\nu} $ の値は、それぞれ $ -2.528664(1) \times 10^{-10} \, \text{s}^{-2} $ および $ 0.531(4) \times 10^{-21} \, \text{s}^{-3} $ であり、以前の推定値よりも精度が向上した。
• 本研究は、以前の結果 $ n \approx 0.031 $ が、計測不能なタイミングノイズによってバイアスを受けていたことを確認し、真の減速指数はより高い値であるが、依然として異常低値であると結論づけた。
• 著者らは、将来的な高分解能観測（例：Chandra）が、パルサーのフラックスとネビュラおよび背景の混入を分離するのに役立つ可能性があり、長期的なフラックス安定性およびタイミング測定に影響を及える可能性があると提言した。