计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] 简谈Gaussian里找过渡态的关键词opt=TS和QST2、3

  [复制链接 Copy URL]

147

帖子

4

威望

1575

eV
积分
1802

Level 5 (御坂)

科研泥瓦匠

16#
发表于 Post on 2020-6-18 21:25:31 | 只看该作者 Only view this author
本帖最后由 欢乐多 于 2020-6-19 20:33 编辑
sobereva 发表于 2020-6-17 11:16
从TS跑IRC,直接连通哪个就把哪个当这一步的反应物
谢谢老师!
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

147

帖子

4

威望

1575

eV
积分
1802

Level 5 (御坂)

科研泥瓦匠

17#
发表于 Post on 2020-6-19 16:45:49 | 只看该作者 Only view this author
本帖最后由 欢乐多 于 2020-6-19 20:13 编辑
sobereva 发表于 2020-6-17 11:16
从TS跑IRC,直接连通哪个就把哪个当这一步的反应物

老师,您好,
图中E即是我要的产物。从IRC1跑的结果再次优化,得到C;从IRC2跑的结果,再次优化得到D,但是C与D的构象不一样,并且C构象也不能优化到D构象去,我的理解是:
1,以上TS1、C、D、TS2等找的都对。图中C与D构象可以互变,若以C为TS1和TS2的中间体时,TS2的能垒为E4;若以D为中间体时,TS2的能垒为E3,TS1的生成能垒为E1,
2,如果1理解错误,可能就是http://sobereva.com/400博文中提到的“(2)优化的时候由于走得步子太大,或者Hessian矩阵不准确,导致优化到了其它邻近的极小点去。......先把IRC跑得尽量长一些,使两端结构尽可能接近极小点时再做几何优化”,就是说TS1、TS2找对了,但是极小点C或D找的不对。
老师,您看理解是否有问题?

屏幕快照 2020-06-19 下午4.20.22.png (24 KB, 下载次数 Times of downloads: 88)

屏幕快照 2020-06-19 下午4.20.22.png
终日寻春不见春
芒鞋踏破岭头云
归来偶把梅花嗅
春在枝头已十分

6万

帖子

99

威望

6万

eV
积分
125116

管理员

公社社长

18#
 楼主 Author| 发表于 Post on 2020-6-21 03:32:33 | 只看该作者 Only view this author
欢乐多 发表于 2020-6-19 16:45
老师,您好,
图中E即是我要的产物。从IRC1跑的结果再次优化,得到C;从IRC2跑的结果,再次优化得到D, ...

当成C和D之间构象变化也需要经历过渡态,作为整个反应中的组成部分
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

32

帖子

0

威望

640

eV
积分
672

Level 4 (黑子)

19#
发表于 Post on 2026-1-5 23:15:11 | 只看该作者 Only view this author
老师您好,请问qst2方法,要加noeigen关键词吗

6万

帖子

99

威望

6万

eV
积分
125116

管理员

公社社长

20#
 楼主 Author| 发表于 Post on 2026-1-6 00:31:16 | 只看该作者 Only view this author
Dave666 发表于 2026-1-5 23:15
老师您好,请问qst2方法,要加noeigen关键词吗

不需要
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

42

帖子

0

威望

485

eV
积分
527

Level 4 (黑子)

21#
发表于 Post on 2026-1-26 16:27:59 | 只看该作者 Only view this author
老师您好,我在计算一个催化过程,其中有一步中间体IM2转IM3整个过程ΔG=3kJ/mol,根据经验推测应该是一个势垒极低的反应。不知道是不是因为势能面太过平缓,用TS(Berry)算得不到过渡态结构,优化所得经常是没虚频的结构。
由于整个过程已知存在其他决速步,这一步能垒对算反应速率常数没啥影响,但整个机理缺失了这一环怕别人不认可。
请问如果不需要得到这一步准确能垒的话,可以做怎样的一个计算来让别人接受IM2转IM3的过程能垒很低,不影响整体动力学的推导。

6万

帖子

99

威望

6万

eV
积分
125116

管理员

公社社长

22#
 楼主 Author| 发表于 Post on 2026-1-26 18:11:09 | 只看该作者 Only view this author
ginlpein 发表于 2026-1-26 16:27
老师您好,我在计算一个催化过程,其中有一步中间体IM2转IM3整个过程ΔG=3kJ/mol,根据经验推测应该是一个 ...

如果要证明没有过渡态(电子能量意义上的),要么做坐标合适的柔性扫描体现出势能面上就是没过渡态,要么证实以反应物为初始结构时做几何优化能直接优化到产物结构因此之间不隔着过渡态。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

42

帖子

0

威望

485

eV
积分
527

Level 4 (黑子)

23#
发表于 Post on 2026-1-26 21:21:20 | 只看该作者 Only view this author
sobereva 发表于 2026-1-26 18:11
如果要证明没有过渡态(电子能量意义上的),要么做坐标合适的柔性扫描体现出势能面上就是没过渡态,要么 ...

sob老师,这个坐标合适的柔性扫描指的是指键断裂/生成方向进行柔扫吗?
另外,那针对势能面很平缓的情况,如果硬要出一个过渡态的话,有什么方法改进可以比较好的捕捉到不咋明显的过渡态?

6万

帖子

99

威望

6万

eV
积分
125116

管理员

公社社长

24#
 楼主 Author| 发表于 Post on 2026-1-28 03:04:40 | 只看该作者 Only view this author
ginlpein 发表于 2026-1-26 21:21
sob老师,这个坐标合适的柔性扫描指的是指键断裂/生成方向进行柔扫吗?
另外,那针对势能面很平缓的情况 ...

对反应过程中最关键的几何变量做扫描。对于简单断裂一个键的过程来说就是这个键长
恰当摆过渡态初猜(最关键。包括借助柔性扫描),用较准确的Hessian(恰当用recalc=n),尝试较小步长上限。有时提升DFT积分格点精度有益
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-16 20:34 , Processed in 0.230299 second(s), 22 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list