计算化学公社

标题: opt+freq，输出文件结果显示opt收敛四个YES，freq却出现四个NO正常吗？ [打印本页]

作者
Author: shiyanan 时间: 2019-1-15 20:20
标题: opt+freq，输出文件结果显示opt收敛四个YES，freq却出现四个NO正常吗？
本帖最后由 shiyanan 于 2019-1-16 14:42 编辑

请问sob老师，opt+freq优化，输出文件结果显示opt收敛四个YES，freq却出现四个NO正常吗？

作者
Author: 蓝色柠檬草 时间: 2019-1-15 21:24
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://bbs.keinsci.com/thread-633-1-1.html

作者
Author: ymygca 时间: 2019-1-15 23:47
不太常见，但是以前遇到过。你贴个截屏让大家看看具体情况？

作者
Author: feng1051leaf 时间: 2019-1-16 04:45
我对你的体系感兴趣。能放一下几何结构或者输入文件么？

作者
Author: shiyanan 时间: 2019-1-16 14:27

ymygca 发表于 2019-1-15 23:47
不太常见，但是以前遇到过。你贴个截屏让大家看看具体情况？

上传输入文件和输出文件了，谢谢呢

作者
Author: shiyanan 时间: 2019-1-16 14:27

feng1051leaf 发表于 2019-1-16 04:45
我对你的体系感兴趣。能放一下几何结构或者输入文件么？

上传输入文件了

作者
Author: shiyanan 时间: 2019-1-16 14:40

蓝色柠檬草发表于 2019-1-15 21:24
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://bbs.keinsci.com/thread-633-1-1.ht ...

这个帖子里面说前两项应该一致，四个no是不是不太正常C:\Users\dell\Desktop\4

作者
Author: ymygca 时间: 2019-1-17 00:27

shiyanan 发表于 2019-1-16 14:27
上传输入文件和输出文件了，谢谢呢

优化过程没看毛病，是你的iop自定义的这些mix HF/DFT, local/non-local的毛病。你用grep查看每一步的能量，可以看到1-18步都是PBE/HF的能量，但是最后一步频率计算是用的PBE/PBE的能量，差别当然大了。你把第18步的几何结构提出来，用同样的方法（那3个iop）单独做一个频率计算（不要做opt+freq, 单独的freq)，应该就解决了。
grep "SCF Done" trans-pic-s0-pbe38.log
SCF Done:  E(RPBE+HF-PBE) =  -840.642848083    A.U. after 12 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.666107820    A.U. after 11 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.670674391    A.U. after 11 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.674449932    A.U. after 11 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.679344990    A.U. after 12 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.681533426    A.U. after 10 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.683231045    A.U. after 10 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.683611353    A.U. after 10 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.684416858    A.U. after 11 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.684914725    A.U. after 13 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.685657731    A.U. after 10 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686355350    A.U. after 11 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686501490    A.U. after 10 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686555783    A.U. after 9 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686565281    A.U. after 9 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686574177    A.U. after 8 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686574449    A.U. after 7 cycles
SCF Done:  E(RPBE+HF-PBE) =  -840.686574485    A.U. after 6 cycles
SCF Done:  E(RPBE-PBE) =  -840.600461256    A.U. after 12 cycles

作者
Author: shiyanan 时间: 2019-1-17 08:54
本帖最后由 shiyanan 于 2019-1-17 09:47 编辑

ymygca 发表于 2019-1-17 00:27
优化过程没看毛病，是你的iop自定义的这些mix HF/DFT, local/non-local的毛病。你用grep查看每一步的能量 ...

谢谢呢，我试一下~是不是自定义的这些泛函做opt+freq都会出现这种情况呀？您帮忙看一下我这个路径写的对吗？#p freq PBEPBE IOP(3/76=1000003750) IOP(3/77=0625006250) IOP(3/78=1000010000)/6-31G* scrf=(solvent=methanol) geom=(check,step=12) guess=read

作者
Author: ymygca 时间: 2019-1-17 12:34

shiyanan 发表于 2019-1-17 08:54
谢谢呢，我试一下~是不是自定义的这些泛函做opt+freq都会出现这种情况呀？您帮忙看一下我这个路径写的对 ...

我没有用过自定义泛函，不知道呢。看你这些iop，是想用PBE0吗？是的话用PBE1PBE就可以了。结构用最后一步，step=18.
#p freq PBE1PBE/6-31G* scrf=(solvent=methanol) geom=(allcheck,step=18) guess=read

作者
Author: shiyanan 时间: 2019-1-17 15:03
本帖最后由 shiyanan 于 2019-1-17 15:32 编辑

ymygca 发表于 2019-1-17 12:34
我没有用过自定义泛函，不知道呢。看你这些iop，是想用PBE0吗？是的话用PBE1PBE就可以了。结构用最后一步 ...

谢谢呢~这个是PBE38的写法。

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