On the Possibility of a Strong First-Order Phase Transition in Neutron Stars

重力波とX線観測による中性子星データに対し、カイラル有効場理論および摂動論的量子色力学からの理論的制約を併用したベイズ推論を行うことで、本研究は、最も重い中性子星の中心密度を超えた領域で発生する可能性が高い高密度の強い一次相転移を支持する証拠を見出しており、それによって、硬い状態方程式の必要性と漸近的な軟化との両立を実現している。

要約

宇宙が、中性子星（死んだ星の崩壊した核）の中にのみ存在する、謎めいた超高密度の「宇宙の生地（cosmic dough）」で満たされていると想像してみてください。何十年もの間、物理学者たちは、この「生地」をより強く、より強く押しつぶしていくと、それがどのように振る舞うのかを正確に解明しようと試みてきました。

この論文は、一種の高レベルな探偵物語のようなものです。著者たちはある特定の謎を解こうとしています。それは、この「宇宙の生地」が、突発的で激しい方法でその質感を変えるのか（「強い一次相転移」）、それとも単にゆっくりと密度が増し、滑らかになっていくだけなのか？ という謎です。

以下に、簡単な比喩を用いた彼らの調査の全容をまとめます。

1. 謎：宇宙の「質感」

中性子星内部の物質を、ゼリーの塊だと考えてみてください。

• 「滑らかな」理論 (NPT): ゼリーを押しつぶしていくと、どんどん硬くなっていくものの、ずっとゼリーのままの状態であると考える科学者たちがいます。これは滑らかな変化です。
• 「突発的な」理論 (FOPT): 他の科学者たちは、ある圧力に達したとき、ゼリーが突然全く別の状態へと変化する――例えば、一瞬にして岩やガスに変わってしまうようなことが起こると考えています。物理学の用語では、これは「強い一次相転移」と呼ばれます。この論文では、これを、ある一定の密度の範囲において、物質が「跳ね返る」能力を失う（音速がゼロになる）瞬間として定義しています。

2. 手がかり：星の声を聞く

著者たちは、中性子星の中に入って直接確認することはできませんでした。そのため、彼らは主に2つの情報源から手がかりを集める探偵として行動しました。

• 「押しつぶし」テスト (重力波): 2つの中性子星が衝突したとき（GW170817と呼ばれるイベント）、彼らは空間に波紋を送りました。衝突する前に星がどれくらい「押しつぶされたか」は、その内部の「生地」がいかに硬いか、あるいは柔らかいかを知る手がかりとなります。
• 「懐中電灯」による測定 (NICER): NICERという宇宙望遠鏡は、いくつかのパルサー（回転する中性子星）の写真を撮影しました。チームは、これらの星のサイズと重さを測定することで、圧力がかかったときの「生地」の振る舞いをより正確に把握しました。
• 「実験室」のルール: 彼らはまた、2組の理論的なルールも使用しました。
  • 低密度ルール: 原子核を用いた実験（カイラル有効場理論）に基づいています。
  • 高密度ルール: 粒子が極限まで押しつぶされたときの振る舞いを記述する数学（摂動QCD）に基づいています。

3. 調査：デジタル・シミュレーション

著者たちは、「ベイズ推論」と呼ばれる手法を用いた大規模なコンピュータ・シミュレーションを構築しました。これは、膨大な数の異なるシナリオを実行して、どのシナリオが最も手がかりに合致するかを確認する作業だと考えてください。

• 彼らは、生地が滑らかに変化するシナリオ（相転移なし）と、突発的に変化するシナリオ（相転移あり）の2つのグループを作成しました。
• そして、現実世界のデータ（衝突の波と星の測定値）をすべてシミュレーションに投入し、どちらのシナリオのグループがより真実に近いかを検証しました。

4. 判決：「スナップ」は存在するが、隠されている

結果は驚くべき、かつ具体的なものでした。

• 「スナップ」は実在する: データは、突発的な「スナップ（相転移）」が実際に起こるという考えをわずかに支持しています。それは、最後までずっと滑らかなゼリーではありません。
• 「スナップ」は深い場所にある: ここにひねりがあります。この転移は、私たちが容易に観察できる星の外層部では起こりません。データは、この「スップ」が星の非常に深い内部、つまり最も重い中性子星の中心部で起こることを示唆しています。
  • 比喩: 重い金属の球を想像してください。外側は滑らかで硬いです。「スナップ」は、球を非常に強く押しつぶして、その核の部分が別のものに変わってしまう場合にのみ起こります。現在の私たちの観測は、球の外側しか見ていないため、その変化を直接見ることはできないのです。
