In-medium hyperon potentials and the quarkyonic hyperon onset: charged $\Sigma$'s in $\beta$-equilibrium and the neutrino connection

技術要約：イン・メディアム（媒質中）ハイペロン・ポテンシャルとクォーキオニック・ハイペロンの出現

問題提起
$2 M_\odot$（太陽質量2倍）のパルサーの存在とGW170817による潮汐変形能の制約は、核飽和密度（$n_0 = 0.16 \text{ fm}^{-3}$）を超える密度において、高密度物質の状態方程式（EoS）が硬（スティフ）であり続ける必要があることを示している。標準的なハドロンモデルは「ハイペロン・パズル」に直面する。すなわち、中性子の化学ポテンシャルが上昇するにつれ、ハイペロン（例：$\Lambda, \Sigma$）を介してストレンジネスが現れ、EoSを軟化させ、通常、中性子星の最大質量を $0.2\text{--}0.4 M_\odot$ 減少させてしまう問題である。

Fujimoto, Kojo, and McLerran (FKM) によって提案された「クォーキオニック（Quarkyonic）」による解決策は、中性子が低運動量の$d$クォークの位相空間を事前占有することで、低運動量のハイペロンを統計的にブロックするというものである。これにより、ハイペロンの出現閾値は真空条件の $\mu_B = M_Y$ から $\mu_B = 2M_Y - M_N$ へとシフトし、出現を $\approx 5\text{--}6 n_0$ まで遅らせる。オリジナルのFKM「IdylliQ」モデルは、真空質量と自由分散を用いており、実験データやニュートリノ誘起の最終状態相互作用（FSI）によって制約されるイン・メディアム（媒質中）単一粒子ポテンシャル（$U_Y$）を無視していた。本論文は、これらのイン・メディアム・ポテンシャルを用いてIdylliQモデルを「ドレス（装飾）」し、ハイペロン・パズルに対するそれらの定量的影響を決定することで、このギャップに対処する。

手法
著者らは、クォーク状態が運動空間の低運動量内部を充填し、バリオンがフェルミ面付近のシェルとして存在するIdylliQモデルを用いる。本研究では、FKMフレームワークを修正し、ゼロ運動量におけるエネルギー分散に密度依存的な単一粒子ポテンシャル（$U_\Lambda, U_\Sigma, U_N$）を組み込んでいる。

主な手法ステップは以下の通りである：

1. ドレスされた閾値の導出： 著者らは、エネルギー分散への代入として $E_Y(0) \to M_Y + U_Y$ および $E_N(0) \to M_N + U_N$ を用いることで、ハイペロン出現の化学ポテンシャル条件を導出した。
2. 自己整合性： ハイペロン出現密度の計算において中性子平均場（$U_N$）を考慮し、中性子ポテンシャルを相互作用エネルギー密度の汎関数微分として扱うことで、熱力学的整合性を確保した。
3. $\beta$平衡解析： レプトンおよび荷電バリオン（$\Sigma^\pm$）を含むようにモデルを拡張し、$\beta$平衡における電子化学ポテンシャル（$\mu_e$）とバリオン分率を解いた。
4. 数値実装： 自己完結的なC++スイートを用い、ポテンシャルの「ドレス」を適用する前にFKMの真空結果を$<1.2\%$の精度で再現した。本研究では、最大質量を計算するためにトールマン・オッペンハイマー・ヴォルコフ（TOV）方程式を統合し、EoSの軟化を分析している。
5. ベイズ伝播： SBNDおよびDUNEのニュートリノ・プログラム（特に$U_\Lambda$に関する予測精度）から得られる実験精度を、結合ベイズフィットを通じて伝播させ、ハイペロンを含まないコアが存在する確率を評価した。

主要な貢献と結果

1. ハイペロン・ポテンシャルの重み2によるレバレッジ：
   ドレスされたクォーキオニック出現条件は次のように導かれる：
   $$\mu_{\text{onset}}^B = (2M_Y - M_N) + 2U_Y - U_N$$
   従来のハドロン閾値（$\mu_B = M_Y + U_Y$）とは異なり、イン・メディアム・ハイペロン・ポテンシャル $U_Y$ は2の重みで入る。これは、静止状態のハイペロンが放出する$d$クォークの位相空間が、中性子の半分に過ぎないためである。その結果、出現密度は標準的なモデルよりも $U_Y$ に対して2倍敏感である。$U_\Lambda$ の $10 \text{ MeV}$ の変化に対し、出現密度は従来モデルの $\approx 0.15 n_0$ に対し、$\approx 0.3 n_0$ シフトする。

