研究チームは実際の月面地形と月殻厚の変化を、重力波に対する月の応答モデルへ組み込んだ。月の裏側に広がる厚い月殻の高地では、シミュレーション信号が概して強かった。この成果は月震計の設置場所や対象周波数の選定に役立つが、あくまで数値・理論上の予測であり、月で重力波を実際に検出したわけではない。

要点

  • 約2キロメートル格子の2次元スペクトル要素法シミュレーションと3次元固有モード結合解析を組み合わせ、実際の地形と月殻厚の横方向変化が重力波への月の応答をどう変えるかを計算した。
  • 周波数帯で積分した信号は、厚い月殻の地域で全体に約1割強い。増幅は月の裏側の高地で最も一貫し、表側の薄い月殻をもつ海では概して弱い。
  • 一部の極めて狭い周波数帯では、局所的な変位エネルギー比が10以上になる。これは現実的なモデルと球対称・成層構造の月モデルとの比較であり、重力波そのものが実測で10倍になったわけではない。
  • 増幅の主因は不均一な月殻によるモード混合である。重力波がまず四重極振動を励起し、その後エネルギーが混合した高次モードへ再配分される。
  • 設置場所は月殻厚だけでは決められない。盆地による集束、検出器の対象周波数帯、長時間の低雑音観測、まだ十分にモデル化されていない小規模な亀裂も感度を左右する。

重力波が月を通過すると、時空のごくわずかな伸縮が月全体に潮汐力を加え、極めて微弱な自由振動を起こす。Lei Zhang、Han Yan、Xian Chen、Jinhai ZhangはPhysical Review Lettersで、実際の地形と月殻厚の変化を計算に含めると、厚い月殻の地域でシミュレーション上の振動が系統的に強まると報告した。月を天然の共振アンテナとして使う構想をより具体化し、将来の月面重力波検出器の設置場所を考える手掛かりになる。

なぜ月で重力波を「聴く」のか

現在の地上干渉計が主に捉える重力波は数十〜数千ヘルツ、パルサータイミングアレイが探るのはナノヘルツ領域である。その間にある0.01〜1ヘルツの中間周波数帯には稼働中の観測所がまだなく、コンパクト連星の旋回接近、中間質量ブラックホールの合体、超大質量ブラックホールの種、初期宇宙の過程に由来する信号が含まれる可能性がある。LISA、TianQin、Taijiなどの宇宙ミッションはこの空白の一部を狙う。別の道は天体そのものの共振を利用することだ。

月の表側と裏側を並べ、下の細い発光層と厚い発光層で月殻厚の違いを示した図
左は薄い月殻の海が多い表側、右は厚い月殻の高地が多い裏側。下の発光断面は定性的な比較に限られ、南極エイトケン盆地自体の月殻は薄く、その周囲の高地は厚い。

月には海洋も大気もなく、人間活動による地面振動の大半から遠いため、背景環境は地球より静かである。重力波の四重極潮汐力は月のl = 2自由振動を励起できる。十分に高感度な月震計を月面で長期間運用すれば、理論上は対応する表面変位を測れる。しかし実際の月は滑らかで均一な球ではない。表側と裏側で地形が異なり、月殻厚も横方向に大きく変化する。こうした構造が共振信号を作り変える。

2つの手法で実際の月殻をモデル化

研究チームはLunar Reconnaissance OrbiterのLOLA高度データとGRAIL月殻モデルを用い、フンボルト海、雨の海、南極エイトケン盆地を通る大円断面に沿って、格子間隔約2キロメートルの高解像度2次元モデルを構築した。スペクトル要素法シミュレーションでは現実的な構造と球対称・成層構造の月を比較し、約0.2ヘルツまでの表面変位エネルギーを追跡した。別の3次元固有モード摂動計算は、低周波の全球モードがどう結合するかを解析した。2つの手法は扱う尺度が異なるが、比較可能な低周波範囲では月殻厚に関して同じ傾向を示した。

完全な月の表面が、4つの葉をもつシアンの四重極振動から複数の紫色の高次モードへ移り変わり、左には広い時空の波紋がある
重力波はまず主に四重極モードを励起する。月殻の横方向の不均一性が固有モードを混合し、複数の高次振動へエネルギーを再配分する。

ある周波数帯で変位エネルギーを積分すると、厚い月殻の高地におけるシミュレーション信号は球対称の基準モデルより全体に約1割強く、薄い月殻の地域では多くが弱くなった。増幅が最も一貫しているのは月の裏側の高地で、表側の海は通常弱い。約0.1ヘルツ付近の一部の極めて狭い帯域では、南極エイトケン盆地周辺の高地で局所変位エネルギー比が10に達する場合もある。これはモデル間の相対比であり、月が重力波そのものを10倍に増幅するという意味ではない。

背後にあるのはモード混合である。月が完全に対称なら、重力波は主にl = 2の四重極モードと結合する。月殻の横方向の不均一性は元のモードを新たな固有振動へ混合し、エネルギーをl > 2の高次モードに再配分する。計算でl = 2成分だけを残した場合、全球差は約4%だった。高次モードを加えて初めて厚い月殻と薄い月殻の差が明確に強まり、モード混合が主因であることを支持する。

最も厚い場所がすべての周波数で最適とは限らない

広い周波数帯の傾向を主に決めるのは月殻厚だが、特定の周波数では局所地形が答えを変える。約0.06ヘルツ付近では、盆地内のノッチや集束効果により、月殻が薄い場所でも強い信号が生じ得ることをシミュレーションが示した。月の裏側の厚い月殻をもつ高地は有望な出発点だが、実際の月震計の位置は装置の対象周波数帯に合わせ、高解像度の地形・構造モデルで地点ごとに評価する必要がある。

最も強いピークの一部は幅が約0.1ミリヘルツしかない。スペクトルでこれほど狭い構造を分離するには、1回の観測を少なくとも約1万秒、すなわち約2.8時間続ける必要がある。論文は、低雑音の月震データを1日以上連続して収集することを勧めている。将来のアレイはさらに、複数観測点間の構造化された相関を利用し、重力波応答とランダムで空間的に非コヒーレントな月震背景を分離しなければならない。

設置場所の指針から月内部の探査へ

この研究はLunar Gravitational-Wave Antenna(LGWA)やLunar Laser Interferometer for Gravitational-wave Antenna(LILA)などの構想に較正の枠組みを提供する。研究者はアレイや解析手法を設計する前に、場所と周波数ごとの月の応答を予測できる。逆に、既知の天体現象による重力波信号と月面振動を将来正確に対応付けられれば、月内部の構造がモードを変える仕方から、月の3次元内部構造に制約を与えられる可能性がある。これは研究が提案する潜在的応用であり、月内部をすでに測定したということではない。

現段階で、高解像度スペクトル要素法シミュレーションは2次元の大円断面にとどまり、3次元摂動計算も高次モードを打ち切っている。キロメートル〜サブキロメートル規模の不均質、月のメガレゴリスに広がる亀裂、個々の重力波事象の方向と偏波はまだ完全には含まれていない。完全な3次元シミュレーション、現実的な装置雑音、長期の月面データによって、モデル上の増幅が実用的な観測上の利点になるかが決まる。