[论文解读] A Reversible Unwrapping Algorithm for Constant Pressure Molecular Dynamics Simulations

Martin Kulke, Josh V. Vermaas|arXiv (Cornell University)|Nov 23, 2021

Nanopore and Nanochannel Transport Studies参考文献 41被引用 8

一句话总结

本文提出了一种用于恒压分子动力学模拟的混合可逆展开算法，通过考虑周期性盒子体积的变化，校正了现有启发式和位移展开方法中的伪影。该方法采用改进的位移方案，并引入基于盒子尺寸演化的附加校正项，通过保留分子几何结构和正确的扩散系数，显著提升了长时间模拟的准确性。

ABSTRACT

Molecular simulation technologies have afforded researchers a unique look into the nanoscale interactions driving physical processes. However, a limitation for molecular dynamics (MD) simulations is that they must be performed on finite-sized systems in order to map onto computational resources. To minimize artifacts arising from finite-sized simulation systems, it is common practice for MD simulations to be performed with periodic boundary conditions (PBC). However, in order to calculate specific physical properties, such as mean square displacements to calculate diffusion coefficients, continuous particle trajectories where the atomic movements are continuous and do not jump between cell faces are required. In these cases, modifying atomic coordinates through unwrapping schemes are an essential post-processing tool to remove these jumps. Here, two established trajectory unwrapping schemes are applied to 1us wrapped trajectories for a small water box. The existing schemes can result in spurious diffusion coefficients, long bonds within unwrapped molecules, and inconsistent atomic coordinates when coordinates are rewrapped after unwrapping. We determine that prior unwrapping schemes do not account for changing periodic box dimensions, and introduce an additional correction term to the existing displacement unwrapping scheme by von Bulow et al. to correct for these artifacts. We also demonstrate that the resulting algorithm is a hybrid between the existing heuristic and displacement unwrapping schemes. After treatment with this new unwrapping scheme, molecular geometries are correct even after long simulations. In anticipation for longer molecular dynamics trajectories, we develop implementations for this new scheme in multiple PBC handling tools.

研究动机与目标

• 解决现有轨迹展开方法在NPT分子动力学模拟中产生的不准确性。
• 识别出在盒子尺寸波动的长时间模拟中，启发式和位移展开方案会失效。
• 开发一种新展开算法，即使在周期性盒子体积变化的情况下，也能保持正确的分子几何结构和扩散行为。
• 将新算法集成到广泛使用的PBC工具中，以促进分子模拟社区的广泛应用。

提出的方法

• 所提出的方法结合了启发式和位移展开方案的要素，以位移矢量作为展开的主要依据。
• 引入校正项 m(Li+1 − Li)，以考虑时间步之间周期性盒子尺寸的变化。
• 该算法使用前一时刻的展开位置、包裹的位移以及基于盒子尺寸调整的体积波动校正项，计算展开坐标。
• 校正项确保原子位置不仅根据晶胞穿越进行调整，也根据盒子的膨胀或收缩进行调整。
• 该方法已实现在VMD工具 pbctools、fastpbc 和 qwrap 中，以确保与常见分子动力学分析工作流的兼容性。
• 该方法具有可逆性，可实现轨迹的精确重包裹，且不会引起几何畸变。

实验结果

研究问题

• RQ1 在盒子尺寸波动的长时间NPT模拟中，启发式展开方法是否能产生准确的扩散系数？
• RQ2 在盒子尺寸变化的长时间模拟中，位移展开方法能否保持正确的分子几何结构？
• RQ3 盒子尺寸波动如何影响恒压分子动力学中轨迹展开的准确性？
• RQ4 一种考虑盒子尺寸演化的混合展开方案，能否同时提升几何保真度和扩散系数的准确性？
• RQ5 新算法是否具有可逆性，并适合集成到标准分子动力学分析工具中？

主要发现

• 启发式方法在NPT模拟中约100 ns后失效，产生虚假的扩散系数，原因在于盒子尺寸波动。
• 位移方法在重包裹后导致错误的键长和不一致的分子间距离，尤其在长时间模拟中更为明显。
• 所提出的混合方案（公式12）即使在1 µs的模拟后，仍能保持水分子几何结构的正确性，键长接近0.1 nm。
• 该混合方法在单精度浮点数精度范围内保持原子位置准确，RMSD的增加仅源于大原子位置，而非几何畸变。
• 新算法通过引入项 m(Li+1 − Li) 有效校正了盒子尺寸变化，消除了启发式方法中对扩散的高估。
• 该方法已成功实现在VMD的 pbctools、fastpbc 和 qwrap 中，实现了在分子动力学分析中的广泛应用。