计算化学公社
标题:
Gaussian和ORCA计算振动分辨吸收光谱结果只有一个单峰,最大吸收波长与实验值相差较大
[打印本页]
作者Author:
llian
时间:
5 hour ago
标题:
Gaussian和ORCA计算振动分辨吸收光谱结果只有一个单峰,最大吸收波长与实验值相差较大
各位老师好,我最近使用Gaussian计算了振动分辨吸收光谱,分三步计算:
1.基态几何优化和频率计算:#P BMK/6-31G* opt freq scrf(solvent=dmso) nosymm em=gd3
2.激发态几何优化和频率计算:#P BMK/6-31G* opt freq=savenm scrf(solvent=dmso) td(nstates=3,root=1) nosymm em=gd3
3.FC分析:#p BMK/6-31G* freq=(readfc,FC,readFCHT) scrf(solvent=dmso) guess=read geom=check nosymm
最后一步出现报错:
-- To: vibronic fundamental state --
Spectrum progression: 0.00%
-- To: single overtones --
Spectrum progression: 0.00%
-- To: combinations of 2 simultaneously excited modes --
Spectrum progression: 0.00%
ERROR: Low progression after class 2. Total convergence = 0.0%.
The vibronic spectrum will likely be unreliable. Stopping.
然后参考了这篇帖子(
Gaussian16计算振动分辨的电子光谱出错?ERROR: Low progression after class 2. - 量子化学 (Quantum Chemistry) - 计算化学公社
)的回答,用orca计算了最后一步,参考这篇帖子(
ORCA5.0.3利用G16解析Hessian进行ESD计算 - 计算化学公社 - Powered by Discuz!
)将基态和激发态的fchk文件转为hess文件进行计算,得到的结果只有一个单峰,最大吸收波长为650 nm(实验值为460 nm),结果已放附件,想请教一下各位老师这是什么原因导致的?我应该如何调整?
下面是我的计算orca输入文件的参数设置,恳请各位老师指正。
! B3LYP DEF2-SVP RIJCOSX DEF2/J TIGHTSCF ESD(ABS) CPCM(DMSO)
%pal nprocs 32 end
%maxcore 8000
%TDDFT
NROOTS 3
IROOT 1
end
%ESD
GSHESSIAN "/home/liulian/MM-PMR/es/orca2/change_abs.hess"
ESHESSIAN "/home/liulian/MM-PMR/es/orca2/es_change.hess"
STATES 1
HESSFLAG AH
LINEW 50
UNIT NM
end
* xyz 0 1
C -11.41984488 -0.11981282 3.07010141
C -11.82879170 -0.38611002 1.75706741
C -10.89029571 -0.70269983 0.77109360
C -9.52370249 -0.77364709 1.09513675
...............
(, 下载次数 Times of downloads: 0)
上传 Uploaded
点击下载Click to download
(, 下载次数 Times of downloads: 1)
上传 Uploaded
点击下载Click to download
作者Author:
wzkchem5
时间:
5 hour ago
检查是不是hess文件用错了,或原子坐标错了。你计算的K*K精确等于0,也就是基态平衡结构和激发态平衡结构严格一致,应该是其中一个结构错了。
作者Author:
llian
时间:
4 hour ago
本帖最后由 llian 于 2026-2-27 15:04 编辑
wzkchem5 发表于 2026-2-27 12:12
检查是不是hess文件用错了,或原子坐标错了。你计算的K*K精确等于0,也就是基态平衡结构和激发态平衡结构严 ...
老师您好,十分感谢您的回复,我确认了hess文件并且重新计算了一遍,发现K*K值仍是0,然后我对比了高斯计算的基态和激发态的结构,发现相差不大,请问老师问题是否出在高斯计算上面?
欢迎光临 计算化学公社 (http://bbs.keinsci.com/)
Powered by Discuz! X3.3