计算化学公社

标题: 一个由有机体系（140多个原子有C,N,O,H）和过渡金属组成的体系OPT计算该怎么选泛函 [打印本页]

作者
Author: a18236822855 时间: 2020-9-30 08:52
标题: 一个由有机体系（140多个原子有C,N,O,H）和过渡金属组成的体系OPT计算该怎么选泛函
本帖最后由 a18236822855 于 2020-9-30 12:17 编辑

本人首先用了ORCA里的B97-3C尝试计算，结果SCF不收敛，后面用Gaussian用PBE0（PBE1PBE）计算（用了6-311g*的基组和D3BJ校正，8核64G算了接近13个小时仍没有SCF收敛的迹象），由于计算资源有限就停止了计算，在这里想向大家请教一下应该怎么选择一个快速有效的泛函来解决这个问题。这里的有机分子基团和过渡金属的模型是分开的，因为我想计算两者在接近时的能量。

作者
Author: sobereva 时间: 2020-9-30 19:05
先用xtb在GFN-xTB下粗略优化一下，之后再用你之前的级别搞。
初始结构偏离平衡结构较大的情况SCF注定难收敛。如果PBE0-D3(BJ)还是不太容易收敛，如果有G09，用scf=conver=6下再优化一下，最后切换到默认收敛限做最终的优化。

这些处理在以下博文里都明确介绍过，其中还有其它诸多解决办法
解决SCF不收敛问题的方法
http://sobereva.com/61

当前问题和泛函选用没有直接关系。绝对不是越适合当前体系的泛函就越容易SCF收敛。

作者
Author: 喵星大佬 时间: 2020-9-30 19:31
结构太离谱了，自己先画的合理，至少键长键角大差不差，不要有原子靠的太近之类的

作者
Author: a18236822855 时间: 2020-9-30 20:18

sobereva 发表于 2020-9-30 19:05
先用xtb在GFN-xTB下粗略优化一下，之后再用你之前的级别搞。
初始结构偏离平衡结构较大的情况SCF注定难收 ...

我懂了，谢谢sob老师

作者
Author: a18236822855 时间: 2020-9-30 20:20

喵星大佬发表于 2020-9-30 19:31
结构太离谱了，自己先画的合理，至少键长键角大差不差，不要有原子靠的太近之类的

这个应该是因为gv打开的缘故，下边的有机结构是我之前优化过的体系，后面又在上边加上了过渡金属化合物，显示是有点离谱

作者
Author: wzkchem5 时间: 2020-9-30 23:11

a18236822855 发表于 2020-9-30 20:20
这个应该是因为gv打开的缘故，下边的有机结构是我之前优化过的体系，后面又在上边加上了过渡金属化合物， ...

这就是你之前问的那个K2Cr2O7吸附的体系吗？
[Cr2O7]2-的Cr是四面体构型，不是平面四方构型，见https://en.wikipedia.org/wiki/Potassium_dichromate。改成四面体构型应该就好收敛了（前提是确保不要有距离太近的原子，尤其是钾和其他原子的距离要合理）。
另外如果你是研究水溶液里的吸附，那么必须加足够的水分子，否则结果很可能定性错误。

作者
Author: a18236822855 时间: 2020-9-30 23:58

wzkchem5 发表于 2020-9-30 23:11
这就是你之前问的那个K2Cr2O7吸附的体系吗？
[Cr2O7]2-的Cr是四面体构型，不是平面四方构型，见https:// ...

对对，这是我之前问的那个问题的延续
我之前还没意识到这个问题，非常感谢你的认真回复，我会仔细考虑一下您提出的问题的

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