实验8：水分子扩散系数分子动力学模拟

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

2、方程来得到力学量的性质。比如，由N个粒子组成的系综的势能计算函数为：UUVDWUint(1-1)(1-2)Uvdw表示粒子内和粒子之间的VanderWaals相互作用；U表示粒子的内部势能(键角弯曲能，键伸缩能、键扭转能等)；根据经典力学方程，系统中第i个粒子的受力大小为：FiiUijkUXiyiZi那么第i个粒子的加速度可以通过牛顿第二定律得到：aitFitmi(1-3)由丁体系有初始位能，每个粒子有初始位置和速度，那么加速度对时间进行积分，速度对时间积分就可以获得各个任意时刻粒子的速度和位置:d22ridt2dVidtai0viviait(1-4)(1-5)(1-6)00.1.2ririV

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

4、律:曲线式中，D为粒子的扩散系数(diffusionconstant)。因此，当时间很长时，均方位移的斜率即为6D。因此可根据MSD曲线斜率求得粒子的扩散系数。二、实验目的(1) 了解分子动力学方法的原理。(2) 掌握MS软件的基本操作和分子动力学模拟的基本步骤。(3) 掌握粒子扩散系数的计算方法。三、实验内容以水分子的扩散系数计算例，展开详述。1.回顾基本操作。模型构建(1)构建分子。首先构建水分子，并命名为water。同lHap3-Batv-rialfadli#-Im椅lji事d-匚X(2)分配力场电荷。利用分子力学方法对分子结构进行优化，首先需要给构建的分子分配力场，选择discover

5、模块中的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模块'*中的Minimize

6、r,出现如下对话框：MMhnd|SitMriirw?eirMoreCcoefjcncclie、dMaiMmLmrerAtonit|50001厂Optimizece|ImEftE”，w»:IMinimife=Fic>.Melp司PicoTErMinx&TzatLonMethod选择SmartMinimizer,Convergencelevel选择Medium，点击Minimize,开女台进行优化。优化完毕后，在Project中会出现文件名为waterDiscoMin的文件火，文件夹内名为water.xsd的文件即为优化后的分子结构。如下图所示:9-SaitiEEiai耳顷

7、山。-顷in*制”句。£lLiElm%r词L0地lLIIo4lLei4.LjE124E心心Jld*F山0'll比审gfl-jli.L-e«LU/-Setuid目dpranf|zoutK曲-,i1TLaJub如；w登1fcju哆t"Nil定义GroupoModify-Charges-ChargeGroup-Define(需选中原子)Cloul如|EditCKegeGroups|罔四己yiuupikalmlglijr-T球Acharwtolerance|0.1Edtchaig&gioup?hofinmalgE可耽mRen»呢chaioeora

8、os(4) 构建水盒子。点击Amorphous-Construction(Legacy)-Add(选中water.xsd点击Add,然后修改添加水分子的数目Number:300,密度设为1)具体参数设置如下图-Construct-弹出waterACConstr文件夹-water.xtd文件即构建的水盒子。-rTx|Fila：£zitVjcaI口dijyjhialdT口。w£tfal3S.liq5.ftcidjdLts虫二日叩Ju.lpJ号孔s中厂书-I*9期EMM加iM|L职dL通血|-'rehrence-rl匕Hi©.DrpfaousCellCona?t

9、rucrt:ionJ'：1-1-I1-Z*A*'|T1'H"t.LaJcijJMluTiFill£aorphaxuCellCantrnctixm-xLfipMtergcriidCorsln-clLiquidCratalIMkuIo3如SILttdhl-rOtxtlucrJtmofecLle：.IDel国已IcmcrrflfijKJCdllIweNuirbpnlrnrt-cuiailDnx2980|Peri&dictslllp1argstcfensb险finalc«rfomratcnr1a/ccI"Ranpdeltafioo

10、ianirKd心遍rfigjccSIparcrctslAngMrariab;c;Ispotfjo.bzJar-laora-/-iBirebixtc<nliouia«&urdbwnsMh-uctcnREftn®czntijjrstonsfoilDwngcancluciohEuldaEEati匚liquid匚rpEfdlri匕f-diZELignstig|Fl蛛-|队耳Erieigpei/aLahcnFixwWd|可协驻Mqi#,WlisttrjpbedcxiteftAtriDmi3l:ei匕,gpzetipJobClTiIIlI口acDwwWsM的匚网如T心I3

11、Jabdxriaian厂Aulcjiic>匚匚glUE2. 优化体系。选择discover模块卜中的Minimizer，具体参数设置Non-Bond,Applysettingsto选择vdW&Coulomb,Summationmethod选择Groupbased,其它同分子结构优化参数设置。最后点击Minimize，完成优化：water.xsd文件。3. 分子动力学模拟。(1)打开上述优化好的3D结构，从工具栏中选择*'，然后选择Dynamics,或者从菜单上选择Modules-Discover-Dynamics,出现如下对话框：将dynamicstime改为100ps,

12、Save改为Full,Frameoutputevery改为1000,点击Run.计算需要花费几分钟的时间，运行完毕后，出现名为waterDiscoDynamics的文件夹-water.xtd。4. 数据分析。(1)定义set：选中所有水分子-Edit-EditSets-New(Name改为water)-OK(2)提取MSD曲线discover模块尹中的Analysis-Dynamic-Meansquareddisplacement-（参数设置如图所示）-Analyze0DistoverAnalyses+FluctuationDynamicMeansquareddispldcemeritDipolemomentadtocorrelatiarifunclionVelocityeutixorreldionfunctionRotationaltinn&correkticnfunctiortSpacetimecorrelationfunctioneortelationfundionMechanicalApplianalysistc,Log