计算化学公社

 找回密码 Forget password
 注册 Register

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

查看数: 10038 | 评论数: 10 | 收藏 Add to favorites 1
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2019-1-15 20:20

正文摘要:

本帖最后由 shiyanan 于 2019-1-16 14:42 编辑 请问sob老师,opt+freq优化,输出文件结果显示opt收敛四个YES,freq却出现四个NO正常吗?

回复 Reply

shiyanan 发表于 Post on 2019-1-17 15:03:40
本帖最后由 shiyanan 于 2019-1-17 15:32 编辑
ymygca 发表于 2019-1-17 12:34
我没有用过自定义泛函,不知道呢。看你这些iop,是想用PBE0吗?是的话用PBE1PBE就可以了。结构用最后一步 ...

谢谢呢~这个是PBE38的写法。
ymygca 发表于 Post on 2019-1-17 12:34:55
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 发表于 Post on 2019-1-17 08:54:41
本帖最后由 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

ymygca 发表于 Post on 2019-1-17 00:27:50
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 发表于 Post on 2019-1-16 14:40:46
蓝色柠檬草 发表于 2019-1-15 21:24
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://bbs.keinsci.com/thread-633-1-1.ht ...

这个帖子里面说前两项应该一致,四个no是不是不太正常C:\Users\dell\Desktop\4
shiyanan 发表于 Post on 2019-1-16 14:27:33
feng1051leaf 发表于 2019-1-16 04:45
我对你的体系感兴趣。能放一下几何结构或者输入文件么?

上传输入文件了
shiyanan 发表于 Post on 2019-1-16 14:27:06
ymygca 发表于 2019-1-15 23:47
不太常见,但是以前遇到过。你贴个截屏让大家看看具体情况?

上传输入文件和输出文件了,谢谢呢
feng1051leaf 发表于 Post on 2019-1-16 04:45:38
我对你的体系感兴趣。能放一下几何结构或者输入文件么?
ymygca 发表于 Post on 2019-1-15 23:47:24
不太常见,但是以前遇到过。你贴个截屏让大家看看具体情况?
蓝色柠檬草 发表于 Post on 2019-1-15 21:24:58
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://bbs.keinsci.com/thread-633-1-1.html

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
sobereva + 1

查看全部评分 View all ratings

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-24 00:56 , Processed in 0.190256 second(s), 32 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list