Gaussian使用FCHT输出Huang-Rhys因子求助

adreamy

我见帖子基于Gaussian FCHT计算分解重组能和计算Huang-Rhys因子的Python小脚本 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)提供了使用Gaussian计算Huang-Rhys因子的方法，然后我计算了一个纯有机分子（仅有CHN元素，基态S0和S1都经过了频率计算，使用Gaussian 16 A.03）后，使用以下输入文件计算S1结构的HR因子
%oldchk=S1.chk#p geom=allcheck freq(readfc,fc,readfcht)

initial=source=chk final=source=chk spectroscopy=onephotonemission
print=(huangrhys,matrix=JK)

S0.chk
S1.chk

其中当freq(readfc,fcht,readfcht)写成这个关键词时会报错无法输出HR因子，所以改成了 freq(readfc,fc,readfcht)，然后我得到HR因子之后，我用使用Dushin分解重组能和计算Huang-Rhys因子 - 思想家公社的门口：量子化学·分子模拟·二次元 (sobereva.com)帖子中提供的计算公式S_i=λ_i/(h*ν_i)，求出各个频率对应的重组能，发现这样计算出来的重组能(~5400 cm-1)远大于两点法计算的重组能(~2100cm-1),请问各位老师，我是哪里的关键词写错了，或者说是什么问题导致的这两种方法计算得到的重组能差异如此巨大。

cokie

首先，假如你表达的：S1结构的HR因子对应的是S1态的振动弛豫过程的话，你的spectroscopy应该写=onephotonabsorption
另外，如果你的结构里有可以旋转的平面并且旋转幅度比较大的话，有可能出现谐振模型无法描述的情况，此时HR因子就会虚高很多很多，这个时候以谐振模型做重组能拆分的程序基本都会失效（包括高斯、FCCLASSES和MOMAP，其他的没试过）

adreamy

cokie 发表于 2024-11-19 17:36
首先，假如你表达的：S1结构的HR因子对应的是S1态的振动弛豫过程的话，你的spectroscopy应该写=onephotonab ...

感谢您的解答，谢谢您