求助-PBE38优化-频率分析后出现很多个虚频

qian04551

输入文件如下：
%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)
虚频如下添加在了下面。
感谢各位老师解答！

snljty

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

喵星大佬

为啥要用pbe38优化基态结构。。。。这玩意不是几乎只有算激发能的时候比较有用么

红米饭1234

#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)  

试试

sobereva

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

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

sobereva

红米饭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是两码事

kay

优化和频率分析分开做。

qian04551

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

谢谢您的回答，感谢

红米饭1234

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

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

sobereva

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

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

sunhe14

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

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

snljty2

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
[空行]
后面什么都不写

sunhe14

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需要分开算对吧，谢谢