nwchem能量计算时正常输出，在优化时却报错电子数错误

Scarlett-ww

请教大家，计算一个铁体系的分子内单电子转移过程的能垒，使用gaussian、利用片段组合波函数，无论初始设置铁为三价或二价，都只能得到电子转移后的二价铁化合物。于是用nwchem计算重组能，取gaussian优化收敛后的二价铁结构，用nwchem计算单点能时可以正常结束，但是优化却始终报错，指出电子数错误，请问是什么原因呢？opt输入文件如下：
charge  1
GEOMETRY
O     -1.49286600   -1.56724100   -0.00002400
O     -2.02744300    1.05426800    0.00006800
C     -0.45362700   -0.72212200   -0.00001200
C      0.86321900   -1.10335300   -0.00004800
C      1.86328600   -0.09624400   -0.00002900
C      1.58134800    1.26876200    0.00002400
C      0.25886100    1.68304900    0.00005700
C     -0.78577700    0.68645000    0.00004000
Fe    -3.06383800   -0.63508700    0.00000100
H      1.13766400   -2.15147700   -0.00009200
H      2.39323200    1.98783300    0.00003900
O     -0.07172400    2.96009600    0.00010200
H     -1.04281800    3.04794100    0.00011200
C      3.31611000   -0.46724700   -0.00006700
O      3.49864500   -1.78928400   -0.00010400
H      4.45952400   -1.95064000   -0.00009900
O      4.19984600    0.35251600   -0.00001500
END
BASIS
* library 6-31G*
END
DFT
XC B3LYP
mult  6
convergence nolevelshifting
END
TASK DFT ENERGY
task scf optimize



报错信息如下：
------------------------------------------------------------------------
scf: no. of closed-shell electrons is not even!                   0
------------------------------------------------------------------------
------------------------------------------------------------------------
  current input line :
    30: task scf optimize
------------------------------------------------------------------------
------------------------------------------------------------------------
There is an error in the input file
------------------------------------------------------------------------
For more information see the NWChem manual at https://github.com/nwchemgit/nwchem/wiki


For further details see manual section:                                                                                                                                                                                       
0:scf: no. of closed-shell electrons is not even!:Received an Error in Communication
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 3 DUP FROM 0
with errorcode -1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------


全部输出如下文件：

Scarlett-ww

wzkchem5 发表于 2022-5-28 20:59
1. 对任何自旋体系都成立，任何自旋体系的任何碱金属、碱土金属、过渡金属都必须这么做。
2. 先解决溶剂 ...

非常感谢，我先溶剂化试一试。

wzkchem5

Scarlett-ww 发表于 2022-5-28 13:55
谢谢，我还有以下几个问题请教
1、请问需要溶剂分子把Fe的所有空配位点都饱和掉，是针对高自旋体系而言 ...

1. 对任何自旋体系都成立，任何自旋体系的任何碱金属、碱土金属、过渡金属都必须这么做。
2. 先解决溶剂化的问题再看
3. 对，但是不是比较结合能，而是比较结合自由能

Scarlett-ww

wzkchem5 发表于 2022-5-28 20:28
先不论报错的问题，你这个结构严重不合理，必须对Fe做显式溶剂化，也就是用溶剂分子把Fe的所有空配位点都饱 ...

谢谢，我还有以下几个问题请教
1、请问需要溶剂分子把Fe的所有空配位点都饱和掉，是针对高自旋体系而言吗？
2、为什么在task dft energy时没有对电子数报错，但是在task scf optimize时报错呢？官网上有相关问题的解释是可能由于对称性的定义，有更多的原子产生，https://www.nwchem-sw.org/index- ... ast_post/id277.html，但是在我的这个文件中，并没有定义对称性。
3、溶剂为醋酸钠缓冲溶液+水，请问是通过比较醋酸根和水与铁化合物的结合能，选择显性溶剂分子吗？

wzkchem5

先不论报错的问题，你这个结构严重不合理，必须对Fe做显式溶剂化，也就是用溶剂分子把Fe的所有空配位点都饱和掉。否则即使计算成功了，结果也是不能用的。