计算化学公社

标题: 求助Gaussian中使用RedDim=ClearLowFreq=N移除低频振动的输入文件正确写法 [打印本页]

作者
Author: X教授 时间: 2025-11-5 00:03
标题: 求助Gaussian中使用RedDim=ClearLowFreq=N移除低频振动的输入文件正确写法
本帖最后由 X教授于 2025-11-5 00:04 编辑

最近在使用Kalinite发布的基于Gaussian FCHT计算分解重组能和计算Huang-Rhys因子的Python小脚本，http://bbs.keinsci.com/thread-28708-1-1.html，遇到了一个消除低频振动的问题。我想计算S0到S1跃迁的Huang-Rhys因子，在使用Gaussian FCHT计算时会报错，报错信息是：
ERROR: At least one near-null frequency detected in final state. Mode: 1 with frequency: 0.8122, cm^-1
Use "ClearLowFreq" to remove automatically those modes.

我检查了结构信息，是在优化S1的结构时有一个低频振动（0.81 cm-1），想利用报错信息提供的方法消除这个低频，但遇到了问题，尝试了以下几种形式，均提示输入文件错误：
#p freq td b3lyp/6-31g(d) em=gd3bj ClearLowFreq=1
#p freq td b3lyp/6-31g(d) em=gd3bj RedDim=ClearLowFreq=1

随后尝试将RedDim=ClearLowFreq=1或者ClearLowFreq=1置于坐标后面，计算正常进行，然而0.81 cm-1的低频振动仍然存在，且在计算FCHT时仍然提示0.81 cm-1的低频振动需要去掉：
#p freq td b3lyp/6-31G(d)

Title

0 1
coordinates

RedDim=ClearLowFreq=1
查阅了高斯的官方文档，https://gaussian.com/freq/，想请教一下RedDim=ClearLowFreq=n应该如何使用？相应的输入输出文件已上传到附件中。

作者
Author: sobereva 时间: 2025-11-9 03:41
参考北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里“振动电子光谱的计算”一节的相关例子：
(, 下载次数 Times of downloads: 0)

(, 下载次数 Times of downloads: 0)

作者
Author: X教授 时间: 2025-11-12 16:53

sobereva 发表于 2025-11-9 03:41
参考北京科音高级量子化学培训班（http://www.keinsci.com/KAQC）里“振动电子光谱的计算”一节的相关例子 ...

感谢Sob老师，问题已解决，实际上检测发现，RedDim=ClearLowFreq=N，N必须要大于1，因此从2开始就没有问题，写法的话还是写在FCHT计算中加上：
%oldchk=im-s1.chk
%nprocshared=128
%mem=128GB
#p geom=allcheck freq(readfc,fc,readfcht)

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

im.chk
im-s1.chk

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