[論文レビュー] A direct simulation approach for the Poisson-Boltzmann equation using the Random Batch Method

Lei Li, Jianguo Liu+1 more|arXiv (Cornell University)|Apr 12, 2020

Electrostatics and Colloid Interactions参考文献 40被引用数 9

ひとこと要約

本稿では、O(N²)の計算コストをO(N)に削減するため、ランダムバッチ法（RBM）を用いたポアソン＝ボルツマン方程式の直接シミュレーション手法を提案する。硬い球体ポテンシャルの必要性を排除し、高次元および複雑な幾何形状でも効率的なシミュレーションを可能にすることで、イオン溶液中の荷電粒子に対する正確でスケーラブルな平衡解の取得が可能になる。

ABSTRACT

The Poisson-Boltzmann equation is a nonlinear elliptic equation that describes how electronic potential changes in the screening layer when a cell is immersed in an ionic solution. It is the equilibrium state formed by many charged particles interacting with each other through Coulomb forces. We propose a direct simulation approach for the dynamics of the charged particles in the solution governed by the Poisson-Nernst-Planck equation and the Poisson-Boltzmann equation can be solved automatically when the particles reach the equilibrium. The interacting $N$ particle system is simulated in the truncated external domain. Directly simulating the $N$ particle system has an $\mathcal{O}(N^2)$ computational cost in each time step. By making use of the Random Batch Method (RBM) proposed by Jin et al [J. Comput. Phys., 400(1): 108877, 2020], the cost per time step is reduced remarkably to $\mathcal{O}(N)$. Another benefit of RBM is that the hard sphere potential is no longer needed in the simulation. An error analysis is provided for the truncation and numerical experiments are performed to validate the method. This particle method is preferable in two aspects: it is simple and effective in high dimensions, especially in non-symmetric cases or with complicated membrane or molecular surface. Meanwhile, direct simulation may be preferred as people may be interested in the dynamics of the physical process.

研究の動機と目的

• 荷電粒子の粒子ダイナミクスを捉え、自然に平衡状態に到達するポアソン＝ナーンスト＝プランクおよびポアソン＝ボルツマン方程式の直接シミュレーション手法の開発。
• 荷電粒子のN体シミュレーションにおいて、従来のO(N²)の高い計算コストを克服すること。
• 粒子の崩壊を防ぐためにしばしば用いられる硬い球体ポテンシャルの必要性を排除すること。
• 複雑な分子的または膜表面を有するような、高次元および非対称な系においても、効率的かつ正確なシミュレーションを可能にすること。
• 非対称かつ幾何学的に複雑な環境においても効果を発揮する、数値的に安定でスケーラブルなアプローチを提供すること。

提案手法

• 有限外部領域に制限されたN粒子系を、ポアソン＝ナーンスト＝プランク方程式に従って直接シミュレートする。
• ランダムバッチ法（RBM）を用いて粒子間相互作用を近似し、1ステップあたりのペアワイズ相互作用コストをO(N²)からO(N)に削減する。
• RBMは各タイムステップで粒子を小さなバッチに確率的にグループ化し、長距離クーロン力の効率的計算を可能にする。
• 人工的なカットオフや硬い球体ポテンシャルを必要とせず、自然に長距離相互作用を扱える。
• 系が平衡状態に近づくに従って、ポアソン＝ボルツマン方程式を暗黙的に解く。
• 有限領域による切断誤差を解析し、数値実験により本手法の精度と効率性を検証する。

実験結果

リサーチクエスチョン

• RQ1 ランダムバッチ法は、計算コストを削減しつつも精度を保持する形で、ポアソン＝ボルツマン系に効果的に適用可能か？
• RQ2 本手法は、荷電粒子のN体シミュレーションにおいて、硬い球体ポテンシャルの必要性を排除できるか？
• RQ3 本手法は、高次元および非対称な系、特に複雑な分子的または膜表面を有する系において、どのように性能を発揮するか？
• RQ4 領域の切断が、平衡解の精度に及ぼす影響は何か？
• RQ5 直接粒子ダイナミクスのシミュレーションにより、本手法は信頼性を持ってポアソン＝ボルツマン解に収束するか？

主な発見

• ランダムバッチ法を用いることで、1ステップあたりの計算コストをO(N²)からO(N)に削減し、大規模N系のスケーラブルなシミュレーションが可能になった。
• 硬い球体ポテンシャルを必要とせず、ポアソン＝ボルツマン方程式の正確な平衡解が得られた。
• 本手法は、高次元および非対称な幾何形状、特に複雑な分子的または膜表面を有する系に対しても、頑健である。
• 有限領域による切断誤差は制御可能であり、解析的に境界が与えられており、信頼性の高い数値結果を支持する。
• 数値実験により、本手法の有効性および正しい平衡ポテンシャル分布への収束が確認された。
• 粒子ダイナミクスの直接シミュレーションにより、自然に平衡状態でポアソン＝ボルツマン解が得られ、物理的ダイナミクスが保持された。