[发明专利]一种计算分子动力学动态时间步长的方法

权利要求书

1. 一种计算分子动力学动态时间步长的方法，包括有多个分子，其特征在于，具体计算步骤如下：

1）测量所有分子的坐标，计算每两个分子之间的距离                                               ，i和j为分子编号，判断数组的数值是否都大于截断半径r_c，若大于则转入步骤2），若小于则结束；

2）测量所有分子的运动速度，计算每两个分子之间的相对速度；

3）根据步骤1）和步骤2）中所测得的每两个分子之间的距离和速度，计算每两个分子之间运动到截断半径所需要的时间t(n)；

4）在数组t(n)中取出最小值t_min作为下一步的时间步长。

2. 如权利要求1所述的一种计算分子动力学动态时间步长的方法，其特征在于：步骤1）中计算每两个分子之间距离的公式为

，

式中，x、y和z为分子的空间坐标值，i和j为分子编号。

3. 如权利要求1所述的一种计算分子动力学动态时间步长的方法，其特征在于：步骤2）中计算每两个分子之间距离的公式为

式中，为分子在x方向的速度，为分子在y方向的速度，为分子在z方向的速度。

4. 如权利要求1所述的一种计算分子动力学动态时间步长的方法，其特征在于，步骤3）中计算每两个分子之间运动到截断半径所需要的时间t(n)的公式为

t(n)为由组成的数组。