《过渡态、反应路径的计算方法及相关问题》由会员分享，可在线阅读，更多相关《过渡态、反应路径的计算方法及相关问题（21页珍藏版）》请在金锄头文库上搜索。
1、过渡态、反应路径的计算方法及相关问题SoberevaDepartment of Chemistry, University of Science and Technology Beijing, Beijing100083, China前言：本文主要介绍过渡态、反应路径的计算方法，并讨论相关问题。由于这类算法极多， 可以互相组合，限于精力不可能面面俱到展开，所以只介绍常用，或者实用价值有限但有启 发性的方法。文中图片来自相关文献，做了一定修改。由于本文作为帖子发布，文中无法插 入复杂公式，故文中尽量将公式转化为文字描述并加以解释，这样必然不如公式形式严谨， 而且过于复杂的公式只能略过，但我想这样
2、做的好处是更易把握方法的梗概，有兴趣可以进 一步阅读原文了解细节。对于Gaussian中可以实现的方法，文中对其在Gaussian中的使用 进行了一些讨论，希望能纠正一些网上流传的误区。虽然绝大多数人不专门研究计算方法， 其中很多方法也不会用到，但多了解一下对开阔思路是很有好处的。文中指的“反应”包括构象变化、异构化、单分子反应等任何涉及到过渡态的变化过程。“反 应物，与“产物，泛指这些过程的初态和末态。“优化”若未注明，包括优化至极小点和优 化至过渡态。势能面是高维的，但为了直观以及表述方便，文中一般用二维势能面模型来讨 论，应推广至高维情况。限于纯文本格式，向量、矩阵无法加粗表示，但容易自
3、行判断。目录：1. 过渡态2. 过渡态搜索算法基于初猜结构的算法牛顿-拉弗森法(Newton-Raphson,NR)与准牛顿法(quasi-Newton,QN)AH 方法(augmented Hessian)RFO 法(Rational Function Optimization，有理函数优化)P-RFO 法(Partitioned-RFO)。入法(Quadratic Approximation， 二次逼近)TRIM 法(trust-region image minimization，置信区域镜像最小化) 在高斯中的常见问题GDIIS 法(Geometry Direct Inversion i
4、n the Iterative Subspace)梯度模优化(gradient norm minimization)Dimer方法基于反应物与产物结构的算法同步转变方法(synchronous transit,ST)STQN 方法(Combined Synchronous Transit and Quasi-Newton Methods)赝坐标法(pseudo reaction coordinate)DHS 方法(Dewar-Healy-Stewart，亦称 Saddle 方法)与 LTP 方法(Line-Then-Plane)Ridge方法Step-and-Slide 方法Muller-Br
5、own 方法CI-NEB、ANEBA 方法基于反应物结构的算法最缓上升法(least steep ascent,shallowest ascent)本征向量/本征值跟踪法(eigenvector/eigenvalue following,EF。也称 mode walking/mode following/Walking up valleys)ARTn(activation-relaxation technique nouveau)梯度极值法(Gradient extremal,GE)约化梯度跟踪(reduced gradient following,RGF)等势面搜索法(Isopotenial
6、 Searching)球形优化(Sphere optimization)全势能面扫描3. 过渡态相关问题无过渡态的反应途径(barrierless reaction pathways)Hammond-Leffler 假设对称性问题溶剂效应计算过渡态的建议流程4. 内禀反应坐标(intrinsic reaction coordinate,IRC)算法最陡下降法(Steepest descent)IMK 方法(Ishida-Morokuma-Kormornicki)Muller-Brown 方法GS(Gonzalez-Schlegel )方法方法Drag method 方法PEB 方法(plain
7、 elastic band)Elber-Karplus 方法SPW 方法(Self-Penalty Walk)LUP 方法(Locally Updated planes)NEB 方法(Nudged Elastic Band)DNEB 方法(Double Nudged Elastic Band)String方法Simplified String 方法寻找过渡态的chain-of-state方法CI-NEB方法ANEBA 方法(adaptive nudged elastic band approach) chain-of-states方法的一些特点高斯中opt关键字的path=M方法CPK 方法(
8、Conjugate Peak Refinement)1.过渡态过渡态结构指的是势能面上反应路径上的能量最高点，它通过最小能量路径（minimumenergy path,MEP）连接着反应物和产物的结构（如果是多步反应的机理，则这里所指反应物 或产物包括中间体）。对于多分子之间的反应，更确切来讲过渡态结构连接的是它们由无穷 远接近后因为范德华力和静电力形成的复合物结构，以及反应完毕但尚未无限远离时的复合 物结构。确定过渡态有助于了解反应机理，以及通过势垒高度计算反应速率。一般来讲，势 垒小于21kcal/mol就可以在室温下发生。在势能面上，过渡态结构的能量对坐标的一阶导数为0，只有在反应坐标方
9、向上曲率（对坐 标二阶导数）为负，而其它方向上皆为正，是能量面上的一阶鞍点。过渡态结构的能量二阶 导数矩阵（Hessian矩阵）的本征值仅有一个负值，这个负值也就是过渡态拥有唯一虚频的 来源。若将分子振动简化成谐振子模型，这个负值便是频率公式中的力常数，开根号后即得 虚数。分子构象转变、化学反应过程中往往都有过渡态的存在，即这个过程在势能面上的运动往往 都会经历满足上述条件的一点。化学反应的过渡态更确切应当成为“反应过渡态”。需要注 意的是化学反应未必都经历过渡态结构。由于过渡态结构存在时间极短，所以很难通过实验方法获得，直到飞秒脉冲激光光谱的出现 才使检验反应机理为可能。计算化学方法在目前是
10、预测过渡态的最有力武器，尽管计算上仍 有一些困难，比如其附近势能面相对于平衡结构更为平坦得多、低水平方法难以准确描述、 难以预测过渡态结构、缺乏绝对可靠的方法（如优化到能量极小点可用的最陡下降法）等。搜索过渡态的算法一般结合从头算、DFT方法，在半经验、或者小基组条件下，难以像描述 平衡结构一样正确描述过渡态结构，使得计算尺度受到了限制。结合分子力场可以描述构象 变化的过渡态，但不适用描述反应过渡态，因为大部分分子力场的势函数不允许分子拓扑结 构的改变，虽然也有一些力场如ReaxFF可以支持，有的力场还有对应的过渡态原子类型， 但目前来看适用面仍然较窄，而且不够精确，尽管更为快速。注：严格来说
11、，“过渡结构”是指势能面上反应路径上的能量最高点，而“过渡态”是指自 由能面上反应路径上的能量最高点，由于自由能变主要贡献自势能部分，所以多数情况二者 结构近似一致。但随着温度升高，往往熵变的贡献导致自由能面与势能面形状发生明显偏离， 从而导致过渡结构与过渡态明显偏离，两个词就不能混用了。但本文不涉及相关问题，故文 中过渡态、过渡结构一律指势能面上反应路径上的能量最高点。2.过渡态搜索算法基于初猜结构的算法牛顿-拉弗森法(Newton-Raphson,NR)与准牛顿法(quasi-Newton,QN)NR法是寻找函数一阶导数为0(驻点)位置的方法。通过对能量函数的泰勒级数的二阶近似 展开，然后
12、使用稳态条件dE/dr=0，可导出步进公式：下一步的坐标向量=当前坐标向量- 能量一阶导数向量* Hessian矩阵的逆矩阵。在势能面上以NR法最终找到的结果是与初猜 位置Hessian矩阵本征值正负号一致、离初猜结构最近的驻点，由于能量极小点、过渡态和 高阶鞍点的能量一阶导数皆为0，故都可以用NR法寻找。对于纯二次形函数NR法仅需一步即可找到正确位置，而势能面远比之复杂，所以需要反复 走步直至收敛。也因为势能面这个特点，为了改进优化，实际应用中NR法一般还结合线搜 索步(line search)，对于优化至极小点，就是找当前点与NR法算出来的下一点的连线上的 能量极小点作为实际下一步结构；若
13、优化至过渡态，且连线方向主要指向过渡态，则找的是 连线上能量极大点，若主要指向其它方向则找连线的能量极小点，若指向二者程度均等则一 般不做线搜索。由于精确的线搜索很花时间，所以一般只是在连线的当前位置附近计算几个 点的能量，以高阶多项式拟和后取其最小/最大点。NR法每一步需要计算Hessian矩阵并且求其逆，所以十分昂贵。QN法与NR法的走步原理一 样，但Hessian矩阵最初是用低级或经验方法猜出来的，每一步优化中通过当前及前一步的 梯度和坐标对Hessian矩阵逆矩阵逐渐修正。由于只需计算一阶导数，即便Hessian矩阵不 准确造成所需收敛步数增加，但一般仍比NR法速度快得多。QN法泛指基
14、于此原理的一类方 法，常用的是 BFGS(Broyden Fletcher Goldfarb Shanno)，此法对 Hessian 的修正保持其 对称性和正定性，最适合几何优化，但显然不能用于找过渡态。还有 DFP(Davidon-Fletcher-Powell)， MS(Murtagh-Sargent，亦称 symmetric rank 1， SR1)， PSB(Powell-symmetric-Broyden)。也有混合方法，如 Bofill 法是 PSB 和 MS 法对 Hessian 修正量的权重线性组合，比二者独立使用更优，权重系数通过位移、梯度改变量和当前 Hessian计算得到
15、，它对Hessian的修正不强制正定，很适宜搜索过渡态。将NR步进公式放到Hessian本征向量空间下其意义更为明显(此时Hessian为对角矩阵)， 可看出在每个方向上的位移就是这个方向势能的负梯度除以对应的本征值，比如在i方向上 的位移可写为 q(i)=-g(i)/h(i)，在受力越大、越平坦的方向位移越大。每一步实际位移 就是这些方向上位移的矢量和。对于寻找过渡态，因为虚频方向对应Hessian本征值为负， 使位移为受力相反方向，所以NR法在过渡态附近每一步都是使虚频方向能量升高，而在其 它正交的方向朝着能量降低的方向位移，通过这个原理步进到过渡态。若有n个虚频，则 NR法就在n个方向升
16、高能量而其它方向降低能量找到n阶鞍点。由于NR法的这个特点，为找到正确类型的驻点，初猜结构必须在目标结构的二次区域 (quadratic region)内。所谓的二次区域，是指驻点附近保持Hessian矩阵本征值符号不变 的区域，它的形状可以用多变量的二次函数近似描述，例如二维势能面情况下这样的区域可 以用F(x,y)=A*x”2+B*y”2+C*x+D*y+E*x*y来近似描述。对于能量极小点，就是指初猜点在 目标结构附近Hessian矩阵为正定矩阵的范围；对于找过渡态，就需要初猜点在它附近含有 且仅含有一个负本征值的范围内。并且这个范围内不能有其它同类驻点比目标结构距离初猜 结构更近。NR法方便之处是只需要提供一个初猜结构即可，但是由于过渡态二次区域很小(相对于
