[論文レビュー] A Measurement Model for Precision Pulsar Timing

J. M. Cordes, R. M. Shannon|arXiv (Cornell University)|Oct 19, 2010

Advanced Frequency and Time Standards参考文献 5被引用数 49

ひとこと要約

本稿は、精度の高いパルサー時刻測定のための包括的な測定モデルを構築し、タイミング誤差の予算を定量化し、電離層散乱とパルス位相ジッターを主な確率的制限要因として特定した。高周波数観測とクロマティック補正が、重力波検出に必要なナノ秒レベルのタイミング精度を達成するために不可欠であることが示された。高精度な電波望遠鏡と補正技術を用いることで、 dispersion measure が 500 pc cm⁻³ までのパルサーも利用可能になる。

ABSTRACT

This paper describes a comprehensive measurement model for the error budget of pulse arrival times with emphasis on intrinsic pulse jitterand plasma propagation effects (particularly interstellar scattering), which are stochastic in time and have diverse dependences on radio frequency. To reduce their contribution, timing measurements can be made over a range of frequencies that depends on a variety of pulsar and instrumentation-dependent factors that we identify. A salient trend for high signal-to-noise measurements of millisecond pulsars is that time-of-arrival precision is limited either by irreducible interstellar scattering or by pulse-phase jitter caused by variable emission within pulsar magnetospheres. A cap on timing errors implies that pulsars must be confined to low dispersion measures (DMs) and observed at high frequencies. Use of wider bandwidths that increase signal-to-noise ratios will degrade timing precision if nondispersive chromatic effects are not mitigated. The allowable region in the DM-frequency plane depends on how chromatic timing perturbations are addressed. Without mitigation, observations at 1.4~GHz or 5~GHz are restricted to $\DM\lesssim 30$ and $\lesssim 100~\DMu$, respectively. With aggressive mitigation of interstellar scattering and use of large telescopes to provide adequate sensitivity at high frequencies (e.g. Arecibo, FAST, phase 1 of the SKA, and the SKA), pulsars with DMs up to 500~$\DMu$ can be used in precision timing applications. We analyze methods that fit arrival times vs. frequency at a given epoch prior to multi-epoch fitting. While the terms of greatest astrophysical interest are achromatic (e.g. orbital and gravitational wave perturbations), measurements may ultimately be limited by similarly achromatic stochasticity in a pulsar's spin rate.

研究の動機と目的

• 中性子星から望遠鏡に至るまでの全過程に及ぶ影響を考慮した、パルサー時刻測定精度の包括的誤差予算モデルの構築。
• 電離層散乱とパルサー磁気圏内のジッターに起因する確率的タイミング誤差を特定・定量化し、ミリ秒パルサーのタイミング精度を制限する要因を明らかにすること。
• 広帯域測定におけるクロマティックタイミング摂動を低減するための観測および処理戦略の特定。
• 特に重力波検出を目的とした高精度タイミングキャンペーンにおいて、利用可能なパルサーの周波数および分散測定（DM）の上限を確立すること。
• 長期にわたるタイミングキャンペーンの設計を支援するため、タイミング誤差を予測し、最適な観測パrameterを同定すること。

提案手法

• モデルは、制御不能な天体物理的要因（非制御要因）と、機器および伝搬に起因する要因（低減可能な要因）にタイミング誤差を分解し、クロマティックな電離層効果に焦点を当てる。
• 複数周波数にわたるTOAの重み付き最小二乗法フィッティングを用いて、クロマティック散乱に起因する系統的タイミング誤差の解析的表現を導出する。
• 主な方程式には、不完全な散乱モデル化に起因するTOAおよびDM推定値のバイアスを定量化する δt∞ と δDM が含まれる。これらは共分散行列と周波数依存遅延から導出される。
• 誤った散乱のべき乗則指数 X の仮定がもたらす影響を評価し、系統的誤差項 δt∞ ∝ δX × Gt∞ を導出する。ここで Gt∞ は帯域全体にわたる周波数加重積分である。
• 狭帯域と広帯域のフィッティング戦略を比較し、t∞、DM、aC、δX を同時に推定する四パrameterフィッティングを用いてモデル不一致を補正する手法を検討する。
• 実際のパrameter（例：nν = 1024、ν1 = 2 GHz、ν2 = 1 GHz）を用いて、さまざまな観測シナリオにおける定量的誤差境界を計算する。

実験結果

リサーチクエスチョン

• RQ1 ミリ秒パルサーのタイミング精度を制限する主な確率的要因は何か。また、それらは観測周波数および分散測定（DM）にどのように依存するか。
• RQ2 広帯域データを適切に補正せずにフィッティングした場合、特に散乱に起因するクロマティックな電離層伝搬効果が、TOA推定値にどのようなバイアスをもたらすか。
• RQ3 信号対雑音比とクロマティック系誤差のトレードオフを考慮した場合、タイミング誤差を最小化する最適な帯域幅および中心周波数は何か。
• RQ4 散乱のべき乗則指数 X を誤って仮定した場合、推定されるタイミング残差および DM 評価にどのような影響が生じるか。
• RQ5 電離層効果を積極的に補正する場合としない場合とで、異なる観測周波数における高精度タイミングに許容可能な最大分散測定（DM）は何か。

主な発見

• ミリ秒パルサーの高信号対雑音比測定において、タイミング精度は機器雑音ではなく、根拠のない電離層散乱またはパルサー位相ジッターによって制限される。
• クロマティック補正が行われない場合、1.4 GHz および 5 GHz の観測では、それぞれ DM ≤ 30 pc cm⁻³ および DM ≤ 100 pc cm⁻³ に制限される。
• 電離層散乱の積極的補正と大型口径電波望遠鏡（例：アレシボ、FAST、SKAフェーズ1）の使用により、DM が 500 pc cm⁻³ までのパルサーも高精度タイミングに利用可能になる。
• オクターブ帯域（例：1–2 GHz）において、クロマティック散乱に起因する系統的TOA誤差は、1 GHz における aC = 1 μs のとき約 -0.29 μs であり、これに伴い DM バイアスは約 1.15 × 10⁻³.⁶² pc cm⁻³ に相当する。
• 誤った散乱指数 X を仮定すると、系統的TOA誤差 δt∞ ≈ −0.069aCδX が生じ、X のわずかな誤差でも測定可能なバイアスが生じることが示された。
• 最適な帯域幅は、信号対雑音比の向上とクロマティック誤差の増加のトレードオフを考慮したバランスである。上側周波数が固定されている場合、帯域が広がると系統的誤差が増加するが、適切にモデル化されていなければならない。