X線パルサー航法(XNAV) — 宇宙灯台を使った自律測位

この回で学ぶこと

ここまでの回で扱ってきた測距・ドップラー・DDORは、どれも1つの共通した前提の上に成り立っていました。地球(あるいは地球周辺の地上局ネットワーク)を基準点として、そこから電波を送り、跳ね返ってくる、あるいは片道で届く信号を使うという前提です。測距(レンジング)は地球-探査機間の往復伝搬遅延を、コヒーレントドップラーは地上局の水素メーザーを基準にした周波数偏移を、DDORは地球上2局の受信時刻差を使います。原子時計とUSOの回で見た通り、これらの精度の根底には「地上局が持つ超高安定な周波数基準」がありました。

しかしこの前提は、探査機が太陽系の奥深く、あるいは太陽系を離れて恒星間空間へ向かうにつれて、次第に崩れていきます。地球からの距離が数十〜数百天文単位(AU)を超えると、往復の光行時間は数日〜数週間に達し、DSNのアンテナが受信できる信号強度は距離の2乗に反比例して急激に弱くなります(リンク設計を思い出してください)。さらに将来構想されている恒星間探査機(たとえばNASAが検討してきたインターステラー・プローブ構想)や、恒星間空間そのものを長期間漂う探査機では、そもそも地球に電波を送り返すこと自体が、電力・アンテナサイズ・運用コストの制約から現実的でない局面が来ます。

このような「地球を頼れない」状況で、探査機が自分の位置を自律的に、外部からの支援なしに決定する手段として研究されてきたのが、この回で扱うX線パルサー航法(X-ray Pulsar-based Navigation, XNAV)です。地球や地上局の代わりに、宇宙にもとから存在する天体、それも極めて時計として優秀なミリ秒パルサーを「灯台」あるいは「GPS衛星」に見立てて、探査機が受信するX線パルスの到達時刻だけから自分の3次元位置を割り出す、という発想の転換がこの技術の核心です。この回では、パルサーが優れた時計たり得る物理的背景から、GNSS(全球測位衛星システム)の疑似距離方程式との数式的な類似性、そして実際に軌道上で実証されたNASAのSEXTANT/NICER実験までを見ていきます。

直感的導入: 地球を離れたら、何を「基準」にするか

まず素朴な疑問から始めましょう。GPS受信機がスマートフォンの位置を決められるのは、複数のGPS衛星がそれぞれ「自分は今この位置にいて、時刻はこれだ」という電波を送り続けているからです。受信機はその電波が届くまでの時間差から、衛星との距離(正確には後述する疑似距離)を求め、複数衛星分の距離を組み合わせて自分の3次元位置と時計のずれを同時に解きます。

この仕組みが機能するための必須条件は2つあります。(1) 送信元(衛星)の位置と時刻が精密にわかっていること、(2) 送信元が極めて正確な時計を積んでいること。GPS衛星にはルビジウムやセシウムの原子時計が搭載されており、これが位置決定精度の根幹を支えています。

では、地球からもGPS衛星からも電波が届かない深宇宙・恒星間空間にいる探査機は、何を「精密な時計を持つ、位置が既知の送信源」の代わりに使えるでしょうか。ここでXNAVが着目するのが、**ミリ秒パルサー(Millisecond Pulsar, MSP)**と呼ばれる種類の中性子星です。

パルサーは、超新星爆発の後に残る、太陽質量程度の質量が半径10km程度に圧縮された超高密度天体(中性子星)のうち、強い磁場を持ち高速で自転しているものです。その磁極方向から放射されるビーム状の電磁波(電波からX線・ガンマ線まで)が、自転に伴って灯台の光のように周期的に地球(や探査機)の方向を掃過するため、外部の観測者からはパルス状の信号として観測されます。中でもミリ秒パルサーは、自転周期が1〜数十ミリ秒という桁外れの速さで回転しており、単独の中性子星として誕生したのではなく、連星系で伴星から角運動量を「奪って」加速されたと考えられている天体です。この回転の高速性と、中性子星という巨大な慣性モーメントを持つ天体であることが組み合わさり、自転周期の変動(スピンダウン)が極めてゆっくりで滑らかという、時計として理想的な性質を生み出します。

