[论文解读] 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

一句话总结

该论文提出了一种基于随机批次法（RBM）的泊松-玻尔兹曼方程直接模拟方法，将N粒子系统的计算成本从O(N²)降低至每时间步O(N)。通过消除对硬球势能的需求，并实现在高维空间和复杂几何结构中的高效模拟，该方法能够实现带电粒子在电解质溶液中平衡态解的精确且可扩展的模拟。

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）来近似粒子间相互作用，将成对相互作用的计算成本从O(N²)降低至每时间步O(N)。
• RBM在每个时间步随机将粒子分组为小批次，从而实现对长程库仑力的高效计算。
• 该方法无需人工截断或硬球势能即可自然处理长程相互作用。
• 随着系统向平衡态演化，泊松-玻尔兹曼方程通过隐式方式求解。
• 分析了有限区域截断带来的截断误差，并通过数值实验验证了该方法的准确性和效率。

实验结果

研究问题

• RQ1 随机批次法能否有效应用于泊松-玻尔兹曼系统，在保持精度的同时降低计算成本？
• RQ2 该方法是否能消除在带电粒子N体模拟中对硬球势能的需求？
• RQ3 在具有复杂分子或膜表面的高维和非对称系统中，该方法表现如何？
• RQ4 区域截断对平衡态解的精度有何影响？
• RQ5 该方法能否通过直接粒子动力学模拟可靠地收敛至泊松-玻尔兹曼解？

主要发现

• 通过随机批次法，每时间步的计算成本从O(N²)降低至O(N)，实现了大规模N系统模拟的可扩展性。
• 该方法在无需硬球势能的情况下，实现了泊松-玻尔兹曼方程的精确平衡态解。
• 该方法在高维和非对称几何结构中表现稳健，包括具有复杂分子或膜表面的系统。
• 有限区域中的截断误差可控且具有解析上界，支持可靠的数值结果。
• 数值实验证实了该方法的有效性及其向正确平衡电势分布的收敛性。
• 通过直接模拟粒子动力学，自然地得到了平衡态下的泊松-玻尔兹曼解，同时保持了物理动力学的完整性。