计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: sobereva
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法

  [复制链接 Copy URL]

38

帖子

0

威望

265

eV
积分
303

Level 3 能力者

136#
发表于 Post on 2021-4-25 10:47:45 | 只看该作者 Only view this author
sobereva 发表于 2021-4-25 00:10
这种问题应当单独发主贴

好的,sob老师,不好意思,下次会注意!

38

帖子

0

威望

265

eV
积分
303

Level 3 能力者

137#
发表于 Post on 2021-4-25 10:49:57 | 只看该作者 Only view this author
wzkchem5 发表于 2021-4-25 01:38
支持在CCSD(T)级别做激发态(EOM-CCSD(T))的程序很少,高斯是不支持的,这就是你报错的原因。
如果你要 ...

好的,太感谢老师了!谢谢!

1632

帖子

4

威望

4088

eV
积分
5800

Level 6 (一方通行)

喵星人

138#
发表于 Post on 2021-4-29 12:29:31 | 只看该作者 Only view this author
计算吸收光谱,体系总体是个阴离子,但是,吸收光的部分应该是个阳离子,应该是非共轭链后面的磺酸阴离子集中带负电荷,这种体系计算激发能的时候有必要加弥散么

5万

帖子

99

威望

5万

eV
积分
112351

管理员

公社社长

139#
 楼主 Author| 发表于 Post on 2021-4-29 15:00:24 | 只看该作者 Only view this author
喵星大佬 发表于 2021-4-29 12:29
计算吸收光谱,体系总体是个阴离子,但是,吸收光的部分应该是个阳离子,应该是非共轭链后面的磺酸阴离子集 ...

不用。为了稳妥的话,也可以只在磺酸阴离子部分加弥散。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

44

帖子

0

威望

1171

eV
积分
1215

Level 4 (黑子)

140#
发表于 Post on 2021-4-29 21:20:26 | 只看该作者 Only view this author
拜读这篇帖子和博文多次,又恭敬学习了其他多篇博文,今天突然彻悟了,非常感谢您

44

帖子

0

威望

1171

eV
积分
1215

Level 4 (黑子)

141#
发表于 Post on 2021-4-30 21:43:02 | 只看该作者 Only view this author
sob老师您好,今天再次学习时有一个疑惑,在优化 T1 结构时,自旋多重度设置成 3,这时需不需要加 “U”,在开壳层下进行优化?

1632

帖子

4

威望

4088

eV
积分
5800

Level 6 (一方通行)

喵星人

142#
发表于 Post on 2021-4-30 22:56:48 | 只看该作者 Only view this author
rendong 发表于 2021-4-30 21:43
sob老师您好,今天再次学习时有一个疑惑,在优化 T1 结构时,自旋多重度设置成 3,这时需不需要加 “U”, ...

3的时候已经默认是U了

44

帖子

0

威望

1171

eV
积分
1215

Level 4 (黑子)

143#
发表于 Post on 2021-5-1 23:51:17 | 只看该作者 Only view this author
喵星大佬 发表于 2021-4-30 22:56
3的时候已经默认是U了

原来是这样,谢谢您释疑,

12

帖子

0

威望

362

eV
积分
374

Level 3 能力者

144#
发表于 Post on 2021-6-2 22:14:56 | 只看该作者 Only view this author
Sob老师,咨询一个关于T1磷光激发能的计算
按方法一计算:
1 优化T1结构:# B3LYP/6-31G* opt,电荷和自旋多重度为0 3。
2 取上一步优化得到的结构,用# B3LYP/6-31G* TD(triplet),电荷和自旋多重度为0 1。第一激发态的激发能就是磷光发射能了。

我计算了两个配合物,第一步用UDFT优化T1,而后进行方法1的第二步,为什么第一个配合物的第一激发态是<S**2>=2,第二个配合物的第一激发态是<S**2>=0.112。我看有的帖子说这是由于第一配合物的参考态是RDFT,而第二个配合物的参考态是UDFT。想知道为什么会出现这种原因。怎么解决,是否可以让第二个配合物算出来的第一激发态的自旋多重度为3,还是说对于UDFT的基函数算出来的磷光,只能自己去指认三重态,找发射能。

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

