计算化学公社

标题: 求助-PBE38优化-频率分析后出现很多个虚频 [打印本页]

作者
Author: qian04551 时间: 2021-12-20 19:19
标题: 求助-PBE38优化-频率分析后出现很多个虚频
输入文件如下：
%mem=5GB
%chk=NDI-ID-RR-PBE38-opt-freq1.chk
%nprocshared=6
#p opt=(calcfc,tight) freq pbepbe/6-31+g(d,p) IOp(3/76=1000003750) IOp(3/77=0625006250) scrf=(pcm,solvent=dichloromethane)
虚频如下添加在了下面。
感谢各位老师解答！

作者
Author: snljty 时间: 2021-12-20 19:57
Gaussian的IOp只对第一步起作用，对第二步无效，而Gaussian的Opt和Freq是两步。所以相当于你在PBE38下优化的结构用PBE算的频率，结果没有意义。自己优化完以后读上一步几何结构和波函数做频率分析即可。可以使用Link1一个文件完成。你注意检查输出文件第二步的提示，肯定有相关内容的。

作者
Author: 喵星大佬 时间: 2021-12-20 20:52
为啥要用pbe38优化基态结构。。。。这玩意不是几乎只有算激发能的时候比较有用么

作者
Author: 红米饭1234 时间: 2021-12-20 21:20
#p opt freq pbe1pbe1/6-31+g(d,p) scrf=(iefpcm,solvent=dichloromethane)

#p opt(calcall,gdiis) freq pbe1pbe1/6-31+g(d,p) scrf=(iefpcm,solvent=dichloromethane)

试试

作者
Author: sobereva 时间: 2021-12-20 21:26

喵星大佬发表于 2021-12-20 20:52
为啥要用pbe38优化基态结构。。。。这玩意不是几乎只有算激发能的时候比较有用么

也就适合（相对于B3LYP、PBE0而言）大共轭体系，但PBE38优化基态结构几乎没有任何文章进行过测试。所以这么用严重不推荐

作者
Author: sobereva 时间: 2021-12-20 21:26

红米饭1234 发表于 2021-12-20 21:20
#p opt freq pbe1pbe1/6-31+g(d,p) scrf=(iefpcm,solvent=dichloromethane)

#p opt(calcall,gdiis) f ...

这是用的PBE0，跟PBE38是两码事

作者
Author: kay 时间: 2021-12-21 08:32
优化和频率分析分开做。

作者
Author: qian04551 时间: 2021-12-21 08:55

snljty 发表于 2021-12-20 19:57
Gaussian的IOp只对第一步起作用，对第二步无效，而Gaussian的Opt和Freq是两步。所以相当于你在PBE38下优化 ...

谢谢您的回答，感谢

作者
Author: 红米饭1234 时间: 2021-12-21 09:10

sobereva 发表于 2021-12-20 21:26
这是用的PBE0，跟PBE38是两码事

我就是想说啊，让他用PBE0试试

作者
Author: sobereva 时间: 2021-12-21 12:08

红米饭1234 发表于 2021-12-21 09:10
我就是想说啊，让他用PBE0试试

LZ问的是当前泛函下的计算问题，而换了泛函结果就变了，在不换泛函的前提下LZ的问题是可以直接解决的，即2L说的做法

作者
Author: sunhe14 时间: 2024-7-2 15:51

kay 发表于 2021-12-21 08:32
优化和频率分析分开做。

您好，再问一个，我要用pbe38计算激发态频率的话，用的是激发态优化完的几何坐标，再算频率是写 td freq 还是直接写freq,写td freq 的话会不会也要分两步，怎么分，求问，谢谢

作者
Author: snljty2 时间: 2024-7-2 16:07
本帖最后由 snljty2 于 2024-7-2 16:08 编辑

sunhe14 发表于 2024-7-2 15:51
您好，再问一个，我要用pbe38计算激发态频率的话，用的是激发态优化完的几何坐标，再算频率是写 td freq ...

%oldchk=上次opt的chk.chk
%chk=新的chk
#P Freq TD=(Root=上次优化得到态的序号，NStates=上次计算态的个数) 泛函/ChkBasis 其他关键词比如IOp Guess=Read Geom=AllCheck
[空行]
后面什么都不写

作者
Author: sunhe14 时间: 2024-7-2 16:36

snljty2 发表于 2024-7-2 16:07
%oldchk=上次opt的chk.chk
%chk=新的chk
#P Freq TD=(Root=上次优化得到态的序号，NStates=上次计算态的 ...

