MD原理与Lammps实例.ppt
* Laboratory of Space Materials Science and Technology 分子动力学模拟与 LAMMPS 软件的使用 李留辉 2016.05 分子动力学模拟基础 MD势函数 LAMMPS软件简介 LAMMPS运行实例 主要内容 MD:根据描述原子间相互作用的势函数, 求解经典牛顿运动方程,得到系统确定的运动轨迹 1. 分子动力学模拟基础 原子间作用势 MD计算流程图 Newton 运动方程的数值求解:有限差分 1、Verlet、 Velocity-Verlet、Leap-Frog算法 2、Gear预测-校正算法 微正则系综 (NVE) —— 孤立体系 正则系综 (NVT) —— 温度恒定 等温等压系综 (NPT) —— 温度和压力恒定 MD中的系综 (Ensemble) 温度标定: ~ MD中的结构表征 1、双体分布函数 (Pair Distribution Function) 局域原子数密度与平均数密度之比 配位数 3、键对分析技术 (Pair Analysis) 区分不同的团簇结构 2、结构因子S(q) q:倒易空间矢量 原子扩散 表面、界面、缺陷 晶体生长、 固态相变 熔化、凝固、非晶形成 力学性质 声子、热力学性质 密度、膨胀系数 焓、比热 表面张力、粘度、扩散系数 熔点、熔化热 MD的应用 热物理性质 2. MD 势函数 对势:Lennard-Jones、Morse 多体势[1] EAM、Finnis-Sinclair, Extended FS SMA-TB、Smoothed SMA-TB Johnson、Mishin、Zhou …… 考虑角度效应 MEAM / 2NN-MEAM Stillinger-Weber、Tersoff Bond-Order potential 对势项 多体项 电子密度: [1] J. H. Li, et al, Physics Reports, 2008, 455:1-134 Lennard-Jones势: e:势阱深度;