g09计算OH自由基振动分辨的电子光谱出错求助

唐妹妹 · 发表于 Post on 2020-3-21 19:26:47

各位老师好！请问我用g09D版本计算OH自由基振动分辨的吸收谱，尝试了m06-2x,wB97XD,cam-B3LYP,PBE0,BMK,B3LYP泛函，也尝试了6-31g*,6-311g**,def2TZVP基组，输出最后均显示如下： Data for initial state taken from current calculation. Data for final state taken from checkpoint file "OH-S1-td-m062x.chk"
Normal modes recovered from file.

Using excited electronic state number  1.

Initial state structure is set in Eckart orientation.
Final state structure is superposed to it.

DSYEV returned Info=          5 IAlg= 2 N=    6 NDim=    6 NE2=       26.
Diagonalization in DiagDN failed.
Error termination via Lnk1e in /home/lzhpc/Downloads/g09/l718.exe at Wed Mar 18 17:07:23 2020.
Job cpu time:    0 days  0 hours  0 minutes 33.4 seconds.
File lengths (MBytes):  RWF=    5 Int=    0 D2E=    0 Chk=    1 Scr=    1
S0和S1态输出见附件。请问大神，还有什么尝试解决方法吗？

sobereva · 发表于 Post on 2020-3-22 14:16:15

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素（仔细看http://bbs.keinsci.com/thread-9348-1-1.html），我已把你的不恰当标题“g09计算OH自由基振动分辨的电子光谱”改了，以后务必注意

这和泛函和基组不会有直接关系。
如果用你同样的步骤计算水分子都没问题，8成是程序bug或者有一些特殊的坑，只有开发者清楚，最好发邮件问官方客服。

唐妹妹 · 发表于 Post on 2020-3-22 19:35:17

sobereva 发表于 2020-3-22 14:16
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出“提问”、“求助”要素（ ...

谢谢Sob老师，以后一定注意。计算您anisole那个例子是没有问题的。计算水分子出现ERROR: The Franck-Condon factor corresponding to the overlap integral between both vibrational ground states is too small: |<0'|0">|^2 = 0.21354E-05。这时用readFCHT读取ForceFCCalc关键词强行进行计算，可是报错
Data for initial state taken from current calculation.
Data for final state taken from checkpoint file "ForceFCCalc"
Error termination in NtrErr:
ntran open failure returned to fopen.
这是该读取.chk文件的位置却读到了ForceFCCalc，可确实是按Sob老师说的顺序写的输入文件，上传一下输入文件请大神指教。
另外，在g09官方手册上没有搜到ForceFCCalc关键词，请问哪里找它的说明呢？
谢谢各位老师的指导！

sobereva · 发表于 Post on 2020-3-22 20:02:44

唐妹妹发表于 2020-3-22 19:35
谢谢Sob老师，以后一定注意。计算您anisole那个例子是没有问题的。计算水分子出现ERROR: The Franck-Cond ...

ForceFCCalc那个是Gaussian开发者Barone在他自己的文档Vibrationally-resolved electronic spectra in GAUSSIAN 09（Google一搜就有）里写的。

唐妹妹 · 发表于 Post on 2020-3-23 09:12:22

sobereva 发表于 2020-3-22 20:02
ForceFCCalc那个是Gaussian开发者Barone在他自己的文档Vibrationally-resolved electronic spectra in GA ...

谢谢Sob老师，我再试试。

唐妹妹 · 发表于 Post on 2020-3-26 11:53:51

唐妹妹发表于 2020-3-23 09:12
谢谢Sob老师，我再试试。

Sob老师，您好！这个错误在g16里不会再出现。如您所说，振动分辨的电子光谱的计算那篇博文，比较适合Gaussian09，而对于Gaussian16，做法不是最好的。但我不知道g16里更好的做法是什么？所以又在g16测试了一下anisole那个例子，按g09的做法依然可以得到博文里的结果。于是我就按g09的做法在g16里算OH，用PBE0,cam-B3LYP算发射谱，光谱横坐标能量很多负值，请问Sob老师，这个负值时怎么回事？如何解决呢？附件为PBE0泛函计算结果。

OH-S1-pbe0.out (129.67 KB, 下载次数 Times of downloads: 4)

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] g09计算OH自由基振动分辨的电子光谱出错求助