清华分子动力学模拟入门.pptVIP

分子间势能及相互作用 一些气体的参数 kB=1.38x10?23(J/K): Boltzmann常数 Neon Argon Krypon Xenon Nitrogen ?(nm) 0.275 0.3405 0.360 0.410 0.370 ?/kB(K) 36 119.8 171 221 95 分子间势能及相互作用 相互作用 标量形式： 直角坐标： 至此，各粒子间相互作用已知，可进行模拟了 模拟的数学方法 Euler法和Euler-Cromer方法？ 不能用：不能保持总能量守恒 Verlet算法：速度形式 模拟的数学方法 Leap-frog算法： x的截断误差为 ，v的截断误差 模拟的数学方法 边界条件 模拟能力限制，不能模拟大量分子，只能模拟有限空间中的有限个分子：有限空间?边界 固体（刚性）边界条件 不仅仅有分子间的相互作 用，还引入了壁面的作用 分子量大时，壁面作用可 忽略不计 模拟的数学方法 取 ~ ， 前比值为0.2~0.01。取前值，模拟粗糙；取后值，模拟计算量太大 处理方法：使用周期性边界条件 周期性边界条件 ￥ ￥ ￥ ￥ 模拟的数学方法 两个不同粒子在x或y方向上的最大分离距离为a/2 最小像约定：两粒子分离距离最大分离距离，相互作用力可以忽略，而加入其中像