パルサーはどうやって宇宙の時を刻むのか

宇宙の彼方には、驚くべき規則正しさで時を刻み続ける天体が存在します。それが「パルサー」です。パルサーとは、太陽より重い恒星が生涯の最後に超新星爆発を起こした後に残る、極めて超高密度な「中性子星」です。強大な磁極から放たれる指向性の高い電波ビームが、高速自転に伴って地球の視線方向を周期的に横切る様子から「宇宙の灯台」と呼ばれます。

天文学者はパルスが地球に届く「パルス到着時刻（TOA）」を極限の精度で計測し、予測モデルとのわずかなズレ（タイミング残差）を解析することで、太陽系外惑星の発見、一般相対性理論の検証、さらには宇宙を満たすナノヘルツ重力波の検出まで成し遂げてきました。パルサーはなぜこれほど正確に時を刻むことができるのか。その物理的仕組みから観測補正の技術、深宇宙航法への応用までを解き明かします。

---

1. 宇宙の灯台が生まれる仕組み——中性子星の超高速自転と放射ビーム

パルサーが時を刻む時計の「振り子」は、天体自身の自転です。その主役である中性子星は、太陽の8倍以上の質量を持つ恒星が核融合を終え、超新星爆発を起こした後に中心核が激しく圧縮されて形成されます。

中性子星の直径はおよそ20キロメートル（半径約10キロメートル）と都市の市街地サイズしかありませんが、太陽の約1.4倍から2倍もの質量が凝縮しています。原子核に匹敵するほどの高密度の物質が詰まった状態にあり、NASAの解説資料によれば、ティースプーン1杯分の質量だけでも約1000万トン規模——山岳一つ分に匹敵します。

この天体が猛烈な速度で自転する理由は「角運動量保存の法則」です。フィギュアスケート選手が腕を胸元に引き込むと回転が急加速するように、もともと数十万キロメートル以上あった恒星がわずか直径20キロメートルへと一気に収縮するため、毎秒数回から数百回という超高速で自転する天体となります。収縮に伴う磁場の増幅は強い磁場を説明する一つの要素ですが、その起源や詳細は研究されています。

中性子星において、自転の中心となる「自転軸」と、磁力線が吹き出す「磁軸」は一般に傾いて交差しています。超強大な磁場が超高速で回転すると猛烈な誘導電場が発生し、表面から引き抜かれた荷電粒子が光速近くまで加速されます。加速粒子は磁力線に沿って運動しながら、磁極を中心とする円錐状（コーン状）の領域に向けて指向性の極めて高い電波やX線のビームを放射します。

中性子星が自転すると、放射ビームは宇宙空間を円を描くように掃過します。このビームが地球の視線方向を周期的に横切る幾何学的配置にある場合、地球の観測者には自転周期と完全に一致したパルス信号として捉えられます。これが「灯台モデル」です。裏を返せば、ビームが地球の視線方向を横切らない配置の中性子星は、どれほど激しく回転していても地球からはパルサーとして観測できません。

また、パルサーの自転は永久不変ではありません。真空の宇宙空間に浮かぶ巨大な慣性モーメントの塊として極めて安定して回転し続けるものの、電磁波（磁気双極子放射）の放出や荷電粒子の宇宙空間への放出によって自転の運動エネルギーを外部へ散逸させています。さらに天体内部の構造変化による微小な影響なども相まって、自転速度はごくわずかずつ遅くなっていきます。この現象を「スピンダウン」と呼びます。

---

2. 若いパルサーとミリ秒パルサー——桁違いの安定性を生むリサイクル進化

パルサーの性質は一様ではありません。自転周期と減速のペースにより、「通常の若いパルサー」と「ミリ秒パルサー」に大別されます。

