计算化学公社

标题: 关于freq=anharmonic计算的频率 [打印本页]

作者
Author:
幻七熏    时间: 2016-6-14 22:04
标题: 关于freq=anharmonic计算的频率
大家好,我计算的体系含Co,C,O;体系不算大,实验的光谱感觉可能会有合频,所以想用freq=anharmonic计算一下。使用的方法是B3LYP,基组对Co是lanl2dz,C和O是DZP基组。计算显示小于1000cm-1的峰的非谐性频率异乎寻常的大,这样子它们的合频也就不准确了,请问是什么原因呢?非常感谢赐教。

作者
Author:
幻七熏    时间: 2016-6-14 22:06
截图没截好,第一列频率是谐振的,第二列是非谐性的
作者
Author:
sobereva    时间: 2016-6-14 22:30
没准是数值不稳定之类。优化时候用tight标准、int=ultrafine结合calcall再看看非谐振结果有无改变。也可以试试HF的结果。
作者
Author:
幻七熏    时间: 2016-6-15 09:15
sobereva 发表于 2016-6-14 22:30
没准是数值不稳定之类。优化时候用tight标准、int=ultrafine结合calcall再看看非谐振结果有无改变。也可以 ...

好的,非常感谢
作者
Author:
幻七熏    时间: 2016-7-12 20:19
前段没时间,现在尝试了,正如sob老师所说,的确是数值问题,用了opt=tignt和int=ultrafine成功解决,再次表示感谢!
作者
Author:
huangtf    时间: 2017-1-17 19:03
请问sob老师、liyuanhe老师和大家,输入文件中freq=anharmonic,在log文件中怎么看非简谐频率的强度?谢谢!
作者
Author:
sobereva    时间: 2017-1-17 19:20
huangtf 发表于 2017-1-17 19:03
请问sob老师、liyuanhe老师和大家,输入文件中freq=anharmonic,在log文件中怎么看非简谐频率的强度?谢谢 ...

找这种地方。I是强度
(, 下载次数 Times of downloads: 84)

作者
Author:
huangtf    时间: 2017-1-17 20:24
本帖最后由 huangtf 于 2017-1-17 20:25 编辑
sobereva 发表于 2017-1-17 19:20
找这种地方。I是强度


sob老师,您看我的log里面是这样的,没有I,怎么办呢?
作者
Author:
sobereva    时间: 2017-1-17 20:36
huangtf 发表于 2017-1-17 20:24
sob老师,您看我的log里面是这样的,没有I,怎么办呢?

用最新版本
作者
Author:
huangtf    时间: 2017-1-17 20:45
sobereva 发表于 2017-1-17 20:36
用最新版本

sob老师,您说的“用最新版本”指的是用高斯的最新版本吗?我用的g09。
作者
Author:
sobereva    时间: 2017-1-17 20:54
huangtf 发表于 2017-1-17 20:45
sob老师,您说的“用最新版本”指的是用高斯的最新版本吗?我用的g09。

D.01及以后
作者
Author:
huangtf    时间: 2017-1-17 20:58
sobereva 发表于 2017-1-17 20:54
D.01及以后

好的,sob老师,谢谢您!学生明白了。
作者
Author:
xiaobogaga    时间: 2019-6-22 11:57
sobereva 发表于 2016-6-14 22:30
没准是数值不稳定之类。优化时候用tight标准、int=ultrafine结合calcall再看看非谐振结果有无改变。也可以 ...

sobereva老师,您好!
我最近也在计算分子的非简谐频率,我查了高斯手册,高斯默认的是GVPT2方法,如果我想用VPT2方法,应该怎么输入呢?
我输入的是 # B3LYP/6-311+G(3df) opt=calcall Freq=(Anharm,ReadAnharm) 然后在坐标后空格一行,输入PT2Model=VPT2 计算结束后发现,跟GVPT2方法相比,anharmonic频率没有区别。
作者
Author:
sobereva    时间: 2019-6-22 15:34
xiaobogaga 发表于 2019-6-22 11:57
sobereva老师,您好!
我最近也在计算分子的非简谐频率,我查了高斯手册,高斯默认的是GVPT2方法,如果 ...

就是你说的这么输入。没差别可能是对当前体系本身就没差别或者差异小到可忽略不计

作者
Author:
xiaobogaga    时间: 2019-6-22 16:36
sobereva 发表于 2019-6-22 15:34
就是你说的这么输入。没差别可能是对当前体系本身就没差别或者差异小到可忽略不计

谢谢sob老师!
作者
Author:
413    时间: 2020-3-6 15:00
借图搭车请教一下。
第一列是谐性的频率、第二列的非谐性的频率。
请教一下,每一行中的谐性和非谐性是一一对应的关系吗?非谐性从上到下是按照什么规则排列的?另,我们咋看每个非谐性频率对应的哪个原子的啥振动模式?

作者
Author:
sobereva    时间: 2020-3-7 11:27
413 发表于 2020-3-6 15:00
借图搭车请教一下。
第一列是谐性的频率、第二列的非谐性的频率。
请教一下,每一行中的谐性和非谐性是一 ...

是。
gview里就可以看,gv6的频率观看窗口会显示各个模式的谐振和非谐振频率




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3