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

简化用了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
 
由于电子激发任务类型很多,此程序不一定兼容所有情况,碰见不兼容时请自行修改程序。

评论已关闭