计算化学公社

标题: 用opt=calcall，进行优化后的情况请教 [打印本页]

作者
Author: 北极天99 时间: 2015-5-4 15:21
标题: 用opt=calcall，进行优化后的情况请教
有两个问题需要请教各位老师：
1.我把物质用opt=calcall这个方式进行优化，最后依然出现两个NO，是不是就证明这个构型不存在？
2.对于其他一些构型我也用opt=calcall，这样优化了之后，out文件里，可以看到震动频率（但是我的输入文件里只是做了优化，并没有写上freq），这样正常吗？

作者
Author: sobereva 时间: 2015-5-5 01:20
1 这么说是错误的。可能只是你用的初猜或者优化方式不合理导致优化没完成，这是数值上的问题，不能从化学意义上说明什么。
并且要弄清楚到底是什么情况，是最后根本没收敛，还是说后两个判据是NO但是受力已小于判据100倍而被认为收敛了。

2 可能高斯根据calcall最后一步的Hessian自动进行了振动分析（因为没有额外耗时），此时的振动频率是没问题的。

作者
Author: 北极天99 时间: 2015-5-5 10:20

sobereva 发表于 2015-5-5 01:20
1 这么说是错误的。可能只是你用的初猜或者优化方式不合理导致优化没完成，这是数值上的问题，不能从化学意 ...

谢谢sob老师。我计算的出现两个No是最后两个判据微NO，而且结果中也出现了虚频。 (, 下载次数 Times of downloads: 30)

作者
Author: 北极天99 时间: 2015-5-5 10:29
补充一下哈，初猜结构是先用B3LYP优化过一次的，不过当时也是计算正常结束，结果中后两个判据是NO但是受力已小于判据100倍而被认为收敛了。这次我用的方法是# M062X/6-311+G(d,p) opt=calcall，结果还是后两个判据是NO但是受力已小于判据100倍而被认为收敛了，同时出现了虚频。我优化的是丙醇二聚体。谢谢sob老师。

作者
Author: sobereva 时间: 2015-5-5 12:11
用此关键词，看看最后有无虚频、最后一次输出的四个收敛判断标准各是多少。
# M062X/6-311+G(d,p) int=ultrafine opt(tight,calcall) freq

作者
Author: 北极天99 时间: 2015-5-5 14:59

sobereva 发表于 2015-5-5 12:11
用此关键词，看看最后有无虚频、最后一次输出的四个收敛判断标准各是多少。
# M062X/6-311+G(d,p) int=ult ...

好的，我试试。谢谢Sob老师。

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