Nonradial oscillations of realistic anisotropic neutron stars: Polar modes

技術要約：現実的な異方性中性子星の非放射振動：極モード

問題提起
中性子星（NS）における非放射振動の研究は、アステロセイモロジーを通じて高密度物質の内部構造を探るための重要なツールとなっている。これまでの研究では等方的な恒星流体が広範に分析されてきたが、超強力な磁場、超流動、あるいはエキゾチックな相の状態といった物理過程は、圧力の異方性（接線方向の圧力が径方向の圧力と異なる状態）を誘発する可能性がある。近年の研究では平衡モデルに異方性を組み込んでいるものの、線形摂動下における異方性の整合的な取り扱いは依然として困難である。具体的には、標準的な異方性処方（例えばBowers-LiangモデルやHorvatモデル）は平衡状態における異方性のみを定義しており、異方性ベクトルが摂動に対してどのように応答するかを指定していない。これにより、摂動問題が未決定の状態となっている。先行研究では、多くの場合、異方性ベクトルのオイラー的摂動が$(\bar{u}, \bar{k})$平面内にある（すなわち$\theta$成分が消失する）と仮定しており、これは流体の変位から一つの力学的な自由度を実質的に除去している。本論文は、この未決定性を解決する、異方性中性子星の極摂動に対する一貫した完全相対論的取り扱いの必要性に対処するものである。

手法
著者らは、異方性圧力を持つ静的かつ球対称な中性子星を研究するために、完全な一般相対論的枠組みを採用している。背景時空は、異方的な流体エネルギー・運動量テンソルによって生成される線素によって記述される。本研究は、計量および流体の線形一次摂動に焦点を当てている。

中心的な手法上の革新は、径方向単位ベクトル $k^\alpha$ のラグランジュ的摂動の取り扱いである。$\delta k^\alpha$ が$(\bar{u}, \bar{k})$平面内にあると仮定する（$\delta k^\theta = 0$とする）のではなく、著者らは、ラグランジュ的摂動 $\Delta k^\alpha$ が$(\bar{u}, \bar{k})$平面内にあるという制約を課している。この選択は、異方性が物質に結びついており、流体と共に移流するという物理的仮定に基づいている。この制約は、追加の力学的自由度を導入し、横方向の変位成分 $V(r,t)$ を代数的に拘束されるものではなく、力学的なものにする。

系は、線形化されたアインシュタイン場の方程式と、摂動されたエネルギー・運動量保存則によって支配される。著者らは、2つの標準的な異方性処方を用いている：

1. Bowers-Liangモデル: $\bar{\sigma}_{BL} = -\lambda_{BL}(\bar{\epsilon} + 3\bar{P})(\bar{\epsilon} + \bar{P})r^2/B$
2. Horvatモデル: $\bar{\sigma}_{H} = -\lambda_{H}(1 - B)\bar{P}$

核子（SLy4）、ハイパーオン（GM1Y6）、およびハイブリッド（QHC21）物質を表す3つの現実的な状態方程式（EOS）が使用されている。数値実装には、極摂動に関する6つの結合一次微分方程式の積分、内部解と外部のZerilli方程式とのマッチング、およびニュートン・ラフソン法を用いた複素準固有振動数（$\omega = \omega_R + i\omega_I$）の算出が含まれる。

主な貢献

1. 一貫した摂動形式: 著者らは、異方性ベクトルのラグランジュ的摂動を定義することにより、前著における軸対称摂動の研究を極セクターへと拡張し、摂動方程式の整合的な閉鎖を与えた。これにより、異方性流体の剪断弾性率に関連する新しい力学的自由度が導入される。
2. 未決定系の解決: 正規化条件と直交条件から $k^\alpha$ の摂動を導出することで、著者らは、$\delta k^\theta = 0$ と恣意的に設定していた従来の研究における曖昧さを解決した。
3. 準普遍関係の導出: 著者らは、無次元のfモード振動成分（$M\omega_R$ および $M\omega_I$）と恒星のコンパクトネス（$C = M/R$）の間の準普遍関係を導出し、フィッティングを行った。これらの関係は異方性パラメータ $\lambda$ によってパラメータ化されており、特定のEOSにはほとんど依存しないことが示された。

結果

• 振動数と減衰の傾向: 3つのEOSおよび両方の異方性モデルにおいて、fモードの振動数（$\nu$）は恒星質量とともに増加し（1から3 kHzの範囲）、減衰時間（$\tau$）は減少する（0.5から1.25 sの範囲）。
• 異方性が振動数に与える影響: 異方性を増大させる（接線方向の圧力が径方向の圧力を上回る）と、一般に振動数は低下する。この効果はBowers-Liangモデルにおいてより顕著である。
• 異方性が減衰時間に与える影響: 減衰時間への影響は、異方性の処方に依存する。Horvatモデルでは、異方性の増大は減衰時間の減少をもたらす。対照的に、Bowers-Liangモデルでは、異方性の増大は減衰時間の増加をもたらす。
• 先行研究との比較: $\delta k^\theta = 0$ を仮定する形式（例：Lauら）と比較して、著者らの手法は、小さな異方性パラメータにおいては同様の結果を与える。しかし、大きな異方性と高い恒星質量においては差異が顕著になり、新しい形式は、わずかに低い振動数と長い減衰時間を予測する。減衰時間は、大きな異方性パラメータに対して最大20%異なる。
• 多項式近似: 著者らは、$M\omega_R(C, \lambda)$ および $M\omega_I(C, \lambda)$ の関係に対する多項式近似を提供している。これらの近似は、数値データに対して10%以内の精度を達成しており、係数はEOSには依存しないが、異方性モデルには依存する。

意義
本論文は、異方性の整合的な摂動処理が、中性子星の振動を完全に特徴付けるために不可欠であると主張している。新しい力学的自由度の導入は、特に大きな異方性において、極準固有モードの挙動を大きく変化させる。導出された準普遍関係は、異方性パラメータによってパラメータ化されており、簡便な現象論的枠組みを提供する。著者らは、これらの関係が将来のアステロセイモロジー観測において、中性子星の圧力異方性を制約するために利用できる可能性があると示唆している。本研究は、摂動の処方の選択が単なる技術的な詳細ではなく、異方性星の重力波シグネチャーに対して観測可能な影響を与えることを確立している。