计算化学公社

标题: Gaussian和ORCA计算振动分辨吸收光谱结果只有一个单峰，最大吸收波长与实验值相差较大 [打印本页]

作者
Author: llian 时间: 2026-2-27 11:59
标题: 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) (, 下载次数 Times of downloads: 1)

作者
Author: wzkchem5 时间: 2026-2-27 12:12
检查是不是hess文件用错了，或原子坐标错了。你计算的K*K精确等于0，也就是基态平衡结构和激发态平衡结构严格一致，应该是其中一个结构错了。

作者
Author: llian 时间: 2026-2-27 13:45
本帖最后由 llian 于 2026-2-27 15:04 编辑

wzkchem5 发表于 2026-2-27 12:12
检查是不是hess文件用错了，或原子坐标错了。你计算的K*K精确等于0，也就是基态平衡结构和激发态平衡结构严 ...

老师您好，十分感谢您的回复，我确认了hess文件并且重新计算了一遍，发现K*K值仍是0，然后我对比了高斯计算的基态和激发态的结构，发现相差不大，请问老师问题是否出在高斯计算上面？

作者
Author: wzkchem5 时间: 2026-2-28 12:13

llian 发表于 2026-2-27 13:45
老师您好，十分感谢您的回复，我确认了hess文件并且重新计算了一遍，发现K*K值仍是0，然后我对比了高斯计 ...

结构相差不大会导致K*K不大，但不会导致K*K精确为0。
检查hess文件转换格式的步骤是否正确，用记事本打开hess文件确认其中分子坐标是否正确，确认ORCA输入文件里面给的那个坐标正确。如果还是找不到原因，从orca说明书里找个算例重复一下，然后一步一步地认真对照检查你现在的操作流程是否有问题

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