计算化学公社

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

[ORCA] ORCA优化过程中能量变化

[复制链接 Copy URL]

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

最近用ORCA优化一个结构,最后把单点能提出来看了看,发现很奇怪,在能量较低的时候没有收敛,反而在能量升高一点的时候收敛了,这是不是就可以说明这个最终结构不是最优结构呢?


6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2016-4-20 18:13:05 | 只看该作者 Only view this author
可以这么说。
ORCA的优化不像高斯那么稳健
北京科音自然科学研究中心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

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2016-4-20 20:10:03 | 只看该作者 Only view this author
sobereva 发表于 2016-4-20 18:13
可以这么说。
ORCA的优化不像高斯那么稳健

sob老师,那针对这种情况有没有好的办法解决呢?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

4#
发表于 Post on 2016-4-20 20:13:14 | 只看该作者 Only view this author
yaochuang 发表于 2016-4-20 20:10
sob老师,那针对这种情况有没有好的办法解决呢?

取最低能量结构,减小步长上限再优化看看。
北京科音自然科学研究中心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

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2016-4-21 10:59:46 | 只看该作者 Only view this author
sobereva 发表于 2016-4-20 20:13
取最低能量结构,减小步长上限再优化看看。

谢谢sob老师,我尝试了,将maxstep 设置的更小,并且用了更加严格的收敛限,但是最终的结果还是一样,能量会升高然后收敛,而且我比较了两种方法得到的结构,基本完全一样(求了一下RMSD为零)。而且有时候orca在判断收敛时会忽略掉能量的变化项,所以会不会是orca软件本身在定义收敛时对能量的要求就没有那么严格呢?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

6#
发表于 Post on 2016-4-21 17:38:14 | 只看该作者 Only view this author
yaochuang 发表于 2016-4-21 10:59
谢谢sob老师,我尝试了,将maxstep 设置的更小,并且用了更加严格的收敛限,但是最终的结果还是一样,能 ...

关键词怎么写的?
北京科音自然科学研究中心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

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2016-4-21 18:49:58 | 只看该作者 Only view this author
sobereva 发表于 2016-4-21 17:38
关键词怎么写的?

关键词是这样的:
!b3lyp/g def2-svp rijcosx def2-svp/C def2-svp/J tightopt   
!pal6                                                         
%maxcore 1000   
%geom
   maxstep 0.1              
end

最初用的是:
!b3lyp/g def2-svp rijcosx def2-svp/C def2-svp/J opt   
!pal6                                                         
%maxcore 1000   

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

8#
发表于 Post on 2016-4-21 19:28:44 | 只看该作者 Only view this author
不要写def2-svp/C def2-svp/J,应当写def2-SVP/JK
你也看看轨迹,考察具体结构是怎么变化的
北京科音自然科学研究中心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

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

9#
 楼主 Author| 发表于 Post on 2016-4-22 14:54:03 | 只看该作者 Only view this author
本帖最后由 yaochuang 于 2016-4-22 15:48 编辑
sobereva 发表于 2016-4-21 19:28
不要写def2-svp/C def2-svp/J,应当写def2-SVP/JK
你也看看轨迹,考察具体结构是怎么变化的

sob老师,我用def2-SVP/JK代替了原来的辅助基组,结果还是一样。并测试了两个结构的RMSD结果是0.007 埃(这应该是由辅助基组改变造成的)。
从轨迹上看就是基团间的小幅的调整,没有看出可疑的地方。
对了,我还想问一下对于def2-svp/C def2-svp/J和def2-svp/jk这两种辅助基组有没有优劣比较呢?因为记得在哪儿好像看到过说分别指定Coulomb和correlation辅助基组会更好一些,不知道是不是这样?

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

10#
发表于 Post on 2016-4-22 18:19:31 | 只看该作者 Only view this author
yaochuang 发表于 2016-4-22 14:54
sob老师,我用def2-SVP/JK代替了原来的辅助基组,结果还是一样。并测试了两个结构的RMSD结果是0.007 埃( ...

那就在优化中用精确计算的Hessian
/C是给后HF计算用的,DFT计算无用。
/J是给纯泛函用的,杂化泛函应当用/JK。
北京科音自然科学研究中心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

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

11#
 楼主 Author| 发表于 Post on 2016-4-24 08:30:31 | 只看该作者 Only view this author
sobereva 发表于 2016-4-22 18:19
那就在优化中用精确计算的Hessian
/C是给后HF计算用的,DFT计算无用。
/J是给纯泛函用的,杂化泛函应当 ...

sob老师,我添加了每一步都计算Hess矩阵,但是结果还是这样:


这是我的输入:
!b3lyp/g def2-svp rijcosx def2-svp/JK tightopt  tightscf
!pal6                                                         
%maxcore 1000   
%geom
   maxstep 0.1
   Recalc_Hess 1
end

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

12#
发表于 Post on 2016-4-24 09:28:15 | 只看该作者 Only view this author
yaochuang 发表于 2016-4-24 08:30
sob老师,我添加了每一步都计算Hess矩阵,但是结果还是这样:

那你索性还是先用高斯搞搞吧,要不然虽然ORCA更快,翻来覆去折腾时间也都浪费了,还不如直接用高斯
北京科音自然科学研究中心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

147

帖子

0

威望

2416

eV
积分
2563

Level 5 (御坂)

13#
 楼主 Author| 发表于 Post on 2016-4-24 10:53:54 | 只看该作者 Only view this author
sobereva 发表于 2016-4-24 09:28
那你索性还是先用高斯搞搞吧,要不然虽然ORCA更快,翻来覆去折腾时间也都浪费了,还不如直接用高斯

嗯,好的!谢谢sob老师耐心的解答。

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 17:48 , Processed in 0.172283 second(s), 23 queries , Gzip On.

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