[論文レビュー] Molecular Dynamics Simulations of Active Matter using LAMMPS

C. S. Dias|arXiv (Cornell University)|Feb 20, 2021

Gas Dynamics and Kinetic Theory被引用数 5

ひとこと要約

本論文は、LAMMPSを用いたアクティブマターのシミュレーションに関する実践的ガイドを提示し、カスタムfixを用いてアクティブブラウン運動粒子（ABP）を実装することに焦点を当てている。Langevinダイナミクスの統合、Verletリストによる力の計算、および粒子の向きに基づく持続的自己駆動力を持つ新しい'langevin/active' fixの拡張を詳細に説明しており、シミュレーションで球形から拡散的移動への遷移と効果的拡散係数を回復することで検証されている。

ABSTRACT

LAMMPS is a widely popular classical Molecular Dynamics package. It was designed for materials modeling but it is well prepared for simulations in Soft Matter. The use packages like LAMMPS has advantages and disadvantages. The main advantage is the optimization of the methods, mainly for parallel computing. The main disadvantage is that, due to the complexity of the code, it ha long learning curve. One purpose of these notes is to shorten that curve for researchers are starting to use LAMMPS to simulate soft active matter. In these notes, we first discuss some Molecular Dynamics methods implemented in LAMMPS. We present an hands-on introduction to first-time users and we finish with an advanced hands-on section, where we implement and test Active Brownian particles simulations.

研究の動機と目的

• ソフトアクティブマターのシミュレーションにおいてLAMMPSに初めて触れる研究者による学習曲線の低減。
• LAMMPSフレームワーク内でのアクティブブラウン運動粒子ダイナミクスの実用的で段階的な実装の提供。
• LAMMPSのモジュラー・アーキテクチャを正しく使用して、カスタム力モデルによるシミュレーション機能の拡張方法の示唆。
• 終端速度と平均二乗変位（MSD）を含む物理的ベンチマークを通じた実装の検証。

提案手法

• メモリ使用量の最小化と並列効率の向上を目的とした、半ステップ速度更新を伴う速度Verlet積分法の実装。
• 「スキン」距離を用いたVerletリストの使用により、近接リストの更新を最適化し、冗長な力計算を削減。
• LAMMPSにおけるLangevinダイナミクスの適応に「langevin」fixを用い、カスタム「langevin/active」fixによりアクティブ力成分を追加。
• クaternionを用いた粒子形状データを活用し、LAMMPSの「fix」システムを拡張して、方向依存の自己駆動力を持つ力を実装。
• クaternionと「angmom」オプションを介して、回転拡散と並進運動を結びつけることで、回転と並進のダイナミクスを統合。
• 還元単位、ゼロ温度での終端速度測定、および平均二乗変位（MSD）解析を用いた実装の検証。

実験結果

リサーチクエスチョン

• RQ1 LAMMPSを用いて自己駆動粒子を有するアクティブマター系を効果的にシミュレートする方法は何か？
• RQ2 クイーンタニオンに基づく粒子の向きを用いて、LAMMPSでカスタムfixとアクティブ力の正しく実装する方法は何か？
• RQ3 ゼロ温度下でのアクティブ粒子の終端速度は、印加力に対してどのようにスケーリングするか？
• RQ4 アクティブブラウン運動粒子における平均二乗変位の時間依存的挙動は何か？そして、球形から拡散的移動への遷移が観察されるか？
• RQ5 並進拡散と回転拡散を含む理論的式を用いて、アクティブ粒子の効果的拡散係数を回復できるか？

主な発見

• ゼロ温度下でのアクティブ粒子の終端速度は、理論で予測されるように、$ v = \frac{\tau_t}{m}F $ の関係に従い、印加力に線形に比例する。
• 平均二乗変位（MSD）曲線は、短時間では明確な球形領域、長時間では拡散的領域を示し、アクティブ粒子のダイナミクスが確認された。
• 効果的拡散係数 $ D_{\text{eff}} = D_t + \frac{1}{6}\frac{v^2}{D_r} $ がシミュレーションで回復され、理論モデルの妥当性が検証された。
• 回転拡散係数 $ D_r $ と並進係数 $ D_t $ の間には、$ \frac{D_t}{D_r} = \alpha \frac{I}{m} $ というスケーリング関係があり、ここで $ \alpha = \frac{\tau_t}{\tau_r} $ である。これは既知の流体力学的関係と整合的である。
• カスタムの「langevin/active」fixは、LAMMPSの既存フレームワークと正常に統合され、持続的自己駆動を持つアクティブな楕円体の安定的かつ正確なシミュレーションを可能にした。
• 小規模系での体系的テストを経て、大規模かつ並列化されたシミュレーションへのスケーリングが成功裏に検証された。