计算化学公社

 找回密码 Forget password
 注册 Register
Views: 31547|回复 Reply: 19
打印 Print 上一主题 Last thread 下一主题 Next thread

[新手求助] 有关过渡态的疑惑

[复制链接 Copy URL]

286

帖子

0

威望

1134

eV
积分
1420

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
#
第一种情况:如图1. 过渡态的能量竟然比产物还低,确认了过渡态的振动方式是对的,也只有一个虚频,请问该怎么解释?可以说这步反应没有过渡态吗?

第二种情况:如图2.在拉大断键距离固定优化直到新建生成的过程中,能量一直在降低,这种情况是不是不存在过渡态,还需要再寻找过渡态吗?

谢谢~

QQ图片20150708095110.png (7.86 KB, 下载次数 Times of downloads: 98)

图1

图1

QQ图片20151015101241.png (9.75 KB, 下载次数 Times of downloads: 111)

图2

图2

5万

帖子

99

威望

5万

eV
积分
112428

管理员

公社社长

19#
发表于 Post on 2023-12-8 05:58:38 | 只看该作者 Only view this author
Kamistry 发表于 2023-12-7 10:14
不好意思老师,本来想写所有关键词的,后来改的时候把em=gd3bj删了忘标出来了,优化和单点能用的其实是B3 ...



能量折线图如果用电子能量,过渡态照常标上;如果用自由能,上一个极小点直接连接下一个极小点
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

35

帖子

3

威望

794

eV
积分
889

Level 4 (黑子)

18#
发表于 Post on 2023-12-7 10:14:03 | 只看该作者 Only view this author
sobereva 发表于 2023-12-7 02:41
老生常谈的问题,电子能垒低的时候这极其正常,说明这不是化学意义上的过渡态

B3LYP算能量早过时了

不好意思老师,本来想写所有关键词的,后来改的时候把em=gd3bj删了忘标出来了,优化和单点能用的其实是B3LYP-D3(BJ)。

不是化学意义上的过渡态的意思是这一步实际是无能垒过程吗?在能量折线图上只需直接连接吗?

5万

帖子

99

威望

5万

eV
积分
112428

管理员

公社社长

17#
发表于 Post on 2023-12-7 02:41:18 | 只看该作者 Only view this author
Kamistry 发表于 2023-12-6 19:21
请问老师,我找到一个过渡态,虚频-1015,Infrared 4380,过渡态结构走完IRC,IRC两端结构拿出来优化,优化 ...

老生常谈的问题,电子能垒低的时候这极其正常,说明这不是化学意义上的过渡态

B3LYP算能量早过时了
仔细看
坚持使用B3LYP算有机体系能量的人的下场
http://bbs.keinsci.com/thread-12773-1-1.html

简谈量子化学计算中DFT泛函的选择
http://sobereva.com/272http://bbs.keinsci.com/thread-536-1-1.html

谈谈量子化学研究中什么时候用B3LYP泛函优化几何结构是适当的
http://sobereva.com/557http://bbs.keinsci.com/thread-17899-1-1.html

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

35

帖子

3

威望

794

eV
积分
889

Level 4 (黑子)

16#
发表于 Post on 2023-12-6 19:21:03 | 只看该作者 Only view this author
请问老师,我找到一个过渡态,虚频-1015,Infrared 4380,过渡态结构走完IRC,IRC两端结构拿出来优化,优化出来的结构单点能确实比过渡态低,但加上热力学校正后焓和自由能反而比过渡态高了。请问可能是什么原因呢?涉及到阴离子,优化用的是b3lyp/6-31g(d),单点能用的是b3lyp/ma-TZVP

5万

帖子

99

威望

5万

eV
积分
112428

管理员

公社社长

15#
发表于 Post on 2016-10-9 22:56:23 | 只看该作者 Only view this author
珊珊来迟 发表于 2016-10-9 16:39
恩恩,好的sob老师,我先做IRC看看。
不过sob老师还有一个问题,您说的U(0)是指输出文件中Sum of electr ...

U(0)对应Sum of electronic and zero-point Energies
电子能垒、内能能垒、焓能垒、自由能垒都叫势垒,用哪个取决于拿它讨论什么
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

274

帖子

0

威望

4278

eV
积分
4552

Level 6 (一方通行)

14#
发表于 Post on 2016-10-9 18:54:47 | 只看该作者 Only view this author
lei234 发表于 2015-10-15 16:41
你是说的图1的情况吧?我以为是图2的情况呢,不好意思。。。

图1有了过渡态,做个IRC是很快的。

146

帖子

0

威望

1273

