计算化学公社

标题: CCSD(T)做频率计算红外强度均为零是出于什么问题？ [打印本页]

作者
Author: YYY 时间: 2021-1-6 21:48
标题: CCSD(T)做频率计算红外强度均为零是出于什么问题？
各位老师，我在用CCSD(T)/aug-cc-pvtz做CO2的频率计算时，结果显示红外强度均为0，请问这是出于什么问题？下面是我的输入和输出文。

作者
Author: sobereva 时间: 2021-1-7 00:32
CCSD(T)这种没有解析一阶导数的方法做freq都不能给出红外强度
本质在于红外强度近似为偶极矩对正则坐标的导数，而有一阶解析导数的方法才直接有偶极矩（有限差分求偶极矩是另一码事）

作者
Author: YYY 时间: 2021-1-7 22:13

sobereva 发表于 2021-1-7 00:32
CCSD(T)这种没有解析一阶导数的方法做freq都不能给出红外强度
本质在于红外强度近似为偶极矩对正则坐标的 ...

明白了，谢谢sob老师！

作者
Author: cxh19990601 时间: 2023-1-17 13:34

sobereva 发表于 2021-1-7 00:32
CCSD(T)这种没有解析一阶导数的方法做freq都不能给出红外强度
本质在于红外强度近似为偶极矩对正则坐标的 ...

请问老师CASSCF是不是也无法解析一阶导数，我用它算频率也是算不出红外强度

作者
Author: zjxitcc 时间: 2023-1-17 13:54

cxh19990601 发表于 2023-1-17 13:34
请问老师CASSCF是不是也无法解析一阶导数，我用它算频率也是算不出红外强度

否。CASSCF解析一阶导数和二阶导数没有CCSD(T)的那么难，有不少程序支持。

作者
Author: sobereva 时间: 2023-1-19 08:19

cxh19990601 发表于 2023-1-17 13:34
请问老师CASSCF是不是也无法解析一阶导数，我用它算频率也是算不出红外强度

Gaussian就有CASSCF的二阶解析导数
直接用freq没给出红外强度是程序设计问题。CASSCF用freq=numer基于解析一阶导数做有限差分的方式做振动分析就能给你红外强度

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