[論文レビュー] Generalized Moment Cancellation for Long-Range Electrostatics

Björn Stenqvist, Vidar Aspelin+1 more|arXiv (Cornell University)|Apr 23, 2019

Electrostatics and Colloid Interactions被引用数 2

ひとこと要約

本稿では、カットオフ球内での無限個の電気双極モーメントを厳密にキャンセルする形式的な手法を提案する。従来の手法で用いられる減衰パラメータの必要性が排除され、分子動力学シミュレーションにおいて水および塩水溶液の長距離電気双極モーメントに優れた精度と効率を示す。この手法は、電気双極モーメント間相互作用への応用も可能である。

ABSTRACT

Describing long-ranged electrostatics using short-ranged pair potentials is appealing since the computational complexity scales linearly with the number of particles. The foundation of this approach is to mimic the long-ranged medium response by cancelling electric multipoles within a small cutoff sphere. We propose a rigorous and formally exact new method that cancels up to infinitely many multipole moments and is free of operational damping parameters often required in existing theories. Using molecular dynamics simulations of water with and without added salt, we discuss radial distribution functions, Kirkwood-Buff integrals, dielectrics, diffusion coefficients, and angular correlations in relation to existing electrostatic models. We find that the proposed method is an efficient and accurate alternative for handling long-ranged electrostatics as compared to Ewald summation schemes. The methodology and proposed parameterization is applicable also for dipole-dipole interactions.

研究の動機と目的

• 有限のカットオフ球内での無限階の電気双極モーメントを正確にキャンセルする手法の開発。
• 既存の短距離電気双極モーメントモデルで一般的に用いられる操作的減衰パラメータの必要性を排除すること。
• 分子動力学シミュレーションにおける長距離電気双極モーメントのEwald和への効率的かつ正確な代替手段を提供すること。
• 水および塩水溶液における主要な熱力学的および動的性質の観点から、手法の性能を検証すること。
• 電気双極モーメント間相互作用への応用可能性を拡張すること。

提案手法

• 電気双極モーメントのすべての高次項までを球状カットオフ領域内で正確にキャンセルするための電気的ポテンシャルの形式的展開を採用する。
• 厳密な数学的枠組みを用いて、経験的減衰パラメータを導入せずに長距離寄与を正確にキャンセルすることを保証する。
• カットオフ球内での相互作用ポテンシャルを変更することで、長距離媒体応答を模倣する形で分子動力学シミュレーションに実装する。
• 相互作用を有限のカットオフ半径に制限することにより、系のサイズに比例する線形スケーリングの計算コストを維持する。
• 一貫した多重極キャンセレーションを通じて、電荷-電荷相互作用および電気双極モーメント-電気双極モーメント相互作用の両方を扱えるように形式を一般化する。
• 径方向分布関数、Kirkwood-Buff積分、誘電率、拡散係数、および角度相関関数を、基準となるEwald結果と比較することで、実装の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 短距離電気双極モーメントモデルが、減衰パラメータを一切用いずに、無限階までのすべての多重極モーメントを正確にキャンセル可能か？
• RQ2 水および塩水系において、提案手法がEwald和に比べて精度と効率でどのように優れているか？
• RQ3 誘電率や拡散係数といった主要な熱力学的および動的性質が、どの程度正確に保持されているか？
• RQ4 電気双極モーメント-電気双極モーメント相互作用を、同程度の精度と効率で扱えるように一般化可能か？
• RQ5 減衰パラメータの欠如が、異なる系にわたるモデルの転送性およびロバスト性を向上させるか？

主な発見

• 提案手法はカットオフ球内においてすべての多重極モーメントを正確にキャンセルし、減衰パラメータの必要性を排除した。
• Ewald和と比較して、径方向分布関数、Kirkwood-Buff積分、誘電率を高い精度で再現した。
• 拡散係数および角度相関関数の計算結果はEwald結果と良好に一致し、動的精度の妥当性を確認した。
• 計算コストは系のサイズに比例して線形に増加し、標準的な短距離モデルと同等の効率性を示した。
• 本手法は電気双極モーメント-電気双極モーメント相互作用に対しても直接適用可能であり、電荷-電荷相互作用を超えた応用性を有する。
• 数値的安定性と転送性に優れた、パラメータフリーのEwald和の代替手段を提供した。