程序中有一片段
程序中有一片段
% chk=test180
(井字号)p rhf/gen ecp 5d 7f test geom=modela
Formaldehyde 5d 7f energy with ECP
0,1
c h h o
C O
1p-31g*
F 1 1.00
0.1000000000D+01 0.1000000000D+01
* * * *
我刚接触不甚了解.
%chk=test180
这一行指定了一个读取和存储计算中必要的数据的文件test180.
#p rhf/genecp 5d 7f test geom=modela
这一行指定计算的任务:p指详细输出计算的结果;计算方法为rhf,也就是闭壳层的自洽场计算;genecp是一个关键字,用来自定义基组和有效芯势(effective core potential),因此任务行里没有指定基组;5d和7f表示自定义基组使用纯d和纯f函数(可以类比于真实亚层里的5个d轨道和7个f轨道);test用来禁止在完成计算后自动存档;geom=modela中geom表示指定计算体系的几何构型,modela表示读取模型建立器的原子间的连接信息来构造Z矩阵.
Formaldehyde 5d 7f energy with ECP
这是任务的标题栏,对实际计算没影响.
0,1
0表示体系的总电荷为0;1表示体系多重度为1.
c h h o
表示体系有C、H、H、O4个原子.因为Z矩阵由模型建立器的信息构造,所以分子说明里的几何构型部分没出现Z矩阵.
从“C 0”起的部分与genecp关键字对应,用来定义基组:第一行里C表示一下的定义用于C(碳)原子,还可以同时指定自定义基组其它的应用范围,最后用0(数字,不是字母O)结束;第二行l(字母,不是数字1)p-31g*是一个基组的名称,表示在这个基组的基础上自定义新基组;第三行里F表示添加自定义的F壳层,1指用1个高斯函数来表示这个壳层,1.00是换算因子f;下一行的2个数字0.1000000000D+01分别指高斯函数的指数α和收缩因子d.f、α、d3个参数可以表明通过高斯函数近似壳层轨道的方式φ(r)=∑d·exp(-α·f^2·r^2).