分子动力学模拟计算水分子扩散系数.pptx

计算材料学实验分子动力学模拟计算水分子扩散系数1、实验目的（1）了解分子动力学方法的原理。（2）熟练MS软件构建分子和体相的基本操作。（3）掌握分子动力学模拟的基本步骤。（4）掌握粒子扩散系数的计算方法。2、实验原理（1）分子动力学模拟（moleculardynamicssimulation，简称MD）方法首先是由Alder和Wainwright提出的，现已逐渐成为预测系统特性、验证理论和改进模型的计算工具。MD方法的基本思想是把物质看成由原子和分子组成的粒子系统，从该体系的某一假定的位能模型出发，并假定体系粒子的运动遵循经典力学或量子力学描述的规律，若已知粒子的所有受力作用，则可以求解出运动方程而得到系统中全体粒子在相空间中的轨道，然后统计得到系统的热力学参数、结构和输运特性等。也就是由体系的微观性质来求算其宏观性质，属于微观尺度的模拟技术。2、实验原理（2）均方根位移：分子动力学计算系统中原子由起始位置不停移动，每一瞬间位置皆不同。以表示时间t时粒子i的位置。粒子位移平方的平均值称为均方根位移（MeanSquareDisplacement=MSD）：（3）粒子扩散系数：根据爱因斯坦扩散定律：式中D为粒子的扩散系数，当时间很长时，均方根位移曲线的斜率就是6D，因此可以通过分子动力学计算MSD曲线斜率得到粒子的扩散系数。3、实验内容（1）回顾创建工程和构建分子的基本操作；（2）构建水体系并进行分子动力学模拟；（3）计算水分子的扩散系数。3.1回顾创建工程和构建分子的基本操作（1）创建工程（2）构建分子3.2构建纯水体系并进行分子动力学模拟。（1）构建水分子，命名：water。3.2构建纯水体系并进行分子动力学模拟。（2）构型优化：点击discover模块中setup，在Energy中选择力场compass；在Typing中点击Calculate，对分子分配电荷和立场参数。点击discover模块中Minimizer，选择SmartMinimizer方法，精度选择Medium，点击Minimize，开始进行优化。3.2构建纯水体系并进行分子动力学模拟。（3）定义Group：选中优化后的整个水分子，点击Modify-Charges选择ChargeGroup，然后点击Define。3.2构建纯水体系并进行分子动力学模拟。（4）构建水盒子：打开water.xsd文件，点击Amorphous-Construction(Legacy)中Add选项Number：300，密度设为1，点击setup选择Usegroupbasedcutoffs，其它具体参数设置如下图，点击Construct弹出waterACConstr文件夹中water.xtd文件即纯水体系模型。将其名字修改为box.xtd3.2构建纯水体系并进行分子动力学模拟。（5）优化水盒子：点击discover模块中setup中将non-bond的vdwCoulomb设置为groupbased，在点击discover中Minimizer，选择SmartMinimizer方法，精度选择Medium，点击Minimize，完成优化，box.xsd文件即优化后的体系。（6）分子动力学模拟:打开上述优化过的模型，点击discover中Dynamics系综选择NVT；温度定为289K；模拟步长1fs模拟时间50ps，每隔1000针记录一次数据。点击Run对体系进行分子动力学模拟，弹出boxDiscoverDynamics文件夹，其中box.xtd文件即体系模拟的轨迹文件。3.2构建纯水体系并进行分子动力学模拟。（7）计算扩散系数：1.定义set：选中所有水分子在工具栏中点击Edit-EditSets中New选项-OK。3.2构建纯水体系并进行分子动力学模拟。（8）计算扩散系数：2.提出水分子的MSD曲线：点击discover模块中Analysis选择Dynamic-Meansquareddisplacement-点击define-Addtolist再点击Analyze。得到MSD曲线如左图所示：3.2构建纯水体系并进行分子动力学模拟。（8）计算扩散系数：3.MSD曲线图中包含水分子在沿X、Y、Z轴方向以及总的四条均方根位移。取总的MSD曲线（在图上右击-DeleteGraph选中要删掉的曲线然后点击OK），拷贝到Exce中，求其斜率，根据MSD与扩散系数的关系式算得水分子的扩散系数。3.2作业模拟H2O在缓蚀剂膜A、B中的扩散系数，判断三种缓蚀剂的缓蚀效率。B-R:A-缓蚀剂waterA(7)B(19)扩散系数（100水分子）2.130.38280.7647????模拟步骤：构建水分子结构，优化结构，定义chargegrou