找过渡态的基本经验

1、交流找过渡态的一些经验已有10人参与gmy1990（金币+5）:谢谢分享！2011-06-0622:14:36zhou2009:编辑内容2011-06-2209:54zhou2009:标题咼亮2011-06-2210:43小红豆:提升帖子2011-07-1418:44小红豆:提升帖子2011-07-2019:24小红豆:提升帖子2011-07-2219:33最近经常见虫友问关于过渡态怎么找的问题，正好我最近比较郁闷就给大家写个简单的教程，仅限新手入门。过渡态简单点说就是化学反应中瞬间形成的高自由能的不稳定复合物，形象点说就是就是反应势能面上的一个鞍点，所谓的找过渡态也就是要找到这个反应势能面上

2、的这个鞍点，这里以DA反应作为例子简单说下我用gaussian找过渡态的经验。Gaussian里面找过渡态主要有3种方法：TS、QST2以及QST3。这里套用gaussian手册里的说明。TS是进行过渡态而不是局部最小值的优化计算。QST2是使用STQN方法寻找过渡态结构。这个选项需要输入反应物和产物分子结构，先后通过两组连续的标题和分子说明部分定义。注意在这两个结构中的原子顺序必需一致TS不能和QST2合用。QST3是使用STQN方法寻找过渡态结构。这个选项需要输入反应物，产物和最初的过渡态结构，先后由三组连续的标题和分子说明部分定义。注意在这三个结构中的原子顺序必需一致。TS不能和QST3

3、合用。接下来贴出几个例子来说明下TS、QST2和QST3的输入文件该怎么写%chk=c2h4_c4h6_ts%nproc=1%mem=200mb#opt=(ts,calcfc)freq=noramanb3lyp/6-31g*TitleCardRequired01C-0.379602001.410261000.50972400H-0.064094001.040254001.48028400H-0.265895002.480672000.40089100C-1.260203000.70558300-0.28508700H-1.846358001.22289600-1.04415700C-1.260

4、23800-0.70552500-0.28508700H-1.84641800-1.22280900-1.04415900C-0.37967300-1.410247000.50972300H-0.06414500-1.040256001.48028200H-0.26601500-2.480663000.40088700C1.456482000.69084700-0.25401700H1.984024001.246920000.51079100H1.292838001.24312300-1.17152100C1.45645300-0.69091000-0.25401800H1.29278400-

5、1.24317800-1.17152200H1.98396500-1.247007000.51079200qst2%chk=c2h4_c4h6_qst2%nproc=1%mem=200mb#opt=(calcfc,qst2)freq=noramanb3lyp/6-31g*TitleCardRequired-3.177919061.13289796-1.01187682-2.644755310.20519304-1.01187682H-4.247919061.13289796-1.01187682C-2.502644752.30787525-1.01187682H-3.035808503.235

6、58018-1.01187682C-0.962644752.30787525-1.01187682H-0.429481003.23558018-1.01187682C-0.280401351.13693077-1.01187682H-0.808043820.20607449-1.01187682H0.789579771.14328807-1.01187682C-2.597764130.988942791.10133233H-3.140310630.066925591.12200992H-3.126998871.918463751.07302387C-1.242608730.983479681.

7、11087410H-0.700062231.905496891.09019651H-0.713373980.053958731.13918256TitleCardRequired01C1.414324000.094601000.31826800H1.423700000.107754001.42899000H2.474456000.127801000.00402400C0.668660001.30146700-0.16992700H1.264389002.14276400-0.50621500C-0.668907001.30134300-0.16992700H-1.264793002.14252

8、900-0.50621700C-1.414346000.094338000.31826600H-1.423727000.107494001.42898900H-2.474484000.127340000.00402000C0.77040400-1.21422100-0.17552700H1.14229400-2.052881000.44176200H1.12544800-1.41421500-1.20543600C-0.77017000-1.21437000-0.17552800H-1.12517200-1.41443600-1.20543700H-1.14189800-2.053102000

9、.44176200qst3%chk=c2h4_c4h6_qst3%nproc=1%mem=200mb#opt=(calcall,qst3,noeigentest)freq=noramanb3lyp/6-31g*TitleCardRequired01C-3.177919061.13289796-1.01187682H-2.644755310.20519304-1.01187682H-4.247919061.13289796-1.01187682C-2.502644752.30787525-1.01187682H-3.035808503.23558018-1.01187682C-0.9626447

10、52.30787525-1.01187682H-0.429481003.23558018-1.01187682C-0.280401351.13693077-1.01187682H-0.808043820.20607449-1.01187682H0.789579771.14328807-1.01187682C-2.597764130.988942791.10133233H-3.140310630.066925591.12200992H-3.126998871.918463751.07302387C-1.242608730.983479681.11087410H-0.700062231.90549

