Constraints on a fifth force from the stellar orbits around the central supermassive black hole of the Milky Way

この論文は、銀河中心の恒星 S2 の軌道観測データをマルコフ連鎖モンテカルロ法で解析し、ヤウカ重力モデルに基づく第五の力のパラメータ（強度と範囲）に制約を与え、その結果が一般相対性理論の予測や既存の研究と整合的であることを示した。

要約

この論文は、銀河の中心にある「巨大なブラックホール」の周りを回る星の動きを詳しく調べることで、「第五の力（フィフス・フォース）」という、まだ見えない不思議な力が存在するかどうかを探った研究です。

まるで探偵が、犯人（第五の力）の足跡を追いかけようとしているような物語です。以下に、難しい専門用語を噛み砕いて、身近な例え話で解説します。

1. 物語の舞台：銀河の中心と「S2」という星

私たちの銀河（天の川銀河）の中心には、太陽の 400 万倍もの質量を持つ**「超大質量ブラックホール」が鎮座しています。その周りを、「S2」という星**が猛スピードで公転しています。

• 通常の話（一般相対性理論）： アインシュタインの理論では、この星の動きは「重力」という目に見えない紐で引っ張られているだけで、完璧に予測できます。
• 今回の疑問： でも、もしそこに**「第五の力」**という、新しい目に見えない力が働いていたらどうなる？その力が星の軌道を少しだけ変えていないか？

2. 第五の力とは？「風」の例え

第五の力は、重力とは少し違う性質を持っています。

• 重力： 物を引き寄せる力（例：磁石の N 極と S 極）。
• 第五の力： 物を**「押す」力**（反重力のようなもの）や、距離によって強さが変わる不思議な力です。

これを**「風」**に例えてみましょう。

• 星がブラックホールを回る際、通常は「重力」という強いロープで引っ張られています。
• もし第五の力があれば、それは**「星を押し流す風」**のようになります。
  • 風の強さ（δ）： 風がどれくらい強いのか。
  • 風の届く範囲（λ）： 風が吹く距離がどれくらいか（数百キロ先までか、それとも銀河全体か）。

3. 研究の方法：3 つの「風のシナリオ」

研究者たちは、この「風」がどれくらいの範囲で吹いているかを 3 つのシナリオで試しました。

1. ケース 1（小さな風）： 星の軌道の**「内側」**だけ吹いている（数百 AU 程度）。
2. ケース 2（中くらいの風）： 星の軌道と**「同じくらい」**の範囲で吹いている（約 1000 AU）。
3. ケース 3（大きな風）： 星の軌道よりも**「ずっと広い」**範囲で吹いている（数千 AU）。

彼らは、コンピュータで「もしこの風が吹いていたら、星はどんな軌道を描くか？」をシミュレーションし、実際の観測データ（S2 星の動き）と照らし合わせました。

4. 見つかった答え：風は「あるかもしれない」が、弱すぎる

結果はどうだったでしょうか？

• 風の強さ（δ）：

  • 風が届く範囲（λ）が広くなるほど、風は強くなる傾向があることがわかりました。
  • 具体的な強さは、ケース 1 で約 0.005、ケース 3 で約 0.15 でした。
  • しかし！ これらの値には**「誤差（不確実性）」が非常に大きく**含まれていました。つまり、「風があるかもしれないし、ないかもしれない（重力だけかもしれない）」という状態です。

• 他の研究との一致：

  • この結果は、他の研究者たち（GRAVITY コラボレーションなど）が異なる方法で出した結果とよく一致していました。これは、「第五の力」のモデル（風の吹き方）をどう変えても、結論はあまり変わらないことを示しています。

5. 重要な発見：アインシュタインの理論は「まだ勝っている」

研究の最後には、**「アインシュタインの理論（一般相対性理論）の予測と、観測結果のズレ」**についてもチェックしました。

• 星の軌道は、ブラックホールの周りを回るたびに少しずれていきます（歳差運動）。
• このズレの大きさを測る指標（fSP というパラメータ）を調べたところ、「アインシュタインの理論の予測」と「観測結果」は、誤差の範囲内で完全に一致していました。
• つまり、「第五の力」が星の動きに大きな影響を与えている証拠は見つかりませんでした。

