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

sobereva

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

simpIOp940: A code to simplify the output file of CIS/TDDFT task of Gaussian with IOp(9/40=x)

文/Sobereva@北京科音

First release: 2019-Feb-19  Last update: 2023-Jun-13

简介

《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对，最后一列是系数。

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

sobereva

更新了此程序，支持了开壳层TDDFT的情况

sobereva

更新了一下此程序和此文，新版本会输出largest_pair.txt文件，包含各个激发态贡献最大的MO对，在基于MO判断激发主要特征时更方便。