中国矿业大学分子模拟课件第六章.pdfVIP

分 子 模 拟 牛继南 njn0516@cumt.edu.cn 2011.3 第六章 经典力学和统计力学回顾 6.1 牛顿运动力学 • 按照牛顿运动定律，粒子所受的力等于此 粒子的质量乘以加速度，即： • 如果以动量表示则为： 牛顿第二定律 • 如果质量和时间无关，则牛顿第二定律可 以写为： r dr r r r r dp d(mv) dv dt r = = m = m = mr= F= ma dt dt dt dt • 式中 为粒子的位置矢量。若 为位置的函数 则上式为一个二次微 分方程。如果给出粒子的初始位置，则可 由此方程解得粒子的位置与时间关系： x(t), y(t), z(t). • 以简谐振子为例： • 设x0为谐振子的自然长度，则根据胡克定律 ，弹簧所系质点所受的力为： • 按照微分形式可改写为： • 如果令 则上式变为： • 此微分方程的解为： • 式中 为振子的振动频率。 • 可将该式化为： • ϕ即我们通常所说的相角。 • 由此可以看出弹簧振子相对于平衡位置的 位移随着时间按三角函数关系变化，即作 简谐振动。 • 以在二维平面内运动的粒子为例，粒子与 原点间为库伦引力： • K为库伦常数，r为原点和粒子间的距离， 即 • 可得： • 此二式中，含x的微分方程中含有变量y， 含y的微分方程中带有变量x，它们之间相 互耦合，不易求解，所以必须将变量转换 除去耦合，才能求解。 • 通常采用极坐标： 1 mx= K cosθ 2 • 所以得到r，θ的运动方程式： r 1 my= K sinθ