计算化学公社

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

[新手求助] 有机分子三重态激发态计算求助

[复制链接 Copy URL]

41

帖子

0

威望

741

eV
积分
782

Level 4 (黑子)

老师好,刚接触一个涉及激发态的课题,需要计算有机分子三重态时激发态的能量,第一次算激发态有很多不懂的地方。在计算上文献提到:“We studied the energetics of 3a(singlet) and 3a* (triplet) computationally using the B3LYP/6-311+G(2d,p) method” ,不太明白3a的singlet和3a*的triplet优化计算怎么处理的
3a的singlet优化是 # opt  B3LYP/6- 311+G(2d,p) freq  电荷0 自旋多重度1
3a*的triplet优化是# opt uB3LYP/6- 311+G(2d,p) td freq    电荷0 自旋多重度3
如上设置关键词合适吗?3a*的triplet优化B3LTP能加U吗?

其次,文献提到“Thecomputed energy gap (58.5 kcal mol−1) was lower than thetriplet energy of Ir*[dF(CF3)ppy]2(dtbbpy)+ (60.1 kcalmol−1), suggesting that energy transfer from the catalyst tothe substrate was feasible” energy gap应该是指上面3a的singlet和3a*的triplet的差,是不是3a的singlet优化时需要加td才能和3a*的triplet相比较?

1593258751(1).png (13.36 KB, 下载次数 Times of downloads: 45)

1593258751(1).png

3

帖子

0

威望

11

eV
积分
14

Level 1 能力者

22#
发表于 Post on 2025-11-24 10:40:10 | 只看该作者 Only view this author
wzkchem5 发表于 2025-11-24 10:21
计算磷光发射能,需要在优化的T1结构下计算T1能量和S0能量之差。
如果你优化T1的方法本身就能给出T1能量 ...

谢谢老师

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

21#
发表于 Post on 2025-11-24 10:21:07 | 只看该作者 Only view this author
lzy5260 发表于 2025-11-23 15:53
谢谢老师
我理解是你的那一步只是优化了T1

计算磷光发射能,需要在优化的T1结构下计算T1能量和S0能量之差。
如果你优化T1的方法本身就能给出T1能量和S0能量之差,那么只需跑优化T1这一个计算,一个计算管两件事。至于准不准是另一码事
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
本团队长期招收研究生,有意者可私信联系

3

帖子

0

威望

11

eV
积分
14

Level 1 能力者

20#
发表于 Post on 2025-11-23 15:53:09 | 只看该作者 Only view this author
wzkchem5 发表于 2025-11-23 12:01
只做# B3LYP/6-31G* opt(自旋多重度3)确实得不到磷光发射能。为了得到磷光发射能,用# B3LYP/6-31G* TD ...

谢谢老师
我理解是你的那一步只是优化了T1

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

19#
发表于 Post on 2025-11-23 12:01:17 | 只看该作者 Only view this author
lzy5260 发表于 2025-11-23 11:36
请问你写的那个进行T1的结构优化了吗?还是说把下面2步合成一步了
1 优化T1结构:# B3LYP/6-31G* opt, ...

只做# B3LYP/6-31G* opt(自旋多重度3)确实得不到磷光发射能。为了得到磷光发射能,用# B3LYP/6-31G* TD(triplet)(自旋多重度1)是一种做法,也可以用# B3LYP/6-31G* opt(自旋多重度3)的能量减去# B3LYP/6-31G* opt(自旋多重度1)的能量。当泛函趋于精确时,两者结果趋于等价,但对于一般的近似泛函两者不等价。此外后一种做法还有一个缺点:难以得到磷光速率。
除此之外还有第三种方法:# B3LYP/6-31G* opt TD(triplet),一个输入文件搞定,缺点是耗时长,结构优化收敛可能困难,且结果未必比前两种更准。
至于什么情况下用哪种方法更好,参见https://mattermodeling.stackexch ... -triplets/9357#9357
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
本团队长期招收研究生,有意者可私信联系

3

帖子

0

威望

11

eV
积分
14

Level 1 能力者

18#
发表于 Post on 2025-11-23 11:36:10 | 只看该作者 Only view this author
喵星大佬 发表于 2020-6-27 21:12
优化激发态结构不需要写td,只要你自旋写的不是1程序会默认用udft所以写不写ub3lyp无所谓
不用加弥散函数 ...

请问你写的那个进行T1的结构优化了吗?还是说把下面2步合成一步了
1 优化T1结构:# B3LYP/6-31G* opt,电荷和自旋多重度为0 3。
2 取上一步优化得到的结构,用# B3LYP/6-31G* TD(triplet),电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。
我刚学不太懂,望解答下,谢谢

7

帖子

0

威望

255

eV
积分
262

Level 3 能力者

17#
发表于 Post on 2023-6-1 21:53:20 | 只看该作者 Only view this author
wzkchem5 发表于 2023-6-1 16:56
这是因为双杂化泛函TDDFT计算是先求解Casida方程,再加doubles correction。双杂化泛函TDDFT的激发能虽然 ...

好的我明白了,我仔细查看了out文件,找到了相关的部分。感谢老师您的解答。

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

16#
发表于 Post on 2023-6-1 16:56:36 | 只看该作者 Only view this author
symbolhone 发表于 2023-6-1 02:29
wzkchem5老师,您好。我最近在用双杂化泛函计算TADF分子ΔEST,我是先用gaussian做的opt,然后再用ORCA做 ...

