计算化学公社

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

[Gaussian/gview] 这种情况是否存在过渡态

[复制链接 Copy URL]

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

优化了一个结构,大概推测了一下产物与反应物的能量差,减小了186kcal/mol,想找过渡态,初猜了一个结构,做完频率分析,有一个虚频,在-600左右,振动方式也是对了,就用了guess=read,geom=check,做TS优化,但优化结束后没有虚频了!!!因为之前一直没有得到稳定的反应物结构,只要一优化就跑到产物了,放出的热量也比较大,所以想是不是没有过渡态,但初猜的确实只有一个虚频啊,感觉很像过渡态啊。。。。求高人指点!!

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

23#
发表于 Post on 2015-10-10 17:05:04 | 只看该作者 Only view this author
lei234 发表于 2015-10-10 14:27
我纠结的是不知道该不该加...最后审稿人问个为什么要加D3BJ,我该怎么回答了

涉及到弱相互作用就加,免费改善精度
北京科音自然科学研究中心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

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

22#
 楼主 Author| 发表于 Post on 2015-10-10 14:27:49 | 只看该作者 Only view this author
sobereva 发表于 2015-10-10 11:03
降低那个程度是正常的。
要加D3的话所有计算从头到尾都得加,否则能量没有可比性。

我纠结的是不知道该不该加...最后审稿人问个为什么要加D3BJ,我该怎么回答了

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

21#
发表于 Post on 2015-10-10 11:03:01 | 只看该作者 Only view this author
lei234 发表于 2015-10-10 10:54
我看了您的这篇文章,试着用您推荐的B3LYP-GD3BJ优化了一下,得到的结构有点儿区别,见图。输入文件命令 ...

降低那个程度是正常的。
要加D3的话所有计算从头到尾都得加,否则能量没有可比性。
北京科音自然科学研究中心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

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

20#
 楼主 Author| 发表于 Post on 2015-10-10 10:54:19 | 只看该作者 Only view this author
sobereva 发表于 2015-10-10 10:31
原来读哪个还读哪个。

R6Disp这个就是色散校正能,高斯内部在算的时候先把这一项纳入到核互斥项里。

我看了您的这篇文章,试着用您推荐的B3LYP-GD3BJ优化了一下,得到的结构有点儿区别,见图。输入文件命令只加了EmpiricalDispersion=GD3BJ,版本是D.01。加了GD3BJ的能量(Sum of electronic and thermal Free Energies)比没有加的降低了-76.56kcal/mol, 我不知道我这个体系接下来的计算是否要加GD3BJ,请老师指教!  不胜感激!

B3LYP.png (34.71 KB, 下载次数 Times of downloads: 50)

B3LYP

B3LYP

B3LYP-GD3BJ.png (39.03 KB, 下载次数 Times of downloads: 51)

B3LYP-GD3BJ

B3LYP-GD3BJ

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

19#
发表于 Post on 2015-10-10 10:31:48 | 只看该作者 Only view this author
lei234 发表于 2015-10-10 08:46
请问加了EmpiricalDispersion=GD3BJ后,能量看哪一个?我原来读的是这一项Sum of electronic and thermal ...

原来读哪个还读哪个。

R6Disp这个就是色散校正能,高斯内部在算的时候先把这一项纳入到核互斥项里。

GD2过时了。一般都用GD3BJ,仅对于明尼苏达系列泛函如M062X才使用GD3。这方面在此文有详述
DFT-D色散校正的使用
http://sobereva.com/210
北京科音自然科学研究中心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

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2015-10-10 08:46:46 | 只看该作者 Only view this author
liyuanhe211 发表于 2015-9-30 09:58
如果是我的话会把上面那个分子单独优化,然后按照优化后的结果锁定C-O键长对这个体系优化。以及会把初猜 ...

请问加了EmpiricalDispersion=GD3BJ后,能量看哪一个?我原来读的是这一项Sum of electronic and thermal Free Energies+Zero-point correction,现在还是读这个吗?这两个能量要管吗?
R6Disp:  Grimme-D3(BJ) Dispersion energy=       -0.1641080015 Hartrees.
Nuclear repulsion after empirical dispersion term =     4632.9845362230 Hartrees.
还有,什么时候用GD2,什么时候用GD3,还有GD3BJ?
谢谢~

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

17#
 楼主 Author| 发表于 Post on 2015-10-6 11:15:44 | 只看该作者 Only view this author
sobereva 发表于 2015-10-6 11:13
再从结构合理性上琢磨吧。

在B3LYP找的过渡态上继续用M062X refine的时候用calcall结合小步长上限, ...

