《分子动力学方法》ppt课件

《《分子动力学方法》ppt课件》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、分子动力学方法MolecularDynamicsSimulationQing-YuZhangStateKeyLaboratoryforMaterialsModificationbyLaser,IonandElectronBeams微观尺度材料设计分子动力学分子动力学是在原子、分子水平上求解多体问题的重要的计算机模拟方法，可以预测纳米尺度上的材料动力学特性。通过求解所有粒子的运动方程，分子动力学方法可以用于模拟与原子运动路径相关的基本过程。在分子动力学中，粒子的运动行为是通过经典的Newton运动方程所描述。微观尺度材料设计分子动力学原子体系的运动方程Lagrangian方程

2、定义Lagrangian函数为L=K-V则运动的Lagrangian方程为微观尺度材料设计分子动力学原子体系的运动方程Lagrangian方程微观尺度材料设计分子动力学原子体系的运动方程微观尺度材料设计分子动力学有限差分方法-预测校正法预测校正法是分子动力学模拟中的常用算法之一，其基本思想是Taylor展开：微观尺度材料设计分子动力学有限差分方法-预测校正法根据新的原子位子rp，通过运动方程可以获得校正后的ac(t+t)。定义预测误差微观尺度材料设计分子动力学有限差分方法-Gear预测校正因子定义一组矢量：微观尺度材料设计分子动力学有限差分方法-Gear预

3、测校正因子微观尺度材料设计分子动力学有限差分方法-Gear预测校正因子微观尺度材料设计分子动力学有限差分方法-Gear预测校正因子对于一阶运动方程Valuesc0c1c2c3c4c535/1211/243/813/41/65251/720111/121/31/24695/288125/2435/725/481/120微观尺度材料设计分子动力学有限差分方法-Gear预测校正因子对于二阶运动方程Valuesc0c1c2c3c4c5301141/65/611/3519/1203/411/21/1263/20251/360111/181/61/60微观尺度材料设计分子动力学有限

4、差分方法-Gear预测校正因子对于二阶运动方程Valuesc0c1c2c3c4c5301141/65/611/3519/903/411/21/1263/16251/360111/181/61/60微观尺度材料设计分子动力学有限差分方法-Verlet算法微观尺度材料设计分子动力学有限差分方法-速度Verlet算法微观尺度材料设计分子动力学有限差分方法-leap-frog算法微观尺度材料设计分子动力学有限差分方法-Verlet算法比较微观尺度材料设计分子动力学力场的截断在分子动力学中，出于计算上的考虑，力场的截断是必须的，即在某一范围内力场是有效的，因此会导致一些计算上

5、的困难。势函数直接截断：微观尺度材料设计分子动力学力场的截断力场连续的势函数截断：微观尺度材料设计分子动力学力场的截断微观尺度材料设计分子动力学近邻表-Verlet近邻表RListUpdateintervalTimeN=256N=500No2.602.702.903.103.433.50-5.7812.5026.3243.4883.33100.003.332.242.172.282.472.89-10.004.934.554.514.79-5.86微观尺度材料设计分子动力学近邻表-网格近邻表l=L/M>rcNc=N/Mdftime=3dNNc/0.5N(N-1)微观尺度

6、材料设计分子动力学近邻表-网格近邻表微观尺度材料设计分子动力学等温分子动力学-随机方法采用随机数重新标定系统的速度动量空间新旧状态的几率比为微观尺度材料设计分子动力学等温分子动力学-扩展系统方法增加一个热池，允许热流在系统和热池之间交换。微观尺度材料设计分子动力学等温分子动力学-约束方法利用速度调整因子约束系统的温度。微观尺度材料设计分子动力学等温分子动力学-约束方法微观尺度材料设计分子动力学等温分子动力学-约束方法微观尺度材料设计分子动力学等压分子动力学-扩展系统方法增加一个“活塞”，“活塞”具有质量Q，体积V。微观尺度材料设计分子动力学等压分子动力学

7、-约束方法微观尺度材料设计分子动力学等压分子动力学-约束方法微观尺度材料设计分子动力学分子动力学演示微观尺度材料设计分子动力学分子动力学的应用-薄膜生长Evolutionofsubstratetemperaturesatvariousincidentenergy微观尺度材料设计分子动力学分子动力学的应用-薄膜生长0.1eV10.0eVMorphologiesofAu/Au(001)filmsatvariousincidenten