模拟实验：水分子扩散系数

《计算材料学》实验讲义实验二：分子动力学模拟-水分子扩散系数一、前言分子动力学模拟的基本思想是将物质看成是原子和分子组成的粒子系统（many-bodysystems），设置初始位能模型，通过分析粒子的受力状况，计算粒子的牛顿运动方程，得到粒子的空间运动轨迹，可以求得复杂体系的热力学参数以及结构和动力学性质。分子动力学模拟的理论是统计力学中的各态历经假说(ErgodicHypothesis)，即保守力学系统从任意初态开始运动，只要时间足够长，它将经过相空间能量曲面上的一切微观运动状态，系统力学量的系综平均等效力学量的时间平均，因此可以通过计算系综的经典运动方程来得到力学量的性质。比如，由N个粒子组成的系综的势能计算函数为：intUUUVDW(1-1)VDWU表示粒子内和粒子之间的VanderWaals相互作用；intU表示粒子的内部势能（键角弯曲能，键伸缩能、键扭转能等）；根据经典力学方程，系统中第i个粒子的受力大小为：UkzjyixUFiiiii(1-2)那么第i个粒子的加速度可以通过牛顿第二定律得到：iiimtFta(1-3)由于体系有初始位能，每个粒子有初始位置和速度，那么加速度对时间进行积分，速度对时间积分就可以获得各个任意时刻粒子的速度和位置：iiiavdtdrdtd22(1-4)tavviii0(1-5)20021tatvrriiii(1-6)ir和v分别是系统中粒子t时刻的位置和速度，0ir和0iv分别是系统中粒子初始时刻的位置和速度。依据各态历经假说，可获得任意物理量Q的系综平均，因此得到体系的相关性质：dttrQttQQtt01lim(1-7)分子动力学模拟能够计算体系的能量，粒子间的相互作用，角动量，角度以及二面角分布，剪切粘度，结构参数，压力参数，热力学参数，弹性性质，动力学性质等。均方位移则是表征粒子动力特性的重要参量，它是粒子位移平方的平均值（meansquaredisplacement,MSD）表达式如下：2()()(0)MSDRtrtr式中，括号表示平均值，()irt表示时间t时粒子i的位置。根据爱因斯坦的扩散定律：式中，D为粒子的扩散系数（diffusionconstant）。因此，当时间很长时，均方位移曲线的斜率即为6D。因此可根据MSD曲线斜率求得粒子的扩散系数。二、实验目的（1）了解分子动力学方法的原理。（2）掌握MS软件的基本操作和分子动力学模拟的基本步骤。（3）掌握粒子扩散系数的计算方法。三、实验内容以水分子的扩散系数计算例，展开详述。1.回顾基本操作。2.模型构建。（1）构建分子。首先构建水分子，并命名为water。（2）分配力场电荷。利用分子力学方法对分子结构进行优化，首先需要给构建的分子分配力场，选择discover模块中的setup，出现如下对话框：选择Energy，确定Forcefield选择为compass；选择Typing，点击Calculate，给定分子每个原子力场参数，即分子内及分子间的相互作用势；点击Automation，确认相关选项均选择No；选择Non-Bond，Applysettingsto选择vdW﹠Coulomb，Summationmethod选择Atombased，Quality选择Medium，点击More，cutoffdistance选择9.5Å；点击JobControl选择MyComputer。（3）优化分子。选择discover模块中的Minimizer，出现如下对话框：Method选择SmartMinimizer，Convergencelevel选择Medium，点击Minimize，开始进行优化。优化完毕后，在Project中会出现文件名为waterDiscoMin的文件夹，文件夹内名为water.xsd的文件即为优化后的分子结构。如下图所示：（4）定义Group。Modify-Charges-ChargeGroup-Define(需选中原子)（5）构建水盒子。点击Amorphous-Construction(Legacy)-Add(选中water.xsd点击Add，然后修改添加水分子的数目Number：300，密度设为1)具体参数设置如下图-Construct-弹出waterACConstr文件夹-water.xtd文件即构建的水盒子。3.优化体系。选择discover模块中的Minimizer，具体参数设置Non-Bond，Applysettingsto选择vdW﹠Coulomb，Summationmethod选择Groupbased，其它同分子结构优化参数设置。最后点击Minimize，完成优化：water.xsd文件。4.分子动力学模拟。（1）打开上述优化好的3D结构，从工具栏中选择，然后选择Dynamics，或者从菜单上选择Modules-Discover-Dynamics，出现如下对话框：将dynamicstime改为100ps，Save改为Full，Frameoutputevery改为1000，点击Run.计算需要花费几分钟的时间，运行完毕后，出现名为waterDiscoDynamics的文件夹-water.xtd。5.数据分析。（1）定义set：选中所有水分子-Edit-EditSets-New（Name改为water）-OK（2）提取MSD曲线discover模块中的Analysis-Dynamic-Meansquareddisplacement-（参数设置如图所示）-AnalyzeMSD曲线如下图：分别是沿X、Y、Z轴方向和总的均方根位移。取总的MSD（在图上右击-DeleteGraph选中要删掉的曲线然后点击OK），拷贝到Exce中，求其斜率从而得到水分子的扩散系数。根据斜率得到水分子的扩散系数：2.13*10-9m2/s，实验值：2.0-2.5*10-9m2/s020406080100020406080100120140r2nmt/psy=1.275x+0.515