计算化学公社

标题: gaussianview里显示Failure reading anharmonic vibrational data [打印本页]

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

作者
Author: sunwukong 时间: 2024-5-22 23:39
定位到错误位置，enter多空一行可能就正常了。gauss bug吧

作者
Author: sobereva 时间: 2024-5-23 03:03
是freq=anharm不是freq=naharm。freq=Anharmonic和这个完全等价

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

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

北京科音基础（中级）量子化学培训班（http://www.keinsci.com/workshop/KBQC_content.html）讲非谐振计算的时候专门说了：
(, 下载次数 Times of downloads: 30)

作者
Author: xyz0101 时间: 2024-5-23 11:09

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

确保GaussView和Gaussian版 ...

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

作者
Author: sobereva 时间: 2024-5-24 04:40

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

正确

作者
Author: xyz0101 时间: 2024-5-24 16:28

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

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

作者
Author: xyz0101 时间: 2024-5-24 16:28
如下

作者
Author: sobereva 时间: 2024-5-25 04:00

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

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

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

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

作者
Author: xyz0101 时间: 2024-5-25 12:06

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文件里有多处非简谐振动频率的强度是星号。

作者
Author: sobereva 时间: 2024-5-26 04:54

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

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

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