本帖最后由 shiyanan 于 2019-1-17 15:32 编辑 ymygca 发表于 2019-1-17 12:34 谢谢呢~这个是PBE38的写法。 |
shiyanan 发表于 2019-1-17 08:54 我没有用过自定义泛函,不知道呢。看你这些iop,是想用PBE0吗?是的话用PBE1PBE就可以了。结构用最后一步,step=18. #p freq PBE1PBE/6-31G* scrf=(solvent=methanol) geom=(allcheck,step=18) guess=read |
本帖最后由 shiyanan 于 2019-1-17 09:47 编辑 ymygca 发表于 2019-1-17 00:27 谢谢呢,我试一下~是不是自定义的这些泛函做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 |
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 |
蓝色柠檬草 发表于 2019-1-15 21:24 这个帖子里面说前两项应该一致,四个no是不是不太正常C:\Users\dell\Desktop\4 |
feng1051leaf 发表于 2019-1-16 04:45 上传输入文件了 |
ymygca 发表于 2019-1-15 23:47 上传输入文件和输出文件了,谢谢呢 |
| 我对你的体系感兴趣。能放一下几何结构或者输入文件么? |
| 不太常见,但是以前遇到过。你贴个截屏让大家看看具体情况? |
|
Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法 http://bbs.keinsci.com/thread-633-1-1.html |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 |
手机版 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.