计算化学公社

标题: Gaussian TDDFT激发态绝热/透热态能量追踪脚本 [打印本页]

作者
Author: exity 时间: 2022-1-4 19:17
标题: Gaussian TDDFT激发态绝热/透热态能量追踪脚本
本帖最后由 exity 于 2022-1-5 15:43 编辑

在高斯里使用TDDFT进行激发态优化，假设使用nstates=5,root=1参数，正常结束任务时对out文件使用该脚本，即可得到S1-S5的能量-优化步数txt文件，直接全选复制黏贴到excel里即可作图，无需分列操作。
(, 下载次数 Times of downloads: 13)

(, 下载次数 Times of downloads: 18)

(, 下载次数 Times of downloads: 12)

用法：放在纯英文，无空格路径的目标文件下赋予权限使用，或加入PATH使用

(, 下载次数 Times of downloads: 50)

作者
Author: 冰释之川 时间: 2022-1-5 09:32
本帖最后由冰释之川于 2022-1-5 09:37 编辑

能否支持对任务还处于运行中时相应的输出文件进行抓取？方便随时跟踪。
既然用shell脚本了，可以考虑用gnuplot之类的画图并且直接导出图片，这样可以省得再打开excel作图了
现成的例子：http://blog.chinaunix.net/uid-24774106-id-3957889.html

作者
Author: exity 时间: 2022-1-5 13:20

冰释之川发表于 2022-1-5 09:32
能否支持对任务还处于运行中时相应的输出文件进行抓取？方便随时跟踪。
既然用shell脚本了，可以考虑用gnu ...

S1态优化不出来，折腾到这里已经无心继续了 T_T

作者
Author: 喵星大佬 时间: 2022-1-5 13:54

exity 发表于 2022-1-5 13:20
S1态优化不出来，折腾到这里已经无心继续了 T_T

也有可能你这个分子就没有S1态

作者
Author: exity 时间: 2022-1-5 15:31

喵星大佬发表于 2022-1-5 13:54
也有可能你这个分子就没有S1态

一般有机分子都应该有S1吧

作者
Author: 喵星大佬 时间: 2022-1-5 16:01

exity 发表于 2022-1-5 15:31
一般有机分子都应该有S1吧

不一定，你看激发能到最后大概多少，如果很小的话可能所谓S1极小点就是CI，也就是没有S1极小点，那么如果有荧光的话可能不是S1发射的

作者
Author: exity 时间: 2022-1-5 16:03

喵星大佬发表于 2022-1-5 16:01
不一定，你看激发能到最后大概多少，如果很小的话可能所谓S1极小点就是CI，也就是没有S1极小点，那么如果 ...

至今没有有效收敛

作者
Author: 喵星大佬 时间: 2022-1-5 16:06

exity 发表于 2022-1-5 16:03
至今没有有效收敛

我看你最后几步激发能越来越低，可能就是在朝着CI方向去，到CI点附近普通DFT就会出现各种问题

作者
Author: exity 时间: 2022-1-5 17:01
本帖最后由 exity 于 2022-1-5 17:22 编辑

喵星大佬发表于 2022-1-5 16:06
我看你最后几步激发能越来越低，可能就是在朝着CI方向去，到CI点附近普通DFT就会出现各种问题

大佬有什么好的解决方案推荐不？CASSCF行不？T_T

作者
Author: snljty 时间: 2022-1-5 17:10
本帖最后由 snljty 于 2022-1-5 17:11 编辑

exity 发表于 2022-1-5 17:01
大佬有什么好的解决方案推荐不？CASSCF行不？还有什么是CI，簇耦合吗？T_T

conical intersection，圆锥交叉点。可以用多参考方法描述，用spin-flip TD-DFT等方法至少也可以描述。可以参考https://mp.weixin.qq.com/s/UY3Ixz6CO86XAeVGtuXn4w

作者
Author: wzkchem5 时间: 2022-1-5 17:40

exity 发表于 2022-1-5 08:31
一般有机分子都应该有S1吧

所有单重态分子都有S1态势能面，但S1态势能面未必有平衡结构。所以研究垂直激发一定有S1态，但研究绝热激发或荧光就是另一码事了。

作者
Author: wzkchem5 时间: 2022-1-5 17:41

exity 发表于 2022-1-5 10:01
大佬有什么好的解决方案推荐不？CASSCF行不？T_T

CASSCF是一种解决方案（但是如果需要可发表的结果，必须加动态相关），此外SF-TDDFT也可以考虑

作者
Author: Mikasa 时间: 2022-1-5 18:02
没有diabatization哪来的diabatic states(透热态)

作者
Author: exity 时间: 2022-1-5 18:41

wzkchem5 发表于 2022-1-5 17:41
CASSCF是一种解决方案（但是如果需要可发表的结果，必须加动态相关），此外SF-TDDFT也可以考虑

我目前可用的软件包有 G09 G16 和 ORCA，他们可以做吗？我搜索了一下看到一个GAMESS的帖子，是不是只有GAMESS能做？

作者
Author: snljty 时间: 2022-1-5 19:14

exity 发表于 2022-1-5 18:41
我目前可用的软件包有 G09 G16 和 ORCA，他们可以做吗？我搜索了一下看到一个GAMESS的帖子，是不是只有GA ...

ORCA就可以，不管是NEVPT2还是spin-flip TD-DFT都可以。前面给你的回复sob老师正在审核。

作者
Author: exity 时间: 2022-1-5 20:09

snljty 发表于 2022-1-5 19:14
ORCA就可以，不管是NEVPT2还是spin-flip TD-DFT都可以。前面给你的回复sob老师正在审核。

谢谢，我正在尝试用SF-TDDFT来做这个工作

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