11、6891.09019651H-0.713373980.053958731.13918256TitleCardRequired01C1.414324000.094601000.31826800H1.423700000.107754001.42899000H2.474456000.127801000.00402400C0.668660001.30146700-0.16992700H1.264389002.14276400-0.50621500C-0.668907001.30134300-0.16992700H-1.264793002.14252900-0.50621700C-1.414346000

12、.094338000.31826600H-1.423727000.107494001.42898900H-2.474484000.127340000.00402000C0.77040400-1.21422100-0.17552700H1.14229400-2.052881000.44176200H1.12544800-1.41421500-1.20543600C-0.77017000-1.21437000-0.17552800H-1.12517200-1.41443600-1.20543700H-1.14189800-2.053102000.44176200TitleCardRequired0

13、1C0.39335100-1.407522000.50945100H0.07415700-1.040745001.48001500H0.28923500-2.478840000.39989000C1.26702700-0.69422900-0.28512700H1.85846600-1.20567200-1.04409200C1.253468000.71680900-0.28478600H1.834993001.23987800-1.04346600C0.365749001.413091000.50971100H0.054145001.040536001.48056200H0.24248700

14、2.482465000.40066900C-1.45041500-0.70335800-0.25427100H-1.97290400-1.264326000.51044000H-1.28115600-1.25385900-1.17177300C-1.462407000.67835900-0.25370700H-1.304679001.23254800-1.17112400H-1.995375001.229107000.51125200这里需要注意下反应物的构型以及产物的构型都必须是优化过的。QST3最后一个构型是自己猜测的过渡态的构型，所以这里需要一定的化学功底。另外使用QST2以及QST3时

15、要注意你所假定的反应物与产物中间应该没有别的中间体并且只有一个过渡态，否则不保证你能优化出来什么东西所以我比较推荐直接用TS进行过渡态优化，这当然需要涉及到自己猜过渡态的一个过程，需要对一些常见的化学反应机理有个细致的了解才能很好的猜过渡态。选定用哪一种方法后接着就是如何获取力常数，是使用rcfc,calcfc还是calcall由于这一步直接关系着整个体系的计算量所以需要慎重选择，如果你是读上一个计算chk文件进行过渡态优化，并且上一个计算做了频率分析，那毫无疑问要加上rcfc，这可以节省大量的时间。如果从头开始找过渡态的话还是calcfc的好，也就是在第一步计算下力常数，至于calcall个

16、人推荐不到万不得已不要使用，虽然有时它也可以最快的找到过渡态，例如我曾经试过加calcall优化15步找到过渡态，而用calcfc花了47步才找到过渡态，但仍不推荐使用-！剩下还有一种方法找过渡态主反应坐标法，该方法通过指定体系中主要参加反应的键的原子通过对其初始Hessian矩阵中该坐标的行和列计算数值二阶导数的方法来获得过渡态，示例如下。%chk=c2h4_c4h6_ts_d%mem=200mb%nproc=1#opt=(ts,modredundant)freq=noramanb3lyp/6-31g*TitleCardRequired01C-0.264421000.424036001.41

17、432400H-1.298998000.828422001.42370100H0.041311000.344168002.47445600C0.618365001.380883000.66866500H1.229978002.049311001.26439800C0.618365001.38088300-0.66866500H1.229978002.04931100-1.26439800C-0.264421000.42403600-1.41432400HT.298998000.82842200-1.42370100H0.041311000.34416800-2.47445600C-0.2644

18、2100-1.616745000.77039900H-1.13801600-2.183518001.14228500H0.62859000-2.167419001.12544200C-0.26442100-1.61674500-0.77039900H0.62859000-2.16741900-1.12544200H-1.13801600-2.18351800-1.14228500B814DB111D这里有一些我个人用TS找过渡态的经验：优先推荐两种方法，第一步都是将主要参加反应的键长延长0.2-0.5埃或者延长20%-40%，然后固定该键进行构象优化，然后再用TS+calcfc找过渡态或者用主

19、反应坐标法找过渡态，一般情况下这样都能找到相应的过渡态，其中主反应坐标法由于不用算力常数在有些条件下会稍微快一点。这样还没有找到的话可以看一下输出文件中每一步的能量以及构象，然后从中选取最可能是过渡态的那个构象缩小步长（maxstep=7oriop（1/8=7）继续找过渡态。如果仍然没有找到可以对你猜测的构象进行频率分析，然后查看有没有你所要的虚频，如果有的话用GV最下面的manualdisplacement调节震动，选中你需要的虚频向左调0.2-0.4或者你不要的虚频向右调0.2-0.4以微调构象，然后接着套用上一种方法的后半部分继续找过渡态。或者从产物构象SCAN参加反应的键长，逐渐拉长键长