[論文レビュー] Tidal distortion and disruption of rubble-pile bodies revisited -- Soft-sphere discrete element analyses

Yun Zhang, Patrick Michel|arXiv (Cornell University)|Jun 24, 2020

Astro and Planetary Science参考文献 50被引用数 7

ひとこと要約

本研究では、pkdgrav N体コードを用いて、柔ららかく球形の離散要素モデル（SSDEM）を用い、地球との重力的接近を伴うルバーパイル小惑星のシミュレーションを実施し、衝突条件と内部構造に応じて、変形、質量の放出、破壊の3つの明確な結果が得られた。破壊は、従来の予測よりもはるかに近接した距離および低い速度で発生し、中程度の接近では極めて細長い形状（最大1:6のアスペクト比）が形成可能であることが判明した。また、内部構造と粒子の配置が、耐性および形状再構築のダイナミクスに強く影響することを確認した。

ABSTRACT

During close approaches to planets or stars, the morphological and dynamical properties of rubble-pile small bodies can be modified, and some may catastrophically break up. This phenomenon is of particular interest for understanding the evolution and population of small bodies, and for making predictions for future encounters. Previous numerical explorations typically used methods that do not adequately represent the nature of rubble piles. The encounter outcomes and influence factors are still poorly constrained. Based on recent advances in modeling rubble piles, we aim to provide a better understanding of the tidal encounter processes through soft-sphere discrete element modeling (SSDEM) and to establish a database of encounter outcomes and the dependencies on encounter conditions and rubble-pile properties. We performed thousands of numerical simulations using the SSDEM implemented in the N-body code pkdgrav to study the dynamical evolution of rubble piles during close encounters with the Earth. The effects of encounter conditions, material strength, arrangement, and resolution of constituent particles are explored. Three typical tidal encounter outcomes are classified, i.e., deformation, mass shedding, and disruption, ranging from mild modifications to severe damages of the progenitor. The encounter speed and distance required for causing disruption events are much smaller than those predicted by previous studies, indicating a smaller creation rate of tidally disrupted small body populations. Extremely elongated fragments with axis ratios ~1:6 can be formed by moderate encounters. Our analyses of the spin-shape evolution of the largest remnants reveal reshaping mechanisms of rubble piles in response to tidal forces, consistent with stable rubble-pile configurations derived by continuum theory. A case study for SL9 suggests a low bulk density (0.2-0.3 g/cc) for its progenitor.

研究の動機と目的

• 従来の連続体モデルや剛体粒子モデルとは異なり、より現実的な数値モデルを用いることで、ルバーパイル天体における潮汐破壊の理解を深めること。
• 接近パラメータ（近点距離、無限遠方での速度）および内部特性（材料強度、粒子の配置、解像度）が潮汐的結果に与える影響を調査すること。
• 物理的および動的条件に依存する、潮汐接近の結果とその依存関係を包括的に記録したデータベースを構築すること。
• 球形粒子モデルが、不規則なルバーパイル構造のダイナミクスをどれだけ的確に再現できるかを評価すること。
• シューケーマー＝リーヴィー9の母天体に関する制約を提示し、体積密度が0.2–0.3 g/cm³である可能性を示唆すること。

提案手法

• ルバーパイル天体内の球形粒子間の重力的および接触力の相互作用を、pkdgrav N体コードに組み込んだ柔ららかく球形の離散要素モデル（SSDEM）を用いてシミュレーションした。
• ヘキサゴナル密着.packing（HCP）およびランダム密着.packing（RCP）を用いた、異なる内部構造を持つルバーパイルをモデル化し、粒子数は247～1000の範囲で変化させた。
• 粒子の回転抵抗および摩擦接触力を組み込み、細粒凝集体の現実的な機械的挙動を再現した。
• 接近パラメータ（近点距離、無限遠方での速度、衝突パラメータ）を変化させ、広範な潮汐的状態を探索した。
• 最大の破片のスピン・形状の進化を追跡し、潮汐応力下での形状再構築メカニズムを分析した。
• 連続体理論（Holsapple, 2004, 2007）および従来の多面体ベースのシミュレーション（Movshovitz et al., 2012）と比較して妥当性を検証し、一貫性があることを確認した。

実験結果

リサーチクエスチョン

• RQ1 地球に近接する際の潮汐力が、どのようにルバーパイル小惑星に変形、質量の放出、または完全な破壊を引き起こすか？
• RQ2 破壊の閾値が、接近距離、速度、および粒子の配置や強度といった内部構造にどのように依存するか？
• RQ3 異なる内部構造（HCP 対 RCP）が、潮汐接近中の耐性および形態的進化にどれほど影響を与えるか？
• RQ4 多面体モデルと比較して、柔ららかく球形の離散要素モデルが、不規則な形状を持つルバーパイル天体のダイナミクスをどれほど的確に捉えられるか？
• RQ5 潮汐接近の結果から、シューケーマー＝リーヴィー9のような既知の天体の体積密度および内部特性にどのような制約を課せられるか？

主な発見

• 破壊は、従来の静的連続体モデルが予測したよりもはるかに近接した近点距離および低い接近速度で発生するため、潮汐的に破壊される小天体の集団形成レートが低くなる可能性がある。
• 変形、質量の放出、および破壊の3つの明確な潮汐的結果が特定され、その結果は接近条件および内部構造に強く依存していた。
• 中程度の潮汐接近では、特にRCP構造において、軸比が最大約1:6にまで細長くなる破片が形成可能である。
• RCPモデルでは材料強度および接近パラメータに高い感受性を示したが、HCPモデルは正確な内部構造と粒子数に強く依存し、細かい構造ほど破壊に対して耐性が高かった。
• 破片のスピン・形状進化の分析から、潮汐力が内部再配置を引き起こし、連続体理論で予測される安定な構造に一致する形状再構築が進行することが明らかになった。
• シューケーマー＝リーヴィー9の事例研究から、母天体の体積密度が0.2–0.3 g/cm³である可能性が示唆され、これはルバーパイル構造と整合的である。