eV
积分
1419

Level 4 (黑子)

13#
发表于 Post on 2016-10-9 16:39:11 | 只看该作者 Only view this author
sobereva 发表于 2016-10-9 16:26
做IRC确认确实连同反应物再说
如果能连通,电子能垒较小时TS的U(0)比反应物更低也是常有的事。

恩恩,好的sob老师,我先做IRC看看。
不过sob老师还有一个问题,您说的U(0)是指输出文件中Sum of electronic and thermal Energies对应的值么?
势垒的话看的到底是要看以下中的哪一个对应的值?
Sum of electronic and zero-point Energies
Sum of electronic and thermal Energies
Sum of electronic and thermal Enthalpies
Sum of electronic and thermal Free Energies
我理解的是第一个
不过从我的结果来看自由能这一项是过渡态能量比反应物和生成物要高

5万

帖子

99

威望

5万

eV
积分
112428

管理员

公社社长

12#
发表于 Post on 2016-10-9 16:26:39 | 只看该作者 Only view this author
珊珊来迟 发表于 2016-10-9 16:14
sob老师我在算一个质子转移的反应的过渡态发现过渡态的能量(零点校正后的电子总能量)比反应物能量低是 ...

做IRC确认确实连同反应物再说
如果能连通,电子能垒较小时TS的U(0)比反应物更低也是常有的事。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

146

帖子

0

威望

1273

eV
积分
1419

Level 4 (黑子)

11#
发表于 Post on 2016-10-9 16:14:16 | 只看该作者 Only view this author
sobereva 发表于 2015-10-15 12:32
1 是指电子能量吧?是的话,过渡态能量比反应物、产物还低必然说明你漏掉了中间体。你从过渡态做IRC(然后 ...

sob老师我在算一个质子转移的反应的过渡态发现过渡态的能量(零点校正后的电子总能量)比反应物能量低是什么原因呢?数值是-5.71Kcal/mol 我得到的过渡态虚频对应的振动模式是H原子在相应的两个原子之间振动,不过还没有对过渡态做IRC

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

10#
发表于 Post on 2015-10-15 17:30:23 | 只看该作者 Only view this author
lei234 发表于 2015-10-15 17:08
电子能变化如图,也是一直在降。。

你这不叫IRC曲线

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

9#
发表于 Post on 2015-10-15 17:29:54 | 只看该作者 Only view this author
lei234 发表于 2015-10-15 16:41
你是说的图1的情况吧?我以为是图2的情况呢,不好意思。。。

我没说任何一个图的情况,我只说IRC。IRC只考虑过渡态的输入。自己查一下IRC是什么意思、过渡态该怎么算。

286

帖子

0

威望

1134

eV
积分
1420

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2015-10-15 17:08:35 | 只看该作者 Only view this author
sobereva 发表于 2015-10-15 14:49
你先别管自由能,先做个IRC曲线,用电子能量,看到底是否TS电子能量比你说的“反应物”和“产物”能量 ...

电子能变化如图,也是一直在降。。

QQ图片20151015170626.png (13.37 KB, 下载次数 Times of downloads: 108)

HF

HF

286

帖子

0

威望

1134

eV
积分
1420

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2015-10-15 16:41:54 | 只看该作者 Only view this author
liyuanhe211 发表于 2015-10-15 16:17
IRC没做呢怎么知道你找到的是你这个基元反应的过渡态?
做IRC就是 # IRC

你是说的图1的情况吧?我以为是图2的情况呢,不好意思。。。

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

6#
发表于 Post on 2015-10-15 16:17:36 | 只看该作者 Only view this author
lei234 发表于 2015-10-15 14:50
才在找过渡态的嘛,怎么做IRC。能量一直在降,都不知道怎么找了。。。原来都是拉大断键的两个原子进行固 ...

IRC没做呢怎么知道你找到的是你这个基元反应的过渡态?
做IRC就是 # IRC

286

帖子

0

威望

1134

eV
积分
1420

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2015-10-15 14:50:12 | 只看该作者 Only view this author
liyuanhe211 发表于 2015-10-15 12:51
键长增大能量降低(如果没做错)显然至少说明你的反应物(产物)不是最优的过渡态经过IRC验证了吗?

才在找过渡态的嘛,怎么做IRC。能量一直在降,都不知道怎么找了。。。原来都是拉大断键的两个原子进行固定优化后能量有个山峰的嘛,这倒好,一直往下降,都不知道怎么哪个是过渡态的结构。。。

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

GMT+8, 2024-11-25 22:30 , Processed in 0.198965 second(s), 24 queries , Gzip On.

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