Numerical calculations of neutron star mountains supported by crustal lattice pressure

Brussels-Montreal 核エネルギー密度汎関数を用いた一貫した計算により、磁場を帯びた降着中性子星の熱勾配によって支持される「山」の形成を詳細に検討した結果、これらが低質量 X 線連星の自転平衡を支配するには小さすぎることが示された。

要約

この論文は、「中性子星（ちゅうせいしん）」という宇宙の超巨大な星が、なぜ「山」を作ってしまうのか、そしてその山が重力波（じゅうりょくは）という宇宙のさざなみを生み出すのかについて、新しい視点から詳しく調べた研究です。

少し難しい専門用語を、身近な例え話に置き換えて解説しましょう。

1. 舞台設定：宇宙の「超硬いパンケーキ」

まず、中性子星とは何か想像してみてください。太陽の質量を、東京ドームくらい（半径 10km 程度）の小さな空間にギュッと押し込んだような星です。
この星の表面（地殻）は、**「世界で最も硬いパンケーキ」**のようなものです。

• 硬さ： 鉄よりも 100 億倍も硬いと言われています。
• 特徴： この硬い殻（地殻）が、星の回転によって歪められると、表面に小さな「山」ができます。

この「山」ができて星が回転すると、重力波という「時空の波」が生まれます。この波を捉えられれば、宇宙の謎が解けるかもしれません。

2. 過去の考え方：「電子のキャッチボール」

これまでの研究では、この「山」ができる主な原因は、**「電子のキャッチボール」**だと思われていました。

• 仕組み： 星の表面に物質が降り積もると、圧力が高まり、原子核が「電子を飲み込む（電子捕獲）」反応が起きます。
• 山ができる理由： 星の温度にムラがあると、電子が飲み込まれる場所（深さ）がズレてしまいます。これが「波打つ層」を作り出し、結果として星の形が歪んで「山」ができる、という考え方でした。

しかし、最近の計算では、この「電子捕獲」が起きる場所や回数が、昔考えられていたよりもはるかに少ないことがわかってきました。つまり、「電子キャッチボール」だけでは、大きな山は作れないかもしれないのです。

3. 新しい発見：「熱いパンケーキの膨らみ」

この論文の著者たちは、新しいアイデアを提案しました。
電子のキャッチボールではなく、**「熱いパンケーキが膨らむ」**という現象に注目したのです。

• 仕組み： 中性子星の表面は、内部の磁場（じこば）の影響で、場所によって温度にムラができます（片側が熱く、片側が冷たい）。
• 山ができる理由： 硬いパンケーキ（中性子星の地殻）は、**「熱い部分は膨らみ、冷たい部分は縮む」**性質を持っています。
  • 熱い場所の原子の結晶（レゴブロックのようなもの）が少し膨らむと、その部分だけ盛り上がります。
  • これを**「結晶格子圧力（けっしょうこうしあつりょく）」**と呼びます。

著者たちは、この「熱による膨らみ」が、実際にどれくらいの「山」を作れるのか、最新の計算機を使って詳しくシミュレーションしました。

4. 結果：「山は小さすぎる」

結論から言うと、**「残念ながら、この仕組みで作られる山は、重力波を捉えるには小さすぎる」**という結果になりました。

• 大きさ： 計算された「山」の高さは、非常に微細です。
• 比較： 重力波を捉えるには、もっと大きな「山」（理論上の限界値）が必要ですが、今回の計算では、その100 万分の 1程度の大きさしかありませんでした。
• 意味： 現在の重力波検出器（LIGO など）では、この小さな「熱いパンケーキの膨らみ」が作る重力波を見つけるのは、まだ無理そうです。

5. なぜこの研究は重要なのか？

「小さくて見つからないなら、意味がないのでは？」と思うかもしれません。しかし、この研究には大きな意義があります。

1. 「最大限の山」の限界を知る：
   宇宙には「電子捕獲」や「熱膨らみ」など、山を作る仕組みがいくつかあります。どれがどれくらい効果があるかを正確に計算することで、「中性子星が作れる最大の山はどれくらいか」という**「天井（限界値）」**をより正確に知ることができます。
