求助！结构依然震荡！已黔驴技穷。。

lei234 · 发表于 Post on 2015-10-15 20:22:01

我知道结构震荡是老掉牙的问题了，但我尝试了sob老师博文中的办法，依然不收敛，还是在震荡，试了noDIIS,conver=6,Xqc,iop(1/8=7)、=3甚至1。过渡态用了gdiis,ts,nrscal,noeigen，用过readfc也用了calcfc。。。。实在没办法了，在此向各位求助，小女子不甚感激！结构如图。（醇分子中的O-H键断裂，生成醛，H与下面的OH结合生成H2O）

lao7 · 发表于 Post on 2015-10-15 21:11:41

关闭对称性试试

liyuanhe211 · 发表于 Post on 2015-10-15 21:20:26

过渡态震荡尝试换个（精度差别不特别大的）基组/方法来试一下，如果成功，以输出为下一步的输入，用你想要的水平重算一下。有时有效。

遇到过加上GD3BJ或者GD3有帮助的

sobereva · 发表于 Post on 2015-10-15 21:30:14

没有正确领会博文
noDIIS,conver=6,Xqc这根本没意义，这是处理SCF不收敛的，跟几何优化毫无关系。而且SCF收敛限设宽了可能更不容易收敛
1/8应当改用maxstep来设，而且光这么用没用，必须加notrust才有用。=3~5就足够了，再小没意义。
readfc跟解决不收敛根本没直接关系
calcfc只是前几步管事，很快就又没用了，要一直管用就得calcall
nrscal没意义，本来就是默认的
而且还可以放宽到loose收敛限，对于极小点优化loose往往已经够了

博文的方法一定有效，只是你用的不对。常用的解决优化不收敛的组合在文末都写明了
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

lei234 · 发表于 Post on 2015-10-16 08:19:04

sobereva 发表于 2015-10-15 21:30
没有正确领会博文
noDIIS,conver=6,Xqc这根本没意义，这是处理SCF不收敛的，跟几何优化毫无关系。而且SCF ...

是的，其实还是没有理解这几个命令的真正含义，只是看到什么加什么罢了。。。谢谢老师的教导！

lei234 · 发表于 Post on 2015-10-16 08:20:07

liyuanhe211 发表于 2015-10-15 21:20
过渡态震荡尝试换个（精度差别不特别大的）基组/方法来试一下，如果成功，以输出为下一步的输入，用你想要 ...

好的。谢谢您！

lei234 · 发表于 Post on 2015-10-16 08:23:33

lao7 发表于 2015-10-15 21:11
关闭对称性试试

行

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

[Gaussian/gview] 求助！结构依然震荡！已黔驴技穷。。

评分 Rate

浏览过的版块