分子动力学模拟-精品文档.docx
分子动力学模拟分子动力学模拟分子动力学就是一门结合物理,数学与化学的综合技术。分子动力学就是一套分子模拟方法,该方法主要就是依靠牛顿力学来模拟分子体系的运动,以在由分子体系的不同状态构成的系统中抽取样本,进而计算体系的构型积分,并以构型积分的结果为基础进一步计算体系的热力学量与其她宏观性质。这门技术的发展进程就是:1980年:恒压条件下的动力学方法(Andersen方法、Parrinello-Rahman法)1983年:非平衡态动力学方法(GillanandDixon)1984年:恒温条件下的动力学方法(能势方法)1985年:第一原理分子动力学法(?法)1991年:巨正则系综的分子动力学方法(CaginandPettit)、最新的巨正则系综,即为组成系综的系统与一温度为T、化学势为的很大的热源、粒子源相接触,此时系统不仅同热源有能量交换,而且能够同粒子源有粒子的交换,最后到达平衡,这种系综称巨正则系综。进行分子动力学模拟的第一步就是确定起始构型,一个能量较低的起始构型就是进行分子模拟的基础,一般分子的其实构型主要就是来自实验数据或量子化学计算。在确定起始构型之后要赋予构成分子的各个原子速度,这一速度就是根据玻尔兹曼分布随机生成,由于速度的分布符合玻尔兹曼统计,因而在这个阶段,体系的温度就是恒定的。另外,在随机生成各个原子的运动速度之后须进行调整,使得体系总体在各个方向上的动量之与为零,即保证体系没有平动位移。由上一步确定的分子组建平衡相,在构建平衡相的时候会对构型、温度等参数加以监控。进入生产相之后体系中的分子与分子中的原子开场根据初始速度运动,能够想象其间会发生吸引、排挤乃至碰撞,这时就根据牛顿力学与预先给定的粒子间互相作用势来对各个例子的运动轨迹进行计算,在这个经过中,体系总能量不变,但分子内部势能与动能不断互相转化,进而体系的温度也不断变化,在整个经过中,体系会遍历势能面上的各个点,计算的样本正就是在这个经过中抽取的。用抽样所得体系的各个状态计算当时体系的势能,进而计算构型积分。作用势的选择与动力学计算的关系极为密切,选择不同的作用势,体系的势能面会有不同的形状,动力学计算所得的分子运动与分子内部运动的轨迹也会不同,进而影响到抽样的结果与抽样结果的势能计算,在计算宏观体积与微观成分关系的时候主要采用刚球模型的二体势,计算系统能量,熵等关系时早期多采用Lennard-Jones、morse势等双体势模型,对于金属计算,主要采用morse势,但就是由于通过实验拟合的对势容易导致柯西关系,与实验不符,因而在后来的模拟中有人提出采用EAM等多体势模型,或者采用第一性原理计算结果通过一定的物理方法来拟合二体势函数。但就是对于二体势模型,多体势往往缺乏明确的表达式,参量很多,模拟收敛速度很慢,给应用带来很大困难,因而在一般应用中,通过第一性原理计算结果拟合势函数的L-J,morse等势模型的应用仍非常广泛。分子动力学计算的基本思想就是赋予分子体系初始运动状态之后,利用分子的自然运动在相空间中抽取样本进行统计计算,时间步长就就是抽样的间隔,因此时间步长的选取对动力学模拟非常重要。太长的时间步长会造成分子间的剧烈碰撞,体系数据溢出;太短的时间步长会降低模拟经过搜索相空间的能力,因而一般选取的时间步长为体系各个自由度中最短运动周期的特别之一。但就是通常情况下,体系各自由度中运动周期最短的就是各个化学键的振动,而这种运动对计算某些宏观性质并不产生影响,因而就产生了屏蔽分子内部振动或其她无关运动的约束动力学,约束动力学能够有效地增长分子动力学模拟时间步长,提高搜索相空间的能力。下面就是做模拟的一般性步骤,详细的步骤与经过依靠于确定的系统或者软件,但这不影响我们把它当做一个入门指南:1首先我们需要对我们所要模拟的系统做一个简单的评估,三个问题就是我们必需要明确的:做什么(whattodo)为什么做(whytodo)怎么做(howtodo)2选择适宜的模拟工具,大前提就是它能够实现您所感兴趣的目的,这需要您非常慎重的查阅文献,瞧瞧别人用这个工具都做了些什么,有没有与您相关的,千万不要做到一半才发现原来这个工具根本就不能实现您所感兴趣的idea,切记！考虑1:软件的选择,这通常与软件主流使用的力场有关,而软件本身就具有一定的偏向性,比方讲,做蛋白体系,Gromacs,Amber,Namd均可;做DNA,RNA体系,首选肯定就是Amber;做界面体系,DI_POLY比拟强大,另外做材料体系,Lammps会就是一个不错的选择。考虑2:力场的选择。力场就是来描绘体系中最小单元间的互相作用的,就是用量化等方法计算拟合后生成的经历式,有人会嫌它粗糙,但就是它确确实实给我们模拟大系统提供了可能,只能讲关注的切入点不同罢了。常见的有三类力场:全原子力场,联合力场,粗粒化力场;当然还有所谓的第一代,第二代,第三代力场的讲法,这里就不逐一列举了。再次提醒注意:必须选择合适我们所关注体系与我们所感兴趣的性质及现象的力场。3通过实验数据或者就是某些工具得到体系内的每一个分子的初始构造坐标文件,之后,我们需要按我们的想法把这些分子根据一定的规则或就是随机的排列在一起,进而得到整个系统的初始结果,这也就是我们模拟的输入文件。4构造输入文件得到了,我们还需要力场参数输入文件,也就就是针对我们系统的力场文件,这通常就是由所选用的力场决定,比方键参数与非键参数等势能函数的输入参数。5体系的大小通常由您所选用的box大小决定,我们必须对可行性与合理性做出评估,进而确定体系的大小,这依靠于详细的体系,这里不细讲了。6由于初始构象可能会存在两个原子挨的太近的情况(称之为badconnect),所以需要在正式模拟开场的第一步进行体系能量最小化,比拟常用的能量最小化有两种,最速下降法与共轭梯度法,最速下降法就是快速移除体系内应力的好方法,但就是接近能量极小点时收敛比拟慢,而共轭梯度法在能量极小点附近收敛相对效率高一些,所以我们一般做能量最小化都就是在最速下降法优化完之后再用共轭梯度法优化,这样做能有效的保证后续模拟的进行。7以平衡态模拟为例,您需要设置适当的模拟参数,并且保证这些参数设置与力场的产生相一致,举个简单的例子,gromos力场就是用的范德华势双截断来定范德华参数的,若您也用gromos力场的话也应该用双截断来处理范德华互相作用。常见的模拟思路就是,先在NVT下约束住您的溶质(剂)做限制性模拟,这就是一个升温的经过,当温度到达您的设定后,接着做NPT模拟,此经过将调整体系的压强进而使体系密度收敛。经过一段时间的平衡模拟,在确定系统弛豫已经完全消除之后,就能够开场取数据了,怎样判定体系到达平衡,这个问题就是比拟技术性的问题,简单的讲能够通过下面几种方式,一,瞧能量(势能,动能与总能)就能否收敛;二,瞧系统的压强,密度等等就能否收敛;三瞧系统的RMSD就能否到达您能接受的范围,等等。8运行足够长时间的模拟以确定我们所感兴趣的现象或就是性质能够被观测到,并且务必确保此现象的可重复性。9数据拿到手后,很容易通过一些可视化软件得到轨迹动画,但这并不能拿来发文章。真正的工作才刚刚开场分析数据,您所感兴趣的现象或性质只就是外表,隐含在它们之中的机理才就是文章中的主题。分子动力学能够用于NPT,NVE,NVT等系综的计算,就是一种基于牛顿力学确定论的热力学计算方法,与蒙特卡洛法相比在宏观性质计算上具有更高的准确度与有效性,能够广泛应用于物理,化学,生物,材料,医学等各个领域。另外,在实际应用中,经常把分子动力学方法与蒙特卡罗法联合使用。