つまりXNAVの発想は、「GPS衛星のように、既知の位置にあって、精密な時計として振る舞う送信源からの信号を複数受信し、そこから自分の位置を逆算する」という点でGPSと数学的にまったく同じ構造を持ちながら、送信源が人工衛星ではなく、太陽系の外にある既知の天体そのものであるという点で、根本的に異なる自律性を実現します。パルサーは探査機がどこにいようと、地球からの支援なしに常に「そこにある」からです。

数式による定式化

ミリ秒パルサーの周期安定度

パルサーからのパルス到達周期(トポセントリックに観測される周期)を $P$ とすると、多くのミリ秒パルサーではスピンダウン(自転の減速)率 $\dot P$ が極めて小さく、代表的な天体では

$$\dot P \sim 10^{-19} \text{–} 10^{-21} \ \text{s/s}$$

程度のオーダーに達します。これがどれほど安定かを、原子時計とUSOの回で導入したアラン偏差 $\sigma_y(\tau)$ の枠組みで捉え直してみましょう。長期的な観測(数年〜十年規模)にわたって蓄積されたタイミングデータを使うと、代表的なミリ秒パルサー(たとえば PSR B1937+21 や PSR J0437−4715)の長期周波数安定度は

$$\sigma_y(\tau) \sim 10^{-14} \text{–} 10^{-15} \quad (\tau \sim 1\text{年程度})$$

という値に達することが知られています。これは、DSN地上局の水素メーザーが $\tau \sim 1000$ 秒程度の積分時間で達成する $\sigma_y \sim 10^{-15}$ に匹敵し、長期の積分時間で比較すると、地上の最良の原子時計群(セシウム原子時計のアンサンブルなど)にも匹敵、あるいは凌駕するという評価すらされています。実際、パルサータイミング観測データを多数集めて逆に「パルサー時系(Pulsar Time Scale)」を構築し、地上の原子時系の長期的な検証に使うという研究さえ行われています。

この安定度は原子時計の量子遷移由来の安定度とはまったく異なる物理起源を持ちます。中性子星は太陽質量程度の質量を持つ巨大な剛体回転子であり、その回転運動量 $L = I\omega$ (慣性モーメント $I \sim 10^{38}\ \text{kg}\cdot\text{m}^2$)は、磁気双極子放射やパルサー風によるごくわずかなエネルギー損失でしか変化しません。この「巨大な慣性を持つ天体そのものが時計として振る舞う」という点が、電子回路や原子遷移に依存する人工の時計とは本質的に異なるスケールの安定性を生み出しています。

到達時刻の予測モデル(タイミングモデル)

XNAVを成立させるには、各パルサーについて「いつ、どの位置からパルスが届くはずか」を精密に予測できるモデルが不可欠です。これは電波天文学で確立されたパルサータイミングの技術をそのまま流用します。太陽系の重心(Solar System Barycenter, SSB)を基準にした、$n$ 番目のパルスの到来時刻 $t_n$ を、周期 $P$、周期変化率 $\dot P$ を使って多項式展開すると、

$$t_n = t_0 + n P + \frac{1}{2} n^2 P \dot P + O(\dot{\dot P})$$

という形になります。実際の観測データ解析(パルサータイミング解析)では、これを逆に使い、観測された到達時刻列から $P$, $\dot P$ などのパラメータを最小二乗的にフィットして「タイミングモデル」を確立します。ひとたび地上の電波望遠鏡群(グリーンバンク望遠鏡、パークス望遠鏡、アレシボなど)による長期観測でこのモデルが確立されれば、それは公開されたパルサーカタログ(たとえばATNFパルサーカタログ)として、任意の探査機がいつでも参照できる「既知の時計」になります。ここが人工衛星の軌道暦(エフェメリス)と同じ役割を果たす部分です。

GNSS疑似距離方程式とのアナロジー

ここからがXNAVの核心である、位置決定の数式です。まず測距・ドップラー・DDORの統合の回で扱った観測方程式の考え方を思い出してください。GNSS(GPS等)の測位原理は、複数の衛星 $i$ から送られてくる信号の送信時刻と受信時刻の差から、疑似距離(pseudorange) $\rho_i$ を構成することにあります。

$$\rho_i = |\mathbf{r}_{\text{sat},i} - \mathbf{r}_{\text{u}}| + c\,\delta t_{\text{u}} + \epsilon_i$$

