Breakdown of hydrodynamic universality in causal first-order viscous neutron star oscillations

技術的要約：因果的一次式粘性中性子星振動における流体力学的普遍性の崩壊

問題提起
粘性中性子星の流体力学的モデリングには、特定の相対論的流体力学の理論を選択する必要がある。異なる定式化（例：異なる流体力学的フレーム）は、長波長かつ近平衡状態の極限において同一の物理的予測を与えることが期待されるが、これは「流体力学的普遍性（hydrodynamic universality）」と呼ばれる特性である。しかし、近年の研究は、十分に粘性の高い中性子星においては、この普遍性が失われる可能性を示唆している。具体的には、因果律の制約や有限サイズ効果によって、系が普遍的な領域に到達できなくなる場合がある。Denicolらによる先行研究では、径方向振動におけるイスラエル・スタス（IS）理論と非因果的なナビエ・ストークス（NS）理論の間で、この崩壊が示された。しかし、この現象が二次の理論に特有のものなのか、それともベムフィカ・ディスコンジ・ノロナ・コブトゥン（BDNK）フレームワークのような広範な因果的一次式理論にまで及ぶのかは不明であった。本論文では、現実的な中性子星の条件下において、BDNK理論における広範な因果的一次式フレームにおいて、流体力学的普遍性が崩壊するかどうかを調査する。

手法
著者らは、BDNK理論における中性子星の径方向振動を研究するために、摂動論的枠組みを用いている。

• 背景および構成関係： 本研究では、ゼロの化学ポテンシャルを持つバロトロピックなBDNK流体に結合した、球対称時空を利用している。構成関係には、体積粘性（$\zeta$）、剪断粘性（$\eta$）、および緩和時間（$\tau_{e,p,q}$）が含まれる。背景となる星は、多項式状態方程式を用いたトルマン・オッペンハイマー・ヴォルコフ（TOV）方程式を用いてモデル化されている。
• 数値的手法： 径方向の摂動方程式を周波数領域で解いている。著者らは、第二次の精度を持つ中心差分法と、非物理的なグリッドモードを回避するための第一のアインシュタイン方程式に対する空間的なクランク・ニコルソン法を用いて、系を離散化している。これにより、二次固有値問題 $Q(\mu)U = 0$ （ここで $\mu = \omega/\omega_{PF}$ は完全流体の固有値で正規化された周波数）が得られる。
• ソルバー： 固有値系を解くために、ニュートン反復法（粘性によるモードの追跡用）とQZアルゴリズム（非流体力学的モード用）の2つの手法を用いている。
• 因果性とパラメータ空間： 解析は、因果的な伝播のための十分条件によって定義されるパラメータ空間の「因果的領域」に焦点を当てている。著者らは、緩和時間と粘性をパラメータ化し、3次元のパラメータ空間 $(\tau_e, \tau_p, \tau_q)$ における「レイ（光線）」を定義している。特に、原点（$\tau \to 0$）に最も近いフレームを表す「最小因果レイ」を分析している。
• 診断指標： 普遍性の崩壊を定量化するために、以下の分数偏差を定義している：
  • $\Delta_{EFT} = |(\omega_{BDDK} - \omega_{NS}) / (\omega_{NS} - \omega_{PF})|$：粘性によるシフトに対する、因果的なBDDK周波数と非因果的なナビエ・ストークス（NS）周波数の偏差を測定する。$\Delta_{EFT} \ll 1$ は、系が普遍的な有効場理論（EFT）領域にあることを示す。
  • $\Delta_{NS} = |(\omega_{BDDK} - \omega_{NS}) / \omega_{NS}|$：NSからBDDKへのシフトを測定する。
  • 普遍性の崩壊は、すべての因果的フレームにおいて $\Delta_{EFT} \geq 0.1$ となったときに診断される。

主な貢献および結果

1. 一次式理論における普遍性の崩壊： 著者らは、BDDK理論における広範な因果的一次式フレームにおいて、流体力学的普遍性が崩壊することを実証した。体積粘性が増加するにつれ、因果的パラメータ空間の領域は、原点（EFT領域が存在する場所）から、より大きな緩和時間の方へと押しやられる。
2. 臨界粘性閾値： 崩壊は、体積粘性が $\zeta \sim 10^{30}–10^{31} \text{ g/cm/s}$ の範囲で発生する。これらの値において、すべての因果的フレームは一次式のEFT領域の外側に位置する（$\zeta_c = 10^{31} \text{ g/cm/s}$ のとき $\Delta_{EFT} \gtrsim 0.85$）。この範囲は、中性子星の体積粘性の従来の見積もり、およびGW170817から導かれた上限値と一致している。
3. 振動モードの構造変化： 普遍性の崩壊は、振動モード構造の著しい変化を伴う：
   • 非流体力学的モードが、流体力学的モードと同等のタイムスケールで減衰する。
   • 第一高調波が基本モードよりも長寿命になることがある。
   • 非物理的な不安定性： 十分に高い粘性（$\zeta_c = 10^{31} \text{ g/cm/s}$）の場合、平坦な時空において理論が線形安定であるにもかかわらず、サブミリ秒のタイムスケールを持つ非物理的な不安定性が現れる。
4. 重力崩壊閾値の不変性： 普遍性の崩壊と非物理的なモードの出現にもかかわらず、著者らは、BDDK理論において重力崩壊の無粘性閾値が変化しないことを見出した。周波数領域で境界条件が適切に処理されている限り、不安定性の開始を示すゼロ周波数モードは、完全流体の場合と同じ中心密度で存続する。
5. フレーム依存性： 本研究は、EFT領域が失われる一方で、粘性が増大するにつれて、物理的予測（固有値）が因果的なフレームの選択に大きく依存し得ることを強調している。

意義および主張
本論文は、BDDKを含む因果的一次式相対論的流体力学の理論が、十分に粘性の高い中性子星のモデリングには不十分である可能性があるという証拠を提示している。流体力学的普遍性の崩壊は、異なる因果的定式化が物理的予測について一致しなくなることを意味し、理論が非物理的な不安定性を許容する領域に入る可能性があることを示唆している。

著者らは、BDDK理論は因果的かつ安定するように設計されているが、それらの性質はあくまでその妥当な領域（EFT領域）内でのみ保証されるものであると強調している。粘性が十分に高く、システムがこの領域の外へ押し出されると、理論の予測能力、特にモード構造や安定性に関する能力が損なわれる。これらの結果は、高い体積粘性を伴う天体物理学的シナリオ（連星中性子星のインスパイラルや合体後のフェーズに関連する可能性がある）において、二次式理論または平衡から遠い記述が必要になる可能性を示唆している。著者らは、これらの結果が自然界においてこれらの非物理的な不安定性が実際に現れることを裏付けるものではなく、一次式因果的流体力学のこのような極端な条件下における適用限界を示すものであることを明示している。