求助：为什么同一分子结构，添加溶剂模型时 SCF 无法收敛，但不加溶剂计算成功

2329568177 · 发表于 Post on 3 day ago

本帖最后由 2329568177 于 2026-5-9 14:36 编辑

想请教各位老师：我对同一个结构进行优化时，加入 SMD 水溶剂后 SCF 收敛失败，去掉溶剂模型做气相计算则可以顺利收敛，想请问出现这种差异的原因是什么？
以及是否可以参考sob老师在解决SCF不收敛问题的方法这篇博文中提到的方法解决：(16) 对于溶剂下的计算若SCF不收敛，也可以试试其它溶剂下或者真空下或者其它溶剂模型下能否收敛，如果可以，将收敛的波函数当做原溶剂下计算的初猜。
         第一次                                                                               第二次
! OPT B3LYP D3BJ def2-SVP RIJCOSX                         ! OPT B3LYP D3BJ def2-SVP RIJCOSX
%pal nprocs 16 end                                                 %pal nprocs 16 end
%MaxCore 3200                                                       %MaxCore 3200
%cpcm                                                                   * xyz 0 1
  smd true                                                                坐标
  SMDsolvent "water"                                              *
end
* xyz 0 1
坐标
*
第一次输出日志

第二次输出日志

而且一直有这个提示，是否对运算有影响：mca_base_component_repository_open: unable to open mca_btl_openib: libosmcomp.so.3: cannot open shared object file: No such file or directory (ignored)

iamjjLin · 发表于 Post on 3 day ago

换pcm溶剂模型试试，smd溶剂模型不收敛的概率挺高的。

2329568177 · 发表于 Post on 3 day ago

iamjjLin 发表于 2026-5-9 15:49
换pcm溶剂模型试试，smd溶剂模型不收敛的概率挺高的。

老师，我是直接用目前能量最低的构想换成pcm去跑吗

iamjjLin · 发表于 Post on 3 day ago

2329568177 发表于 2026-5-9 15:51
老师，我是直接用目前能量最低的构想换成pcm去跑吗

不加溶剂模型优化后的结构就行

2329568177 · 发表于 Post on 3 day ago

iamjjLin 发表于 2026-5-9 15:57
不加溶剂模型优化后的结构就行

谢谢老师

2329568177 · 发表于 Post on 3 day ago

iamjjLin 发表于 2026-5-9 15:57
不加溶剂模型优化后的结构就行

老师，我想请教对于中等体系（255原子）做几何优化和振动分析，def2-svp够用吗，有必要用到def2-TZVP吗

陈AG · 发表于 Post on 3 day ago

2329568177 发表于 2026-5-9 16:32
老师，我想请教对于中等体系（255原子）做几何优化和振动分析，def2-svp够用吗，有必要用到def2-TZVP吗

足够了

2329568177 · 发表于 Post on 3 day ago

陈AG 发表于 2026-5-9 16:44
足够了

谢谢老师

sobereva · 发表于 Post on 前天 15:34

2329568177 发表于 2026-5-9 16:32
老师，我想请教对于中等体系（255原子）做几何优化和振动分析，def2-svp够用吗，有必要用到def2-TZVP吗

def2-TZVP对这么大体系做振动分析想都别想，本来ORCA振动分析就慢（这方面明显不如Gaussian），Gaussian用def2-SVP跑这么大体系振动分析都相当耗时。255个原子完全是很大体系而非中等体系

前面的回复有误。SMD下往往因为势能面噪音导致几何优化比IEFPCM或CPCM更难收敛，但对于SCF收敛并没有SMD更困难一说。更何况ORCA里的SMD的非极性部分本来就是CPCM，而非极性部分只影响总能量而和SCF收敛完全无关。

ORCA里SCF不收敛该怎么解决怎么解决、该试什么试什么，带不带溶剂模型的情况都一样（PS：北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）专门讲ORCA的一节花了10页左右幻灯片专门讲ORCA怎么对付SCF不收敛）

2329568177 · 发表于 Post on 前天 15:37

sobereva 发表于 2026-5-10 15:34
def2-TZVP对这么大体系做振动分析想都别想，本来ORCA振动分析就慢（远不如Gaussian），Gaussian用def2-SV ...

谢谢老师，我用def2-svp做的优化几何，但是MAX step是no，请问老师这算收敛完成吗，可以做振动分析吗

sobereva · 发表于 Post on 前天 15:39

2329568177 发表于 2026-5-10 15:37
谢谢老师，我用def2-svp做的优化几何，但是MAX step是no，请问老师这算收敛完成吗，可以做振动分析吗

算。离收敛限相差甚微了

2329568177 · 发表于 Post on 前天 15:49

sobereva 发表于 2026-5-10 15:39
算。离收敛限相差甚微了

谢谢老师，那我就做振动分析了，几何优化用了7个半小时，振动分析会慢很多吗（只能用免费的orca哈哈哈），而且输入文件是xxx.xyz文件做振动分析会自动保留.hess文件吗
! Freq B3LYP D3BJ def2-SVP RIJCOSX
%pal nprocs 96 end
%MaxCore 3200
%cpcm
smd true
SMDsolvent "water"
end
* xyz 3 1

sobereva · 发表于 Post on 前天 16:11

2329568177 发表于 2026-5-10 15:49
谢谢老师，那我就做振动分析了，几何优化用了7个半小时，振动分析会慢很多吗（只能用免费的orca哈哈哈） ...

跟几何优化耗时相比如何，直接取决于你之前优化了多少步收敛。
振动分析你就硬着头皮计算就完了，算不动的话可以尝试改成基于纯泛函的B97-3c做opt freq，耗时更低

2329568177 · 发表于 Post on 前天 17:09

sobereva 发表于 2026-5-10 16:11
跟几何优化耗时相比如何，直接取决于你之前优化了多少步收敛。
振动分析你就硬着头皮计算就完了，算不动 ...

老师您好，我又重新学习了您的博文（http://sobereva.com/soft/Sobtop，例1、例2），我只需要构建非标准氨基酸残基的力场并进行MD模拟，在几何优化完成后，用 def2-TZVP 基组做单点计算，得到优化后的 xyz 结构和 gbw 波函数文件，再通过 Sobtop 结合 Multiwfn 生成力场文件即可，根本不需要再做振动分析。

sobereva · 发表于 Post on 前天 22:15

2329568177 发表于 2026-5-10 17:09
老师您好，我又重新学习了您的博文（http://sobereva.com/soft/Sobtop，例1、例2），我只需要构建非标准 ...

这种体系完全不需要做振动分析，根据原子类型匹配的bond/angle/dihedral参数通常在AMBER力场里一般都有现成的。就算没有现成的，bond和angle都可以猜个大小差不多的力常数凑合，柔性二面角需要势能面扫描拟合成周期性扭转势，这也不依靠振动分析

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] 求助：为什么同一分子结构，添加溶剂模型时 SCF 无法收敛，但不加溶剂计算成功

浏览过的版块