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

《分子力场及其在分子动力学模拟中的应用》由会员分享，可在线阅读，更多相关《分子力场及其在分子动力学模拟中的应用（9页珍藏版）》请在装配图网上搜索。

1、分子力场及其在分子动力学模拟中的应用分子模拟的概述自从20 世纪量子力学的快速发展后，几乎有关分子的一切性质，如结构、 构想、偶极矩、电离能、电子亲和力、电子密度等，皆可由量子力学计算获得。 计算与实验结果往往想当吻合，并且可由分析计算的结果得到一些实验无法获得 的资料，有助于对实际问题的了解。与实验相比较，利用计算机计算研究化学有下列几项优点：（1）成本降低； （2）增加安全性；（3）可研究极快速的反应或变化；（4）得到较佳的准确度； （5）增进对问题的了解。基于这些原因，分子的量子力学计算子 1970年后逐渐 受到重视。利用计算先行了解分子的特性，一成为合成化学家和药物设计学家所 依赖的重

2、要方法。化学家们借此可设计出最佳的反应途径，预测合成的可能性， 并评估所欲合成的分子的适用性，节省许多时间和避免材料的浪费。以欧美的许 多大型药厂为例，在采用计算以前，合成新药的成功率约为 17%-20%，但自从 1980 年后，由于在合成前先利用计算预测，其成功率已提高到 50%-60%。图一 为1955年美国化学会数据库（ACS database）所作的统计图。图中的纵坐标为引 用计算机计算程序报告所占的比例，横坐标为年份。由图中可清楚看出计算受重视的程度逐年增加。图一 美国化学会所发表的计算机计算在化学报告中的比例分子动力模拟（MD）,是时下最广泛为人米用的计算庞大复杂系统的方法。自 1

3、970 年起，由于分子力学的发展迅速，人们又系统地建立了许多适用于生化 分子体系、聚合物、金属与非金属材料的力场，使得计算复杂体系的结构与一些热力学与光谱性质的能力及准确性大为提升。分子动力模拟是应用这些力场及根据牛顿运动力学原理所发展的计算方法。此方法的优点为准确性高，可同时获得 系统的动态与热力学统计资料，并可广泛地使用与各种系统及各类特性的探讨。二 力场力场可以看作是势能面的经验表达式，是分子动力学模拟的基础。2.1 力场简述分子的总能量为动能与势能的和，分子的势能通常可表示为简单的几何坐标 的函数。例如，可将双原子分子 AB 的振动势能表示为 A 与 B 间键长的函数，即：2-1)式中

4、，k为弹性常数；r为键长；N为AB的平衡键长。这样以简单数学形式 表示的势能函数称为力场。复杂分子的总势能一般可以为各类型势能的和，这些 类型包括：总势能=非键结势能+键伸缩势能+键角弯曲势能+二面角扭曲势能+离平面振 动能+库伦静电势能势能项习惯以符号表示为：1；儿八 5 U, Ut，:1）非键结势能Unb一般在分子力场中，若 A、B 两原子属于同一分子但其间隔多余两个连接 的化学键（如A-C-C-B）,或两原子分属两个不同的分子，则原子对间存在非键结 作用力。2）键的伸缩项 分子中互相键结的原子形成的化学键。化学键的键长并非维持恒定，而是于其平衡值附近呈小幅度的震荡。描述此种作用的势能项称

5、为键伸缩项。3）键角弯曲项几 分子中连续键结的三原子形成键角。与键的伸缩一样，这些键角并非维持恒定不变，而是于其平衡位置附近呈小幅度的震荡。描述此种作用的势能为键角弯 曲项。4）二面角扭曲项九 分子中连续键结的四个原子，形成二面角。一般分子中的二面角较为松软，易于扭动。描述二面角扭转的势能项称为二面角扭曲项。5）离平面振动项， 分子中有些部分的原子有共平面的倾向，碳原子与氧原子的平衡位置位于共同平面。这些原子有形成共平面的倾向，通常共平面的四个原子的中心原子离平面小幅度振动。描述此种振动的势能项称为离平面振动项。6) 库伦作用项 离子或分子中的原子带有部分电荷，则这些带电荷的粒子间存在静电吸引

