计算化学公社

标题: nwchem能量计算时正常输出，在优化时却报错电子数错误 [打印本页]

作者
Author: Scarlett-ww 时间: 2022-5-28 19:58
标题: nwchem能量计算时正常输出，在优化时却报错电子数错误
本帖最后由 Scarlett-ww 于 2022-5-28 19:59 编辑

请教大家，计算一个铁体系的分子内单电子转移过程的能垒，使用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.
--------------------------------------------------------------------------

全部输出如下文件：

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

作者
Author: Scarlett-ww 时间: 2022-5-28 20:55

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、溶剂为醋酸钠缓冲溶液+水，请问是通过比较醋酸根和水与铁化合物的结合能，选择显性溶剂分子吗？

作者
Author: wzkchem5 时间: 2022-5-28 20:59

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

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

作者
Author: Scarlett-ww 时间: 2022-5-28 21:42

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)