超新星爆発で誕生して間もない若いパルサーは、自転周期が約0.1〜数秒程度で強大な磁場を持ちます。電磁波放射によるエネルギー散逸が大きいため減速が比較的早く、スピンダウン率（自転周期の時間変化率 \dot{P}）は典型的に1秒間に約 10^{-15} 秒（秒/秒）程度です。

一方、精密計時で主役となるのが「ミリ秒パルサー（MSP）」です。自転周期はわずか1〜10ミリ秒（毎秒100〜700回以上回転）と極めて短く、しかもスピンダウン率 \dot{P} はおよそ 10^{-18}\text{ s/s} 以下、極めて安定した天体では 10^{-20} \sim 10^{-21}\text{ s/s} に達します。若いパルサーと比べて1万倍から100万倍も自転速度が落ちにくい抜群の安定性を誇ります。

この安定性を生み出すのが、連星系における「リサイクル過程」です。中性子星が通常の恒星とペアを組む連星系では、自転が遅くなって一度活動を停止した中性子星に対し、寿命を迎えて膨張した伴星からガスが流れ込みます（質量降着）。ガスが形成する降着円盤から物質が落下する際、軌道の角運動量が中性子星へと継続的に注入され、何千万年もかけて再び猛烈な勢いで再加速（スピンアップ）されます。この過程で磁場は若いパルサーの1万分の1程度に弱まるため、電磁波放射によるエネルギー散逸が大幅に抑えられ、超高速自転を長く維持できるリサイクル中性子星が完成します。

天文学では、自転周期 P とその変化率 \dot{P} から、一定の減速則や初期周期に関する仮定を置いた「特性年齢（\tau）」という代表的な時間尺度を見積もります。

\tau = \frac{P}{2\dot{P}}

特性年齢は、パルサーの減速の緩やかさを示す目安です。実際に自転周期が2倍になるまでの時間や、必ずしも実年齢を直接示すものではありません。「数億年に1秒しか狂わない」という時計の誤差とも別です。

さらに、ミリ秒パルサーの計時を支える要因が「グリッチ（突発的な自転加速）」の稀少さです。中性子星内部の超流動体と固体地殻との回転差の歪みが限界に達すると、星震とともに内部の角運動量が地殻へ伝達され自転速度が跳ね上がります。かにパルサーなどの若いパルサーではグリッチが数年〜十数年ごとに頻発しますが、ミリ秒パルサーは過去の物質降着や長い進化で歪みが緩和されており、発生頻度は若いパルサーに比べてはるかに稀です（J0613−0200で微小なグリッチが観測された実例はあるため皆無ではありませんが、観測上は極めて例外的です）。

ミリ秒パルサーと地上の原子時計との関係についても正確な理解が必要です。「パルサーは原子時計を超える究極の基準器であり、原子時計を校正するための時計である」という見方は科学的に適切ではありません。

原子時計は、到着時刻を測る地上の基準です。パルサーの自転は人工時計とは異なる物理過程なので、長期の変化を照合する手がかりになります。ただし、電波の到着時刻を測定・校正するには地上の時刻系と観測装置を使います。「計測全体が原子時計から完全に独立する」という意味ではありません。

---

3. 到着時刻（TOA）の計測と補正技術——太陽系重心（SSB）への換算

パルサーを時計として機能させるには、電波信号から正確な到着時刻を導き出し、地球環境や太陽系に起因する様々なズレを極限まで補正する必要があります。

受信した個々のパルスは弱いため、自転周期に合わせて多数回の信号を重ね、平均的なパルス形を求めます。その基準位相から観測所での到着時刻（Site TOA）を推定します。この時刻の記録には地上の時計を使い、その後に太陽系重心への補正や伝搬遅延の補正を行います。

