计算化学公社

 找回密码 Forget password
 注册 Register
Views: 2855|回复 Reply: 15

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

[复制链接 Copy URL]

312

帖子

1

威望

3106

eV
积分
3438

Level 5 (御坂)

发表于 Post on 2022-1-4 19:17:42 | 显示全部楼层 Show all |阅读模式 Reading model
本帖最后由 exity 于 2022-1-5 15:43 编辑

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

2.png

3.png

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

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

评分 Rate

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

查看全部评分 View all ratings

1061

帖子

16

威望

5780

eV
积分
7161

Level 6 (一方通行)

計算化学の社畜

发表于 Post on 2022-1-5 09:32:09 | 显示全部楼层 Show all
本帖最后由 冰释之川 于 2022-1-5 09:37 编辑

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

312

帖子

1

威望

3106

eV
积分
3438

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-1-5 13:20:34 | 显示全部楼层 Show all
冰释之川 发表于 2022-1-5 09:32
能否支持对任务还处于运行中时相应的输出文件进行抓取?方便随时跟踪。
既然用shell脚本了,可以考虑用gnu ...

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

1464

帖子

1

威望

2591

eV
积分
4075

Level 6 (一方通行)

喵星人

发表于 Post on 2022-1-5 13:54:01 | 显示全部楼层 Show all
exity 发表于 2022-1-5 13:20
S1态优化不出来,折腾到这里已经无心继续了 T_T

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

312

帖子

1

威望

3106

eV
积分
3438

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-1-5 15:31:53 | 显示全部楼层 Show all
喵星大佬 发表于 2022-1-5 13:54
也有可能你这个分子就没有S1态

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

1464

帖子

1

威望

2591

eV
积分
4075

Level 6 (一方通行)

喵星人

发表于 Post on 2022-1-5 16:01:40 | 显示全部楼层 Show all
exity 发表于 2022-1-5 15:31
一般有机分子都应该有S1吧

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

312

帖子

1

威望

3106

eV
积分
3438

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-1-5 16:03:43 | 显示全部楼层 Show all
喵星大佬 发表于 2022-1-5 16:01
不一定,你看激发能到最后大概多少,如果很小的话可能所谓S1极小点就是CI,也就是没有S1极小点,那么如果 ...

至今没有有效收敛

1464

帖子

1

威望

2591

eV
积分
4075

Level 6 (一方通行)

喵星人

发表于 Post on 2022-1-5 16:06:28 | 显示全部楼层 Show all
exity 发表于 2022-1-5 16:03
至今没有有效收敛

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

312

帖子

1

威望

3106

eV
积分
3438

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-1-5 17:01:55 | 显示全部楼层 Show all
本帖最后由 exity 于 2022-1-5 17:22 编辑
喵星大佬 发表于 2022-1-5 16:06
我看你最后几步激发能越来越低,可能就是在朝着CI方向去,到CI点附近普通DFT就会出现各种问题

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

1188

帖子

5

威望

2758

eV
积分
4046

Level 6 (一方通行)

发表于 Post on 2022-1-5 17:10:19 | 显示全部楼层 Show all
本帖最后由 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

6672

帖子

0

威望

4002

eV
积分
10674

Level 6 (一方通行)

发表于 Post on 2022-1-5 17:40:34 | 显示全部楼层 Show all
exity 发表于 2022-1-5 08:31
一般有机分子都应该有S1吧

所有单重态分子都有S1态势能面,但S1态势能面未必有平衡结构。所以研究垂直激发一定有S1态,但研究绝热激发或荧光就是另一码事了。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

6672

帖子

0

威望

4002

eV
积分
10674

Level 6 (一方通行)

发表于 Post on 2022-1-5 17:41:08 | 显示全部楼层 Show all
exity 发表于 2022-1-5 10:01
大佬有什么好的解决方案推荐不?CASSCF行不?T_T

CASSCF是一种解决方案(但是如果需要可发表的结果,必须加动态相关),此外SF-TDDFT也可以考虑
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

214

帖子

0

威望

4830

eV
积分
5044

Level 6 (一方通行)

跳跳猪

发表于 Post on 2022-1-5 18:02:14 | 显示全部楼层 Show all
没有diabatization哪来的diabatic states(透热态)
流年似水,浮生如梦。

312

帖子

1

威望

3106

eV
积分
3438

Level 5 (御坂)

 楼主 Author| 发表于 Post on 2022-1-5 18:41:39 | 显示全部楼层 Show all
wzkchem5 发表于 2022-1-5 17:41
CASSCF是一种解决方案(但是如果需要可发表的结果,必须加动态相关),此外SF-TDDFT也可以考虑

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

1188

帖子

5

威望

2758

eV
积分
4046

Level 6 (一方通行)

发表于 Post on 2022-1-5 19:14:53 | 显示全部楼层 Show all
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, 2023-2-2 23:41 , Processed in 0.656461 second(s), 25 queries .

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