实验8：水分子扩散系数分子动力学模拟-(2).pdf

《实验8：水分子扩散系数分子动力学模拟-(2).pdf》由会员分享，可在线阅读，更多相关《实验8：水分子扩散系数分子动力学模拟-(2).pdf（8页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、计算材料学实验讲义 实验二：分子动力学模拟-水分子扩散系数 一、前言 分子动力学模拟的基本思想是将物质看成是原子和分子组成的粒子系统（many-body systems），设置初始位能模型，通过分析粒子的受力状况，计算粒子的牛顿运动方程，得到粒子的空间运动轨迹，可以求得复杂体系的热力学参数以及结构和动力学性质。分子动力学模拟的理论是统计力学中的各态历经假说(Ergodic Hypothesis)，即保守力学系统从任意初态开始运动，只要时间足够长，它将经过相空间能量曲面上的一切微观运动状态，系统力学量的系综平均等效力学量的时间平均，因此可以通过计算系综的经典运动方程来得到力学量的性质。比如，由N

2、个粒子组成的系综的势能计算函数为：intUUUVDW (1-1)VDWU表示粒子内和粒子之间的Van der Waals相互作用；intU表示粒子的内部势能（键角弯曲能，键伸缩能、键扭转能等）；根据经典力学方程，系统中第i个粒子的受力大小为：UkzjyixUFiiiii (1-2)那么第i个粒子的加速度可以通过牛顿第二定律得到：iiimtFta (1-3)由于体系有初始位能，每个粒子有初始位置和速度，那么加速度对时间进行积分，速度对时间积分就可以获得各个任意时刻粒子的速度和位置：iiiavdtdrdtd22 (1-4)tavviii0 (1-5)20021tatvrriiii (1-6)ir和

3、v分别是系统中粒子t时刻的位置和速度，0ir和0iv分别是系统中粒子初始时刻的位置和速度。依据各态历经假说，可获得任意物理量 Q的系综平均，因此得到体系的相关性质：dttrQttQQtt01lim (1-7)分子动力学模拟能够计算体系的能量，粒子间的相互作用，角动量，角度以及二面角分布，剪切粘度，结构参数，压力参数，热力学参数，弹性性质，动力学性质等。均方位移则是表征粒子动力特性的重要参量，它是粒子位移平方的平均值（mean square displacement,MSD）表达式如下：2()()(0)MSDR tr tr 式中，括号表示平均值，()ir t表示时间t时粒子i的位置。根据爱因斯坦

4、的扩散定律：式中，D为粒子的扩散系数（diffusion constant）。因此，当时间很长时，均方位移曲线的斜率即为6D。因此可根据MSD曲线斜率求得粒子的扩散系数。二、实验目的（1）了解分子动力学方法的原理。（2）掌握 MS 软件的基本操作和分子动力学模拟的基本步骤。（3）掌握粒子扩散系数的计算方法。三、实验内容 以水分子的扩散系数计算例，展开详述。1.回顾基本操作。2.模型构建。（1）构建分子。首先构建水分子，并命名为water。（2）分配力场电荷。利用分子力学方法对分子结构进行优化，首先需要给构建的分子分配力场，选择 discover 模块中的 setup，出现如下对话框：选择 En

5、ergy，确定 Force field 选择为 compass；选择 Typing，点击 Calculate，给定分子每个原子力场参数，即分子内及分子间的相互作用势；点击 Automation，确认相关选项均选择 No；选择 Non-Bond，Apply settings to 选择 vdWCoulomb，Summation method 选择 Atom based，Quality选择Medium，点击More，cutoff distance选择9.5；点击Job Control选择My Computer。（3）优化分子。选择 discover 模块中的 Minimizer，出现如下对话框：M

6、ethod 选择 Smart Minimizer，Convergence level 选择 Medium，点击 Minimize，开始进行优化。优化完毕后，在 Project 中会出现文件名为 water Disco Min 的文件夹，文件夹内名为water.xsd 的文件即为优化后的分子结构。如下图所示：（4）定义 Group。Modify-Charges-Charge Group-Define(需选中原子)（5）构建水盒子。点击 Amorphous-Construction(Legacy)-Add(选中 water.xsd点击 Add，然后修改添加水分子的数目 Number：300，密度设

7、为 1)具体参数设置如下图-Construct-弹出 water AC Constr 文件夹-water.xtd 文件即构建的水盒子。3.优化体系。选择 discover 模块中的 Minimizer，具体参数设置 Non-Bond，Apply settings to 选择vdWCoulomb，Summation method 选择 Group based，其它同分子结构优化参数设置。最后点击 Minimize，完成优化：water.xsd 文件。4.分子动力学模拟。（1）打开上述优化好的 3D 结构，从工具栏中选择，然后选择 Dynamics，或者从菜单上选择 Modules-Discove

8、r-Dynamics，出现如下对话框：将 dynamics time 改为 100ps，Save 改为 Full，Frame output every 改为 1000，点击 Run.计算需要花费几分钟的时间，运行完毕后，出现名为 water Disco Dynamics 的文件夹-water.xtd。5.数据分析。（1）定义 set：选中所有水分子-Edit-Edit Sets-New（Name 改为 water）-OK（2）提取 MSD 曲线 discover 模块中的 Analysis-Dynamic-Mean squared displacement-（参数设置如图所示）-Analyze MSD 曲线如下图：分别是沿 X、Y、Z 轴方向和总的均方根位移。取总的 MSD（在图上右击-Delete Graph 选中要删掉的曲线然后点击 OK），拷贝到 Exce 中，求其斜率从而得到水分子的扩散系数。根据斜率得到水分子的扩散系数：2.13*10-9m2/s，实验值：2.0-2.5*10-9m2/s 020406080100020406080100120140 nmt/psy=1.275x+0.515