计算化学公社

标题: ORCA计算单个金属离子的单点能时报错：电子数为0 [打印本页]

作者
Author: Kitkim 时间: 2023-11-2 18:01
标题: ORCA计算单个金属离子的单点能时报错：电子数为0
各位大佬，我在用ORCA计算一个络合物体系的结合能的时候，涉及到计算金属离子的单点能，计算时报错为number of electrons (NEL) is zero
以下为inp文件

! PBE D3 RIJCOSX noautostart printbasis nopop PAL8  LARGEPRINT
%maxcore  20480
%pal nprocs 8 end

%basis
NewECP Sn "SDD" end
end

* xyz 4 1
Sn       0    0    0
*

我看ORCA论坛上涉及该问题的解法是直接设置电荷为0或者为负值，但是之前在gaussian上就能够计算这个离子的能量，想问一下各位大佬在orca上这类计算除了改变电荷外有没有更优解

作者
Author: zjxitcc 时间: 2023-11-2 19:31
本帖最后由 zjxitcc 于 2023-11-2 19:32 编辑

Gaussian确实不会报错，但算出来电子能量是0，这个不用算也知道，46个芯电子全被赝化了。用ORCA不用算也是0，这个跟程序没有关系。

如果你真的算+4价的Sn，我个人建议使用def2TZVP赝势基组，它只有28个芯电子被赝化。当然，这样的话各个物种里的Sn元素都得用def2TZVP。

作者
Author: Kitkim 时间: 2023-11-2 20:06

zjxitcc 发表于 2023-11-2 19:31
Gaussian确实不会报错，但算出来电子能量是0，这个不用算也知道，46个芯电子全被赝化了。用ORCA不用算也是0 ...

谢谢大佬，我刚刚确认了一下，之前用gaussian算之所以会有数值，是因为我加入了溶剂化，如果去掉溶剂化的话就会直接报错了。我再考虑一下用def2TZVP，感谢

作者
Author: sobereva 时间: 2023-11-3 03:51
SDD根本不适合主族，连重要的极化函数都没有（除非是Stuttgart官网上带极化函数的版本）
该用什么看下文
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

作者
Author: Kitkim 时间: 2023-11-3 10:21

sobereva 发表于 2023-11-3 03:51
SDD根本不适合主族，连重要的极化函数都没有（除非是Stuttgart官网上带极化函数的版本）
该用什么看下文
...

谢谢卢老师指点，确实忽略这个元素主族的问题，我换def2-TZVPP再试试，感谢

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