ここで $\mathbf{r}_{\text{sat},i}$ は衛星 $i$ の既知位置、$\mathbf{r}_{\text{u}}$ は受信機(ユーザー)の未知位置、$\delta t_{\text{u}}$ は受信機の時計誤差(受信機は原子時計を積んでいないため、この誤差項が必ず未知数として残る)、$\epsilon_i$ は電離層・対流圏遅延などの誤差項です。「疑似(pseudo)」距離と呼ばれるのは、実際の幾何学的距離に加えて、この未知の時計誤差の寄与が混入しているためです。未知数は位置3成分と時計誤差1つの合計4つなので、最低4機の衛星からの疑似距離が必要になります。

XNAVでは、この構造をほぼそのまま踏襲しつつ、送信源を人工衛星からパルサーに置き換えます。$i$ 番目のパルサーの方向単位ベクトルを $\hat{\mathbf{n}}_i$(パルサーは地球から見て極めて遠方にあるため、探査機周辺のスケールでは事実上「無限遠の方向」として扱える平面波近似が成立します)、探査機の太陽系重心座標系での位置を $\mathbf{r}$ とすると、探査機で観測されるパルス到達時刻 $t_{\text{sc}}$ と、太陽系重心で予測されるモデル上の到達時刻 $t_{\text{SSB}}$ との差(バリセントリック補正、光行差補正と呼ばれる項)は、

$$t_{\text{SSB}} - t_{\text{sc}} = \frac{\hat{\mathbf{n}}_i \cdot \mathbf{r}}{c} + \Delta_{\text{rel}}$$

という形で書けます。ここで右辺第1項が幾何学的な光行時間差、$\Delta_{\text{rel}}$ は太陽の重力によるシャピロ遅延や、探査機の運動による相対論的補正項です。この式は、GNSSの疑似距離方程式と数式の構造がぴったり対応していることに注目してください。「衛星までの距離」が「パルサー方向への視線距離の投影 $\hat{\mathbf{n}}_i \cdot \mathbf{r}$」に、「受信機の時計誤差 $c\,\delta t_{\text{u}}$」が「探査機搭載時計(USO)のオフセット $\delta t_{\text{sc}}$」に、それぞれ対応する、という構造です。

$$\underbrace{c\,(t_{\text{SSB}} - t_{\text{sc}})}_{\text{疑似距離に相当}} = \hat{\mathbf{n}}_i \cdot \mathbf{r} + c\,\delta t_{\text{sc}} + c\,\Delta_{\text{rel}} + \epsilon_i$$

未知数はここでも探査機の位置 $\mathbf{r} \in \mathbb{R}^3$(3成分)と探査機搭載時計のオフセット $\delta t_{\text{sc}}$(1成分)の合計4つなので、GNSSとまったく同じ理屈で、少なくとも4つの異なる方向のミリ秒パルサーを同時に観測すれば、探査機は自分の位置と時計誤差を地球からの支援なしに解けることになります。実際には冗長性を持たせるため、より多くのパルサーを使い、測距・ドップラー・DDORの統合の回で見たのと同じ重み付き最小二乗法

$$\Delta \mathbf{x} = (H^T W H)^{-1} H^T W \Delta \mathbf{z}, \qquad \mathbf{x} = \begin{pmatrix} \mathbf{r} \\ \delta t_{\text{sc}} \end{pmatrix}$$

によって観測方程式を線形化し、状態ベクトル $\mathbf{x}$(位置3成分+時計誤差1成分、合計4次元)を逐次推定する、というのがXNAVの位置決定アルゴリズムの標準的な骨格です。ここで観測行列 $H$ の各行は、各パルサー方向の単位ベクトル $\hat{\mathbf{n}}_i^T$ (時計項に対応する列は定数1)から構成され、これも測距・ドップラー・DDORの統合の回でのDOP(Dilution of Precision)の議論がそのまま適用できます。すなわち、観測に使うパルサーが天球上で偏った方向にしかない場合、その方向に直交する成分の位置決定精度が悪化するという幾何学的な制約が、XNAVにも同様に効いてきます。

到達時刻の測定精度と位置決定精度の関係

