版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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参加反应的键长,逐渐拉长键长
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 面食加工过程中的品质监控考核试卷
- 锡冶炼行业产业链分析与发展策略考核试卷
- 风扇在空气净化市场的发展趋势考核试卷
- 鞋类行业创新管理与实践案例分析考核试卷
- 长途客运企业安全生产责任制落实考核试卷
- 颜料在声学材料中的应用考核试卷
- 橡胶制品在电子产品中的功能化设计考核试卷
- 集装箱道路运输创新与发展趋势考核试卷
- 鹅的饲养与饲养环境监测技术发展展望分析考核试卷
- 文化用品行业品牌与企业文化考核试卷
- 八年级物理第一次月考卷(考试版A3)(辽宁专用人教版2024第1~3章第2节)
- 第一次月考 (1-2单元)(月考)- 2024-2025学年五年级上册数学人教版
- 建筑施工安全技术操作规程
- 2024至2030年中国生物质能发电行业市场深度调研及发展前景分析报告
- 财务主管岗位招聘笔试题及解答(某大型集团公司)
- 简易呼吸器的使用-
- 2024年《金融市场基础知识》冲刺复习讲义
- 2024年新青岛版六年级上册(六三制)科学全册知识点
- 5.2 外力作用与地表形态高三地理一轮复习课件
- NB/T 11446-2023煤矿连采连充技术要求
- 2024年度山东省招聘社区工作者试题及答案
评论
0/150
提交评论