计算化学公社

 找回密码 Forget password
 注册 Register
Views: 7184|回复 Reply: 1

[Gaussian/gview] 简化用了IOp(9/40=x)的Gaussian的CIS/TDDFT任务的输出文件的程序:simpIOp940

[复制链接 Copy URL]

4万

帖子

99

威望

4万

eV
积分
89933

管理员

公社社长+计算化学玩家

发表于 Post on 2019-2-19 19:50:24 | 显示全部楼层 Show all |阅读模式 Reading model
简化用了IOp(9/40=x)的Gaussian的CIS/TDDFT任务的输出文件的程序:simpIOp940

文/Sobereva@北京科音  2019-Feb-19


简介

《Multiwfn支持的电子激发分析方法一览》(http://sobereva.com/437)一文全面介绍了Multiwfn支持的丰富的电子激发分析,其中有些功能是依赖于组态系数的。由于Gaussian默认情况下只输出绝对值大于0.1的组态系数,光靠这些对激发态描述的精度不够,因此对于Gaussian用户,用Multiwfn做那些电子激发分析时一般需要用诸如IOp(9/40=x)关键词来输出所有绝对值大于10^-x的组态系数(一般x=4)。但是此时造成的问题就是输出文件里的组态系数特别多,有人觉得通过这样的文件人工查看激发态主要信息不方便。

为解决这个问题,笔者写了个小程序simpIOp940,可在此下载:http://sobereva.com/soft/simpIOp940.rar

此程序可以载入用了IOp(9/40=x)的Gaussian输出文件,然后自行输入一个组态系数绝对值的阈值,程序就会把激发态的信息连同大于阈值的组态系数都输出出来,使得人工查看比较方便。另外,这个程序还会把各个激发态里系数绝对值最大的轨道跃迁输出出来,令指认主要跃迁类型更为方便。

此程序文件包里带.exe后缀的是Windows版可执行文件,不带后缀的是Linux版可执行文件,Fortran源文件也附上了。

例子

例如此程序的压缩包里有一个uracil.gjf,TDDFT计算时用了IOp(9/40=4),算出来的文件时uracil.out,里面列出的组态系数特别多:
Excited State   1:      Singlet-A"     4.7968 eV  258.47 nm  f=0.0001  <S**2>=0.000
       9 -> 30        -0.00071
       9 -> 31        -0.00043
       9 -> 35         0.00028
       9 -> 67        -0.00011
       9 -> 73        -0.00016
       9 -> 76         0.00012
       9 -> 78        -0.00015
       9 -> 79        -0.00020
       9 -> 81        -0.00011
       9 -> 82         0.00023
       9 -> 88        -0.00021
       9 -> 89        -0.00059
...略

启动simpIOp940之后,输入uracil.out的路径,然后输入阈值0.1,当前目录下会产生new.out,这是处理好的结果文件,内容如下

Excited State   1:      Singlet-A"     4.7968 eV  258.47 nm  f=0.0001  <S**2>=0.000
       26 -> 30         0.11588
       28 -> 30         0.67734
       28 -> 31         0.13625

Excited State   2:      Singlet-A'     5.4219 eV  228.67 nm  f=0.1317  <S**2>=0.000
       27 -> 30        -0.13650
       27 -> 31        -0.14060
       29 -> 30         0.67086

Excited State   3:      Singlet-A"     6.0180 eV  206.02 nm  f=0.0000  <S**2>=0.000
       26 -> 30         0.54135
       26 -> 31        -0.20634
       28 -> 30        -0.15424
       28 -> 31         0.36715


当前目录下还出现了largest_pair.txt:
       1      28      30   0.67734
       2      29      30   0.67086
       3      26      30   0.54135

第一列是激发态序号,第2、3列是贡献最大的MO对,最后一列是系数。

由于电子激发任务类型很多,此程序不一定兼容所有情况,碰见不兼容时请自行修改程序。

评分 Rate

参与人数
Participants 5
eV +25 收起 理由
Reason
rendong + 5 谢谢
Adair + 5 好物!
liyuanhe211 + 5
哇哇吐 + 5 我之前都单独再算一遍,谢谢老师
Novice + 5 好物!

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口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!

4万

帖子

99

威望

4万

eV
积分
89933

管理员

公社社长+计算化学玩家

 楼主 Author| 发表于 Post on 2019-9-26 00:27:07 | 显示全部楼层 Show all
更新了一下此程序和此文,新版本会输出largest_pair.txt文件,包含各个激发态贡献最大的MO对,在基于MO判断激发主要特征时更方便。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办最高质量的各种计算化学类培训:初级量子化学培训班基础(中级)量子化学培训班分子动力学与GROMACS培训班量子化学波函数分析与Multiwfn程序培训班,内容介绍以及往届资料购买请点击链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的最佳途径。培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取培训最新消息、避免错过网上最有价值的计算化学文章!
欢迎加入人气最高、水准最高的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人,讨论范畴相同
思想家公社的门口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!

本版积分规则 Credits rule

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

GMT+8, 2023-2-5 10:08 , Processed in 0.182688 second(s), 22 queries .

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