2. 自己整合的中性子ポテンシャルとプロテクション・クリフ（保護の崖）：
   中性子ポテンシャル $U_N$ は、出現密度において重み **$-2$** で入る（化学ポテンシャルの$-U_N$シフトと、$\mu_B(n_B)$曲線の$+U_N$シフトの両方が組み合わさるため）。これにより「プロテクション・クリフ」が生じる。クォーキオニック機構が有効であるためには（出現密度をクォーク飽和密度以上に保つためには）、出現密度における中性子ポテンシャルが$U_N(\text{onset}) \lesssim +96 \text{ MeV}$（$U_\Lambda = -28 \text{ MeV}$ と仮定した場合）を満たす必要がある。中性子の反発が強すぎると、出現密度が押し下げられ、メカニズムが無効になる可能性がある。

3. $\beta$平衡における $\Sigma^-$ の反転と禁止：
   $\beta$平衡において、$\Sigma^-$ (dds) の閾値は純粋なレプトン・セクターの条件へと崩壊する：
   $$\mu_e \ge (M_{\Sigma^-} - M_N) + U_\Sigma - U_N \approx 258 \text{ MeV} + U_\Sigma - U_N$$
   $2 M_\odot$ のコア内部では、電子化学ポテンシャル $\mu_e$ は（相互作用を含むモデルであっても）この閾値に達せず（$\approx 225 \text{ MeV}$ 以下に留まる）、結果として $\Sigma^-$ は**事実上禁止（バン）**される。これによりハイペロンの順序が逆転し、ブロックされていない $\Sigma^+$ (uus) が最も安価な $\Sigma$ 種となる。この禁止は、$U_N - U_\Sigma \lesssim +35 \text{ MeV}$ である限り、EoSの硬さの変動に対して頑健である。

4. EoS軟化の抑制：
   FKMアンサッツのファミリー内では、ハイペロンによる最大質量減少は $< 0.05 M_\odot$（現実的な相互作用を持つ星では $< 0.025 M_\odot$）に抑えられる。これは、従来のハドロンモデルが予測する軟化と比較して、4〜8分の1の小ささである。論文では、この上限は検討された特定のモデルファミリーにおける「天井（シーリング）」であるが、$u$クォーク飽和密度以上のハイペロンを省略しているという点では、クォーキオニック物質全般に対する「床（フロア）」であると述べている。

5. ベイズ制約と「ハイペロン・フリー」コア：
   SBND+DUNEのFSI精度（$\delta U_\Lambda \approx 5.6 \text{ MeV}$）をベイズフィットを通じて伝播させると、低密度アンカー $U_\Lambda(n_0)$ は特定されるが、超飽和領域のターンオーバー・パラメータ $c_\Lambda$ はニュートリノ・データでは制約されない。重イオン実験の事前分布（プライア）を用いると、ハイペロンを含まないコアの確率は $P \approx 0.90$ となる。著者らは、これは直接的な組成測定ではなく、「事前分布に支配された記述」であることを強調している。

意義と主張
本論文は、イン・メディアム・ポテンシャルによるドレスが、単にFKMモデルを修正するだけでなく、その現象論を根本的に変えるものであると主張している。最も重要な知見は、微分可能な観測量である。すなわち、最大質量に対するハイペロン・ポテンシャルの感度（$dM_{\text{max}}/dU_Y$）は、平均場ハドロンシナリオよりもクォーキオニック・シナリオの方が一桁小さい。これは、ハイペロン・パズルの解決策を識別するための具体的な判別法を提供する。

著者らは、電荷を持つ $\Sigma$ ハイペロンやレプトンを考慮しても、中性子物質における $\Sigma$ の反発が測定値に従う限り、クォーキオニック・メカニズムは頑健であると断言している。本研究は、ハイペロン・パズルがクォーキオニック物質においては位相空間のブロッキングによって統計的に解決されることを確立しており、イン・メディアム・ポテンシャルは、出現密度を精密に決定する極めて敏感なレバーとして機能するものの、必ずしも $2 M_\odot$ の星を支える能力を破壊するものではない。結論として、コアの組成を決定する決定的な量は、依然としてハイペロン・核子の相互作用の超飽和挙動（$c_\Lambda$）であり、現在のニュートリノ・データではこれを制約できないとしている。