分子力场及其在分子动力学模拟中的应用.docx

分子力场及其在分子动力学模拟中的应用一分子模拟的概述自从20世纪量子力学的快速进展后，几乎有关分子的一切性质，如结构、构想、偶极矩、电离能、电子亲和力、电子密度等，皆可由量子力学计算获得。计算与试验结果往往想当吻合，并且可由分析计算的结果得到一些试验无法获得的资料，有助于对实际问题的了解。与试验相比较，采用计算机计算讨论化学有下列几项优点：（1）成本降低;（2）增加平安性；（3）可讨论极快速的反应或变化；（4）得到较佳的精确度；（5）增进对问题的了解。基于这些缘由，分子的量子力学计算子1970年后渐渐受到重视。采用计算先行了解分子的特性，一成为合成化学家和药物设计学家所依靠的重要方法。化学家们借此可设计出最佳的反应途径，猜测合成的可能性，并评估所欲合成的分子的适用性，节约很多时间和避开材料的铺张。以欧美的很多大型药厂为例，在采纳计算以前，合成新药的胜利率约为17%-20%,但自从1980年后，由于在合成前先采用计算猜测，其胜利率已提高到50%-60%°图一为1955年美国化学会数据库（ACSdatabase）所作的统计图。图中的纵坐标为引用计算机计算程序报告所占的比例，横坐标为年份。由图中可清晰看出计算受重视的程度逐年增加。分子动力模拟（MD）,是时下最广泛为人采纳的计算浩大简单系统的方法。自1970年起，由于分子力学的进展快速，人们又系统地建立了很多适用于生化分子体系、聚合物、金属与非金属材料的力场，使得计算简单体系的结构与一些热力学与光谱性质的力量及精确性大为提升。分子动力模拟是应用这些力场及依据牛顿运动力学原理所进展的计算方法。此方法的优点为精确性高，可同时获得系统的动态与热力学统计资料，并可广泛地使用与各种系统及各类特性的研讨。二力场力场可以看作是势能面的阅历表达式，是分子动力学模拟的基础。2.1力场简述分子的总能量为动能与势能的和，分子的势能通常可表示为简洁的几何坐标的函数。例如，可将双原子分子AB的振动势能表示为A与B间键长的函数，即：U（r）=2（r-r0）2（2-1）式中，k为弹性常数；r为键长；r°为AB的平衡键长。这样以简洁数学形式表示的势能函数称为力场。简单分子的总势能一般可以为各类型势能的和，这些类型包括：总势能:非键结势能+键伸缩势能+键角弯曲势能十二面角扭曲势能+离平面振动能+库伦静电势能势能项习惯以符号表示为：U=Unb+%+%+%+%+Uel1）非键结势能Unb一般在分子力场中，若A、B两原子属于同一分子但其间隔多余两个连接的化学键（如A-C-C-B）,或两原子分属两个不同的分子，则原子对间存在非键结作用力。2）键的伸缩项/分子中相互键结的原子形成的化学键。化学键的键长并非维持恒定，而是于其平衡值四周呈小幅度的震荡。描述此种作用的势能项称为键伸缩项。3）键角弯曲项为分子中连续键结的三原子形成键角。与键的伸缩一样，这些键角并非维持恒定不变，而是于其平衡位置四周呈小幅度的震荡。描述此种作用的势能为键角弯曲项。4）二面角扭曲项分子中连续键结的四个原子，形成二面角。一般分子中的二面角较为松软，易于扭动。描述二面角扭转的势能项称为二面角扭曲项。5）离平面振动项力分子中有些部分的原子有共平面的倾向，碳原子与氧原子的平衡位置位于共同平面。这些原子有形成共平面的倾向，通常共平面的四个原子的中心原子离平面小幅度振动。描述此种振动的势能项称为离平面振动项。6)库伦作用项”,离子或分子中的原子带有部分电荷，则这些带电荷的粒子间存在静电吸引或排斥作用。描述静电作用的势能项为昆仑作用项。图二综合显示分子中势能的各种作用项。Ub图二综合显示分子中势能的各种作用项2.2力场作用项的一般式计算非键结作用，通常将原子视为位于其原子核坐标的一点。一般力场中最常见的非键结势能形式为lennard-jones(LJ)势能。此种势能又称为126势能,其数学式为U(r)=4疝巴-(Wrr式中，r为原子对间的距离；£与。为势能参数，因原子的种类而异。键伸缩势能项之一般形式为简谐振动，即：Ub=-r；)22i式中，尤为键伸缩的弹力常数；心及分别表示第i个键的键长及其平衡键长。弹力常数愈大，振动愈快，振动频率愈大。为了提高计算的精准性，有的出了简谐振动项外，亦增加了非简谐振动项。键角弯曲项的一般形式为键角的简谐振荡：4=：(4-£)22i式中，的为键角弯曲的弹力常数；。.及小分别表示第i个键角及其平衡键角的角度。二面角扭曲项的一般形式为：%=1Zq(l+cos¢)+/(1-COS20)+匕(1+COS30)2i式中，、及匕为二面角扭曲项的弹力常数，。为二面角的角度。平面振动项的一般形式为：=g)式中，内为离平面振动项的弹力常数；Z定义为离平面振动的角度。库伦作用项的形式为：式中，及为分子中第i个离子与第j个离子所带的电荷；为距离；D为有效介电常数。三.分子动力模拟计算的原理分子动力模拟计算，简称为MD计算，自1966年起进展至今约四十年。应用特别广泛，随着计算机的快速进展，越来越多地被应用于各种体系的模拟。目前已有很多成熟的商用软件，使用便利，图形清晰，功能较强。但是，计算的成败很大程度上取决于力场的适用性、计算速度的快慢、计算方法的争取性以及起始结构的合理性。分子动力学的基本原理，即为采用牛顿运动定律。先由分子中各分子位置计算系统的势能，再由(1)(2)式计算系统中各原子所受的力及加速度，然后在式(3)中令t=t,t表示一特别短的时间间隔，则可得到经过t后各分子的位置及速度重复以上步骤，由新的位置计算分子的势能，计算各原子所受的力及加速度，猜测再经过t后各分子的位置及速度如此反复循环，可得到各时间下系统中分子运动的位置，速度及加速度资料。四.热力学特性的计算通过分子模拟获得系统中粒子与速度的轨迹后，通过分析这些轨迹可以获得各种的性质和参数(1)计算系统动能的平均值可以得到系统的温度：T一处1f式中，K为动能，T为温度，f为系统的自由度。(2)若计算的(N,V,E)系统中含有N个原子，且原子间的作用为成对加成，即系统的势能为：fN=UijijJ>i(3)系统的压力为：P=PklST/净后P为系统的密度，吃为原子i与j之间的作用力，亏为其距离向量。(4)动能的统计涨落定义为：附=卜、-阴由此物理量，可计算系统的定容比热容CV：23413%)2N126J四.统计力学与模拟计算的关系系综各种特性的平均值称为系综平均值。若特性X为连续的值，贝U：<x>=Ixw(x)dxJ-co式中，积分的范围为全部可能的X值，W(x)dx称为X的概率配分函数。W(x)dx的意义为于范围XTX+dx间X的概率。一个系统的特性X亦可能与时间相关，x(t)o由时间所得的平均值为：X=Iimx(t)dt式中，积分所表示全部时间的范围。依据各态历经假设，系综的平均值与长时间的平均值应相等，即：X=<X>依此假说，平均值可由很多系统的特性计算，或是计算单独一个系统长时间的平均值。模拟计算系统的特性所依据的即是统计学的原理。参考文献：1陈正隆等，分子模拟的理论与实践，2007