しかし、観測局のある地球は自転しながら太陽の周りを秒速約30キロメートルで公転しているため、望遠鏡とパルサーの距離は刻々と変化します。そこで地上の局所到着時刻を、太陽系全体の質量中心である「太陽系重心（SSB: Solar System Barycenter）」を基準とする慣性座標系時刻（太陽系重心力学時: TDBなど）へと数学的に換算します。なお太陽系重心自体も銀河中心に対して秒速200キロメートル以上で公転しており、SSBは絶対静止点ではなく、太陽系内の運動を統一して扱うための慣性基準系です。

地上時刻からSSB時刻への換算には、次の3大補正が不可欠です。

1. 幾何学的遅延（ローマー遅延）

地球が太陽を公転することで生じるパルサーとの距離変動です。地球軌道の半径（約1億5000万キロメートル）を行き来するため、電波の到達時間差は最大で約500秒（約8.3分）に達します。これをナノ秒精度で補正するには、NASAジェット推進研究所（JPL）などが開発する「太陽系天体暦（惑星暦）」に基づき地球の位置を精密計算します。ただし惑星暦における木星など巨大惑星の質量や軌道のわずかな不確実性は、後述の重力波探索において系統ノイズとなり得るため、現在も精密な検証が続けられています。

1. 一般相対論的遅延（アインシュタイン遅延とシャピロ遅延）

一般相対性理論が予言する時空の歪みの補正です。地球の公転軌道が楕円であるため、太陽重力ポテンシャルの変化と速度変化によって時計の進み方が周期的に変わる「アインシュタイン遅延」が生じます。また電波が大質量天体の重力場近傍を通過する際に伝播時間が遅れる「シャピロ遅延」も補正します。太陽近傍をかすめる場合の遅延は最大約120マイクロ秒（120 \mu\text{s}）に達する一方、木星近傍通過時の遅延は最大約200ナノ秒（200 ns）程度であり、天体質量に応じた桁違いのスケール差を精密にモデル化します。

1. 星間プラズマによる周波数分散遅延

星間空間に広がる希薄な電離ガス（フリー電子）中を電波が伝播する際、低周波ほど伝播速度が遅れ、高周波より遅れて地球に届きます。この遅延は周波数の2乗に反比例し、視線方向の電子密度を積分した「分散尺度（DM）」に比例します。この影響を差し引くため、複数の周波数帯（例えば1.4 GHzと2 GHzなど）で同時に観測する「多周波観測」を行い、周波数ごとの到着時間差からDMを逆算して補正します。 ただし、星間プラズマは乱流として常に揺れ動いており、地球の公転やパルサーの固有運動によって電波が通過する領域が変わるため、DMの値自体が時間変動します（DM variation）。この変動ノイズを完全に消去することは原理的に困難であり、ミリ秒パルサーの到達時刻精度を制約する主要な要因となっています。

これらの補正とパルサー自身の自転減速モデルを適用した上で、実際の観測到着時刻と予測時刻の差を算出したものが「タイミング残差」です。ミリ秒パルサーで達成されるタイミング残差の標準偏差（RMS）は典型的にサブマイクロ秒級（約0.1〜1マイクロ秒）であり、最良の天体では数十ナノ秒（数十 ns）という極限の精度に達します。

---

4. 到着時刻のずれから何を調べるか――単一天体と多数天体

前掲の模式図は、パルサーのビームが見える仕組み、望遠鏡での受信、観測時刻に補正を加える工程の入口を示します。この節では、その後に得られる『予測時刻との差』をどう読むかを文章で整理します。図にタイミング残差の観測グラフや重力波の検出結果は載っていません。

回転する放射ビームが視線を横切る灯台モデルと、受信時刻を予測と比較する解析は別の段階です。まず地上の時計で受信時刻を記録し、必要な補正を施してからモデルとの差を調べます。

到着時刻の解析では、観測所と太陽系重心の位置関係、相対論的な遅れ、星間プラズマの影響などを補正します。図は処理の概略を示すもので、実際の積算波形、各補正量の値、解析結果は描いていません。

補正後の到着時刻と予測時刻との差をタイミング残差といいます。その変化を一つの天体について時間に沿って調べる方法と、多数の天体の残差を比較する方法では、問いが異なります。

