[論文レビュー] Neighbor List Collision-Driven Molecular Dynamics Simulation for Nonspherical Particles. I. Algorithmic Details II. Applications to Ellipses and Ellipsoids

Aleksandar Donev, Salvatore Torquato+1 more|ArXiv.org|May 17, 2004

Scientific Research and Discoveries参考文献 42被引用数 53

ひとこと要約

本稿では、非球形粒子のための衝突駆動型分子動力学アルゴリズムを提示する。新しい部分更新近隣リスト（NNL）とバウンディングスフィアコンプレックス（BSC）を用いて、硬い楕円形および楕円体粒子の効率的かつ高密度なシミュレーションを実現する。この手法により、ジャムドパッキング、力の鎖、結晶の融解といった現象を、従来の手法よりも計算効率を向上させつつ正確に再現可能となる。

ABSTRACT

In the first part of a series of two papers, we present in considerable detail a collision-driven molecular dynamics algorithm for a system of nonspherical particles, within a parallelepiped simulation domain, under both periodic or hard-wall boundary conditions. The algorithm extends previous event-driven molecular dynamics algorithms for spheres. We present a novel partial-update near-neighbor list (NNL) algorithm that is superior to previous algorithms at high densities, without compromising the correctness of the algorithm. This efficiency of the algorithm is further increased for systems of very aspherical particles by using bounding sphere complexes (BSC). In the second part of this series of papers we apply the algorithm presented in the first part of this series of papers to systems of hard ellipses and ellipsoids. The theoretical machinery needed to treat such particles, including the overlap potentials, is developed in full detail. We describe an algorithm for predicting the time of collision for two moving ellipses or ellipsoids. We present performance results for our implementation of the algorithm. The practical utility of the algorithm is demonstrated by presenting several interesting physical applications, including the generation of jammed packings inside spherical containers, the study of contact force chains in jammed packings, and melting the densest-known equilibrium crystals of prolate spheroids.

研究の動機と目的

• 楕円形や楕円体のような非球形粒子のための、堅牢かつ効率的な分子動力学アルゴリズムの開発。
• 非球形粒子の高密度なシミュレーションにおける計算ボトルネックを、最適化された近隣リスト管理によって克服すること。
• 異方的粒子系におけるジャミング、力の鎖、結晶の融解といった物理現象の正確なシミュレーションを可能にすること。
• 球形粒子にとどまらないイベント駆動型MDを、楕円体形状の幾何的重なり検出を組み込むことで拡張すること。
• 高密度な非球形粒子パッキングを研究するためのスケーラブルで正確かつ高性能なフレームワークを提供すること。

提案手法

• 計算コストを低減しつつアルゴリズムの正しさを保持するため、新しい部分更新近隣リスト（NNL）アルゴリズムを導入。
• 高非球形粒子の衝突検出を高速化するために、入れ子の球体で形状を近似するバウンディングスフィアコンプレックス（BSC）を用いる。
• 幾何的重なり条件の解析的解を用いて、移動中の楕円形および楕円体粒子間の衝突時刻を予測する。
• 平行六面体シミュレーション領域内で、周期的境界条件およびハードウォール境界条件の両方をサポートする。
• 効率性を向上させるために、衝突イベントでのみシステムを進める衝突駆動型統合スキームを採用。
• 硬い楕円形および楕円体粒子の重なりポテンシャルと幾何的制約の理論的基盤を厳密に構築。

実験結果

リサーチクエスチョン

• RQ1 イベント駆動型分子動力学を、楕円形や楕円体のような非球形・異方的粒子に効率的に拡張するにはどうすればよいか？
• RQ2 非球形粒子の高密度系において、近隣リストを最適に維持する戦略は何か？
• RQ3 バウンディングスフィアコンプレックス（BSC）は、非常に細長い粒子の衝突検出においてどのように性能を向上させるか？
• RQ4 このアルゴリズムを用いて、ジャミングや結晶の融解といった物理現象をどれほど信頼性高くシミュレートできるか？
• RQ5 部分更新NNLアルゴリズムは、計算コストを削減しつつも、正確性をどの程度維持できるか？

主な発見

• 部分更新NNLアルゴリズムは、従来の手法と比較して高密度領域において著しく性能向上を示し、正しさを損なわない。
• バウンディングスフィアコンプレックス（BSC）の導入により、高非球形粒子のシミュレーションにおいて詳細な幾何的チェックの回数を削減し、効率性が向上した。
• 球状の容器内に楕円形および楕円体粒子のジャムドパッキングを効果的に生成でき、実用的価値を示した。
• ジャムドパッキング内の接触力の鎖を可視化・分析し、異方的力伝達パターンを明らかにした。
• プロライトスフェロイドの最も高い密度に達する平衡結晶を、成功裏に融解させ、本手法が相転移をシミュレートできる能力を検証した。
• パフォーマンスベンチマークにより、特に形状の異方性が強い系において、本アルゴリズムのスケーラビリティと効率性が確認された。