6. まとめ：まだ謎は解けていない

この論文の結論を一言で言うと、以下のようになります。

「銀河の中心には、もしかしたら『第五の力』という新しい風が吹いているかもしれない。でも、今の観測技術では、その風が『重力というロープ』に比べてどれくらい強いのか、はっきりと区別できないほど弱いか、あるいは存在しない可能性が高い。

もっと精密な望遠鏡や観測技術を使って、さらに小さな『風の揺らぎ』を見つけ出す必要があるよ！」

今後の展望：
将来的に、より高性能な観測装置（GRAVITY 装置の改良など）を使えば、この「風の強さ」をより正確に測れるようになるでしょう。もしそれができれば、宇宙の謎（ダークマターやダークエネルギー）を解く大きな手がかりになるかもしれません。

---

簡単な比喩のまとめ：

• ブラックホール ＝ 巨大な回転するテーブル。
• S2 星 ＝ テーブルの上を滑るビー玉。
• 重力 ＝ ビー玉をテーブルの中心に引き寄せる磁石。
• 第五の力 ＝ ビー玉を横から押す「風」。
• この研究 ＝ 「ビー玉の動きを見て、風が吹いているかどうかを調べたが、風は弱すぎて、磁石の力だけだと言わざるを得なかった（でも、もっと強い風が吹いていないか、もっと精密に調べる必要がある）」という話です。

技術概要

論文の技術的サマリー：銀河中心の恒星軌道から第五の力を制約する

1. 研究の背景と課題

本論文は、我々の銀河系中心部（GC）にある超大質量ブラックホール（Sgr A*）の周囲を公転する恒星（特に S2 星）の軌道データを用いて、「第五の力」の存在とその特性を制約することを目的としています。

• 第五の力とは: 銀河スケールでのニュートン重力を補完し、暗黒物質（DM）の仮説なしに銀河の回転曲線を説明したり、大規模スケールでダークエネルギー（DE）の効果を模倣したりする可能性のある反発力（反重力）です。
• 理論的枠組み: 大質量重力理論や拡張重力理論（ETG）の弱場極限において、ニュートン重力ポテンシャルに指数関数的なヤウカ（Yukawa）型の補正項が現れると予測されます。
  • ポテンシャル式：$\Phi(r) = -\frac{G_\infty M}{r} (1 + \delta e^{-r/\lambda})$
  • ここで、$\delta$ は第五の力の強度（結合定数）、$\lambda$ はその到達距離（範囲）を表します。
• 課題: 既存の研究では異なるヤウカ型ポテンシャルが用いられており、モデル依存性が懸念されていました。また、S2 星の軌道における一般相対性理論（GR）からの微小な乖離が、第五の力によって説明できるかどうかの検証が十分に行われていませんでした。

2. 研究方法

著者らは、以下の手法を組み合わせて分析を行いました。

2.1 拡張 PPN 形式を用いた軌道シミュレーション

• 運動方程式: 一般相対性理論（GR）の標準的なパラメータ化された後ニュートン（PPN）形式を拡張し、ヤウカ重力における運動方程式（式 5）を導出・使用しました。これにより、第五の力が存在する場合の S2 星の真の軌道をシミュレートします。
• 観測データとの比較: 得られたシミュレーション軌道を、S2 星の実際の天測観測データ（位置データ：[56] 参照）に適合させました。

2.2 マルコフ連鎖モンテカルロ（MCMC）解析

• パラメータ推定: 第五の力の強度 $\delta$ と範囲 $\lambda$、および S2 星の軌道要素（半長径 $a$、離心率 $e$、傾斜角 $i$ など）と Sgr A* の質量 $M$ を同時に推定するために MCMC 法（emcee パッケージ）を使用しました。
• 事前分布（Prior）の検討:
  • 一様事前分布（Uniform Priors）: パラメータの広い範囲を探索するために使用。
  • ガウス事前分布（Gaussian Priors）: 既存の観測結果に基づいたより狭い範囲を設定し、結果の感度（事前分布への依存性）を分析するために使用。
