[論文レビュー] Small-N Collisional Dynamics V: Beyond the Realm of Not-So-Small-N

Carlos Barrera, Nathan W. C. Leigh+3 more|arXiv (Cornell University)|Oct 30, 2020

Stellar, planetary, and galactic studies参考文献 2被引用数 4

ひとこと要約

本研究は、N ≈ 10 から N ≈ 1000 の粒子数の範囲で、直接 N-体シミュレーションと解析的モデルを比較することにより、密度の高い星団における星の衝突率を予測するための平均自由行程近似の妥当性を検証する。モデルは Leigh ら (2017) に基づくが、シミュレーションと通常1標準偏差以内で良好な一致を示し、類似した粒子種と中程度の N に対してはその頑健性が確認されるが、リラクゼーション時間の影響や考慮されていない団子の進化過程のため、N が高くなると精度が低下する。

ABSTRACT

Direct collisions between finite-sized particles occur commonly in many areas of astrophysics. Such collisions are typically mediated by chaotic, bound gravitational interactions involving small numbers of particles. An important application is stellar collisions, which occur commonly in dense star clusters, and their relevance for the formation of various types of stellar exotica. In this paper, we return to our study of the collision rates and probabilities during small-number chaotic gravitational interactions ($N$ $\lesssim$ 10), moving beyond the small-number particle limit and into the realm of larger particle numbers ($N$ $\gtrsim$ 10$^3$) to test the extent of validity of our analytic model as a function of the particle properties and the number of interacting particles. This is done using direct $N$-body simulations of stellar collisions in dense star clusters, by varying the relative numbers of particles with different particle masses and radii. We compute the predicted rate of collisions using the mean free path approximation, adopting the point-particle limit and using the sticky-star approximation as our collision criterion. We evaluate its efficacy in the regime where gravitational-focusing is important by comparing the theoretical rates to numerical simulations. Using the tools developed in previous papers in this series, in particular Collision Rate Diagrams, we illustrate that our predicted and simulated rates are in excellent agreement, typically consistent with each other to within one standard deviation.

研究の動機と目的

• 小N領域 (N ≤ 10) を超えて大N領域 (N ≳ 10^3) における、密度の高い星団における星の衝突率に対する平均自由行程近似の妥当性を検証すること。
• Leigh ら (2017) の解析的モデルが、異なる粒子質量および半径分布における衝突率をどの程度正しく予測できるかを評価すること。
• 小N解析的モデルの仮定が、リラクゼーション時間の影響および団子の進化に起因して破綻し始める臨界粒子数を特定すること。
• 重力的焦点効果が顕著な領域において、スティッキーサーモン近似および平均自由行程形式の有効性を評価すること。
• モンテカルロ的星団進化シミュレーションにおいて、計算コストの高い N-体積分器に代わる解析的衝突確率モデルを提供すること。

提案手法

• 初期条件にプラマーパターンを用い、粒子質量および半径を変化させた N ≈ 1000 粒子の星系に対する直接 N-体シミュレーションを実施する。
• 衝突基準としてスティッキーサーモン近似を適用し、接触した粒子は完全に融合するとみなす。
• 理論的衝突率を平均自由行程近似を用いて計算する：R_coll ≈ nσ_collv で、σ_coll は重力的焦点効果に基づいて導出する。
• 衝突率図を用いて、異なる粒子種および空間領域における予測値とシミュレーション値の比較を可視化する。
• Leigh ら (2017) の解析的予測とシミュレーションによる衝突回数を比較し、1標準偏差内での一致を評価する。
• 質量比、サイズ比、総粒子数を変化させることで、さまざまな天体物理学的条件下でのモデルの頑健性を調査する。

実験結果

リサーチクエスチョン

• RQ1 平均自由行程近似は、大N星系 (N ≳ 1000) において直接 N-体シミュレーションと比較して、どの程度衝突率を正しく予測できるか？
• RQ2 小N解析的モデル (Leigh ら, 2017) がリラクゼーション時間の影響および団子の進化に起因して破綻し始める粒子数は何か？
• RQ3 粒子種間の質量比および半径比が、解析的衝突率モデルの精度にどのように影響するか？
• RQ4 重力的焦点効果は、密度の高い星団において解析的予測とシミュレーションの一致度にどの程度影響を与えるか？
• RQ5 解析的モデルは、モンテカルロ的星団進化コードにおける小数の重力積分器に信頼性を持って置き換えられるか？

主な発見

• Leigh ら (2017) の解析的モデルは、N-体シミュレーションと良好な一致を示し、予測値とシミュレーション値は通常1標準偏差以内に収束する。
• 等質量粒子の場合、コア領域の衝突に限定してはモデルがシミュレーション値の約5倍のオーダーで予測値を出し、全団子体積にわたる統合を行うとこのずれは是正される。
• 同種粒子（質量比およびサイズ比が小さい）に対してモデルは最も良好に機能する。これは、質量分離やシミュレーションで観測される団子の進化効果をモデルが考慮していないためである。
• N が高くなる（≳1000）と性能が低下するが、これはリラクゼーション時間が衝突時間よりも長くなるためであり、小Nモデルにおける時間平均仮定が成立しなくなるためと推定される。
• 個々の少数体相互作用（N = 3–6）の間における衝突確率を予測するという点で、モデルは有効であるため、モンテカルロ的星団進化コードにおける N-体積分器の置き換えに有望である。
• 今後の研究では、小Nからシェルベースのモンテカルロモデリングに移行する必要が生じる臨界 N を特定する予定である。