• なぜこれが重要なのか: この発見は一つのパズルを解きました。「滑らかな」理論は、中性子星が崩壊せずにこれほど重くなれる理由を説明するのに苦労しますが、「突発的な」理論は通常、星を重さを支えられないほど柔らかくしてしまいます。この「スナップ」を（星の外形にあまり影響を与えない）深い中心部に配置することで、著者たちは、高密度の物理法則に従いつつ、重い星を維持できる方法を見出したのです。

5. これが将来何を意味するか

この論文は、現在の観測技術ではこの「スナップ」を見ることはできないものの、それはおそらく私たちの手の届かないところに存在していると結論づけています。

• 「双子星」の神話: この研究によれば、この「スナップ」によって、一部の人々が予想していたような「双子星」（重さは同じだがサイズが異なる2つの星）は、おそらく生まれないとのことです。
• 次なる手がかり: この転移を実際に「見る」ためには、中性子星の衝突の直後の様子を観察する必要があります。2つの星が合体するとき、そこには、安定した星では決して到達できないほど深い密度を持つ、超高密度の残骸が一時的に形成されます。将来の検出器が、これらの衝突の「鳴り響き」を聴き取ることができれば、この宇宙の生地が「スナップ」する音をついに捉えられるかもしれません。

要約すると、 著者たちは星のデータを用いて、宇宙で最も高密度な物質の「レシピ」を推測しました。彼らは、物質が最も重い星の深部で劇的な変化を起こす可能性が高いことを見出しました。それは、星が崩壊するのを防ぎつつ、物理学の法則を満たすための、宇宙の秘密なのです。

技術概要

技術要約：中性子星における強一次相転移の可能性について

問題提起
核飽和密度（$n_0 \approx 0.16 \text{ fm}^{-3}$）の数倍に達する高密度領域における強相互作用物質の相構造は、核物理学および天体物理学における中心的な未解決問題である。格子QCDおよび重イオン衝突実験により、高温・低バリオン化学ポテンシャルにおけるハドロン物質からクォーク・グルーオン・プラズマへの転移は解析的なクロスオーバーであることが確立されているが、低温・高バリオン密度の領域における性質は、第一原理に基づく格子計算における符号問題のため、未知である。具体的には、冷たい高密度QCD物質が、有限の密度区間にわたって音速が消失する（$c_s^2 = 0$）強一次相転移（FOPT）を起こすかどうかが不明である。これまでの研究では、ベイズモデル選択を用いてFOPT仮説と相転移なし（NPT）仮説を直接比較した例は一般的ではなく、また、FPTの開始密度を安定な中性子星内で既に探査されている密度に限定している場合も多かった。

手法
著者らは、冷たく$\beta$平衡状態にある中性子星物質に対して、非パラメトリックなガウス過程（GP）を用いた非パラメトリックな方程式の状態（EOS）の構成を行い、ベイズ推論を実施している。その手法は以下の通りである：

1. EOSの構築: 因果律（$0 \le c_s^2 \le 1$）を強制するために、補助変数 $\phi(n) \equiv -\ln[1/c_s^2(n) - 1]$ に対してGP回帰を用いて、音速の二乗 $c_s^2(n)$ をモデル化する。
   • NPT仮説: GPは全密度範囲にわたって滑らかな $c_s^2(n)$ を生成する。
   • FOPT仮説: モデルは一様分布から2つの密度 $n_S$ と $n_E$ をサンプリングする。音速は $n_S \le n \le n_E$ においてゼロ（$c_s^2=0$）に設定され、転移の両側で独立したGP実現が行われる。
2. 制約条件: 推論には複数のデータソースを組み合わせている：
   • 観測データ: 連星中性子星合体GW170817からの潮汐変形能、およびPSR J0740+6620、PSR J0030+0451、PSR J0437−4715、PSR J0614−3329に関するNICERによる質量・半径測定。
   • 理論的制約: $1.5 n_0$ 以下でのカイラル有効場理論（ChEFT）による制約、および高密度における摂動論的QCD（pQCD）による制約。
