我试图弄清楚如何处理耦合。我正在尝试对聚合物链进行建模。我必须处理相邻角之间的角度、二面角和振动。我能够分别计算所有三个:我通过计算力来创建一个加速度矩阵。加速度矩阵的行=维度,列=粒子。矩阵中的每个向量代表真实空间中的一个单体。然后,我将加速度矩阵乘以时间步长的平方并将其添加到初始位置。这让我知道了单体是如何进化的。但是,我不确定如何将这三个结合起来。由于大多数程序按顺序执行代码,我将如何解决这个问题?有没有通用的方法?我正在考虑两种方法:计算包括所有加速度及其耦合的有效加速度,并在获得所有加速度后校正加速度。还有其他方法吗?任何建议或参考将不胜感激。
建模耦合:振动,旋转等
计算科学
分子动力学
2021-12-03 18:39:56
1个回答
用于更新分子动力学中位置和速度的著名且精确的算法是 Verlet 算法及其变体(速度和跨越式)。您可以将 Verlet 算法视为类似牛顿的积分器,将速度或跳跃式算法视为 Runge-Kutta 版本。
你会发现数学,在任何分子模拟教科书中都很简单,你会发现它非常接近你正在做的事情(这是一个简短的版本http://www.fisica.uniud.it/~ercolessi/md /md/node21.html)。
只要您的时间步长比系统的较短周期(大致)短十倍左右,它们就计算您的加速度矢量将原子中的所有力求和是没有问题的。否则,您将错过部分动态,或者由于位置漂移而完全不连贯。
其它你可能感兴趣的问题