2. 将来の探査への道しるべ：
   もし将来、重力波検出器がもっと敏感になったら、この小さな「熱い山」を見つける日が来るかもしれません。そのために、正確な計算データが必要なのです。
3. 他の可能性への示唆：
   この研究は、低質量 X 線連星（LMXB）という特定の星に焦点を当てましたが、もっと強い磁場を持つ「マグネター」や、超高速で物質を吸い込む「超輝性 X 線源（ULX）」では、もっと大きな山ができる可能性が残されています。

まとめ

この論文は、「中性子星の表面が熱くなって膨らむことで山ができる」という新しい仮説を、最新の計算で検証したものです。

結果は「残念ながら、今の技術では小さすぎて見つけられない」でしたが、「宇宙の硬いパンケーキが、どのように歪むのか」というメカニズムを、これまでで最も正確に解明したという点で、重力波天文学の重要な一歩となりました。

まるで、「宇宙の山がどれくらい高いか」を測るための、より正確な定規を作ったような研究だと言えます。

技術概要

論文の技術的サマリー：「格子圧力に支えられた中性子星の山（マウンテン）の数値計算」

論文タイトル: Numerical calculations of neutron star mountains supported by crustal lattice pressure
著者: T. J. Hutchins, D. I. Jones
掲載誌: MNRAS (2026 年予稿)

1. 研究の背景と問題提起

連星 X 線源（LMXB）に存在する中性子星の自転速度は、重力波（GW）によるエネルギー放射によって平衡状態（トルクバランス）に達している可能性があります。この重力波放射を駆動するメカニズムとして、中性子星の固体地殻に生じる質量非対称性、すなわち「山（マウンテン）」の形成が注目されています。

これまでの研究では、以下の 2 つの点に課題がありました。

1. 非対称性の発生源の分離: 温度勾配による弾性ひずみ（熱的マウンテン）の形成過程において、温度非対称性の発生源と、それによって生じる弾性ひずみのメカニズムが別々に扱われてきた。
2. 電子捕獲モデルの限界: 従来の主流モデル（UCB モデル）は、温度依存性を持つ電子捕獲反応（キャプチャレイヤーのシフト）を質量非対称性の源としていた。しかし、最新の核物質の状態方程式（EOS）計算では、電子捕獲反応が予測より遥かに少なく、かつ閾値エネルギーが低いため、UCB モデルが予測するほど大きな楕円率（$\varepsilon$）が得られない可能性が指摘されていた。

本研究は、これらの課題を解決し、**「磁場誘起の異方性熱伝導」と「地殻格子の熱的圧力」**を組み合わせた、自己整合的な熱弾性マウンテンの形成シミュレーションを行うことを目的としています。

2. 手法とモデル

本研究では、以下のステップで自己整合的な計算を行いました。

2.1 状態方程式（EOS）の統一

中性子星の全領域（外層地殻、内層地殻、コア）を熱力学的に整合的に記述する、ブリュッセル・モントリオール核エネルギー密度汎関数（BSk19, BSk20, BSk21）を使用しました。これにより、従来のように異なる EOS をつなぎ合わせることで生じる不整合や数値的不安定性を回避しました。

• 数値計算の安定化のため、離散的な EOS データを解析的な関数（Potekhin et al. 2013 の形式を拡張）で近似し、密度の不連続点（核種変化界面）における数値積分の収束問題を解決しました。

2.2 熱構造と温度非対称性の計算

• 背景熱構造: 定常的な降着を行う中性子星の熱構造を、深層地殻加熱（電子捕獲反応など）と浅層地殻加熱（未解明の加熱源）を考慮して計算しました。
• 温度非対称性の生成: 内部磁場（地殻内またはコア内）による相対論的電子の異方性熱伝導を考慮し、温度勾配（$\delta T$）を計算しました。磁場強度は地殻で $10^{12}$ G、コアで $10^8$ G を仮定しました。

2.3 質量非対称性（山）の生成メカニズム

従来の「電子捕獲レイヤーのシフト」に代わり、**「地殻格子圧力（Crustal Lattice Pressure）」**を主要なメカニズムとして採用しました。