3. 密度範囲: GPは密度 $n_L$ で終了する。主要な解析では、ChEFTとpQCDの領域間の全窓をカバーするため $n_L = 25 n_0$ を使用し、比較のための二次解析として $n_L = 12 n_0$ を用いる。
4. 統計的枠組み: 階層的ベイズ法を用いて、事後分布および2つの仮説を比較するためのベイズ因子（$B_{NPT}^{FOPT}$）を計算する。解析では、「FOPT-in」（転移が最も重い安定中性子星の内部で起こる、$n_S < n_c$）と「FOPT-out」（$n_S \ge n_c$）を区別する。

主な貢献と結果

• FOPTの証拠: $n_L = 25 n_0$ の場合、データはNPT仮説に対してFOPT仮説を支持する中程度の証拠を示しており、ベイズ因子は $B_{NPT}^{FOPT} = 3.2$ である。
• 転移の場所: FOPT仮説内において、データは「FOPT-out」シナリオ（$B_{in}^{out} = 5.1$）を強く支持しており、これは相転移の開始（$n_S$）が最も重い中性子星の中心密度（$n_c$）よりも高い位置にあることを示唆している。
  • 推定された開始密度の中央値は $n_S \approx 7.35 n_0$ （68%信用区間は $5.30$–$11.85 n_0$）であり、最大質量星の中心密度は$n_c \approx 6.13 n_0$ である。
  • 解析を $n_L = 12 n_0$ に制限した場合、FOPTへの証拠は逸話的なレベル（$B_{NPT}^{FOPT} = 2.8$）まで弱まるが、転移が $n_c$ より上で起こるという傾向（$B_{in}^{out} = 2.6$）は維持される。
• 恒星観測量:
  • 標準的な質量: 転移は安定な恒星の内部より高い密度で起こると推定されるため、標準的な質量の観測量（半径 $R_{1.4}$ および潮汐変形能 $\Lambda_{1.4}$）は、FOPTとNPTの仮説間で68%信用レベルにおいて統計的に区別できない。
  • 最大質量: FOPT仮説は、最大非回転中性子星質量（$M_{TOV}$）のわずかな増加を許容する。$n_L = 25 n_0$ のとき、$M_{TOV}$ は $2.07^{+0.10}_{-0.08} M_\odot$ (NPT) から $2.15^{+0.13}_{-0.11} M_\odot$ (FOPT) へとシフトする。
  • ツイン・スター: 「ツイン・スター」解（離れた枝）の事後確率は $\le 0.1\%$ と算出され、このシグネチャーは強く否定される。
• 音速と共形性: FOPT仮説の下では、音速 $c_s^2$ は $4 n_0$ 付近でピークに達し、転移プラトーを通じてゼロに低下した後、pQCD限界に近づく前に小さなピークへと再上昇する。トレース・アノマリー $\Delta$ は、FOPTシナリオ下であっても、中性子星内部の物質は非共形的であり、強く結合していることを示している。

意義と主張
本論文は、主に最も重い安定な中性子星の中心密度よりも高い密度で発生すると推定されるFOPTが、物理的な緊張関係を自然に解決することを主張している。すなわち、$\sim 2 M_\odot$ の中性子星を支えるための硬いEOSの必要性と、漸近的に高密度になるpQCDが好む軟化との間の矛盾である。軟化（$c_s^2=0$ のプラトー）を安定な恒星の内部の外側に配置することで、モデルは重い恒星に必要な剛性を維持しつつ、高密度のpQCD制約を満たすことができる。

著者らは、現在の安定な中性子星の観測からは、この特定のタイプのFOPTを直接検出することはできないと結論付けている。代わりに、直接的なテストは、$n_c$ を超える密度に一時的に到達する系、例えば連星中性子星合体の過渡的な残骸などを対象とすべきであると示唆している。次世代の検出器による合体後の信号の重力波観測が、この結論を検証するための自然な候補として特定されている。本論文はFOPTの存在を決定的に証明したと主張しているのではなく、広い密度範囲にわたってテストした結果、現在のデータによって中程度に支持される一つの仮説として提示しているものである。