请教结构优化中溶剂用二氯甲烷不收敛的问题

lr8616

请教各位老师，我现在计算的一个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，果然有虚频

RAL

scf=conver=6对计算受力来说精度不够，应该让高斯用默认值

lr8616

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

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

fatpig

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

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

lr8616

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

谢谢指点，我试试看

lr8616

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

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

RAL

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

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

wjc404

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

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

lr8616

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

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