Discrepancies between Chandrasekhar's theory of relaxation and $N$-body simulations

問題の記述
二体緩和は、特に球状星団において、恒星力学における基本的なプロセスである。チャンドラセカールによる軌道平均理論（1941年）は、長らく標準的なモデルとして用いられてきた。先行研究（Aarseth et al. 1974; Theuns 1996; Tep et al. 2022、以下T22）は、概して理論とN体シミュレーションの一致を見出してきたが、依然として不一致が存在する。具体的には、T22は、等方的なプラマーモデルにおいて、チャンドラセカールの理論が分布関数の変化率（$\partial F/\partial t$）を約1.4倍過大評価していることを報告しており、この不一致は半質量半径の内部で最も顕著であった。さらに、Theuns (1996) による初期の研究は、特定の粒子エネルギーにおいて理論が拡散を過小評価している可能性を示唆しており、これらを整合させる必要があった。また、異方性を持つクラスターにおける緩和率の振幅に関する理論の予測精度についても、未解決の課題として残っている。

手法
著者らは、更新された手法とデータを用いて、チャンドラセカール理論とN体シミュレーションの比較を再検討している。

1. グローバルなコア崩壊の比較: 本研究は、現代的なN体データ（Pavlík & Šubr 2018による $N=10^4$ および $5\times10^4$）と、非確率論的な有限要素法によるフォッカー・プランク方程式の解（Takahashi 1995）を用いることで、Aarseth et al. (1974) によるグローバルな比較を更新している。著者らは、従来の導出では無視されがちであった非主要項を含む、Hénon (1975) に由来する特定のクーロン対数引数（$\Lambda = 0.15N$）を採用している。
2. 局所的な緩和率の測定: $\partial F/\partial t$ における具体的な不一致（T22が報告したもの）に対処するため、著者らは様々な異方性パラメータ（$q$）を持つ異方的なプラマー球（Dejonghe 1987）のファミリーを解析している。
   • 固定された時間間隔 $T$ に対する単純な有限差分法を用いる代わりに（これは有限$N$のゆらぎやビン分割の選択による大きな誤差を導入する）、作用空間における分布関数（DF）の初期時間進化に対する多項式フィッティング手法を採用している。
   • 最初の100〜1000ヘノン単位（HU）におけるDFの進化に対してフィッティングを行うことで、統計的ノイズやパラメータ選択の影響を最小限に抑えた、より堅牢な $\partial F/\partial t$ の推定値を導出している。
   • また、ゆらぎの影響を受けにくい二次的な力学的量、すなわち重力ポテンシャルの時間微分（$\partial \psi/\partial t$）およびコア半径の時間微分（$\dot{R}_c$）についても検証を行っている。

主な貢献と結果

• グローバルな不一致の解消: 現代的なN体コア崩壊時間と有限要素法によるフォッカー・プランク解、および修正されたクーロン対数を組み合わせることで、著者らは、等方的なプラマークラスターにおいてチャンドラセカールの理論とN体シミュレーションが数パーセント以内の範囲で一致することを見出した（$t_{cc, NB} / t_{cc, FP} \approx 1.02$）。これは、先行研究で見られた1.5倍の不一致とは対照的である。
• 局所的な振幅のミスマッチの低減: 多項式フィッティングを用いた $\partial F/\partial t$ の精緻な測定により、以前に報告されていたミスマッチは大幅に減少した。等方的な場合（$q=0$）において、理論的な率とN体の率のグローバルな比率は約1.0へと低下し、T22が主張した1.4倍の過大評価は事実上解消された。
• 異方性と位置への依存性: 等方的なケースにおける一致の改善にもかかわらず、残留するミスマッチが存在し、それは二つの要因に依存している：
  1. 異方性: 初期的な接線方向の異方性が強まる（$q$ がより負になる）につれて、理論的な率とN体の率の比が増大する（例：$q=-30$ で約2.2に達する）。著者らは、これが異方的な系におけるクーロン対数の減少によって定性的に説明できる可能性があると考えているが、これだけでは不一致を完全には解決できていない。
  2. 半径方向の位置: ミスマッチはクラスターの中心領域で最も顕著である。これは、中心部では大きなインパクトパラメータが抑制されるためクーロン対数が小さくなるべきであるという示唆と定性的に一致しているが、理論はこれを定量的に説明できていない。
• 諸量の整合性: コア半径や重力ポテンシャルの進化を解析した際にも、接線方向の異方性の増大に伴ってミスマッチが増大するという同様の傾向が見られた。これは、不一致がDFの測定手法のみによるアーティファクトではないことを裏付けている。
• 等方化のタイムスケール: 著者らは、異方的なクラスターが（内側領域において）最初の半質量緩和時間内に急速に等方化することを指摘している。したがって、初期の異方性に起因するミスマッチは、理論が予測する長期的なグローバルな進化に対して限定的な影響しか及ぼさない可能性がある。

意義と主張
本論文は、N体の測定、特に時間微分の抽出における多項式フィッティングの使用を通じた注意深い再評価により、等方的なクラスターにおけるチャンドラセカール理論とシミュレーションとの間の長年の振幅のミスマッチが大部分解決されると主張している。

しかし、著者らは、異方的なクラスターおよびすべてのクラスターの中心部においては、残留するミスマッチが残っていると控えめに結論づけている。彼らは、異方性への依存性は（クーロン対数を通じて）定性的に理解可能であるが、半径方向の依存性は、標準的な軌道平均チャンドラセカール理論の範囲を超えた、空間的な不均一性と集団効果を考慮することが、不一致を完全に解消するために不可避であることを示唆していると主張している。本論文は、チャンドラセカールに代わる新しい理論を提案するものではなく、今後の研究において、これらの残された差異を体系的に説明するために、不均質なランダウ理論の拡張や集団効果（バレスク・レナード方程式による）を含めるべきであることを示唆している。