欢迎来到三一办公! | 帮助中心 三一办公31ppt.com(应用文档模板下载平台)
三一办公
全部分类
  • 办公文档>
  • PPT模板>
  • 建筑/施工/环境>
  • 毕业设计>
  • 工程图纸>
  • 教育教学>
  • 素材源码>
  • 生活休闲>
  • 临时分类>
  • ImageVerifierCode 换一换
    首页 三一办公 > 资源分类 > DOC文档下载  

    实验8水分子扩散系数分子动力学模拟_2.doc

    • 资源ID:1132012       资源大小:458.50KB        全文页数:8页
    • 资源格式: DOC        下载积分:10金币
    快捷下载 游客一键下载
    会员登录下载
    三方登录下载: 微信开放平台登录 QQ登录  
    下载资源需要10金币
    邮箱/手机:
    温馨提示:
    用户名和密码都是您填写的邮箱或者手机号,方便查询和重复下载(系统自动生成)
    支付方式: 支付宝    微信支付   
    验证码:   换一换

    加入VIP免费专享
     
    账号:
    密码:
    验证码:   换一换
      忘记密码?
        
    友情提示
    2、PDF文件下载后,可能会被浏览器默认打开,此种情况可以点击浏览器菜单,保存网页到桌面,就可以正常下载了。
    3、本站不支持迅雷下载,请使用电脑自带的IE浏览器,或者360浏览器、谷歌浏览器下载即可。
    4、本站资源下载后的文档和图纸-无水印,预览文档经过压缩,下载后原文更清晰。
    5、试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓。

    实验8水分子扩散系数分子动力学模拟_2.doc

    计算材料学实验讲义实验二:分子动力学模拟-水分子扩散系数一、前言分子动力学模拟的根本思想是将物质看成是原子和分子组成的粒子系统many-body systems,设置初始位能模型,通过分析粒子的受力状况,计算粒子的牛顿运动方程,得到粒子的空间运动轨迹,可以求得复杂体系的热力学参数以与结构和动力学性质。分子动力学模拟的理论是统计力学中的各态历经假说(Ergodic Hypothesis),即保守力学系统从任意初态开始运动,只要时间足够长,它将经过相空间能量曲面上的一切微观运动状态,系统力学量的系综平均等效力学量的时间平均,因此可以通过计算系综的经典运动方程来得到力学量的性质。比如,由N个粒子组成的系综的势能计算函数为: (1-1)表示粒子内和粒子之间的Van der Waals相互作用;表示粒子的内部势能键角弯曲能,键伸缩能、键扭转能等;根据经典力学方程,系统中第i个粒子的受力大小为: (1-2)那么第i个粒子的加速度可以通过牛顿第二定律得到: (1-3)由于体系有初始位能,每个粒子有初始位置和速度,那么加速度对时间进展积分,速度对时间积分就可以获得各个任意时刻粒子的速度和位置: (1-4) (1-5) (1-6)和分别是系统中粒子t时刻的位置和速度,和分别是系统中粒子初始时刻的位置和速度。依据各态历经假说,可获得任意物理量Q的系综平均,因此得到体系的相关性质: (1-7)分子动力学模拟能够计算体系的能量,粒子间的相互作用,角动量,角度以与二面角分布,剪切粘度,结构参数,压力参数,热力学参数,弹性性质,动力学性质等。均方位移如此是表征粒子动力特性的重要参量,它是粒子位移平方的平均值mean square displacement, MSD表达式如下:式中,括号表示平均值,表示时间t时粒子i的位置。 根据爱因斯坦的扩散定律:式中,D为粒子的扩散系数diffusion constant。因此,当时间很长时,均方位移曲线的斜率即为6D。因此可根据MSD曲线斜率求得粒子的扩散系数。二、实验目的1了解分子动力学方法的原理。2掌握MS软件的根本操作和分子动力学模拟的根本步骤。3掌握粒子扩散系数的计算方法。三、实验内容以水分子的扩散系数计算例,展开详述。1. 回顾根本操作。2. 模型构建。1构建分子。首先构建水分子,并命名为water。2分配力场电荷。利用分子力学方法对分子结构进展优化,首先需要给构建的分子分配力场,选择discover模块中的setup,出现如下对话框:选择Energy,确定Force field选择为pass;选择Typing,点击Calculate,给定分子每个原子力场参数,即分子内与分子间的相互作用势;点击Automation,确认相关选项均选择No;选择Non-Bond,Apply settings to选择vdWÅ;点击Job Control选择My puter。3优化分子。选择discover模块中的Minimizer,出现如下对话框:Method选择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,密度设为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-Discover-Dynamics,出现如下对话框:将dynamics time改为100ps,Save改为Full,Frame output every改为1000,点击Run.计算需要花费几分钟的时间,运行完毕后,出现名为water Disco Dynamics的文件夹-water.xtd。5. 数据分析。1定义set:选中所有水分子-Edit-Edit Sets-NewName改为water-OK2提取MSD曲线discover模块中的Analysis-Dynamic-Mean squared displacement-参数设置如下列图-AnalyzeMSD曲线如如下图:分别是沿X、Y、Z轴方向和总的均方根位移。取总的MSD在图上右击-Delete Graph选中要删掉的曲线然后点击OK,拷贝到Exce中,求其斜率从而得到水分子的扩散系数。根据斜率得到水分子的扩散系数:2.13*10-9m2/s,实验值:2.0-2.5*10-9m2/s

    注意事项

    本文(实验8水分子扩散系数分子动力学模拟_2.doc)为本站会员(李司机)主动上传,三一办公仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知三一办公(点击联系客服),我们立即给予删除!

    温馨提示:如果因为网速或其他原因下载失败请重新下载,重复下载不扣分。




    备案号:宁ICP备20000045号-2

    经营许可证:宁B2-20210002

    宁公网安备 64010402000987号

    三一办公
    收起
    展开