[论文解读] Elastic Collision Based Dynamic Partitioning Scheme for Hybrid Simulations

Björn Kirchhoff, Elvar Ö. Jónsson+3 more|arXiv (Cornell University)|Mar 30, 2021

Theoretical and Computational Physics参考文献 64被引用 8

一句话总结

该论文提出SAFIRES，一种用于混合分子模拟的新型动态分区方案，通过可移动边界处的弹性碰撞来保持能量和动量守恒，从而实现精确的系综统计。该方法通过动态调整时间步长以检测碰撞，相较于静态或近似边界方法（如FIRES），可消除粒子穿插现象，保持热力学一致性，且与无约束模拟的偏差极小。

ABSTRACT

The scattering-adapted flexible inner region ensemble separator (SAFIRES) is a partitioning scheme designed to divide a simulation cell into two regions to be treated with different computational methodologies. SAFIRES prevents particles from crossing between regions and resolves boundary events through elastic collisions of the particles mediated by the boundary, conserving energy and momenta. A multiple-time-step propagation algorithm is introduced where the time step is scaled automatically to identify the moment a collision occurs. If the length of the time step is kept constant, the new propagator reduces to a regular algorithm for Langevin dynamics, and to the velocity Verlet algorithm for classical dynamics if the friction coefficient is set to zero. SAFIRES constitutes the exact limit of the premise behind boundary-based methods such as FIRES, BEST, and BCC which take advantage of the indistinguishability of molecules on opposite sides of the separator. It gives correct average ensemble statistics despite the introduction of an ensemble separator. SAFIRES is tested in simulations where the molecules on the two sides are treated in the same way, for a Lennard-Jones (LJ) liquid and a LJ liquid in contact with a surface, as well as for liquid modelling simulations using the TIP4P force field. Simulations using SAFIRES are shown to reproduce the unconstrained reference simulations without significant deviations.

研究动机与目标

• 开发一种基于边界的分区方法，确保在对内层和外层区域应用不同计算方法的混合模拟中，维持正确的系综统计。
• 消除区域间的粒子穿插，避免自适应耦合方案中常见的不连续性和能量漂移。
• 通过在界面处强制实现精确的弹性碰撞，改进现有边界方法（如FIRES、BEST和BCC），以保持动量和能量守恒。
• 实现动态、时间步长自适应的传播，使其在极限情况下退化为标准的朗之万或Verlet动力学。
• 在多个系统中验证该方法，包括Lennard-Jones液体和TIP4P水，结果与无约束模拟高度一致。

提出的方法

• SAFIRES使用可移动边界将模拟单元划分为内层和外层区域，该边界根据内层区域中最外侧粒子的位置进行扩展或收缩。
• 该方法在粒子穿越边界时强制执行弹性碰撞，通过基于旋转的碰撞模型同时守恒能量和动量。
• 引入多时间步长推进器，通过动态调整时间步长以精确检测碰撞发生的时刻。
• 碰撞求解采用旋转矩阵将相对速度矢量对齐至法线方向，随后基于约化质量进行动量交换。
• 该算法集成朗之万或速度Verlet推进器，时间步长缩放机制在摩擦系数为零或设为零时退化为标准动力学。
• 提供伪代码框架，包括冲突检测、时间步长外推以及基于向量代数和力传播的碰撞求解。

实验结果

研究问题

• RQ1 基于弹性碰撞的动态边界方法是否能在不发生区域间粒子转移的情况下，维持混合模拟中的正确系综统计？
• RQ2 与固定时间步长或近似边界方法（如FIRES）相比，SAFIRES中的动态时间步长调整在多大程度上提升了精度？
• RQ3 SAFIRES在Lennard-Jones液体和TIP4P水系统中，能否准确再现无约束模拟的径向分布函数（g(r)）？
• RQ4 弹性碰撞模型是否能防止粒子在边界处产生虚假累积，这一问题在FIRES和硬球边界方法中已被广泛报道？
• RQ5 SAFIRES是否能在极限情况下无缝退化为标准的朗之万或Verlet动力学，确保与现有模拟框架的兼容性？

主要发现

• SAFIRES模拟对密度ρ = 1.374 g cm⁻³、温度T = 94.4 K的Lennard-Jones液体，其径向分布函数（g(r)）与无约束参考模拟结果一致，无显著偏差。
• 对于较低密度的Lennard-Jones液体（ρ = 1.100 g cm⁻³），SAFIRES保持了准确的g(r)和边界位置概率分布，与参考模拟高度吻合。
• SAFIRES中的边界位置概率分布比FIRES更宽，表明其界面更具物理真实性，避免了人为的粒子累积。
• 当摩擦系数设为0.05时，该方法成功退化为标准朗之万动力学；当摩擦系数为零时，退化为速度Verlet动力学。
• 动态时间步长调整确保了碰撞的高精度检测，无需使用缓冲区或能量校正方案。
• 弹性碰撞模型通过强制实现瞬时动量和能量守恒，有效防止了粒子在边界处的虚假累积，解决了FIRES和硬球边界方法中的已知缺陷。