计算化学公社

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

[Gaussian/gview] Gaussian TDDFT激发态绝热/透热态能量追踪脚本

[复制链接 Copy URL]

393

帖子

1

威望

5143

eV
积分
5556

Level 6 (一方通行)

本帖最后由 exity 于 2022-1-5 15:43 编辑

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






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

states_tracker.sh (999 Bytes, 下载次数 Times of downloads: 50)


评分 Rate

参与人数
Participants 1
eV +8 收起 理由
Reason
sobereva + 8

查看全部评分 View all ratings

1102

帖子

18

威望

6703

eV
积分
8165

Level 6 (一方通行)

計算化学の社畜

2#
发表于 Post on 2022-1-5 09:32:09 | 只看该作者 Only view this author
本帖最后由 冰释之川 于 2022-1-5 09:37 编辑

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

393

帖子

1

威望

5143

eV
积分
5556

Level 6 (一方通行)

3#
 楼主 Author| 发表于 Post on 2022-1-5 13:20:34 | 只看该作者 Only view this author
冰释之川 发表于 2022-1-5 09:32
能否支持对任务还处于运行中时相应的输出文件进行抓取?方便随时跟踪。
既然用shell脚本了,可以考虑用gnu ...

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

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

4#
发表于 Post on 2022-1-5 13:54:01 | 只看该作者 Only view this author
exity 发表于 2022-1-5 13:20
S1态优化不出来,折腾到这里已经无心继续了 T_T

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

393

帖子

1

威望

5143

eV
积分
5556

Level 6 (一方通行)

5#
 楼主 Author| 发表于 Post on 2022-1-5 15:31:53 | 只看该作者 Only view this author
喵星大佬 发表于 2022-1-5 13:54
也有可能你这个分子就没有S1态

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

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

6#
发表于 Post on 2022-1-5 16:01:40 | 只看该作者 Only view this author
exity 发表于 2022-1-5 15:31
一般有机分子都应该有S1吧

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

393

帖子

1

威望

5143

eV
积分
5556

Level 6 (一方通行)

7#
 楼主 Author| 发表于 Post on 2022-1-5 16:03:43 | 只看该作者 Only view this author
喵星大佬 发表于 2022-1-5 16:01
不一定,你看激发能到最后大概多少,如果很小的话可能所谓S1极小点就是CI,也就是没有S1极小点,那么如果 ...

至今没有有效收敛

1665

帖子

5

威望

4788

eV
积分
6553

Level 6 (一方通行)

喵星人

8#
发表于 Post on 2022-1-5 16:06:28 | 只看该作者 Only view this author
exity 发表于 2022-1-5 16:03
至今没有有效收敛

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

393

帖子

1

威望

5143

eV
积分
5556

Level 6 (一方通行)

9#
 楼主 Author| 发表于 Post on 2022-1-5 17:01:55 | 只看该作者 Only view this author
本帖最后由 exity 于 2022-1-5 17:22 编辑
喵星大佬 发表于 2022-1-5 16:06
我看你最后几步激发能越来越低,可能就是在朝着CI方向去,到CI点附近普通DFT就会出现各种问题

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

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

10#
发表于 Post on 2022-1-5 17:10:19 | 只看该作者 Only view this author
本帖最后由 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

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

11#
发表于 Post on 2022-1-5 17:40:34 | 只看该作者 Only view this author
exity 发表于 2022-1-5 08:31
一般有机分子都应该有S1吧

所有单重态分子都有S1态势能面,但S1态势能面未必有平衡结构。所以研究垂直激发一定有S1态,但研究绝热激发或荧光就是另一码事了。
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
本团队长期招收研究生,有意者可私信联系

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

12#
发表于 Post on 2022-1-5 17:41:08 | 只看该作者 Only view this author
exity 发表于 2022-1-5 10:01
大佬有什么好的解决方案推荐不?CASSCF行不?T_T

CASSCF是一种解决方案(但是如果需要可发表的结果,必须加动态相关),此外SF-TDDFT也可以考虑
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
本团队长期招收研究生,有意者可私信联系

221

帖子

0

威望

6214

eV
积分
6435

Level 6 (一方通行)

跳跳猪

13#
发表于 Post on 2022-1-5 18:02:14 | 只看该作者 Only view this author
没有diabatization哪来的diabatic states(透热态)
流年似水,浮生如梦。

393

帖子

1

威望

5143

eV
积分
5556

Level 6 (一方通行)

14#
 楼主 Author| 发表于 Post on 2022-1-5 18:41:39 | 只看该作者 Only view this author
wzkchem5 发表于 2022-1-5 17:41
CASSCF是一种解决方案(但是如果需要可发表的结果,必须加动态相关),此外SF-TDDFT也可以考虑

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

1187

帖子

5

威望

2876

eV
积分
4163

Level 6 (一方通行)

15#
发表于 Post on 2022-1-5 19:14:53 | 只看该作者 Only view this author
exity 发表于 2022-1-5 18:41
我目前可用的软件包有 G09 G16 和 ORCA,他们可以做吗?我搜索了一下看到一个GAMESS的帖子,是不是只有GA ...

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

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 04:54 , Processed in 0.277424 second(s), 24 queries , Gzip On.

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