知道了,谢谢您。

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

16#
发表于 Post on 2015-10-6 11:13:20 | 只看该作者 Only view this author
lei234 发表于 2015-10-6 10:27
用您说的M062X/6-311G**试了,还是跑得远远的了,而且用M062X/6-311G**试了下原来用B3LYP/6-31g*找到的过 ...


再从结构合理性上琢磨吧。

在B3LYP找的过渡态上继续用M062X refine的时候用calcall结合小步长上限,以保证可靠和避免一下跑飞了。
北京科音自然科学研究中心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

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

15#
 楼主 Author| 发表于 Post on 2015-10-6 10:27:37 | 只看该作者 Only view this author
sobereva 发表于 2015-10-3 15:44
按理说应该形成氢键不该跑那么远
用M062X/6-311G**试试,Ti还用那个

用您说的M062X/6-311G**试了,还是跑得远远的了,而且用M062X/6-311G**试了下原来用B3LYP/6-31g*找到的过渡态,然后,没有了!过渡态不见了!!怎么办?

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

14#
发表于 Post on 2015-10-3 15:44:49 | 只看该作者 Only view this author
lei234 发表于 2015-10-3 10:31
图1中是IRC跑完的一个反应物结构,H(11)-O(60)的键长是1.90571,用这个结构优化,上面的分子就跑了,如图 ...

按理说应该形成氢键不该跑那么远
用M062X/6-311G**试试,Ti还用那个
北京科音自然科学研究中心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

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

13#
 楼主 Author| 发表于 Post on 2015-10-3 10:31:21 | 只看该作者 Only view this author
sobereva 发表于 2015-10-2 14:37
跑得远远的是多远?如果特别远,说明当前用的计算级别不合适,没法描述好弱相互作用,否则应当得到反应 ...

图1中是IRC跑完的一个反应物结构,H(11)-O(60)的键长是1.90571,用这个结构优化,上面的分子就跑了,如图2。目前用的方法是B3LYP,基组C H O是6-31g*,Ti 是lanl2dz,请问如果需要换方法的话,应该考虑哪几个方法,请指教~

irc-rea.png (58.73 KB, 下载次数 Times of downloads: 53)

图1

图1

irc-opt.png (52.84 KB, 下载次数 Times of downloads: 65)

图2

图2

6万

帖子

99

威望

5万

eV
积分
120160

管理员

公社社长

12#
发表于 Post on 2015-10-2 14:37:08 | 只看该作者 Only view this author
lei234 发表于 2015-10-2 11:31
用了您说的方法,得到了一个虚频的过渡态,但是振动变了,得到的产物应该不是之前的,现在也得到了我想要 ...


跑得远远的是多远?如果特别远,说明当前用的计算级别不合适,没法描述好弱相互作用,否则应当得到反应复合物。

IRC跑得特别长的时候两端就能通向反应物和产物,但与之精确结构有所偏离,理应将两端的结构作为初猜再优化反应物和产物精确结构。所以直接拿IRC两端的结构来算单点与过渡态结构的单点求差值是不适合拿来算势垒的(如果IRC本来就没跑完整,则这样得到的势垒误差会很大)。
北京科音自然科学研究中心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

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

11#
 楼主 Author| 发表于 Post on 2015-10-2 11:31:23 | 只看该作者 Only view this author
sobereva 发表于 2015-10-1 04:29
GD3是用零阻尼形式的DFT-D3校正,见
DFT-D色散校正的使用
http://sobereva.com/210

用了您说的方法,得到了一个虚频的过渡态,但是振动变了,得到的产物应该不是之前的,现在也得到了我想要的过渡态(自己不知道怎么瞎整给整出来了)。又出现了一个问题,就是跑完irc再优化的时候,反应物分子就跑了,解离了,跑得远远的,我想能不能用最后irc得到的反应物和产物的结构直接做个单点计算,不知道这种方法是否可行?就是过渡态是优化得到的,但反应物和产物是跑完irc的单点能?这样似乎也不能直接计算能垒,如果给过渡态最后也算个单点,这样可以吗?所有的能量都是单点能...总觉得不对劲...请您多多指教~谢谢~

286

帖子

0

威望

1136

eV
积分
1422

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2015-10-1 15:57:36 | 只看该作者 Only view this author
sobereva 发表于 2015-10-1 04:29
GD3是用零阻尼形式的DFT-D3校正,见
DFT-D色散校正的使用
http://sobereva.com/210

谢谢您的回复~

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

GMT+8, 2025-8-16 15:57 , Processed in 0.187609 second(s), 24 queries , Gzip On.

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