计算化学公社

标题: 请教结构优化中溶剂用二氯甲烷不收敛的问题 [打印本页]

作者
Author: lr8616 时间: 2023-2-9 09:26
标题: 请教结构优化中溶剂用二氯甲烷不收敛的问题
本帖最后由 lr8616 于 2023-2-13 09:34 编辑

请教各位老师，我现在计算的一个160个原子、带4个正电荷的结构，溶剂用的二氯甲烷，结构优化的时候三个YES，第三个Maximum Displacement总是No，在0.03左右波动，例如：
      Item             Value    Threshold  Converged?
Maximum Force          0.002398    0.002500    YES
RMS    Force          0.000309    0.001667    YES
Maximum Displacement    0.023583    0.010000    NO
RMS    Displacement    0.004722    0.006667    YES
Predicted change in Energy=-7.561671D-05

根据老师帖子中提到的方法，在真空和DMSO中做了优化都可以收敛，用这两个的收敛的波函数作为初猜，二氯甲烷溶剂还是一样Maximum Displacement No。现在的输入文件如下：

#p opt=(loose,gdiis,maxstep=2,notrust) b3lyp/6-31g(d) scrf=(solvent=chloroform,pcm) guess=read em=gd3bj iop(7/127=-99) scf=(noincfock,novaracc,conver=6)

折腾了一个星期了还没搞定，不知道该怎么试了，还请各位老师指点一二，谢谢！！！

【新补充】根据2楼和4楼的建议，去掉了scf=conver=6，用默认的，其它关键词不变，但还是Maximum Displacement不收敛，最小的时候大概在0.013，试着计算了一下freq，果然有虚频

作者
Author: RAL 时间: 2023-2-9 09:39
scf=conver=6对计算受力来说精度不够，应该让高斯用默认值

作者
Author: lr8616 时间: 2023-2-9 09:47

RAL 发表于 2023-2-9 09:39
scf=conver=6对计算受力来说精度不够，应该让高斯用默认值

我是想先初步优化一下，结构收敛之后再换默认优化+freq

作者
Author: fatpig 时间: 2023-2-9 09:56

lr8616 发表于 2023-2-9 09:47
我是想先初步优化一下，结构收敛之后再换默认优化+freq

初步loose优化的话，这个结果已经可以用于进行下一步了，没必要纠结...

作者
Author: lr8616 时间: 2023-2-9 10:36

fatpig 发表于 2023-2-9 09:56
初步loose优化的话，这个结果已经可以用于进行下一步了，没必要纠结...

谢谢指点，我试试看

作者
Author: lr8616 时间: 2023-2-13 09:35

RAL 发表于 2023-2-9 09:39
scf=conver=6对计算受力来说精度不够，应该让高斯用默认值

去掉了还是不收敛，不知道还有啥办法呀？

作者
Author: RAL 时间: 2023-2-13 09:58

lr8616 发表于 2023-2-13 09:35
去掉了还是不收敛，不知道还有啥办法呀？

四楼的意思是这个结构可以做正常收敛限的opt+freq了，你不会直接做freq了吧

作者
Author: wjc404 时间: 2023-2-13 10:04
本帖最后由 wjc404 于 2023-2-13 10:19 编辑

lr8616 发表于 2023-2-13 09:35
去掉了还是不收敛，不知道还有啥办法呀？

http://sobereva.com/164这里面提到的recalc选项尝试了没？如果recalc甚至calcall都不能收敛，就要仔细看一下分子结构中哪些基团在优化时一直在显著变动位置，模型搭建是否有不合理之处。
iop(7/127=-99)不建议用。如果发现scf收敛（不是opt收敛）比预期要困难，或者所研究体系不是闭壳层，做优化前最好看一下波函数稳定性（输入文件中关键词去掉opt的内容，然后加上stable=opt）。

作者
Author: lr8616 时间: 2023-2-13 10:18

wjc404 发表于 2023-2-13 10:04
http://sobereva.com/164这里面提到的recalc选项尝试了没？如果recalc甚至calcall都不能收敛，就要仔细看 ...

是一个和C60复合的结构，优化到现在的程度每一步结构变化都非常非常小了。还没有试recalc，不是太懂这对opt不收敛的情况帮助大吗？之前遇到opt四个YES，但是freq有NO的时候用recalc效果比较好。

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