opt+freq，输出文件结果显示opt收敛四个YES，freq却出现四个NO正常吗？

shiyanan

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

蓝色柠檬草

Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://bbs.keinsci.com/thread-633-1-1.html

ymygca

不太常见，但是以前遇到过。你贴个截屏让大家看看具体情况？

feng1051leaf

我对你的体系感兴趣。能放一下几何结构或者输入文件么？

shiyanan

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

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

shiyanan

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

上传输入文件了

shiyanan

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

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

ymygca

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

shiyanan

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

ymygca

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

shiyanan

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

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