多体系统—分子动力学Multi-bodySystemMolecularDynamics引言物质基本构成—分子、原子•在分子、原子这个微观水平上来考察物质:多体世界•查清楚微观世界,宏观就清楚了从微观考虑问题的现实可行性从微观考虑问题的必要性•物性的观测性参数:热传导、温度、压力、粘性、...…微观处理的前提•已知微观粒子间的相互作用假设•分子为球,惰性,分子间的作用只取决于分子间的距离分子动力学(MolecularDynamics,MD)MD的应用领域:物理、化学、生物、材料等MD方法能实时将分子的动态行为显示到计算机屏幕上,便于直观了解体系在一定条件下的演变过程MD含温度与时间,因此还可得到如材料的玻璃化转变温度、热容、晶体结晶过程、输送过程、膨胀过程、动态弛豫(relax)以及体系在外场作用下的变化过程等水和离子在微小硅孔中的运动聚乙烯的结晶MD的基本原理用牛顿经典力学计算许多分子在相空间中的轨迹•求解系统中的分子或原子间作用势能和系统外加约束共同作用的分子或原子的牛顿方程。•模拟系统随时间推进的微观过程。•通过统计方法得到系统的平衡参数或输运性质•计算程序较为复杂,占用较多内存MD的主要步骤选取要研究的系统及其边界,选取系统内粒子间的作用势能模型设定系统中粒子的初始位置和初始动量建立模拟算法,计算粒子间作用力及各粒子的速度和位置当体系达到平衡后,依据相关的统计公式,获得各宏观参数和输运性质分子间势能及相互作用N个粒子系统的总势能NjiijNNNNrVrVrVrVrVrVrVV11,-22311312)()()()()()()(dVrVr刚球模型()0rdVrrd()()dVrr斥力力心点模型drV-()()rdVrdrdrSoutherland模型分子间势能及相互作用Lennard-Jones势能能量尺度;长度尺度为方便,时常归一化:记;6124)(rrrV612/1/14)(rrrVVV/rr/612114)(rrrVV(r)F(r)排斥力V(r),F(r)r吸引力分子间势能及相互作用一些气体的参数kB=1.38x1023(J/K):Boltzmann常数NeonArgonKryponXenonNitrogen(nm)0.2750.34050.3600.4100.370/kB(K)36119.817122195分子间势能及相互作用相互作用•标量形式:•直角坐标:至此,各粒子间相互作用已知,可进行模拟了)(rVrmf()(),yxffVrxVrymrrmrr)(rVmF模拟的数学方法Euler法和Euler-Cromer方法?•不能用:不能保持总能量守恒•Verlet算法:速度形式12121212111212112nnnnnnnnnnnnnnxxvtatvvaatvvatvvat模拟的数学方法Leap-frog算法:x的截断误差为,v的截断误差1/21/211/21/22nnnnnntnnnvvtaxxtvvva4()Ot2()Ot模拟的数学方法边界条件•模拟能力限制,不能模拟大量分子,只能模拟有限空间中的有限个分子:有限空间边界•固体(刚性)边界条件不仅仅有分子间的相互作用,还引入了壁面的作用分子量大时,壁面作用可忽略不计233611aaaN和壁面作用分子数壁面积总分子数体积3aN总分子数壁面积和壁面作用分子数模拟的数学方法取~,前比值为0.2~0.01。取前值,模拟粗糙;取后值,模拟计算量太大处理方法:使用周期性边界条件•周期性边界条件210N610a模拟的数学方法两个不同粒子在x或y方向上的最大分离距离为a/2•最小像约定:两粒子分离距离>最大分离距离,相互作用力可以忽略,而加入其中像粒子之一相互作用力来考虑?2a2/aa模拟的数学方法考虑的粒子总数不变初始条件•随机初始条件给法之一要求大小:条件一:规则给法条件二:随机给法max||Vv)1random2(maxVv模拟的数学方法random:随机数产生函数,产生(0,1)之间的随机数。方向(按球坐标给法):分量:xyz))5.0random(2arccos()]5.0random(2[arccos)0.5random(sign180180;1800coss...