计算化学公社

标题: 求助使用高斯进行优化时频繁报错的解决方法 [打印本页]

作者
Author: yyllxx 时间: 2021-7-13 10:53
标题: 求助使用高斯进行优化时频繁报错的解决方法
本帖最后由 yyllxx 于 2021-7-13 11:26 编辑

请教各位老师，在使用M06-2X方法对七分子团簇进行结构优化的过程中频繁报错，已经进行了两次计算，多次报错，以l103错误为主，第一次优化到200多步感觉无果停止结算（如图1，优化结果文件如1），第二次优化到500多步震荡感觉无果停止结算（如图2，优化结果文件如2）。图3是第三次计算过程中的第一次报错，老师这个错误还是按l103的方法处理吗？对这个分子团簇使用B3LYP-D3方法优化与频率计算过程非常顺利，为什么M06方法频繁报错甚至无果呢？
抱歉，文件1和文件2因为太大，上传不了，不知道我是否表达清楚我的问题与困惑，如果表达不清楚，我随时回复各位老师的疑问，谢谢你们！另外，我大概说一些我课题内容，我是研究太赫兹光谱的，最终研究的是弱相互作用，我用的方法是静电势分析+范德华势分析，贴几张已发表论文的图（图4,5,6,7）。图7是回答审稿人提问的为什么不使用其他计算方法给的结果，M06确实效果不太好，但是他可以解释B3LYP-D3解释不了的吸收峰，所以我还是想用M06方法计算现在的物质

作者
Author: yyllxx 时间: 2021-7-13 10:54
本帖最后由 yyllxx 于 2021-7-13 11:31 编辑

因为现在很多文章还是解释不了太赫兹吸收峰与弱相互作用的关系，所以我想通过不同的计算方法（计算方法的侧重点不一样）来初步确认吸收峰与弱相互作用的关系，比如图6中B3LYP计算的1.16THz处的吸收峰没有出现在加色散校正的B3LYP的光谱中，所以我初步判定，1.19处的实验吸收峰，静电作用占主导。类似这样的想法，最后通过其他方法确认这个吸收峰确实是静电作用引起的，我执着与M06也是这样的原因，M06可以很好地描述范德华作用，那么是否可以解释一些范德华作用引起的吸收峰呢

作者
Author: qlyfe 时间: 2021-7-13 11:09
看sob老师的量子化学计算中帮助几何优化收敛的常用方法：http://sobereva.com/164
也不知道你的输入关键词是什么样的，用g09 m06-2x泛函首先考虑加上int=ultrafine
还是震荡不收敛的话再考虑使用opt=calcfc或者calcall吧

作者
Author: yyllxx 时间: 2021-7-13 11:35
本帖最后由 yyllxx 于 2021-7-13 11:38 编辑

我现在课题出现的问题是不管哪个方法计算出的光谱都与实验光谱有差大的差别，也没有很大的优势体现出实验光谱中的某个吸收峰确实有某种作用引起。请教各位老师，有什么方法提高光谱吻合度吗，并且能初步确认吸收峰的来源。图中是现在的课题计算的结果

作者
Author: wzkchem5 时间: 2021-7-13 15:05

yyllxx 发表于 2021-7-13 04:35
我现在课题出现的问题是不管哪个方法计算出的光谱都与实验光谱有差大的差别，也没有很大的优势体现出实验 ...

搞不好是因为非谐性效应，谐振近似做不准

作者
Author: sobereva 时间: 2021-7-13 16:43
分清楚M06和M06-2X。不要用前者

用M06-2X不是好选择，计算慢，而且对积分格点敏感，对这种柔性体系几何优化收敛费劲还容易出虚频。

说清楚用的什么基组

尝试把团簇再取大一圈。或者用CP2K算周期性体系的振动分析

作者
Author: yyllxx 时间: 2021-7-13 20:41

wzkchem5 发表于 2021-7-13 15:05
搞不好是因为非谐性效应，谐振近似做不准

谢谢您的思路，我应该考虑频率校正因子的。

作者
Author: yyllxx 时间: 2021-7-13 20:42

qlyfe 发表于 2021-7-13 11:09
看sob老师的量子化学计算中帮助几何优化收敛的常用方法：http://sobereva.com/164
也不知道你的输入关键词 ...

好嘞，谢谢你呀，我加关键词再试试

作者
Author: wzkchem5 时间: 2021-7-13 21:48

yyllxx 发表于 2021-7-13 13:41
谢谢您的思路，我应该考虑频率校正因子的。

这个应该不是单纯一个频率校正因子能解决的，频率校正因子只会整体缩放光谱的x轴，不会改变光谱的形状。可能需要考虑AIMD之类的方法

作者
Author: yyllxx 时间: 2021-7-15 11:40

wzkchem5 发表于 2021-7-13 21:48
这个应该不是单纯一个频率校正因子能解决的，频率校正因子只会整体缩放光谱的x轴，不会改变光谱的形状。 ...

非常感谢您的耐心解释，我多了解一下其他方法，比如AIMD等

欢迎光临计算化学公社 (http://bbs.keinsci.com/)