145#
发表于 Post on 2021-6-2 22:31:02 | 只看该作者 Only view this author
Mance 发表于 2021-6-2 15:14
Sob老师,咨询一个关于T1磷光激发能的计算
按方法一计算:
1 优化T1结构:# B3LYP/6-31G* opt,电荷和自 ...

要自己根据<S**2>找第一个三重态,不能盲目取第一个激发态
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

12

帖子

0

威望

362

eV
积分
374

Level 3 能力者

146#
发表于 Post on 2021-6-3 15:49:50 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-2 22:31
要自己根据找第一个三重态,不能盲目取第一个激发态

好的,谢谢老师的解答,我还想问一下,同样用UDFT优化完T1后,用TD(Triplet) 自旋多重度为1,计算磷光,是什么原因导致在第二步计算的时候,一个配合物为开壳层,一个配合物为闭壳层

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

147#
发表于 Post on 2021-6-3 16:46:08 | 只看该作者 Only view this author
Mance 发表于 2021-6-3 08:49
好的,谢谢老师的解答,我还想问一下,同样用UDFT优化完T1后,用TD(Triplet) 自旋多重度为1,计算磷光, ...

把完整的输入输出文件放上来
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

12

帖子

0

威望

362

eV
积分
374

Level 3 能力者

148#
发表于 Post on 2021-6-3 22:00:54 | 只看该作者 Only view this author
本帖最后由 Mance 于 2021-6-5 11:52 编辑
wzkchem5 发表于 2021-6-3 16:46
把完整的输入输出文件放上来

请管理员删除

UDFT-T1优化关键词.png (6.46 KB, 下载次数 Times of downloads: 58)

UDFT-T1优化关键词.png

Ru.rar

986.77 KB, 下载次数 Times of downloads: 4

1万

帖子

0

威望

7391

eV
积分
18144

Level 6 (一方通行)

149#
发表于 Post on 2021-6-3 22:22:47 | 只看该作者 Only view this author
Mance 发表于 2021-6-3 15:00
麻烦了老师,T1使用UDFT优化的,磷光计算的输入输出文件已经添加至附件。万分感谢

名字里有mpa的那个输入文件读了初猜波函数,另一个没有读。
如果初猜波函数是UDFT的,即使自旋多重度是0 1,程序也会做UDFT计算而不是RDFT计算。
不过最后收敛的S**2不是0或2,说明这个体系的闭壳层单重态波函数是不稳定的。如果fmp2和mpa是很类似的体系,建议检查一下fmp2那个计算的波函数的稳定性,如果是不稳定的,需要找到稳定的波函数才能计算。
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员

12

帖子

0

威望

362

eV
积分
374

Level 3 能力者

150#
发表于 Post on 2021-6-3 22:56:18 | 只看该作者 Only view this author
wzkchem5 发表于 2021-6-3 22:22
名字里有mpa的那个输入文件读了初猜波函数,另一个没有读。
如果初猜波函数是UDFT的,即使自旋多重度是0 ...

是的,在优化了T1结构后,发现mpa在开始的时候磷光激发能为负值,所以用stable=opt,自旋多重度为1做了单点,然后读取了初猜波函数进行的磷光计算,他的波函数就变成开壳层了。然后,我也发现用stable=opt做单点的时候就变是UDFT。
老师,你说的mpa体系闭壳层单重态波函数不稳定,需要找到稳定波函数,这个如何解决,是重新优化T1几何结构,找最低点吗?还是重S0开始优化。
fmp2是能计算出S**2磷光发射能的,这个是需要重新验证它的波函数稳定性是吗?

本版积分规则 Credits rule

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

GMT+8, 2024-11-23 10:02 , Processed in 0.198129 second(s), 23 queries , Gzip On.

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