小木虫emuch-net清华分子动力学模拟入门.pptx
《小木虫emuch-net清华分子动力学模拟入门.pptx》由会员分享,可在线阅读,更多相关《小木虫emuch-net清华分子动力学模拟入门.pptx(38页珍藏版)》请在咨信网上搜索。
1、引 言物质基本构成分子、原子在分子、原子这个微观水平上来考察物质:多体世界查清楚微观世界,宏观就清楚了从微观考虑问题的现实可行性从微观考虑问题的必要性物性的观测性参数:热传导、温度、压力、粘性、.微观处理的前提已知微观粒子间的相互作用假设分子为球,惰性,分子间的作用只取决于分子间的距离分子动力学(Molecular Dynamics,MD)MD的应用领域:物理、化学、生物、材料等MD方法能实时将分子的动态行为显示到计算机屏幕上,便于直观了解体系在一定条件下的演变过程MD含温度与时间,因此还可得到如材料的玻璃化转变温度、热容、晶体结晶过程、输送过程、膨胀过程、动态弛豫(relax)以及体系在外场
2、作用下的变化过程等水和离子在微小硅孔中的运动聚乙烯的结晶MD的基本原理用牛顿经典力学计算许多分子在相空间中的轨迹求解系统中的分子或原子间作用势能和系统外加约束共同作用的分子或原子的牛顿方程。模拟系统随时间推进的微观过程。通过统计方法得到系统的平衡参数或输运性质计算程序较为复杂,占用较多内存MD的主要步骤选取要研究的系统及其边界,选取系统内粒子间的作用势能模型设定系统中粒子的初始位置和初始动量建立模拟算法,计算粒子间作用力及各粒子的速度和位置当体系达到平衡后,依据相关的统计公式,获得各宏观参数和输运性质分子间势能及相互作用N个粒子系统的总势能刚球模型斥力力心点模型Southerland模型分子间
3、势能及相互作用Lennard-Jones势能 能量尺度;长度尺度 为方便,时常归一化:记 ;分子间势能及相互作用一些气体的参数kB=1.38x1023(J/K):Boltzmann常数NeonArgonKryponXenonNitrogen(nm)0.2750.34050.3600.4100.370/kB(K)36119.817122195分子间势能及相互作用相互作用标量形式:直角坐标:至此,各粒子间相互作用已知,可进行模拟了模拟的数学方法Euler法和Euler-Cromer方法?不能用:不能保持总能量守恒Verlet算法:速度形式模拟的数学方法Leap-frog算法:x的截断误差为 ,v的
4、截断误差模拟的数学方法边界条件模拟能力限制,不能模拟大量分子,只能模拟有限空间中的有限个分子:有限空间边界固体(刚性)边界条件不仅仅有分子间的相互作用,还引入了壁面的作用分子量大时,壁面作用可忽略不计模拟的数学方法取 ,前比值为0.20.01。取前值,模拟粗糙;取后值,模拟计算量太大处理方法:使用周期性边界条件周期性边界条件模拟的数学方法两个不同粒子在x或y方向上的最大分离距离为a/2最小像约定:两粒子分离距离最大分离距离,相互作用力可以忽略,而加入其中像粒子之一相互作用力来考虑模拟的数学方法考虑的粒子总数不变初始条件随机初始条件给法之一 要求 大小:条件一:规则给法条件二:随机给法模拟的数学
5、方法 random:随机数产生函数,产生(0,1)之间的随机数。方向(按球坐标给法):分量:模拟微观量温度 根据统计热力学,平衡态下经典系统的能量中的每一个二次项具有平均值kBT/2,即 注意:上式在系统质心速度为0时适用空间维数粒子个数:取时间平均模拟问题:如何给定系统的初始条件,得到所需要的平衡态温度Teq?解决方法之一:速度标定法 任给初始条件,模拟到平衡,得到系统平衡态温度T。一般TTeq。令 用速度 再模拟直到平衡,若所得温度仍不等于Teq,再进行上述过程给定初始条件:xi,vi计算到平衡态|T-Teq|f=Teq/Tvi=vi f 1/d计算结束YN其它方法:Gaussian热浴法
6、(约束温度调节方法)其基本原理在运动方程中加入“摩擦力”项,并将其与粒子速度联系起来。平衡态时,系统温度不变,因此dEk/dt=0宏观性质的统计系统的势能系统的内能系统的总能 E=Ep+Ek系统的温度模拟热容 定义热容 计算系统在温度T和T+T时的总能ET、ET+T,E:系统总能模拟压强对壁面的压强t时间里作用在单位面积壁上的压力时刻,速度为时刻,速度为刚性壁模拟粒子速度分布 选速度间隔v,模拟nt个时间步,记录在每个速度间隔中的粒子数,最后归一化。模拟气、液状态方程维里定理(Virial Theorem)压强体积粒子i的位置矢量粒子i所受到的其它粒子的合相互作用力温度的模拟可得此项在温度的模
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 小木虫 emuch net 清华 分子 动力学 模拟 入门
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【人****来】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【人****来】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。