%chk=1TPA-jifa-freq-h2o.chk
%nprocshared=32
%mem=10GB
#p freq td pbepbe/6-31+g(d) scrf=iefpcm iop(3/76=1000003750,3/77=0625006250,3/78=1000010000)

Title Card Required

0 1
C                -1.25817300 2.99484600 -0.17404500
C                0.08910800 2.67247600 -0.15219200
C                1.09328100 3.66779400 -0.20423700
C                0.69025100 5.04289000 -0.27707700
C                -0.66169400 5.35777400 -0.29547700
C                -1.61982200 4.35820600 -0.24251300
H                0.40273000 1.64063100 -0.06661700
H                -0.94754100 6.39981400 -0.36088600
H                -2.66525800 4.63589900 -0.28698300
C                -4.28095000 -0.06272600 -0.07264900
C                -3.04196400 -0.29739100 -0.67861700
C                -2.07658900 0.69186900 -0.71483900
C                -2.29454000 1.95186100 -0.14413400
C                -3.53750600 2.17425300 0.46104400
C                -4.51280900 1.19469500 0.49549500
H                -2.83817900 -1.26088700 -1.12793200
H                -1.13626100 0.48080600 -1.20930700
H                -3.74772400 3.12891800 0.92783400
H                -5.46351900 1.40229100 0.96954800
N                -5.26361700 -1.05883100 -0.03718300
C                -6.43595200 -1.21447000 3.47401500
C                -5.59562600 -1.09438600 2.37629800
C                -6.11317100 -1.18166000 1.08279800
C                -7.48039500 -1.40136000 0.90842100
C                -8.31201300 -1.53672900 2.01127100
C                -7.79723700 -1.43987700 3.29959900
H                -6.01964500 -1.14526900 4.47178700
H                -4.53395000 -0.93433900 2.51783100
H                -7.88745000 -1.46602600 -0.09305100
H                -9.37122600 -1.70681300 1.86032100
H                -8.44937800 -1.54010200 4.15821700
C                -5.77675400 -4.18550500 -1.95466100
C                -5.72656500 -3.71849700 -3.26366600
C                -5.52128900 -2.36261100 -3.49566800
C                -5.37504800 -1.48063800 -2.43441300
C                -5.41517000 -1.95020600 -1.12037200
C                -5.61491400 -3.31199100 -0.88840100
H                -5.93054300 -5.23989800 -1.75906100
H                -5.84713100 -4.40329100 -4.09379100
H                -5.48844300 -1.98335800 -4.50993900
H                -5.23041900 -0.42360300 -2.62026000
H                -5.64248600 -3.68201100 0.12899400
O                1.57391100 6.03601000 -0.32906900
H                2.48796600 5.66676900 -0.30052300
C                2.45517400 3.34694500 -0.16858300
C                4.58591800 3.46899600 -0.14032900
C                4.24351600 2.10299500 -0.06873600
O                2.87276900 2.05428100 -0.09033800
H                5.57288800 3.90245000 -0.14994700
C                4.98525700 0.92614300 0.01648300
C                6.41136700 0.95393700 0.05326600
C                4.35270500 -0.35027100 0.07415900
C                7.13395200 -0.20389800 0.13888300
H                6.93242200 1.90277400 0.02697200
C                5.09493400 -1.49739300 0.16787700
H                3.27291300 -0.40867900 0.03160300
C                6.51106500 -1.47709900 0.20316500
H                8.21234100 -0.13109200 0.19219000
H                4.57241200 -2.44514100 0.18245900
C                7.28885400 -2.70106400 0.30061100
C                8.63657500 -2.75045800 -0.11088300
C                6.72816200 -3.88924500 0.81251000
C                9.37550500 -3.91718900 -0.01617200
H                9.10305600 -1.87307500 -0.53983300
C                7.47161400 -5.05300400 0.90726200
H                5.70690600 -3.89435100 1.17086900
C                8.80166000 -5.07929600 0.49410900
H                10.40573200 -3.92342700 -0.35169900
H                7.01322200 -5.94562600 1.31603600
H                9.38160100 -5.99085700 0.56864200
N                3.48607900 4.20495800 -0.19936100

我是这样写的，坐标是激发态优化完的最后一步坐标，nstates=3，root=1，可以不写吧，也就是说这个pbe38,加的iop是可以将td和freq 一起算的是吧，不用像opt和freq需要分开算对吧，谢谢

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