关于含赝势的分子wfn输出文件格式问题

兵库北先生

起因由于在用相关软件读取gaussian输出的wfn文件时，发现读取错误。

究其原因，发现我这个含Pd的分子的wfn输出文件是这样一个格式：

在一堆centre assignmnets, type assignments 和 Exponents之后，MO之前。

还有几个centre assignmnets, type assignments 和 Exponents之类的数据，如下所示（因文件太大，就截取一段）。

……
EXPONENTS  0.3272304D+00 0.1030724D+00 0.8000000D+00 0.8000000D+00 0.8000000D+00
EXPONENTS  0.3406134D+02 0.5123575D+01 0.1164663D+01 0.3272304D+00 0.1030724D+00
EXPONENTS  0.8000000D+00 0.8000000D+00 0.8000000D+00 0.3406134D+02 0.5123575D+01
EXPONENTS  0.1164663D+01 0.3272304D+00 0.1030724D+00 0.8000000D+00 0.8000000D+00
EXPONENTS  0.8000000D+00 0.3406134D+02 0.5123575D+01 0.1164663D+01 0.3272304D+00
EXPONENTS  0.1030724D+00 0.8000000D+00 0.8000000D+00 0.8000000D+00
CENTRE ASSIGNMENTS   22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22 22
CENTRE ASSIGNMENTS   22 22 22 22 22
TYPE ASSIGNMENTS      1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1
TYPE ASSIGNMENTS      1  1  1  1  1
EXPONENTS  0.4459566D+08 0.2277439D+08 0.1163057D+08 0.5939569D+07 0.3033256D+07
EXPONENTS  0.1549042D+07 0.7910743D+06 0.4039907D+06 0.2063125D+06 0.1053609D+06
EXPONENTS  0.5380637D+05 0.2747817D+05 0.1403272D+05 0.7166317D+04 0.3659739D+04
EXPONENTS  0.1868978D+04 0.9544614D+03 0.4874302D+03 0.2489239D+03 0.1271220D+03
EXPONENTS  0.6491944D+02 0.3315346D+02 0.1693102D+02 0.8646436D+01 0.2254994D+01
MO    1     MO 0.0        OCC NO =    2.0000000  ORB. ENERGY =  -77.353055
  0.20317665D-06  0.37920772D-06  0.64720622D-06  0.10058231D-05  0.13472289D-05
  0.13838959D-05 -0.16330708D-05 -0.48392617D-06 -0.77666372D-04  0.45668736D-03
  0.92570547D-03 -0.34574864D-05 -0.36903999D-05 -0.32318248D-05 -0.20408571D-05
-0.71645728D-05 -0.76472141D-05 -0.66969588D-05 -0.42290460D-05 -0.11426885D-04

……

22是我的Pd原子，我在查完basissetexchange之后发现Pd原子的高斯型函数已经在前面一堆centre assignmnets, type assignments 和 Exponents之中都有了，而且相应的PRIMITIVES数量也都对的上。

不知道这几个多出来的几行是什么意思，删去之后可以正确运行相关软件并读入该wfn文件。希望大佬解惑！

兵库北先生

sobereva 发表于 2019-8-20 22:11
估计是EDF字段，原本是wfx文件里才有的，现在可能在wfn文件里也有了。EDF字段在这里有说明
在赝势下做波函 ...

感谢sob老师！

sobereva

估计是EDF字段，原本是wfx文件里才有的，现在可能在wfn文件里也有了。EDF字段在这里有说明
在赝势下做波函数分析的一些说明
http://sobereva.com/156

可以用Multiwfn载入fch文件，用主功能100的子功能2导出成wfn文件，就没有这部分了。