清华分子动力学模拟入门课件.ppt
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《清华分子动力学模拟入门课件.ppt》由会员分享,可在线阅读,更多相关《清华分子动力学模拟入门课件.ppt(39页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、多体系统分子动力学Multi-body SystemMolecular Dynamics引 言?物质基本构成分子、原子? 在分子、原子这个微观水平上来考察物质:多体世界? 查清楚微观世界,宏观就清楚了?从微观考虑问题的现实可行性?从微观考虑问题的必要性? 物性的观测性参数:热传导、温度、压力、粘性、. ?微观处理的前提? 已知微观粒子间的相互作用?假设? 分子为球,惰性,分子间的作用只取决于分子间的距离?分子动力学 (Molecular Dynamics,MD)MD的应用?领域:物理、化学、生物、材料等?MD方法能实时将分子的动态行为显示到计算机屏幕上, 便于直观了解体系在一定条件下的演变过程
2、?MD含温度与时间, 因此还可得到如材料的玻璃化转变温度、热容、晶体结晶过程、输送过程、膨胀过程、动态弛豫(relax)以及体系在外场作用下的变化过程等?水和离子在微小硅孔中的运动?聚乙烯的结晶MD的基本原理?用牛顿经典力学计算许多分子在相空间中的轨迹? 求解系统中的分子或原子间作用势能和系统外加约束共同作用的分子或原子的牛顿方程。? 模拟系统随时间推进的微观过程。? 通过统计方法得到系统的平衡参数或输运性质? 计算程序较为复杂,占用较多内存MD的主要步骤?选取要研究的系统及其边界,选取系统内粒子间的作用势能模型?设定系统中粒子的初始位置和初始动量?建立模拟算法,计算粒子间作用力及各粒子的速度
3、和位置?当体系达到平衡后,依据相关的统计公式,获得各宏观参数和输运性质分子间势能及相互作用?N个粒子系统的总势能?NjiijNNNNrVrVrVrVrVrVrVV11,-22311312)()()()()()()(?dVrVr刚球模型( )0rdV rrd?( )()dV rr?斥力力心点模型dr V-?( )()rdV rdrdr?Southerland模型分子间势能及相互作用?Lennard-Jones势能能量尺度;长度尺度为方便,时常归一化:记;?6124)(rrrV?612/1/14)(?rrrVVV?/rr?/?612114)(rrrV V(r) F(r)排斥力 V(r), F(r)
4、r吸引力分子间势能及相互作用?一些气体的参数kB=1.38x10?23(J/K): Boltzmann常数NeonArgonKryponXenonNitrogen?(nm)0.2750.34050.3600.4100.370?/kB(K)36119.817122195分子间势能及相互作用?相互作用? 标量形式:? 直角坐标:至此,各粒子间相互作用已知,可进行模拟了)(rVrmf?( )( ),yxffV r xV r ymrrmrr? ? ?)(rVmF?模拟的数学方法?Euler法和Euler-Cromer方法?? 不能用:不能保持总能量守恒? Verlet算法:速度形式?121212121
5、11212112nnnnnnnnnnnnnnxxvtatvvaatvvatvvat? ?模拟的数学方法?Leap-frog算法:x的截断误差为,v的截断误差1/21/211/21/22nnnnnntnnnvvtaxxtvvva? ? ?4()Ot?2()Ot?模拟的数学方法?边界条件? 模拟能力限制,不能模拟大量分子,只能模拟有限空间中的有限个分子:有限空间?边界? 固体(刚性)边界条件?不仅仅有分子间的相互作用,还引入了壁面的作用?分子量大时,壁面作用可忽略不计233611aaaN?和壁面作用分子数壁面积总分子数体积3aN?总分子数壁面积和壁面作用分子数?模拟的数学方法取 , 前比值为0.2
6、0.01。取前值,模拟粗糙;取后值,模拟计算量太大处理方法:使用周期性边界条件? 周期性边界条件210?N610? ?a?模拟的数学方法两个不同粒子在x或y方向上的最大分离距离为 a/2? 最小像约定:两粒子分离距离 最大分离距离,相互作用力可以忽略,而加入其中像粒子之一相互作用力来考虑?2a?2/a?a模拟的数学方法考虑的粒子总数不变?初始条件? 随机初始条件给法之一要求大小:条件一:规则给法条件二:随机给法max|Vv ?) 1random2(max? Vv模拟的数学方法random:随机数产生函数,产生 (0,1)之间的随机数。方向(按球坐标给法):分量 :xyz?)5 . 0rando
7、m(2arccos(?)5 . 0random(2arccos)0.5random(sign?180180;1800?cossinvvx?sinsinvvy?cosvvz?模拟?微观量? 温度根据统计热力学,平衡态下经典系统的能量中的每一个二次项具有平均值 kBT/2,即注意:上式在系统质心速度为 0时适用?iiiBvmTNkd2212空间维数粒子个数:取时间平均模拟问题:如何给定系统的初始条件,得到所需要的平衡态温度Teq?解决方法之一:速度标定法任给初始条件,模拟到平衡,得到系统平衡态温度T。一般T?Teq。令用速度再模拟直到平衡,若所得温度仍不等于 Teq,再进行上述过程TTfeq/?1
8、/ diifvv?给定初始条件:xi,vi计算到平衡态|T - Teq| f = Teq/ Tvi= vi f 1/d计算结束YN其它方法:Gaussian热浴法(约束温度调节方法)其基本原理在运动方程中加入“摩擦力”项,并将其与粒子速度联系起来。平衡态时,系统温度不变,因此dEk/dt=0iiifmamv?0iiiva ?2i iiiifvmv?宏观性质的统计?系统的势能?系统的内能?系统的总能E = Ep+Ek?系统的温度1()pijij NEV r? ? ?22kiipEm?21iiiBTmvdNk?模拟? 热容定义热容计算系统在温度T和T+?T时的总能ET、ET +?T,VvTEC?E
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 清华 分子 动力学 模拟 入门 课件
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内