如何对某个原子的单点能进行确定

Opalite

请问大家我在利用计算H原子的单点能的时候出现了问题，我直接把H原子的原始坐标（未opt和freq优化，因为我记得看过说对某个单原子进行结构优化无意义），但是ORCA会进行报错：Error (ORCA_MDCI): Number of processes (6) in parallel calculation exceeds number of pairs (1)。请大家指教。我之前利用MS进行O原子的能量计算的时候（密度泛函理论），是将O2的能量算出来，然后除以2就是O原子的能量，单点能的计算也可以这样吗？请大家解惑

Uus/pMeC6H4-/キ

氢原子只有一个电子，不像其他多电子的体系那样需要多核并行跑高级别电子相关方法得到单点能，参见
http://bbs.keinsci.com/thread-10115-1-1.html
http://bbs.keinsci.com/thread-13859-1-1.html
http://bbs.keinsci.com/thread-38956-1-1.html
http://bbs.keinsci.com/thread-43832-1-1.html
以及ORCA论坛官方说法https://orcaforum.kofo.mpg.de/viewtopic.php?f=18&t=7913。

另外，“将O2的能量算出来，然后除以2就是O原子的能量”听起来不对劲，这相当于强行忽略了1/2 O2 -> O的能量变化，想想双原子分子的键解离曲线就知道不合理。

rugals

孤立氢原子的能级有精确解，随便翻一本量子力学教材都有。

非孤立的氢原子就没有单点能的概念了，至少得去做能量分解。

Opalite

Uus/pMeC6H4-/キ 发表于 2026-3-13 10:43
氢原子只有一个电子，不像其他多电子的体系那样需要多核并行跑高级别电子相关方法得到单点能，参见
http:/ ...

好的谢谢 根据您的指导我算出来了  另外关于那个O原子的能量其实我也是觉得很奇怪的 但是我们之前做金属催化材料都是这样去算 然后去算某个催化剂的氧空位形成能，可能和算法有关。

sobereva

将O2的能量算出来，然后除以2就是O原子的能量

没有前提情况下说某个原子能量，大家一律默认是指真空下的能量，O2能量除以2显然不是孤立存在的单个氧原子的能量。做某些能量数据计算（如计算形成能）时中间利用到的每原子能量是另一码事。