浅谈分子动力学模拟.PPTVIP

MD模拟参数及注意事项 分子动力学模拟参数及注意事项 1.周期性边界条件 2.势函数 3.系综 4.积分步长的选取 5.注意事项 含义：只模拟实际物质中很小一部分，使模拟体系成为无限的具有相同性质的分子体系的，该部分在三维空间中周期性的存在，代表整个体系。 1.周期性边界条件 * * 发展与整个20世纪 * * * * * 浅谈分子动力学模拟 作者：罗万力 专业：道路与铁道工程 指导老师：孔令云 内容 1.MD模拟的发展历程 2.分子动力学模拟方法 3.MD模拟参数及注意事项 4.实例展示 MD模拟的发展历程 经典力学 经典力学的基本定律是牛顿运动定律或与牛顿定律有关且等价的其他力学原理，它是20世纪以前的力学。 两个基本假定 1.假定时间和空间是绝对的，长度和时间间隔的测量与观测者的运动无关，物质间相互作用的传递是瞬时到达的。 2.一切可观测的物理量在原则上可以无限精确地加以测定。 根据经典力学的基本假定可以看出经典力学有何局限性？ 量子力学 量子力学理论形成于20世纪初期，彻底改变了人们对物质组成成分的认识。微观世界里，粒子是不断跳跃的?概率云，它们不只存在一个位置，也不会从点A通过一条单一路径到达点B。 统计力学（统计物理学）是研究原子、分子集合的宏观运动规律的科学，根据对物质微观结构及微观粒子相互作用的认识，用概率统计的方法，对由大量粒子组成的宏观物体的 物理性质及 宏观规律作出微观解释的理论物理学分支。 统计力学 什么分子动力学模拟？ 分子动力学模拟（molecular dynamics）是近年来飞速发展的一种分子模拟方法，它以经典力学、量子力学、统计力学为基础，利用计算机数值求解分子体系经典力学运动方程的方法得到体系的相轨迹，并统计体系的结构特征与性质。 分子动力学简史 1957年：基于刚球势的分子動力学法（Alder and Wainwright） 1964年：利用Lennard-Jone势函数法对液态氩性质的模拟（Rahman） 1971年：模拟具有分子团簇行为的水的性质（Rahman and Stillinger） 1977年：约束动力学方法（Rychaert, Ciccotti Berendsen; van Gunsteren） 1980年：恒压条件下的动力学方法（Andersen法、Parrinello-Rahman法） 1983年：非平衡态动力学方法（Gillan and Dixon） 1984年: 恒温条件下的动力学方法(Berendsen et al.) 1984年：恒温条件下的动力学方法（Nosé-Hoover法） 1985年：第一原理分子動力学法（→Car-Parrinello法） 1991年：巨正则系综的分子动力学方法（Cagin and Pettit） 分子动力学模拟方法 分子动力学方法特征： 1.分子动力学是在原子、分子水平上求解多体问题的重要的计算机模拟方法，可以预测纳米尺度上的材料动力学特性。 2.通过求解所有粒子的运动方程，分子动力学方法可以用于模拟与原子运动路径相关的基本过程。 3.在分子动力学中，粒子的运动行为是通过经典的Newton运动方程所描述。 4.分子动力学方法是确定性方法，一旦初始构型和速度确定了，分子随时间所产生的运动轨迹也就确定了。 MD的理论依据是什么？ 每个粒子的运动都取决于牛顿定律 （F=ma） 分子动力学方法基本原理： 初始条件： 分子动力学----牛顿运动方程式的数值解法 1.解析法（牛顿第二定律） 由于分子体系相互作用非常复杂，难以用解析法求解分子体系运动方程，通常只能采用差分法求解分子体系运动方程的近似解。 2.有限差分法 Verlet算法 优点： 1、精确，坐标精度误差四阶O(Δ4) 2、每次积分只计算一次力 3、时间可逆 缺点： 1、速度有较大误差O(Δ2) 2、轨迹与速度无关，无法与热浴耦联 分子动力学----牛顿运动方程式的数值解法 蛙跳算法 ?. 首先利用当前时刻的加速度，计算半个时间步长后的速度： ?. 计算下一步长时刻的位置： 分子动力学----牛顿运动方程式的数值解法 ?. 计算当前时刻的速度： V a0 v r v r 分子动力学----牛顿运动方程式的数值解法 优点： 1、提高精确度 2、轨迹与速度有关，可与热浴耦联 缺点： 1、速度近似 2、比Verlet算子多花时间 分子动力学----牛顿运动方程式的数值解法 * * 发展与整个20世纪 * * * * *