パルス到達時刻をどれだけ精密に測定できるかが、最終的な位置決定精度を決めます。X線検出器で観測されるパルスは、単一光子レベルでは統計的にゆらぐため、パルスの立ち上がり位相を折り畳み(フォールディング、pulse-folding)によって統計的に推定する必要があります。観測時間 $T_{\text{obs}}$、パルサーからの検出X線光子計数率 $\dot{N}$、パルスの変調度(パルス分率)を考慮すると、到達時刻の推定誤差 $\sigma_{\text{TOA}}$ は、おおまかに

$$\sigma_{\text{TOA}} \propto \frac{W}{\text{SNR}} \propto \frac{P}{\sqrt{\dot N \, T_{\text{obs}}}}$$

という形でスケールします。ここで $W$ はパルス幅、SNRは折り畳み後のパルスプロファイルの信号対雑音比です(この形は整合フィルタやクラメール・ラオ限界の回で見た「タイミング精度は信号のシャープさとSNRで決まる」という一般則の一例になっています)。到達時刻の誤差 $\sigma_{\text{TOA}}$ は、疑似距離方程式の右辺の誤差項 $\epsilon_i/c$ にそのまま伝わるため、位置決定精度は

$$\sigma_{\mathbf{r}} \sim c\,\sigma_{\text{TOA}} \times \text{DOP}$$

という形で見積もれます。実証実験レベルでは $\sigma_{\text{TOA}}$ は数マイクロ秒〜数十マイクロ秒程度が達成されており、これは $c\,\sigma_{\text{TOA}}$ に換算すると数百m〜数kmオーダーの位置決定精度に相当します。GNSSのメートル級精度と比べると見劣りしますが、地上支援なしで恒星間空間でも原理的に機能するという点に、XNAVの独自の価値があります。

実務での使われ方

XNAVは長年理論研究の域を出ませんでしたが、2010年代後半に軌道上での技術実証が行われ、実際に機能することが示されました。

• NASA SEXTANT (Station Explorer for X-ray Timing and Navigation Technology) 実験: 2017年、国際宇宙ステーション(ISS)に搭載されたX線望遠鏡 NICER (Neutron star Interior Composition Explorer) の観測データを使い、複数のミリ秒パルサー(PSR B1821−24、PSR B1937+21 など)からのX線パルスをリアルタイムに観測・処理して、ISSの軌道上の位置を自律的に決定する実証実験を行いました。地上のGPS/軌道暦データと突き合わせた結果、数時間の観測で数km(公表値では概ね10 km以内)の精度で位置決定が可能であることが実証され、XNAVが単なる理論上の可能性ではなく実装可能な技術であることを示す重要なマイルストーンになりました。
• 中国の「慧眼」(HXMT: Hard X-ray Modulation Telescope、Insight衛星) でも同種のパルサー航法実験が行われ、独立に技術実証が報告されています。
• 将来ミッションでの応用構想: 木星・土星より外側を航行する外惑星探査ミッションや、太陽系を離れる恒星間探査機(インターステラー・プローブのような構想)では、往復光行時間が長すぎて地球からのリアルタイムな軌道決定支援が事実上不可能になります。こうしたミッションで、探査機が自律的に自己位置を認識し、姿勢制御や軌道修正の判断材料にする手段として、XNAVは有力な候補技術の1つに位置づけられています。深宇宙探査機が太陽合で地球との通信が困難になる期間の自律運用の補完としても構想されています。
• 現状の技術的課題: (1) X線検出器の感度と有効面積の制約により、十分なSNRのパルスプロファイルを得るには数千秒〜数万秒オーダーの長時間観測が必要になることが多く、リアルタイム性に乏しい。(2) ミリ秒パルサーのタイミングモデル自体が、長期的なパルサーの「グリッチ」(内部構造の急激な変化に由来する周期の跳躍)や、未知の伴星の影響などによって時々刻々わずかに更新され続けており、そのモデル更新を探査機側がどう受け取るか(通信できない状況では搭載カタログの陳腐化が避けられない)という運用上の課題がある。(3) 小型探査機に搭載可能なX線検出器の小型・軽量化と、そのSNR向上を両立させる技術開発が、実用化の鍵として現在も進められている。

演習問題