6、或 排斥作用。描述静电作用的势能项为昆仑作用项。图二综合显示分子中势能的各种作用项。图二综合显示分子中势能的各种作用项2.2力场作用项的一般式计算非键结作用，通常将原子视为位于其原子核坐标的一点。一般力场中最常见的非键结势能形式为lennard-jones (LJ)势能。此种势能又称为126势能,其数学式为U(r) = 4s (-)12 - (-)6rr式中，r为原子对间的距离；与o为势能参数，因原子的种类而异。键伸缩势能项之一般形式为简谐振动，即：U二1工k (r - ro)2b 2 b i ii式中，kb为键伸缩的弹力常数；ri及ri0分别表示第i个键的键长及其平衡键长。弹力常数愈大，振动

7、愈快，振动频率愈大。为了提高计算的精准性，有的 出了简谐振动项外，亦增加了非简谐振动项。键角弯曲项的一般形式为键角的简谐振荡：U 二 1 工 k （9 一90）2929 iii式中，k为键角弯曲的弹力常数；9 .及9 0分别表示第i个键角及其平衡键角的 ei i角度。二 面 角 扭 曲 项 的 一 般 形 式 为 ：U =工V (1 + cos ) + V (1 一 cos 2) + V (1 + cos 3)e 2 i23i式中，1、V2及3为二面角扭曲项的弹力常数，为二面角的角度。平面振动项的一般形式为：式中， 为离平面振动项的弹力常数； 定义为离平面振动的角度 库伦作用项的形式为:Uel

8、工i,jqqi_jDrij式中，qi及qj为分子中第i个离子与第j个离子所带的电荷；rij为距离；D为有 效介电常数。三分子动力模拟计算的原理分子动力模拟计算，简称为 MD 计算，自 1966年起发展至今约四十年。应用 十分广泛，随着计算机的迅速发展，越来越多地被应用于各种体系的模拟。目前 已有很多成熟的商用软件，使用方便，图形清晰，功能较强。但是，计算的成败 很大程度上取决于力场的适用性、计算速度的快慢、计算方法的争取性以及起始 结构的合理性。分子动力学的基本原理，即为利用牛顿运动定律。先由分子中各分子位置计 算系统的势能，再由(1)(2)式计算系统中各原子所受的力及加速度，然后在式（3）中

9、令t= 6 t, 6 t表示一非常短的时间间隔，则可得到经过6 t后各 分子的位置及速度重复以上步骤，由新的位置计算分子的势能，计算各原子所受 的力及加速度，预测再经过6 t后各分子的位置及速度如此反复循环，可得 到各时间下系统中分子运动的位置，速度及加速度资料。四.热力学特性的计算通过分子模拟获得系统中粒子与速度的轨迹后，通过分析这些轨迹可以获得各种 的性质和参数（1）计算系统动能的平均值可以得到系统的温度：式中，K为动能，T为温度，f为系统的自由度。（2）若计算的（N, V，E）系统中含有N个原子，且原子间的作用为成对加成, 即系统的势能为：U =工为Uij i,j j i3）系统的压力为

10、：F C) rij ij ijP为系统的密度，F.为原子i与j之间的作用力，f.为其距离向量。 ijij（4）动能的统计涨落定义为：6）=心）-K；：2由此物理量，可计算系统的定容比热容 CV：K2K ：-23kB2CV3k2T 2B -2N四统计力学与模拟计算的关系系综各种特性的平均值称为系综平均值。若特性 x 为连续的值，则/ + oo= I xw(x)dxoo式中，积分的范围为所有可能的X值，w(x)dx称为x的概率配分函数。w(x )dx的意义为于范围*I h间x的概率。一个系统的特性x亦可能与时间相关，x(t)。由时间所得的平均值为：x= limToo式中，积分所表示所有时间的范围。

11、依据各态历经假设，系综的平均值与长时间 的平均值应相等，即：X 依此假说，平均值可由许多系统的特性计算，或是计算单独一个系统长时间 的平均值。模拟计算系统的特性所依据的即是统计学的原理。参考文献：1陈正隆等，分子模拟的理论与实践，2007仅供个人用于学习、研究；不得用于商业用途。For personal use only in study and research; not for commercial use.Nur fur den persdnlichen fur Studien, Forschung, zu kommerziellen Zwecken verwendet werden.Pour l etude et la recherche uniquement a des fins personnelles; pas a des fins commerciales.TonbKO gn只 nwgen, KOTopbie ucnonb3yQTC只 gn只 o6yqeHU只，uccnegoBaHU叼 u He gonxHN ucnonb3OBaTbc只 b KOMMepqecKux qen只x.以下无正文