计算化学公社

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

[Gaussian/gview] opt freq

[复制链接 Copy URL]

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
请问:我用的相同的基组和方法,初始构型一样,一种是算opt  一种是freq  为什么两种结果优化得到的构型不一样呢      只算opt的出现4个yes
  算opt和freq的只前两个出现yes  这是为什么呢  求指教?

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

15#
发表于 Post on 2014-11-17 10:57:44 | 只看该作者 Only view this author
直接把输出文件压缩一下发上来(如果压缩后>5MB,用网盘)
北京科音自然科学研究中心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

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

14#
 楼主 Author| 发表于 Post on 2014-11-17 10:23:56 | 只看该作者 Only view this author
sobereva 发表于 2014-11-17 09:53
光从你的叙述上我弄不清楚。
freq过程之后肯定会判断一次YES/NO,这和写不写tight没关系。

..skipping...
         Item               Value     Threshold  Converged?
Maximum Force            0.000000     0.000015     YES
RMS     Force            0.000000     0.000010     YES
Maximum Displacement     0.000007     0.000060     YES
RMS     Displacement     0.000001     0.000040     YES
Predicted change in Energy=-9.705411D-14
Optimization completed.


这是最后一个converged   之前这个不是在倒数第二个吗?为什么现在变成最后一个了呢

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

13#
发表于 Post on 2014-11-17 09:53:21 | 只看该作者 Only view this author
小菜鸟求关照 发表于 2014-11-17 09:35
我发现加上opt(tight)以后就没有freq后面的converged了,搜索的converged直接就是complete上面那一个, ...

光从你的叙述上我弄不清楚。
freq过程之后肯定会判断一次YES/NO,这和写不写tight没关系。
北京科音自然科学研究中心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

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

12#
 楼主 Author| 发表于 Post on 2014-11-17 09:35:11 | 只看该作者 Only view this author
sobereva 发表于 2014-11-12 09:23
你所说的倒数第二圈是opt最后一步。正是因为这里全yes了所以优化才收敛了。这步判断是否yes是根据优化时 ...

我发现加上opt(tight)以后就没有freq后面的converged了,搜索的converged直接就是complete上面那一个,这是为什么呢?这样计算的结果准确吗?

93

帖子

0

威望

5308

eV
积分
5401

Level 6 (一方通行)

11#
发表于 Post on 2014-11-12 09:39:40 | 只看该作者 Only view this author
是的 这样可以减少很多计算时间
树欲静而风不止

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

10#
发表于 Post on 2014-11-12 09:35:02 | 只看该作者 Only view this author
lyfchem 发表于 2014-11-12 09:27
哦 如果用opt(calcall)计算相当于opt也是用精确Hessian 但是成本太高 没必要这样 是吗 老师

是的。
实际上,就算最后freq出来不是都yes,但NO的项里离收敛限偏差不太远,也没关系,结构是可用的。
北京科音自然科学研究中心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

93

帖子

0

威望

5308

eV
积分
5401

Level 6 (一方通行)

9#
发表于 Post on 2014-11-12 09:27:50 | 只看该作者 Only view this author
哦 如果用opt(calcall)计算相当于opt也是用精确Hessian 但是成本太高 没必要这样 是吗 老师
树欲静而风不止

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

8#
发表于 Post on 2014-11-12 09:23:13 | 只看该作者 Only view this author
lyfchem 发表于 2014-11-12 09:18
两个YES是因为在opt那一步用的Hessian不是精确的?那为什么往往opt+freq计算的倒数第二圈是4个yes呢?

你所说的倒数第二圈是opt最后一步。正是因为这里全yes了所以优化才收敛了。这步判断是否yes是根据优化时用的近似Hessian判断的,而最后一次是freq之后用精确Hessian判断的,结论不同正因为Hessian用的不同。
北京科音自然科学研究中心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

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2014-11-12 09:22:46 | 只看该作者 Only view this author
lyfchem 发表于 2014-11-12 09:18
两个YES是因为在opt那一步用的Hessian不是精确的?那为什么往往opt+freq计算的倒数第二圈是4个yes呢?

我那种也是这样情况 很想知道为什么

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2014-11-12 09:22:21 | 只看该作者 Only view this author
sobereva 发表于 2014-11-11 21:52
如果你是指第一种情况只写opt,第二种情况写opt freq,那么第二种情况最后只有两个YES是因为在opt那一步用 ...

opt的收敛限设成tight  是写成这样的关键词吗:opt(tight)

145

帖子

0

威望

290

eV
积分
435

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2014-11-12 09:21:43 | 只看该作者 Only view this author
378024572 发表于 2014-11-11 19:31
为什么opt和freq不一起用呢?

因为一起用的时候freq后面出现两个no我想看一下是不是结构的问题

93

帖子

0

威望

5308

eV
积分
5401

Level 6 (一方通行)

4#
发表于 Post on 2014-11-12 09:18:12 | 只看该作者 Only view this author
两个YES是因为在opt那一步用的Hessian不是精确的?那为什么往往opt+freq计算的倒数第二圈是4个yes呢?
树欲静而风不止

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

3#
发表于 Post on 2014-11-11 21:52:22 | 只看该作者 Only view this author
如果你是指第一种情况只写opt,第二种情况写opt freq,那么第二种情况最后只有两个YES是因为在opt那一步用的Hessian不是精确的,而freq的时候则用的Hessian总是精确的。写opt(calcall) freq使得opt时也用精确Hessian才能绝对确保freq最后显示的也是四个yes,当然这样在优化过程中成本很高,故也可以考虑不用calcall而把opt的收敛限设成tight,这样使得优化出的结构更精准,也往往能使freq之后NO出现的可能性更少。
北京科音自然科学研究中心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

15

帖子

0

威望

100

eV
积分
115

Level 2 能力者

2#
发表于 Post on 2014-11-11 19:31:15 | 只看该作者 Only view this author
为什么opt和freq不一起用呢?

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

GMT+8, 2025-8-14 08:15 , Processed in 0.171734 second(s), 21 queries , Gzip On.

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