The effect of matter discreteness on gravitational wave propagation in post-geometrical optics

技術要約：ポスト幾何光学における物質の離散性が重力波伝播に及ぼす影響

問題提起
重力波（GW）の伝播は、通常、波が測地線に従う幾何学的光学近似を用いてモデル化される。この近似は、時空の曲率半径が重力波の波長よりも十分に大きいことを前提としている。これはマクロな物体や平均的な宇宙曲率に対しては成立するが、物質の微視的な離散性を考慮する場合、この前提が崩れる可能性がある。個々の粒子（電子、バリオン、ダークマター）は、その位置において時空曲率の鋭いスパイクを生成する。平均的な曲率は無視できるほど小さいが、粒子の位置における局所的な曲率は、特に粒子がそのコンプトン波長内に局在している場合、非常に大きくなる可能性がある。

標準的な扱いでは、運動方程式におけるリーマンテンソルへの重力波の直接的な結合を無視し、背景メトリックに焦点を当てる。しかし、重力波は電磁波よりも波長が長くエネルギーが低いため、これらの曲率項が伝播、具体的には角直径距離に有意な影響を与える可能性がある。先行研究[8]では、物質の離散性が宇宙マイクロ波背景放射（CMB）への角直径距離を著しく変化させ、特定のダークマター質量範囲を排除する可能性が示唆されている。本論文では、同様の効果が重力波にも適用されるか、また「ポスト幾何光学」近似が曲率スパイクが存在する状況下で依然として有効であるかどうかを調査する。

方法論
著者らは、[8]における電磁波のために開発されたポスト幾何光学近似を重力波に適用する。その手法は以下の通りである：

1. 運動方程式： 線形化アインシュタイン方程式から出発し、逆転トレース・メトリック摂動 $\bar{h}_{\alpha\beta}$ の波動方程式を導出する。電磁波の場合とは異われ、重力波の方程式には、テンソル $M_{\alpha\beta\mu\nu}$ で表される、リーマンテンソル（具体的にはワイルテンソル $C_{\alpha\mu\beta\nu}$ およびリッチテンソル $R_{\alpha\beta}$）への直接的な結合が含まれる。
2. ポスト幾何学的展開： 解を $\bar{h}_{\alpha\beta} = \sum \text{Re}(\bar{h}_{\alpha\beta}^n \epsilon^n e^{iS/\epsilon})$ という級数で展開する。ここで $\epsilon$ は波長と他のスケールの比である。決定的な違いとして、標準的な幾何学的光学とは異なり、曲率項 $M_{\alpha\beta\mu\nu}$ は $\epsilon^{-2}$ のオーダーとして扱われる。つまり、主要な運動項に対して無視できるとは仮定されない。
3. 分散関係： 先行項（$\epsilon^{-2}$）を分析することで、修正された分散関係を導出する。物質内部では、有効なグラビトンの質量平方が負となり（$k^2 = -2\rho/3M_{Pl}^2$）、これは波のベクトルがタイムライク（時間的）になり、重力波が減速することを意味する。
4. 角直径距離の計算： レイ束（ray bundle）の面積膨張率を積分することにより、角直径距離（$D_A$）の補正を計算する。これには、粒子の幅をコンプトン波長としたガウス型密度プロファイルを仮定した上で、密度勾配の二乗（$\partial_r \rho$）および密度自体の項を評価するプロセスが含まれる。
5. 妥当性の検証： 誘起された波面の方向の変化（偏向）が、位相変化と比較して十分に小さいかどうかを確認することで、局所平面波近似の整合性を検証する。特に、急峻な密度勾配が存在する場合に注目する。

主要な貢献と結果

• 修正された分散と質量： 本研究は、物質内部において、重力波がエネルギー密度 $\rho$ に比例する有効質量項を獲得することを裏付けている。分散関係は $E^2 = |\vec{k}|^2 + 2\rho/3M_{Pl}^2$ である。これにより、密度勾配によって駆動される非測地線運動が生じる。
• 角直径距離の補正： 角直径距離への補正は、物質密度とその導関数の二乗に依存する。コンプトン波長内に局在する単一粒子の場合、補正項は $m^6 / (M_{Pl}^2 E^4)$ のようにスケールする。
• 有効領域とコースティック（焦散）の形成：
  • 著者らは、LIGO–Virgo–KAGRAで検出される重力波（周波数 $\sim 10$ Hz）に対して、この近似が粒子質量 $m \lesssim 100$ MeVにおいて有効であることを発見した。
  • しかし、電子の質量（$m_e \approx 0.5$ MeV）の場合、曲率効果があまりに大きく、極めて短距離（$\sim 0.1$ pc）で**コースティックス（集束による特異点）**の形成を招く。
  • コースティックスの形成は、局所平面波近似の破綻を意味する。したがって、曲率効果（およびそれによる距離補正）が大きくなる領域において、ポスト幾何光学近似は適用不可能となる。
• 電磁波との比較：
  • 光子とは異なり、重力波はすべての物質（電子やバリオンを含む電荷を持つ粒子）と相互作用し、それらを通過する。
  • バリオンは現在の重力波検出器の有効領域外であるが、電子はそうではない。電子の累積的な影響は、サブ宇宙論的スケールでの強力な集束とコースティックス形成を引き起こす。
  • 重力波方程式に存在し、電磁波のベクトルポテンシャル方程式には存在しないワイルテンソルは、有効質量と偏光の進化に寄与する。しかし、検討されたガウス型プロファイルにおいては、ワイルの寄与はリッチテンソル／密度勾配項よりも副次的である。
• 偏光： 解析は、偏光ベクトルが平行移動ではなく、円錐双屈折を示すことを示唆している。ただし、現在の近似における階層構造の式は偏光に関して閉じていないため、この効果の定量的な評価は困難である。

意義と結論
本論文は、「ポスト幾何光学近似の興味深い領域」――すなわち、曲率効果が有意であり、角直径距離を変化させるほど大きい領域――は、自己限定的であると結論付けている。すなわち、曲率効果が大きくなると、ちょうどその時に近似が（コースティックス形成を通じて）破綻するのである。

• 限界： 本研究は、局在化した粒子（電子など）からの曲率スパイクが、宇宙論的な距離にわたる重力波伝播に与える影響を、現在の近似を用いて信頼性高く評価することはできないことを明らかにしている。局在化された粒子の曲率スパイクが、局所平面波の仮定を破る（＝近似が適用できなくなる）のは、誘起されたグラビトン質量が波のエネルギーと同程度になるか、あるいは集束が極端になった時である。
• 今後の課題： 著者らは、曲率スパイクの影響を評価し、粒子の局在化スケール（デコヒーレンスや相互作用によって決定される）をより正確に理解するために、新しい手法が必要であると強調している。局在化スケールは極めて重要である。ここで仮定されたコンプトン波長を用いることは、その効果を最大化するものである。もし粒子がより広がっていれば効果は減少するかもしれないが、より局在していれば、近似の破綻はより深刻なものとなる。
• 控えめな主張： 本論文は、特定のダークマターモデルを排除したり、新しい観測的シグネチャーを発見したりすることを目的としたものではない。むしろ、電磁波に対してこのようなシグネチャーを予測するために用いられた理論的枠組み[8]が、重力波に適用した場合、コースティックスの形成によって根本的な妥当性の問題に直面することを実証している。著者らは、有効領域は限定的であり、曲率効果が大きくなる場合には近似が適用できないことを強調している。