一つのパルサーを長く追うと、タイミング残差が時間とともにどう変わるかを調べられます。周期的な変化は伴星などを考える手掛かりになりますが、残差だけで一つの原因を特定できません。ここでいう時間変化は文章で説明しており、前掲の図の波形は実際の残差データではありません。

多数のパルサーを使うタイミングアレイでは、別々の天体の残差を比較し、天球上の角距離に応じた相関を調べます。一つの天体の残差を時間に沿って見る解析とは、横軸にあたる量も目的も異なります。重力波背景の検証で使う理論上の関係は『ヘリングス・ダウンズ関係』と呼ばれますが、前掲の模式図はその理論曲線や観測結果を描いたものではありません。

一つの残差パターンには、天体固有のノイズ、時刻系、太陽系暦、星間プラズマ、観測装置など複数の要因が入り得ます。重力波背景を調べるには、多数のパルサーの残差に角距離依存の相関があるかを統計的に検討し、他の共通ノイズも評価します。

---

5. 時空のゆらぎを捉える——パルサータイミングアレイ（PTA）と重力波

高精度なタイミング残差の解析は、天文学と基礎物理学に画期的な成果をもたらしてきました。

その最初の金字塔が、1992年にアレクサンデル・ヴォルシュチャンらがミリ秒パルサー「PSR B1257+12」で達成した史上初の太陽系外惑星の発見です。パルサーのタイミング残差に現れた約66日と約98日周期のサインカーブ変動から、地球質量級（地球の約4.3倍、約3.9倍など）の複数惑星の存在が突き止められました。主系列星ではなく死んだ星の周りで人類初の系外惑星が見つかったことは、世界に大きな衝撃を与えました。

さらにパルサーは一般相対性理論の究極の検証場となりました。1974年にラッセル・ハルスとジョゼフ・テイラーが発見した連星パルサー「PSR B1913+16」は、中性子星同士が約7.75時間周期で公転する連星系です。30年以上に及ぶ観測の結果、重力波放出に伴うエネルギー損失によって軌道周期が年間76マイクロ秒ずつ短縮していることが判明し、一般相対論の予言と「約0.2%の精度で一致」しました。この重力波の間接証明により、両氏は1993年のノーベル物理学賞を受賞しました。

そして21世紀、この技術は時空のさざ波を直接捉える「パルサータイミングアレイ（PTA）」へと進化しました。

銀河系内に点在する数十個以上の高精度なミリ秒パルサーを10〜20年以上にわたり監視し、地球とパルサーを結ぶ空間を「銀河スケールの重力波望遠鏡」として利用します。地上レーザー干渉計（LIGO等）がPTAより高周波の重力波（ブラックホール合体など）を捉えるのに対し、PTAが標的とするのは波長1光年未満〜数十光年に相当する、周波数「ナノヘルツ帯（10^{-9} \sim 10^{-7}\text{ Hz}、周期約数か月〜数十年）」という途方もなく長波長の重力波です。重力波が銀河を横切ると時空がゆっくり伸縮し、全パルサーの到着時刻にナノ秒単位の相関したズレが生じます。

2023年6月、複数のパルサータイミングアレイ共同研究が、それぞれの観測データや共同解析から、ナノヘルツ帯の重力波背景を支持する結果を報告しました。これは全グループが一つの統合データセットで同一の解析をしたという意味ではありません。

このナノヘルツ重力波の発生源としては、銀河合体に伴い形成された「超巨大ブラックホール連星系」の群れが最も有力視されています。ただし現行データでは超巨大ブラックホール連星と完全に確定したわけではなく、初期宇宙の相転移や宇宙ひもなど宇宙論的起源の可能性も理論的に検証されています。国際PTA（IPTA）による全天データの統合解析を通じて、発生源の正体を確定する研究が今も進められています。