• 3 つのケース設定: 第五の力の範囲 $\lambda$ について、S2 星の軌道に対する相対的な大きさで 3 つのケースを設定しました。
  1. ケース 1: $\lambda \approx$ 数百 AU（S2 星軌道の深部）
  2. ケース 2: $\lambda \approx$ 1,000 AU（S2 星軌道サイズ程度）
  3. ケース 3: $\lambda \approx$ 数千 AU（S2 星軌道より遥かに大きい）

2.3 一般相対性理論からの乖離の検証

• GRAVITY 協力団が 2020 年に測定したシュワルツシルト歳差パラメータ $f_{SP} = 1.10 \pm 0.19$ を用いて、得られた $\delta$ と $\lambda$ の値が GR からの微小な乖離を説明可能か、式 (6) と (8) を用いて整合性を確認しました。

3. 主要な結果

3.1 第五の力の強度と範囲の関係

MCMC 解析の結果、以下の傾向が確認されました。

• 強度の増大: 第五の力の範囲 $\lambda$ が増加するにつれて、その強度 $\delta$ も増加する傾向が見られました。
  • ケース 1 ($\lambda \sim$ 数百 AU): $\delta \sim 0.005$
  • ケース 2 ($\lambda \sim$ 1,000 AU): $\delta \sim 0.02$
  • ケース 3 ($\lambda \sim$ 数千 AU): $\delta \sim 0.15$
• 誤差の減少: $\lambda$ が増加するにつれて、$\delta$ の相対誤差（$\Delta\delta/\delta$）は減少しました。
• モデル独立性: 異なるヤウカ型ポテンシャル（式 2, 3, 4）を用いた既存の研究や、GRAVITY 協力団の結果と比較しても、得られた制約は整合的であり、第五の力の強度はモデルに依存しない可能性が高いことが示唆されました。

3.2 事前分布への感度

• 一様事前分布 vs ガウス事前分布: ガウス事前分布を用いた場合、事後分布がより狭くなり、パラメータの誤差が大幅に減少しました。
• パラメータの依存性: 黒孔質量 $M$ や第五の力パラメータ $\delta, \lambda$ は、事前分布の選択に多少依存する傾向が見られましたが、軌道要素（$a, e, i$ など）は事前分布に依存しないことが確認されました。
• 相関関係: 事後分布の相関行列を分析した結果、$M$ と $\delta$ の間に一定の相関（特にガウス事前分布で約 50%）が見られましたが、$\delta$ と $\lambda$ の間には統計的に有意な線形相関は見られませんでした。ただし、理論的な関係式 (6) や非線形な相関の存在が示唆されました。

3.3 一般相対性理論（GR）との整合性

• 得られた $\delta$ と $\lambda$ の組み合わせは、GRAVITY 協力団が測定した $f_{SP} = 1.10 \pm 0.19$ の値と、誤差範囲内で整合していることが確認されました。
• 現在の観測精度では、GR の予測と第五の力による説明の区別は困難であり、GR のケースを排除することはできません。

4. 結論と意義

4.1 結論

• 銀河中心の S2 星の軌道データを用いることで、第五の力の存在を制約する新たな枠組みが確立されました。
• 範囲 $\lambda$ が大きくなるほど、許容される強度 $\delta$ も大きくなるという傾向が確認されました。
• 得られた制約は、異なる重力理論モデルや既存の観測結果（GRAVITY 協力団など）と矛盾しません。
• しかし、パラメータの誤差が依然として大きいため、第五の力の存在を確定したり、GR を排除したりするには至っていません。

4.2 科学的意義

• 重力理論の検証: 銀河中心という極限環境における重力理論の検証に新たなアプローチを提供しました。
• モデル比較: 異なるヤウカ型ポテンシャルを用いた研究との比較を通じて、第五の力の制約がモデルに依存しない普遍的な性質を持つ可能性を示唆しました。
• 将来の展望: 現在の制約は観測誤差に大きく依存しています。将来、より高精度な天測観測（特に $f_{SP}$ の誤差を $\pm 0.1$ や $\pm 0.05$ まで縮小する観測）が行われれば、第五の力の存在についてより厳密な制約が得られると期待されます。

本論文は、銀河中心の恒星軌道解析を通じて、現代物理学の未解決問題である「第五の力」や「暗黒物質/エネルギー」の正体解明に向けた重要なステップを示すものです。