请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 667|回复: 1

[Multiwfn资源与经验] 使用Multiwfn便利地查看所有激发态中的主要轨道跃迁贡献

[复制链接]

2万

帖子

25

威望

2万

eV
积分
52358

管理员

公社社长

发表于 2020-1-28 04:22:48 | 显示全部楼层 |阅读模式
使用Multiwfn便利地查看所有激发态中的主要轨道跃迁贡献

文/Sobereva@北京科音  2020-Jan-28


之前笔者写过一篇文章《电子激发任务中轨道跃迁贡献的计算》(http://sobereva.com/230),介绍了如何对TDDFT及类似方法计算的激发态计算轨道跃迁的贡献。在很多电子激发的研究文章中都会给出各个激发态中各种主要MO跃迁情况。虽然计算贡献很简单,但当考察的态比较多的时候,一个一个去考察、记录贡献值还是比较费事的,初学者还容易弄错。

从2020-Jan-28更新的Multiwfn开始,在电子激发分析(主功能18)里加入了子功能15,专门用于快速输出所有激发态中的主要的轨道跃迁贡献,使得从MO角度考察各个激发态的特征非常简单容易。此功能支持Gaussian、ORCA、GAMESS-US、Firefly的CIS/TDHF/TDA-DFT/TDDFT/ZINDO的输出文件。Multiwfn可以在其主页http://sobereva.com/multiwfn免费下载。不熟悉Multiwfn者推荐阅读《Multiwfn FAQ》(http://sobereva.com/452)。

下面是个简单的例子。启动Multiwfn后输入以下内容
examples\excit\D-pi-A.out  //Gaussian做标准TDDFT任务的输出文件(其中的IOp(9/40=4)关键词对于当前功能是多余的)
18  //电子激发分析
15  //输出各个激发态中主要的轨道跃迁贡献
然后马上看到以下信息
HOMO index:    56
LUMO index:    57

Only MO transitions with absolute contribution >=  5.0 % are shown below. It co
rresponds to 10 times of "compthres" parameter in settings.ini

#   1    3.907 eV  Multi= 1: H-4 -> L 81.9%, H-4 -> L+2 12.1%
#   2    4.062 eV  Multi= 1: H -> L 86.0%, H-3 -> L 5.3%
#   3    4.417 eV  Multi= 1: H-6 -> L 85.3%, H-6 -> L+2 11.9%
#   4    4.791 eV  Multi= 1: H-2 -> L 54.5%, H -> L+1 27.6%, H-3 -> L+1 6.4%
#   5    4.887 eV  Multi= 1: H -> L+3 57.3%, H-2 -> L 17.0%, H-1 -> L+2 8.8%, H
-1 -> L 8.0%


由输出可见,比如对于第3激发态,主要特征是HOMO-6 -> LUMO贡献85.3%,其次是HOMO-6 -> LUMO+2贡献11.9%,输出的格式非常明确直观。Multi代表自旋多重度,激发态序号和激发能也都一起给出了。

之后如果你输入y,上面的信息将会被导出到当前目录下的一个文本文件里,便于后处理。

默认情况下,如提示所示,只有MO跃迁贡献大于5%的才会被输出。如果你想提升或者降低阈值,可以修改settings.ini里的compthres,当前功能输出阈值是compthres值的10倍。

此功能输出的信息可以放到你写的文章的补充材料去,目前有不少文章都这么做。

下面是基于一个Gaussian对开壳层体系做TDDFT的输出文件,用此功能输出的信息
#   1    2.274 eV  Multi= ?: Hb-1 -> Lb 96.6%
#   2    2.851 eV  Multi= ?: Hb -> Lb 90.7%, Ha -> La 8.4%
#   3    3.453 eV  Multi= ?: Hb-2 -> Lb 69.7%, Ha -> La+1 27.6%
#   4    4.191 eV  Multi= ?: Ha -> La 85.5%, Hb -> Lb 7.6%
#   5    4.404 eV  Multi= ?: Ha-1 -> La 48.5%, Hb -> Lb+1 40.4%
#   6    5.164 eV  Multi= ?: Hb-3 -> Lb 96.3%
#   7    5.267 eV  Multi= ?: Ha -> La+1 65.4%, Hb-2 -> Lb 22.4%, Hb -> Lb+1 8.8%
#   8    5.686 eV  Multi= ?: Hb-4 -> Lb 89.1%, Hb-1 -> Lb+1 6.8%

由于参考态是开壳层时TDDFT算的激发态有自旋污染,因此Multi后面是问号。当前情况轨道自旋都明确标出了,a、b分别代表alpha和beta。比如基态到第5激发态中alpha的HOMO-1到alpha的LUMO跃迁贡献了48.5%。

如果你是ORCA用户,当前功能显得更有意义了。ORCA里有个令很多人讨厌的习俗是轨道序号是从0开始记,导致很多人搞错HOMO和LUMO的序号,指使在文章中描述激发态的MO跃迁情况时出现错误。而使用Multiwfn查看主要MO跃迁贡献就可以完全杜绝这个问题。

如果你是GAMESS-US用户,别忘了要把电子激发任务输出文件设成.gms后缀,否则此功能无法正常使用。比如examples\excit\H2CO_TDDFT_GAMESS.gms就是GAMESS-US对甲醛做TDDFT的输出文件,可以直接给Multiwfn作为输入文件使用当前功能。

评分

参与人数 2eV +6 收起 理由
liyuanhe211 + 1 GJ!之前还特意自己写了一个
zsu007 + 5 赞!

查看全部评分

北京科音自然科学研究中心http://www.keinsci.com  致力于计算化学的发展和传播,长期开办最高水准的各种量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训,是提升计算化学研究水平的最佳选择。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯!培训相关信息见《北京科音办的培训班FAQ》(http://bbs.keinsci.com/thread-5098-1-1.html)。
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群:思想家公社QQ群1号:18616395,2号:466017436。合计约6000人。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一定会被拒绝加入。
思想家公社的门口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!

21

帖子

0

威望

414

eV
积分
435

Level 3 能力者

发表于 2020-1-28 11:15:16 | 显示全部楼层
谢谢社长。当前的输出结果中跃迁能量部分是以ev单位给出的,如果同时加上nm单位的能量值会更方便,因为光谱实验数据中多用nm。当然这不是什么大问题,自己看看log文件一对照就可以了,或者用excel一转换就行。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-4-3 12:27 , Processed in 0.155431 second(s), 25 queries .

快速回复 返回顶部 返回列表