计算化学公社

标题: 激发态势能面扫描请教 [打印本页]

作者
Author: shiyz 时间: 2019-1-7 20:37
标题: 激发态势能面扫描请教
本帖最后由 shiyz 于 2019-1-7 20:40 编辑

请教诸位大神，要获得如图1所示的势能面扫描图，其中3号线的绝热第一激发态在不同角度下的值是怎么计算的？
我通过固定一个二面角，再计算此二面角下分子的绝热S1态，
%chk=10.chk
%mem=8000MB
%nprocshared=14
#p opt b3lyp/3-21g td(nstates=15,singlets,root=1,direct) IOp(9/40=4)

*******
*******
D 1 2 3 4 10.0 F

但是出现如下的结果，且没有优化完全，改怎么破？
It is also a good rule not to put too much confidence in
experimental results until they have been confirmed by
Theory.
                           -- Sir Arthur Eddington
Error termination request processed by link 9999.
Error termination via Lnk1e in /opt/software/g09D01/g09/l9999.exe at Mon Dec 24 12:09:54 2018.
Job cpu time:    18 days  5 hours 27 minutes 15.8 seconds.
File lengths (MBytes):  RWF= 956 Int=    0 D2E=    0 Chk= 101 Scr=    2

作者
Author: 哇哇吐 时间: 2019-1-7 22:00
我的看法是先不用固定角度，就正常把S1优化出来，然后在S1结构的基础上做柔性扫描(opt=modredundant)即得到3号线，至于你的报错论坛搜一搜就有吧。

作者
Author: sobereva 时间: 2019-1-7 23:27
别在优化时候用IOp(9/40=4)，不仅造成输出文件大，Multiwfn还没法正常利用
甭写没意义的direct
优化root=1设nstates=15完全是浪费时间
3-21G下优化激发态不要指望能得到靠谱结果
冻结时候不能设初值（10.0），仔细看
在Gaussian中做限制性优化的方法
http://sobereva.com/404

几何优化不收敛看
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

作者
Author: shiyz 时间: 2019-1-8 15:43

sobereva 发表于 2019-1-7 23:27
别在优化时候用IOp(9/40=4)，不仅造成输出文件大，Multiwfn还没法正常利用
甭写没意义的direct
优化root= ...

谢谢sobereva老师，我修改一下，再尝试一下

作者
Author: shiyz 时间: 2019-1-28 10:23

哇哇吐发表于 2019-1-7 22:00
我的看法是先不用固定角度，就正常把S1优化出来，然后在S1结构的基础上做柔性扫描(opt=modredundant)即得到 ...

在S1的结构上做scan，得到的是基态的数据吧

作者
Author: 小西瓜123 时间: 2021-7-6 17:15

shiyz 发表于 2019-1-28 10:23
在S1的结构上做scan，得到的是基态的数据吧

请问老师，如果要做激发态的扫描，可以是在把基态扫描的每一个结构做td opt吗？还有什么办法可以做激发态扫描吗？

作者
Author: wzkchem5 时间: 2021-7-6 17:50

小西瓜123 发表于 2021-7-6 10:15
请问老师，如果要做激发态的扫描，可以是在把基态扫描的每一个结构做td opt吗？还有什么办法可以做激发态 ...

一般但凡支持基态扫描和激发态结构优化的程序，都支持激发态扫描。就把算激发态的关键词和做扫描的关键词都写上就行了。
从基态扫描的结构做td opt不等于激发态扫描，除非你把你想扫描的那个内坐标限制住。

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