第四章 分子动力学模拟方法.pptVIP

练习： 推导LJ流体分子间力的表达式(fx, fy, fz及其对比量)： 势能函数形式： 力： 采用对比量： ?=x, y, z LJ分子间的维里项： ?=x, y, z 采用对比量的运动方程形式：(以蛙跳(Leap-frog)算法为例) 采用对比量： 最终得到： 同理得到： 速度的标度(Velocity Scaling)： 根据能量均分原理，可知： 标度因子: 对比量 速度标度: 或 微正则系综MD模拟程序F3讲解（LJ, NVE）： 初始化： READ (*,‘(A)’) TITLE ! 运行作业题目 READ (*,*) NSTEP ! 运行步数 READ (*,*) IPRINT ! 打印步数 READ (*,‘(A)’) CNFILE ! 位型文件 READ (*,*) DENS ! 对比密度 READ (*,*) RTEMP ! 对比温度 READ (*,*) RCUT ! 对比截断半径 READ (*,*) DT ! 对比时间步长 CALL READCN ( CNFILE ) 初始位型： 面心立方 (face-centered cubic, FCC)： 每面中心有一格点 体心立方 (body-c