• 格子圧力は、イオン格子自体の熱的振動に由来する圧力成分であり、温度に依存します。
• 温度非対称性 $\delta T$ が格子圧力 $\delta P_{\text{th}}$ の非対称性を生み出し、これが地殻の弾性変形（ラグランジュ変位 $\xi_i$）を引き起こし、最終的に質量四重極モーメント（$Q_{22}$）および楕円率 $\varepsilon$ を生成します。
• Ushomirsky et al. (2000) が開発した摂動理論の枠組みを流用しつつ、温度摂動の定義（固定された地殻に対する摂動 vs 弾性変形後のラグランジュ摂動）を厳密に再定義し、自己整合的な解を得ました。

3. 主要な結果

3.1 楕円率の推定値

計算された「磁気 - 熱 - 弾性マウンテン」の楕円率 $\varepsilon$ は、以下の条件で最大となりました。

• 最大楕円率: $\varepsilon_{\text{max}} \approx 2.32 \times 10^{-11}$
  • 条件：BSk21 EOS、降着率 $\dot{M} = 4.4 \times 10^{-9} M_\odot/\text{yr}$、浅層加熱 $Q_S = 10$ MeV、地殻磁場 $B = 2 \times 10^{12}$ G。
• 一般的な傾向: 楕円率は降着率や浅層加熱量の増加とともに増大しますが、降着率が極端に高い場合（$\sim 5 \times 10^{-9} M_\odot/\text{yr}$ 以上）、地殻の一部が融解し、寄与する領域が制限されるため、増加が頭打ちになります。

3.2 トルクバランス限界との比較

重力波放射による自転停止（トルクバランス）を説明するために必要な楕円率は、観測された降着率と自転速度から推定される $\varepsilon_{\text{TB}} \sim 10^{-9} \text{--} 10^{-7}$ です。

• 本研究で計算された $\varepsilon \lesssim 10^{-11}$ は、トルクバランス限界よりも2〜6 桁小さい値です。
• したがって、このメカニズム単独では、LMXB 中性子星の自転速度を重力波によって制限することはできないと結論付けられます。

3.3 電子捕獲モデルとの比較

最新の EOS（Gusakov & Chugunov 2020 など）に基づき、電子捕獲レイヤーシフトによる楕円率を再評価しました。

• 従来の UCB モデル（閾値エネルギー 90 MeV を仮定）では $\varepsilon \sim 10^{-7}$ 程度が期待されていましたが、現実的な EOS（閾値エネルギー $\sim 25$ MeV、内層地殻での捕獲レイヤー欠如）では、$\varepsilon \lesssim 10^{-10}$ 程度に減少します。
• 格子圧力メカニズムによる楕円率（$\sim 10^{-11}$）は、現実的な電子捕獲モデルによる値と同程度か、やや小さいレベルです。

4. 結論と意義

1. メカニズムの検証: 磁場誘起の温度非対称性と、地殻格子圧力を組み合わせた自己整合的な計算を初めて行い、このメカニズムが物理的に可能であることを示しました。
2. 観測への示唆: 現在の重力波検出器（LIGO/Virgo/KAGRA）の感度では、LMXB 中性子星からこのメカニズムによる連続重力波を検出することは極めて困難です。検出には、より高感度な将来の検出器（Einstein Telescope や Cosmic Explorer）が必要となる可能性があります。
3. モデルの限界と将来展望:
   • 本研究はニュートン重力近似で行われており、一般相対論的効果（重力ポテンシャルの摂動など）を完全に含んでいません。これを含めると楕円率はさらに小さくなる可能性があります。
   • 降着率が高く、磁場が極めて強い「超輝 X 線源（ULX）」のようなシステムでは、温度非対称性や格子圧力がさらに増大する可能性があり、ULX 中性子星が重力波源としての候補となり得ることを示唆しています（ただし、ULX は自転が遅いため、信号強度は不明瞭です）。

総括:
本研究は、中性子星の「山」形成において、電子捕獲反応に依存しない新しい経路（格子圧力）を詳細に検討しましたが、その結果として生じる質量非対称性は、現在の重力波観測で検出可能なレベル（トルクバランス限界）には達しないことが示されました。これは、中性子星の自転速度制限に重力波が寄与している場合、より大きな変形（最大山限界に近いもの）を生み出す別のメカニズム、あるいはより極端な環境（ULX など）が存在する可能性を示唆しています。