这是因为双杂化泛函TDDFT计算是先求解Casida方程,再加doubles correction。双杂化泛函TDDFT的激发能虽然应该读包括了doubles correction的值,但是是按Casida方程的解的大小顺序排序的,而不是按包括了doubles correction的值的大小顺序排序的。所以如果只看最终的激发能的话,可能不是单调递增的,这是正常现象。此时应该把你的STATE 2叫T1,把你的STATE 1叫T2。
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
本团队长期招收研究生,有意者可私信联系

7

帖子

0

威望

255

eV
积分
262

Level 3 能力者

15#
发表于 Post on 2023-6-1 09:29:13 | 只看该作者 Only view this author
wzkchem5 发表于 2023-4-18 17:15
建议用高精度方法算单点能,例如双杂化泛函、CC2、ADC(2)、DLPNO-STEOM-CCSD、NEVPT2等

wzkchem5老师,您好。我最近在用双杂化泛函计算TADF分子ΔEST,我是先用gaussian做的opt,然后再用ORCA做的双杂话泛函计算,其中一个分子,它的out里显示这样:
TD-DFT/TDA EXCITED STATES (SINGLETS)
STATE  1:  E=   0.099942 au      2.720 eV    21934.7 cm**-1 <S**2> =   0.000000
STATE  2:  E=   0.116849 au      3.180 eV    25645.5 cm**-1 <S**2> =   0.000000
STATE  3:  E=   0.126277 au      3.436 eV    27714.6 cm**-1 <S**2> =   0.000000

TD-DFT/TDA EXCITED STATES (TRIPLETS)
STATE  1:  E=   0.102962 au      2.802 eV    22597.6 cm**-1 <S**2> =   2.000000
STATE  2:  E=   0.096927 au      2.638 eV    21273.0 cm**-1 <S**2> =   2.000000
STATE  3:  E=   0.114811 au      3.124 eV    25198.2 cm**-1 <S**2> =   2.000000

其中STATE  1 (TRIPLETS)的能量是2.802eV, STATE  2 (TRIPLETS)则是2.638eV,这个是什么原因导致的呢,出现这种情况是不是意味着我这个计算任务就需要重新进行呢。

10

帖子

0

威望

161

eV
积分
171

Level 3 能力者

14#
发表于 Post on 2023-4-19 09:21:05 | 只看该作者 Only view this author
wzkchem5 发表于 2023-4-18 17:15
建议用高精度方法算单点能,例如双杂化泛函、CC2、ADC(2)、DLPNO-STEOM-CCSD、NEVPT2等

谢谢老师的建议,我尝试一下。

1万

帖子

0

威望

9739

eV
积分
21935

Level 6 (一方通行)

13#
发表于 Post on 2023-4-18 17:15:29 | 只看该作者 Only view this author
柚子 发表于 2023-4-18 10:09
感谢老师的回答。
这两种方法我都做了,但是两个能量相差较大。研究对象是一个TADF分子,计算得到的S1态 ...

建议用高精度方法算单点能,例如双杂化泛函、CC2、ADC(2)、DLPNO-STEOM-CCSD、NEVPT2等
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
本团队长期招收研究生,有意者可私信联系

10

帖子

0

威望

161

eV
积分
171

Level 3 能力者

12#
发表于 Post on 2023-4-18 17:09:39 | 只看该作者 Only view this author
zjxitcc 发表于 2023-4-18 16:31
有两种计算方式,以高斯为例:
(1)UDFT计算,坐标上面一行的自旋多重度写3,不要搭配TD关键词;
(2 ...

感谢老师的回答。
这两种方法我都做了,但是两个能量相差较大。研究对象是一个TADF分子,计算得到的S1态能量为40kcal/mol,我是以基态单重态的能量为基准,用UDFT得到的能量为33kcal/mol,用TDDFT得到的能量为40kcal/mol。
请问老师这种情况下怎么选择呢?

4279

帖子

4

威望

9473

eV
积分
13832

Level 6 (一方通行)

MOKIT开发者

11#
发表于 Post on 2023-4-18 16:31:27 | 只看该作者 Only view this author
柚子 发表于 2023-4-18 16:10
老师您好,请问如果基态为闭壳层单重态,DFT,若要计算T1态的单点能,自旋多重度是“0 1”还是“0 3”呢 ...

有两种计算方式,以高斯为例:
(1)UDFT计算,坐标上面一行的自旋多重度写3,不要搭配TD关键词;
(2)TDDFT计算,坐标上面一行的自旋多重度写1,搭配TD(triplet)关键词。
自动做多参考态计算的程序MOKIT

10

帖子

0

威望

161

eV
积分
171

Level 3 能力者

10#
发表于 Post on 2023-4-18 16:10:12 | 只看该作者 Only view this author
zjxitcc 发表于 2023-2-11 15:44
(1)是自旋,不是自选;
(2)说清楚要用什么理论方法算,以及基态是否为闭壳层单重态

老师您好,请问如果基态为闭壳层单重态,DFT,若要计算T1态的单点能,自旋多重度是“0 1”还是“0 3”呢?

4279

帖子

4

威望

9473

eV
积分
13832

Level 6 (一方通行)

MOKIT开发者

9#
发表于 Post on 2023-2-11 15:44:12 | 只看该作者 Only view this author
yiwind 发表于 2023-2-11 14:27
我得到了 三重态的结构,现在想计算三重态的单点能,自选多重度是设置 “0 1”还是设置为“0 3”呀?

(1)是自旋,不是自选;
(2)说清楚要用什么理论方法算,以及基态是否为闭壳层单重态
自动做多参考态计算的程序MOKIT

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

GMT+8, 2026-1-24 11:43 , Processed in 0.182733 second(s), 24 queries , Gzip On.

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