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

计算化学公社

 找回密码
 现在注册!
查看: 468|回复: 0

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

[复制链接]

1万

帖子

25

威望

2万

eV
积分
44694

管理员

公社社长

发表于 2019-2-19 19:50:24 | 显示全部楼层 |阅读模式
简化用了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

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

评分

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

查看全部评分

北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。达5000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

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

GMT+8, 2019-8-21 09:48 , Processed in 0.146848 second(s), 25 queries .

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