gaussianview里显示Failure reading anharmonic vibrational data

xyz0101

用高斯计算体系的非简谐振动频率，加入关键词freq=naharm，计算正常结束，将log文件拖到gview显示Failure reading anharmonic vibrational data。
gview看到输出的频率信息和直接加freq一致，请问该如何查看体系的非简谐振动频率？
另一种方法，将关键词写出freq=Anharmonic，用g16运算，计算非常耗时，是为什么呢？

sunwukong

定位到错误位置，enter多空一行可能就正常了。gauss bug吧

sobereva

是freq=anharm不是freq=naharm。freq=Anharmonic和这个完全等价

确保GaussView和Gaussian版本的兼容性，甭用太老的GaussView

GVPT2计算过程中需要算能量对坐标的三阶导数，耗时当然比只需要二阶导数的谐振频率高得多得多

北京科音基础（中级）量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）讲非谐振计算的时候专门说了：

xyz0101

sobereva 发表于 2024-5-23 03:03
是freq=anharm不是freq=naharm。freq=Anharmonic和这个完全等价

确保GaussView和Gaussian版 ...

谢谢老师的解答！
我想要计算体系的非谐振拉曼强度该如何添加关键词呢？我尝试加了freq=(anharm, raman)，这样的计算结果是否合理呢。

sobereva

xyz0101 发表于 2024-5-23 11:09
谢谢老师的解答！
我想要计算体系的非谐振拉曼强度该如何添加关键词呢？我尝试加了freq=(anharm, raman) ...

正确

xyz0101

sobereva 发表于 2024-5-24 04:40
正确

老师，我用到的是GaussView6.0.16，打开加freq=(anharm, raman)关键词，计算完的log文件，还是出现Failure reading anharmonic vibrational data
读取非谐波振动数据失败
在log文件里看这里可以吗？

xyz0101

如下

sobereva

xyz0101 发表于 2024-5-24 16:28
如下

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

当前数据明显异常，出现一堆星号，gview自然无法读取

E(anharm)的数值也明显异常

xyz0101

sobereva 发表于 2024-5-25 04:00
有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员 ...

老师，想知道异常原因是什么呢？以下是我的输入文件，用g09运算正常结束的：
%chk=test.chk
%nprocshared=20
%mem=30000MB
#p freq=(anharm, raman) genecp scrf=(smd,solvent=water) nosymm EmpiricalDispersion=GD3 geom=allcheck pbe1pbe pseudo=read scf=xqc

U 0
SDD
****
O H 0
6-311g(d,p)
****

U
SDD

log文件里有多处非简谐振动频率的强度是星号。

sobereva

xyz0101 发表于 2024-5-25 12:06
老师，想知道异常原因是什么呢？以下是我的输入文件，用g09运算正常结束的：
%chk=test.chk
%nprocshar ...

写了genecp就不要再写pseudo=read
做波函数稳定性测试确保波函数稳定
并且尝试用G16、去掉溶剂模型以简化情况