---

6. 宇宙のGPSへ——X線パルサー航法（XNAV）の現在地と展望

パルサーが刻む高精度な時は、深宇宙探査機の自律的な位置決定を可能にする実践的な技術「X線パルサー航法（XNAV）」へと応用されています。

従来の深宇宙探査機は、地上の大型アンテナ網（ディープスペースネットワーク: DSN）による電波追跡に全面的に依存してきました。しかし遠方へ進むほど通信遅延が数十分〜数時間と増大し、探査機の増加に伴って地上アンテナの運用枠が逼迫するという構造的課題を抱えています。

パルサーを航法の手掛かりにする方法がX線パルサー航法（XNAV）です。宇宙機搭載の基準時計とX線検出器を用いてパルスの到着時刻を測り、事前の地上観測などに基づく到着時刻モデルと比べます。搭載装置の大きさなどの制約から、X線で観測できるパルサーを使う方式が検討されてきました。

ただし、地上GPSのように瞬時に測位できるわけではありません。ミリ秒パルサーから届くX線光子（フォトン）は極めて微弱であり、小型検出器には毎秒数個〜数十個程度しか飛び込んできません。十分な統計精度でTOAを導出するには、一定時間（数時間〜数日程度）にわたる光子積算が必要です。そのため実運用では、軌道力学モデルと「拡張カルマンフィルター」等の状態推定アルゴリズムを組み合わせ、時間をかけて推測位置の誤差を絞り込んでいきます。

この方式を宇宙で実証したのがNASAのSEXTANT実験です。ISSのNICER望遠鏡で複数のX線パルサーを観測し、約2日間のデータを用いて機上で位置を推定しました。GPSとの照合では、多くの推定位置が約5キロメートル以内に収まりました。これは地上の事前観測や到着時刻モデルが不要という意味ではありません。

現在、月周回拠点「ゲートウェイ」を含むシスルナ空間での運用や、火星以遠を目指す無人探査機の自律航法に向け、機器のさらなる小型軽量化と高精度化が進められています。

---

7. よくある疑問と誤解（FAQ）

Q1: パルサーは地上の人工原子時計を置き換える「究極の時計」なのですか？

A: 置き換えるものではありません。地上の原子時計は到着時刻を測る基準です。ミリ秒パルサーの自転は人工時計とは異なる物理過程なので、長い期間の変化を照合する手掛かりになります。ただし、パルサーの時刻測定自体は地上の時計と観測装置に依存し、計測全体が独立しているわけではありません。

Q2: パルサーの電波を受信すれば、その場ですぐ正確に時刻を合わせられますか？

A: すぐには合わせられません。微弱信号を数千〜数万回積算（エポックフォールディング）して波形を取り出し、さらに地球の公転・自転や重力、星間プラズマによる遅延を太陽系重心（SSB）基準で精密補正して初めて高精度な時刻が得られます。

Q3: パルサーの自転速度は未来永劫まったく変わらないのですか？

A: わずかに減速（スピンダウン）しています。電磁波放射などで回転エネルギーを失うためです。P/(2\dot P) で表す特性年齢は、仮定に依存する減速の時間尺度であり、文字どおり自転周期が2倍になる時間ではありません。

---

結び

パルサーが時を刻む仕組みは、天体物理学、一般相対性理論、そして高度な観測補正技術の融合によって成り立っています。

中性子星の規則的な回転は、地上の時計とは異なる物理過程から生じる長期の比較対象です。一方、パルス到着時刻の測定には地上の時刻基準と補正が欠かせません。残差を使った解析も、単一の曲線だけで惑星や重力波を特定するものではありません。

パルサーの規則的な信号は、恒星や連星の研究に加え、ナノヘルツ重力波背景を探る手掛かりにもなっています。X線パルサーを使う自律航法への応用は宇宙で実証され、実用化に向けた研究が進んでいます。

---