计算化学公社

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

[辅助/分析程序] dushin算反应物和产物的重组能距离判断报错求助

[复制链接 Copy URL]

143

帖子

0

威望

1554

eV
积分
1697

Level 5 (御坂)

各位老师,请教关于dushin的问题,我用dushin算反应物和产物的重组能,反应物和产物是沿着irc路径取的没有虚频的结构,在用dushin算的时候出现某两个原子距离太远的问题,我的结构和报错分别如图所示,请问这个距离判断可以修改吗?如果可以,请问应该怎么修改呢?结构与报错如图所示

dushin1.jpg (46.42 KB, 下载次数 Times of downloads: 17)

产物结构

产物结构

dushin2.jpg (45.42 KB, 下载次数 Times of downloads: 15)

反应物结构

反应物结构

dushin3.jpg (90.13 KB, 下载次数 Times of downloads: 19)

dushin报错

dushin报错

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

2#
发表于 Post on 2021-6-22 17:33:07 | 只看该作者 Only view this author
如果结构变化太大,说明dushin的那种计算方法的结果不靠谱,因为势能面不能认为是谐振的。应该改用四点法
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
本团队长期招收研究生,有意者可私信联系

143

帖子

0

威望

1554

eV
积分
1697

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2021-6-22 19:42:10 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-22 17:33
如果结构变化太大,说明dushin的那种计算方法的结果不靠谱,因为势能面不能认为是谐振的。应该改用四点法

谢谢老师,可是我这种情况怎么用四点法呢?两个态之间只有结构的区别,没有电荷和自旋多重度的区别

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

4#
发表于 Post on 2021-6-22 20:39:28 | 只看该作者 Only view this author
小西瓜123 发表于 2021-6-22 12:42
谢谢老师,可是我这种情况怎么用四点法呢?两个态之间只有结构的区别,没有电荷和自旋多重度的区别

你这个不是反应物标的是T1,产物标的是S0吗?还是说实际电子态和这里画的不一样?
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
本团队长期招收研究生,有意者可私信联系

143

帖子

0

威望

1554

eV
积分
1697

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2021-6-22 21:02:39 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-22 20:39
你这个不是反应物标的是T1,产物标的是S0吗?还是说实际电子态和这里画的不一样?

嗯嗯实际电子态和这里画的不一样,产物和反应物的都是03,本来应该反应物是S1的结构的,但是做tddft优化S1的时候,优化出来的结构不符合化学直觉了,不是我所算的反应物了,所以我就只在基态下得到了反应物和产物的结构和频率,命名t1和s0是和dushin里面的对应

1万

帖子

0

威望

9855

eV
积分
22090

Level 6 (一方通行)

6#
发表于 Post on 2021-6-22 22:56:42 | 只看该作者 Only view this author
小西瓜123 发表于 2021-6-22 14:02
嗯嗯实际电子态和这里画的不一样,产物和反应物的都是03,本来应该反应物是S1的结构的,但是做tddft优化S ...

同一个电子态的不同极小值点,不好定义重组能。重组能是按四点法定义的,dushin的做法只是近似计算重组能,以及把重组能拆分成各个振动模式的贡献的一种方法,所以四点法做不了的情况,重组能是没有严格定义的。
一定要定义重组能的话,必须在透热态框架下做。想办法构造两个透热态,使得这两个结构分别处于不同的透热态上,然后才可以定义重组能。
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
本团队长期招收研究生,有意者可私信联系

143

帖子

0

威望

1554

eV
积分
1697

Level 5 (御坂)

7#
 楼主 Author| 发表于 Post on 2021-6-23 09:55:10 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-22 22:56
同一个电子态的不同极小值点,不好定义重组能。重组能是按四点法定义的,dushin的做法只是近似计算重组能 ...

那还是要尽量得到T1的结构哈,谢谢老师

本版积分规则 Credits rule

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

GMT+8, 2026-2-19 02:35 , Processed in 0.173847 second(s), 23 queries , Gzip On.

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