1. GNSSの疑似距離方程式 $\rho_i = |\mathbf{r}_{\text{sat},i} - \mathbf{r}_{\text{u}}| + c\,\delta t_{\text{u}} + \epsilon_i$ と、本文で導いたXNAVの観測方程式 $c(t_{\text{SSB}}-t_{\text{sc}}) = \hat{\mathbf{n}}_i\cdot\mathbf{r} + c\,\delta t_{\text{sc}} + c\,\Delta_{\text{rel}} + \epsilon_i$ を比較し、両者の未知数・既知量の対応関係を表にまとめてください。また、なぜXNAVでは衛星位置ベクトルの差ではなく、方向単位ベクトルとの内積という形になるのか、幾何学的な理由(パルサーまでの距離のスケール)を説明してください。
2. あるミリ秒パルサーの長期アラン偏差が $\sigma_y(\tau=1\text{年}) \sim 3\times10^{-15}$ であったとする。原子時計とUSOの回で扱った水素メーザーの $\sigma_y(\tau=1000\text{s}) \sim 10^{-15}$ と比較したとき、両者を単純に同じ物差しで比べることの妥当性・限界について、積分時間 $\tau$ の違いに注意しながら論じてください。
3. パルス到達時刻の測定誤差が $\sigma_{\text{TOA}} = 10\ \mu\text{s}$ のとき、これに起因する視線方向の距離換算誤差 $c\,\sigma_{\text{TOA}}$ を計算してください。また、DOPが$3$程度のとき、最終的な3次元位置決定誤差の目安を見積もってください。
4. XNAVは最低4つの未知数(位置3成分+時計オフセット1成分)を解くために最低4つのパルサー方向が必要だと説明しました。もし探査機がすでに極めて安定なUSOを搭載しており、事前に地上の原子時計群と同期を取ってから出発したため時計オフセットが無視できるほど小さいとみなせる場合、必要なパルサーの最小数はいくつになるか、そしてそれでもなお複数のパルサーを観測することが望ましい理由を、DOPの観点も含めて説明してください。

まとめと次回予告

この回では、これまでの回で暗黙のうちに前提としてきた「地球(または地上局ネットワーク)を基準にした航法」という枠組みそのものを問い直しました。ミリ秒パルサーが、中性子星という巨大な剛体回転子としての物理的性質ゆえに、原子時計に匹敵する周期安定度を持つこと、そしてその安定な「宇宙の灯台」を複数観測することで、GNSSの疑似距離方程式とほぼ同じ数学的構造の観測方程式を使って、探査機が地球の支援なしに自分の位置を自律的に決定できることを見ました。NASAのSEXTANT/NICER実験によって、この技術はすでに理論から軌道上実証の段階に進んでいますが、X線検出器の感度やタイミングモデルの精度といった課題は今も研究が続いています。

XNAVは「観測データから自分の位置を求める」という問題を扱いましたが、探査機の自律性という観点では、この位置決定の結果を使って軌道そのものをどう推定・伝播させ続けるかという、もう一段上のレイヤーの問題が残っています。次回は、地上のカルマンフィルタによる軌道決定とは異なり、探査機が搭載コンピュータ上でリアルタイムに自分の軌道を推定・更新し続けるオンボード自律軌道決定の考え方に、軽く触れる形で見ていく予定です。

参考文献

• K. C. Gendreau, Z. Arzoumanian, T. Okajima, “The Neutron star Interior Composition Explorer (NICER): design and development,” Proc. SPIE 8443 (2012)
• S. R. Winternitz et al., “X-ray Pulsar Navigation Algorithms and Testbed for SEXTANT,” IEEE Aerospace Conference (2015)
• J. E. Mitchell, S. R. Winternitz et al., “SEXTANT X-ray Pulsar Navigation Demonstration: Initial On-Orbit Results,” AAS/AIAA Space Flight Mechanics Meeting (2018)
• S. I. Sheikh, “The Use of Variable Celestial X-ray Sources for Spacecraft Navigation,” Ph.D. Dissertation, University of Maryland (2005)
• D. N. Matsakis, J. H. Taylor, T. Marshall Eubanks, “A statistic for describing pulsar and clock stabilities,” Astronomy & Astrophysics, 326, 924–928 (1997)
• G. Hobbs et al., “The ATNF Pulsar Catalogue,” Astronomical Journal, 129, 1993–2006 (2005)
• J. H. Yuen (ed.), Deep Space Telecommunications Systems Engineering, JPL Publication 82-76