Gaussian单点能自动相互运算工具enecalc

Gaussian单点能自动相互运算工具enecalc

文/Sobereva(3)  2015-Nov-3


1.0版下载地址:/usr/uploads/file/20151103/20151103190920_48250.rar


1 用途

enecalc程序用来对指定的结构在不同级别下调用Gaussian来计算单点能,然后根据预设的系数进行相互运算,得到最终能量,很适合用来做基组外推计算。

2 压缩包内容

enecalc.exe:编译好的Windows版可执行文件
enecalc:编译好的Linux 64bit版可执行文件
enecalc.f90:源代码。如果自行编译可以诸如这样:gfortran enecalc.f90 -o enecalc
settings.ini:用来设定运行参数
template.gjf:模板文件,里面包含体系结构和内存、内核数等设定,关键词部分留空。

3 实例1:一般情况

例如对某个结构,按照以下外推公式计算能量:
E[QCISD(T)/CBS]=1.4629*E[QCISD(T)/TZ]-0.4629*E[QCISD(T)/DZ]+1.6938*E[MP2/QZ]-2.1567*E[MP2/TZ]+0.4629*E[MP2/DZ]

可以将settings.ini写为
"d:\study\g09w\g09"   // Gaussian的路径,注意用双引号扩住
1   // 设定当前的系统,1为Windows,2为Linux或MacOS
0   // enecalc支持ONIOM计算,如果是ONIOM任务写1,否则为0
5   // 要考虑的计算级别数,可以为无穷多
1.4629 -0.4629 1.6938 -2.1567 0.4629     // 每个级别的能量要乘的系数
QCISD(T)/cc-pVTZ  // 第1个级别计算时用的关键词,后同
QCISD(T)=   // enecalc从输出文件末尾读取能量,这里设定用于定位第1个级别能量的字符串,后同
QCISD(T)/cc-pVDZ
QCISD(T)=
MP2/cc-pVQZ
MP2=
MP2/cc-pVTZ
MP2=
MP2/cc-pVDZ
MP2=

以下是模板文件template.gjf示例
%mem=1500MB
%nproc=4
#

B3LYP/6-31G* opted

0 1
 O                  0.33287540   -1.10635328    1.50332174
 H                  0.33287540   -0.34740028    0.90678174
 H                  0.33287540   -1.86530628    0.90678174


启动enecalc之后,程序先读取当前目录下的settings.ini,然后将第一个计算级别的关键词套入当前目录下的模板文件template.gjf里,生成001.gjf,然后调用Gaussian产生001.out,从中读取能量并显示出来。然后再如此处理第二个级别,直至全部级别都处理完,最后输出根据定义的规则运算出最终总能量。

4 实例2:ONIOM的情况

此例和上例一样进行能量外推计算,但是是基于ONIOM计算。下面是settings.ini文件,注意不需要写用于定位的字符串,因为enecalc直接读取输出文件中的ONIOM: extrapolated energy =后面的值。
"g09"
2
1
5
1.4629 -0.4629 1.6938 -2.1567 0.4629
oniom(CCSD(t)/cc-pVTZ:M062X/TZVP)
oniom(CCSD(t)/cc-pVDZ:M062X/TZVP)
oniom(MP2/cc-pVQZ:M062X/TZVP)
oniom(MP2/cc-pVTZ:M062X/TZVP)
oniom(MP2/cc-pVDZ:M062X/TZVP)

模板文件template.gjf示例:
%mem=7GB
%nproc=4
#

Title Card Required

0 1
 C                0    0.87378638   -1.26213590    0.00000000 L
 H                0    1.23044081   -2.27094591    0.00000000 L
 H                0    1.23045922   -0.75773771   -0.87365150 L
 H                0   -0.19621362   -1.26212272    0.00000000 L
 C                0    1.38712860   -0.53617963    1.25740497 L H 8
 H                0    1.03207226    0.47319346    1.25642745 L
 H                0    1.02885949   -1.03944806    2.13105486 L
 C                0    2.92712600   -0.53863598    1.25881364 H
 H                0    3.28539448   -0.03569410    0.38497547 H
 H                0    3.28218223   -1.54800868    1.26016738 H
 O                0    3.40380116    0.13590181    2.42615231 H
 H                0    4.36379924    0.13453019    2.42749555 H

添加新评论