计算化学公社

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

[ORCA] ORCA做几何优化的问题总是得到一个很扭曲的结构

[复制链接 Copy URL]

39

帖子

0

威望

3203

eV
积分
3242

Level 5 (御坂)

16#
发表于 Post on 2020-3-23 10:43:47 | 只看该作者 Only view this author
我看到过很多共轭分子的工作中,拿DFT优化的系列化合物。共轭片段之间的二面角差个几度到十几度来说明共轭性的问题,连二面角旋转势能面都不给,这基本没有什么实际意义,尤其是B3LYP严重高估共轭片段之间的平面性。

6万

帖子

99

威望

6万

eV
积分
125127

管理员

公社社长

17#
发表于 Post on 2020-3-23 13:46:06 | 只看该作者 Only view this author
ORCA几何优化功能本身做得并没有什么明显问题,虽然SCF和几何优化收敛性肯定不如Gaussian
北京科音自然科学研究中心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

65

帖子

1

威望

589

eV
积分
674

Level 4 (黑子)

18#
 楼主 Author| 发表于 Post on 2020-3-23 23:38:54 | 只看该作者 Only view this author
liyuanhe211 发表于 2020-3-23 01:12
你这个问题有几个因素的影响:
  • 泛函不同(含相同名字的泛函定义不同)

  • 感谢解答,排查了一天,问题找到了,我把 RI去了就正常了,得到了最终稳定的结构。我其实一开始把RI这个因素放在了后面去排查,因为b97-3c应该也是自动开RI的,虽然这两个RI不一样。所以下意识地觉得RI应该没问题。没想到最终是RI导致的。(所以可能是COSX引起的?)
    我还有个小问题就是那么对于这个分子,我后续在单点计算的时候还要不要开RI呢,会有影响吗?

    6万

    帖子

    99

    威望

    6万

    eV
    积分
    125127

    管理员

    公社社长

    19#
    发表于 Post on 2020-3-25 22:06:45 | 只看该作者 Only view this author
    qinsc 发表于 2020-3-23 23:38
    感谢解答,排查了一天,问题找到了,我把 RI去了就正常了,得到了最终稳定的结构。我其实一开始把RI这个 ...

    单点开RI很稳妥
    应当区分RIJK和RIJCOSX,如果之前用的RIJCOSX,加大gridx再试,或者尝试RIJK
    北京科音自然科学研究中心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

    65

    帖子

    1

    威望

    589

    eV
    积分
    674

    Level 4 (黑子)

    20#
     楼主 Author| 发表于 Post on 2020-3-26 21:26:25 | 只看该作者 Only view this author
    sobereva 发表于 2020-3-25 22:06
    单点开RI很稳妥
    应当区分RIJK和RIJCOSX,如果之前用的RIJCOSX,加大gridx再试,或者尝试RIJK

    感谢sob老师的解答,谢谢!
    我尝试过grid4 gridx4 finalgrid5,没有得到预期的结构,我会试一下RIJK看看效果如何,谢谢!

    65

    帖子

    1

    威望

    589

    eV
    积分
    674

    Level 4 (黑子)

    21#
     楼主 Author| 发表于 Post on 2020-4-13 19:39:12 | 只看该作者 Only view this author
    sobereva 发表于 2020-3-25 22:06
    单点开RI很稳妥
    应当区分RIJK和RIJCOSX,如果之前用的RIJCOSX,加大gridx再试,或者尝试RIJK

    sob老师好
    今天在查看手册的时候意外发现了这段话,虽然是在MP2部分的,但和我出现的问题几乎一摸一样,都是平面共轭结构优化会扭曲。手册中提到在几何结构优化的时候即使交换部分积分精度提高到Gridx6,仍然会有扭曲的问题。我有一个小问题,假设我把精度提高到Gridx8甚至9避免了问题,那我相较不开加速还会获得很大程度的时间节省吗?如果节省的不多,那还不如不开加速,避免误差的引入。
    手册建议可以使用RIJONX,我试了一下,是没有问题的,不会出现结构扭曲问题。
    RIJK的话因为我的体系比较大,大概150个原子左右,优化的时候基组在def2-SVP,我看官方手册有个比较表格,在这个基组下,大体系RIJK甚至比default用时还要多一点,所以我感觉还不如不开,不知道我这么想对不对。

    最近我做了很多类似多稠环共轭的体系,几乎都出现了这个问题,所以我觉得这是一个值得注意的问题。
    另外虽然我提问中的体系使用B97-3c可以优化,但在第一次SCF迭代中足足花了123步才收敛,另外的一些类似体系在默认的125次迭代中完全没有收敛迹象,而其他杂化泛函虽然结构不对,但每一次的SCF迭代还是很快的,基本15步左右就可以收敛。不知道是纯泛函的问题还是B97-3c的问题,我有时间用BLYP试一下。

    另外我还有一个小问题,对于多稠环平面共轭体系,进行结构优化的时候,我发现wB97X-D3优化出来的结构单双键更分明一点(单键更长,双键更短),B3LYP优化出来的单双键更不分明(键长更平均)。我在选择泛函的时候对这一类结构哪个优化出来的结构更会贴近真实结构呢,我看您在有机结构优化中推荐PBE0,这个泛函对于这类体系会比wB97X-D3和B3LYP更好吗?

    Annotation 2020-04-13 191053.jpg (378.13 KB, 下载次数 Times of downloads: 29)

    Annotation 2020-04-13 191053.jpg

    6万

    帖子

    99

    威望

    6万

    eV
    积分
    125127

    管理员

    公社社长

    22#
    发表于 Post on 2020-4-15 00:31:04 | 只看该作者 Only view this author
    qinsc 发表于 2020-4-13 19:39
    sob老师好
    今天在查看手册的时候意外发现了这段话,虽然是在MP2部分的,但和我出现的问题几乎一摸一样, ...

    可以实测一下gridx8,应该还是比不用RIJCOSX快。到了这种格点精度,COSX带来的误差也算可以忽略不计了

    多稠环共轭偏离平面弯曲的力常数非常小,8成是导致平面性对积分格点很敏感的主要原因

    对于有机体系,我极少碰到B97-3c出现SCF不收敛的。但B97-3c毕竟是纯泛函,SCF收敛难度必然整体高于杂化泛函。如果你算得动,可以用PBEh-3c,整体精度更好一些,也应该比B97-3c更容易收敛,但也更贵。

    B3LYP算大范围pi共轭众所周知有你提到的单双键均衡化过度的问题。这种情况明显wB97X-D3对结构描述得更好。换句话说,这类情况更适合用HF成分更高的泛函。18碳环体系用B3LYP优化完全失败也体现出这一点,本质大抵在于SIE问题,见此文介绍的文章里的专门的讨论
    一篇最全面、系统的研究新颖独特的18碳环的理论文章
    http://sobereva.com/524http://bbs.keinsci.com/thread-15792-1-1.html
    北京科音自然科学研究中心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

    65

    帖子

    1

    威望

    589

    eV
    积分
    674

    Level 4 (黑子)

    23#
     楼主 Author| 发表于 Post on 2020-4-15 14:39:30 | 只看该作者 Only view this author
    本帖最后由 qinsc 于 2020-4-15 15:05 编辑
    sobereva 发表于 2020-4-15 00:31
    可以实测一下gridx8,应该还是比不用RIJCOSX快。到了这种格点精度,COSX带来的误差也算可以忽略不计了

    ...

    感谢老师百忙之中的回复。
    我看了您18碳环的文章,对于PBE0,B3LYP这一类HF成分比较低的泛函,对于这类体系的处理都不是很理想。正好看到您说也可以用PBEh-3c,我想到PBEh-3c的HF成分达到了42%,会不会描述这类体系还不错,所以我用PBEh-3c优化了一下这个18碳环,最后也是得到了键长交替的结构,不过不如wB97XD那么分明,键长1.234783/1.327069,差值0.092286。能不能说这个PBEh-3c其实还不错(另外我在计算的时候发现PBEh-3c会自动开RIJONX)。

    由于我这个体系wB97XD3结合def2-TZVP完全算不动,所以如果在wB97XD3结合def2-SVP与PBEh-3c之间做选择的话是不是还是wB97XD3好一些?不过同样开RIJONX,对于我这个分子,PBEh-3c比wB97XD3用时大概快了30%。综合来看PBEh-3c是不是也是个不错的选择呢?


    我还有一个小问题,结构优化完之后做一些轨道分析例如LOL分析pi轨道,我做单点的时候用过wB97XD3(直接使用结构优化的结果),B3LYP-D3,PWPB95-D3,B2GP-PLYP-D3,我发现这些给出的LOL图不太相同,后两个双杂泛函会在一些奇怪的地方出现等值面(比如环外杂原子在平行环平面的p轨道而非垂直环平面),而B2GP-PLYP在主要考察区域和其他几个泛函差别比较大,具体表现在有些等值面覆盖明显更大。所以我的小问题是对于这一类轨道分析,有没有哪一类泛函适合呢,或者说需要和结构优化时的泛函保持一致?

    另外就是我用wB97XD3的结果的话,我的HOMO-LUMO gap有些偏大,我使用B3PW91单独计算它的gap会好很多,这是不是一个合理的切换泛函的理由呢。用HF成分高的优化结构好,但是会高估它的HOMO-LUMO gap或是吸收gap,这会不会有些矛盾。

    6万

    帖子

    99

    威望

    6万

    eV
    积分
    125127

    管理员

    公社社长

    24#
    发表于 Post on 2020-4-16 13:38:32 | 只看该作者 Only view this author
    qinsc 发表于 2020-4-15 14:39
    感谢老师百忙之中的回复。
    我看了您18碳环的文章,对于PBE0,B3LYP这一类HF成分比较低的泛函,对于这类 ...

    C18这个体系比较极端。PBEh-3c对于C18还是差了点,不过对于描述普通有机共轭体系完全够了。

    PBEh-3c比wB97X-D3/def2-SVP结果有较大概率更好。而且本身PBEh-3c对应的基组质量就更高。

    我不知道你用双杂化的时候具体是怎么产生波函数的。要注意默认情况下双杂化算完之后,给你的波函数实际上只是双杂化泛函中杂化泛函那一部分的波函数,这个波函数质量几乎一定不如常规的泛函的波函数。如果你需要用真正的双杂化的波函数做波函数分析,需要产生自然轨道,而这需要多花很多时间,还需要额外步骤,而且此时还没法做轨道定域化。所以一般情况下完全没必要用双杂化波函数作分析,费力不讨好。你就用做结构优化时候的级别做波函数分析就够了。

    算HOMO-LUMO gap的时候完全可以切换成B3LYP或B3PW91,这有理有据。
    北京科音自然科学研究中心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

    65

    帖子

    1

    威望

    589

    eV
    积分
    674

    Level 4 (黑子)

    25#
     楼主 Author| 发表于 Post on 2020-4-16 14:50:05 | 只看该作者 Only view this author
    sobereva 发表于 2020-4-16 13:38
    C18这个体系比较极端。PBEh-3c对于C18还是差了点,不过对于描述普通有机共轭体系完全够了。

    PBEh-3c比 ...

    明白了!谢谢老师耐心解答!

    本版积分规则 Credits rule

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

    GMT+8, 2026-2-18 22:07 , Processed in 0.263166 second(s), 22 queries , Gzip On.

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