高斯优化结构震荡严重

来源:量子化学  › 热点前沿问题讨论  › 倒序浏览
  • 回答数

    7

  • 浏览数

    3461

  • 收藏数

    0

作者:疯不觉 发表于 2020-5-11 13:52:28
跳转到指定楼层
震荡严重,优化不出来结构,是因为初始结构不合理吗?如果直接用这个结构跑freq,虚频比较多,但是可以得到想要虚频,该怎么解决?求解答

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册 微信登录

x
分享:
回复

使用道具

该用户从未签到

新手上路

Rank: 1

积分
73
极客币
17
主题
3
帖子
11
注册时间
2020-5-11
在线时间
1 小时
性别
保密
发表于 2020-5-11 14:52:51 | 显示全部楼层
从你这个电子能量来看,分子也不是很大。
那过渡态这个样子大概率是初始结构不够好导致的。要不可能就是该加溶剂模型的没加或者不该加溶剂模型的加了。之类的问题
回复

使用道具 举报

该用户从未签到

管理员

Rank: 9Rank: 9Rank: 9

积分
158
极客币
49
主题
12
帖子
18
注册时间
2020-4-25
在线时间
3 小时
性别
保密
 楼主| 发表于 2020-5-11 14:54:06 | 显示全部楼层
科研小宝 发表于 2020-5-11 14:52
从你这个电子能量来看,分子也不是很大。
那过渡态这个样子大概率是初始结构不够好导致的。要不可能就是该 ...

您好!是加了溶剂模型,溶剂会影响优化和过渡态结果吗?
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
73
极客币
17
主题
3
帖子
11
注册时间
2020-5-11
在线时间
1 小时
性别
保密
发表于 2020-5-11 14:55:00 | 显示全部楼层
疯不觉 发表于 2020-5-11 14:54
您好!是加了溶剂模型,溶剂会影响优化和过渡态结果吗?

溶剂模型肯定会有影响。不过90%以上的过渡态不收敛都是因为初始结构不够好导致的。 可以尝试调整初始结构。或者你先把这个过渡态两侧的中间体找到。根据中间体结构再去猜测过渡态结构。
回复

使用道具 举报

该用户从未签到

管理员

Rank: 9Rank: 9Rank: 9

积分
158
极客币
49
主题
12
帖子
18
注册时间
2020-4-25
在线时间
3 小时
性别
保密
 楼主| 发表于 2020-5-11 14:59:14 | 显示全部楼层
科研小宝 发表于 2020-5-11 14:55
溶剂模型肯定会有影响。不过90%以上的过渡态不收敛都是因为初始结构不够好导致的。 可以尝试调整初始结构 ...

老师,不是过渡态振荡,而是前一步优化就出现了振荡,气相中优化正常,加了溶剂模型,优化不出来。
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
73
极客币
17
主题
3
帖子
11
注册时间
2020-5-11
在线时间
1 小时
性别
保密
发表于 2020-5-11 14:59:51 | 显示全部楼层
疯不觉 发表于 2020-5-11 14:59
老师,不是过渡态振荡,而是前一步优化就出现了振荡,气相中优化正常,加了溶剂模型,优化不出来。
...

这个状况也不少见。可以试试别的溶剂模型看看。
回复

使用道具 举报

该用户从未签到

管理员

Rank: 9Rank: 9Rank: 9

积分
158
极客币
49
主题
12
帖子
18
注册时间
2020-4-25
在线时间
3 小时
性别
保密
 楼主| 发表于 2020-5-11 15:00:16 | 显示全部楼层
科研小宝 发表于 2020-5-11 14:59
这个状况也不少见。可以试试别的溶剂模型看看。

老师,之前的计算都是用的这个模型!如果就想用这个模型优化,还有没有其他办法?出现这种情况还是因为结果不合理吗?
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
73
极客币
17
主题
3
帖子
11
注册时间
2020-5-11
在线时间
1 小时
性别
保密
发表于 2020-5-11 15:01:17 | 显示全部楼层
疯不觉 发表于 2020-5-11 15:00
老师,之前的计算都是用的这个模型!如果就想用这个模型优化,还有没有其他办法?出现这种情况还是因为结 ...

对付几何优化不收敛。而你不想换溶剂模型那就只能用常规手段来一套。

1        OPT=GDIIS (或者OPT=RFO)

高斯的默认优化算法是OPT=GEDIIS ,在优化比较“soft” 比较巨大的体系时候可能不好用。没准换成别的优化算法可以帮助收敛。
2        OPT=CALCALL (或者 OPT=RECALC=3)
OPT=CACALL这个关键词帮助高斯在优化的步骤中进行严格的HESSIAN矩阵计算。这样会帮助收敛。但是慎用,对于原子很多的体系每一步都计算HESSIAN其实就是每一步都算一次频率一样。会让计算总时长翻十倍都不止可能。只有在非常难收敛的时候再考虑使用。OPT=RECALC=3 意思是每三步计算一次HESSIAN,会比CALCALL快很多。G16中支持。
3        OPT=MAXCYC=200
这个就是增加优化步数的上限。这个仅仅对于原子数小于10个的体系使用才有效。因为原子数目比较少的时候,高斯的默认步数太小了 (最大步数=2倍的冗余内坐标数目)
4        OPT(MAXSTEP=10, NOTRUST)   (默认值是30 * 0.01Bohr)
这个关键词是缩短优化过程每一步的最大步长。对于优化中反复出现震荡的体系可能会有效。(查看优化进程是否出现振荡,查看本文GaussView使用部分3.1)NOTRUST的意思是让优化中每一步的步长上限都是10.不写这个的话从第二个点开始步长又回到30.
5        OPT=LOOSE
对于有的很大的体系,就是很难收敛这时候可以释放放宽收敛限。但是如果要进行频率计算。则这个做法无效。因为频率需要很准的结构,一般至少是OPT=TIGHT 甚至要OPT=VERYTIGHT才行
6    以上提到的方法可以各种组合,帮助收敛。比如
        OPT(MAXSTEP=5, NOTRUST, RECALC=3,GDIIS)
  而且如果分子又高对称性也最好打破对称性去算。要不然也容易振荡。
回复

使用道具 举报

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册 微信登录