The crust of dark-matter admixed neutron stars: bulk properties and torsional oscillations

技術要約：ダークマター混入中性子星の地殻：バルク特性とねじれ振動

問題提起
中性子星は、重力捕獲によるダークマター（DM）の蓄積を研究するための自然な実験室として機能する。もしダークマターが中性子星の内部または周囲に蓄積されると、それはダークマター混入中性子星（DANS）を形成する。これまでの研究では、ダークマターが中性子星の全域的な質量-半径関係をどのように変化させるかについて検討されてきたが、固体の地殻およびそれに伴う地震学的シグネチャーに対するダークマターの具体的な影響については、まだ十分に探索されていない。本論文では、重力的に結合したダークマター流体が、バリオン地殻の厚さおよび、マグネター・フレアで観測される準周期振動（QPO）に関連するねじれ振動モードのスペクトルをどのように変化させるかを調査する。

手法
著者らは、静的かつ球対称なDANSをモデル化するために、二流体形式を採用している。この系は、二つの完全流体（バリオン物質とダークマター）で構成され、それらは重力を通じてのみ相互作用する。

• 状態方程式（EoS）： バリオン成分には、核（コア）と地殻を整合的に扱うために、統一NL3$\omega\rho$L55 EoSを用いている。ダークマターは、粒子質量 $m_\chi$ と無次元の自己相互作用強度 $y$ によって特徴付けられる、反発的な自己相互作用を持つフェルミオン流体としてモデル化されている。
• 数値構築： 著者らは、結合されたトルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を解き、平衡配置を生成する。総重力質量（$M_t = 1.4 M_\odot$ および $2.0 M_\odot$）および固定されたダークマター質量分率（$f_\chi = 0.05$）のシーケンスに焦点を当てている。
• 振動解析： ねじれ地殻振動は、相対論的コーリング近似（計量摂動を無視）を用いて導出される。変位関数のマスター波方程式は、地殻-核界面および星表面におけるゼロトルク境界条件を持つ固有値問題として解かれる。
• 剪断弾性率： 地殻の剪断弾性率は、電子スクリーニング効果（弾性率を軟化させる）を含めた場合と含まない場合の両方について、標準的な格子モデルを用いて計算される。

主要な貢献と結果

1. 地殻の厚さの減少：
   ダークマターの存在は、同等の質量を持つ純粋なバリオン中性子星（PNS）と比較して、バリオン地殻の厚さ（$\Delta R$）を減少させる。

• DM分布への依存性： ダークママーがバリオン表面を超えて広がる拡散したハローを形成する場合、減少は無視できる。しかし、ダークマターがバリオン半径内に集中している場合（「ダークコア」を形成する場合）、地殻は著しく薄くなる。
• 規模： $M_t = 1.4 M_\odot$ の場合、重いダークマター粒子（$m_\chi \approx 1$ GeV）のダークコア領域において、地殻の厚さは最大で約12%減少する。$M_t = 2.0 M_\odot$ の場合、減少率は約16%に達する。
• 解析的近似： 著者らは、バリオンのコンパクトネスに基づく地殻の厚さの近似的な解析公式を導出している。これらの公式は、最適なシナリオにおいて数値結果と1%未満の精度で一致する。
• 自己相互作用： 反発的な自己相互作用強度（$y$）を増加させると、ダークマターの分布が外側へと押し出され、星はダークコアからハロー構成へと遷移し、地殻の厚さがPNSの値へと回復する。

2. ねじれ振動周波数：
   地殻の薄層化は、ねじれ振動モード（基本モード $2t_0$ および第一高次モード $2t_1$）の周波数の上昇をもたらす。

• 周波数シフト： DANSの周波数は、比較可能なPNSよりも一般に高い。このシフトは、地殻が最も薄くなるパラメータ空間（コンパクト・ハローおよびダークコア領域）で最も顕著になる。
• 感度： 第一高次モード（$2t_1$）は、基本モードよりもダークマターの存在に対して敏感である。$M_t = 1.4 M_\odot$ の場合、$m_\chi$ が 1 GeV に近づくにつれて、周波数比は基本モードで約6.4%、第一高次モードで約16.4%増加する。

3. 微視物理学との縮退：
   本研究は、ダークマターによる周波数シフトと、地殻の微視物理学（特に電子スクリーニング）の不確実性との間の縮退に対処している。

• 基本モード： $M_t = 1.4 M_\odot$ の場合、ダークマターと電子スクリーニングによる周波数シフトは、全パラメータ空間にわたって縮退している。$M_t = 2.0 M_\odot$ の場合、縮退はダークコア領域（$m_\chi \gtrsim 400$ MeV）においてのみ解かれる。
• 第一高次モード： 第一高次モードの方が、より容易に縮退を解くことができる。ダークマターの分布が十分にコンパクト（コンパクト・ハローまたはダークコア）であれば、両方の質量シーケンスにおいて、$m_\chi \gtrsim 200\text{--}250$ MeV でDANSとPNSのバンドは分離する。

意義と主張
本論文は、ねじれ振動が、大質量中性子星内におけるダークマターコアの検出のための実行可能な地震学的プローブを提供することを主張している。著者らは、もしねじれモードに関連するQPOが測定されれば、その特定の周波数シフトによって、標準的な中性子星とダークマターコアを保持する中性子星を区別できる可能性があると断言している。これは、特に大質量星（$M_t \approx 2 M_\odot$）および重いダークマター粒子において顕著である。

本研究は、地殻がダークマターの空間分布に対する敏感な指標となることを強調している。著者らは、単一のEoSの使用、コーリング近似、および磁場や回転の除外といった限界を認めているが、これらの結果が、より完全なDANSのアステロセイズモロジー（星震学）に向けた基礎的な一歩を提供すると示唆している。彼らは、制約を精緻化するために、今後の研究において追加の微視物理学（例：核パスタ相）や完全な計量摂動を取り入れるべきであると提案している。