计算化学公社

 找回密码 Forget password
 注册 Register
Views: 53|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] Gaussian和ORCA计算振动分辨吸收光谱结果只有一个单峰,最大吸收波长与实验值相差较大

[复制链接 Copy URL]

7

帖子

0

威望

143

eV
积分
150

Level 3 能力者

各位老师好,我最近使用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

...............
gaussian.zip (2.08 MB, 下载次数 Times of downloads: 0) orca.zip (1.03 MB, 下载次数 Times of downloads: 1)

1万

帖子

0

威望

9900

eV
积分
22154

Level 6 (一方通行)

2#
发表于 Post on 4 hour ago | 只看该作者 Only view this author
检查是不是hess文件用错了,或原子坐标错了。你计算的K*K精确等于0,也就是基态平衡结构和激发态平衡结构严格一致,应该是其中一个结构错了。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?user=XW6C6eQAAAAJ
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1133/1776.htm
GitHub:https://github.com/wzkchem5
本团队长期招收研究生,有意者可私信联系

7

帖子

0

威望

143

eV
积分
150

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 3 hour ago | 只看该作者 Only view this author
本帖最后由 llian 于 2026-2-27 15:04 编辑
wzkchem5 发表于 2026-2-27 12:12
检查是不是hess文件用错了,或原子坐标错了。你计算的K*K精确等于0,也就是基态平衡结构和激发态平衡结构严 ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-27 16:48 , Processed in 0.677906 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list