求助怎么算一个金属原子位于真空层中的能量

Nora

我目前在计算单个原子能量方面遇到了问题，我看过一个视频说是单原子在真空层中的化学势与计算出一个晶胞的化学势除以原子数的能量相同，但是我现在计算一个Si原子的能量用两种方法差别比较大，这是什么原因呢？第一个方案的能量大约为0.08eV，第二个方案得到的结果是43.3828eV。而且我发现一个奇怪的现象，第一个方案的原子用vasp算出来是有磁性的，第二个没有磁性。这是为什么，另附INCAR的参数设置图。

wzkchem5

这两个能量显然不一样，因为前者没有考虑cohesive energy而后者考虑了。如果你要算自由原子，显然只能用第一种方法。
不要机械性地理解视频说的话，这个不是说做视频的人不懂计算，而是有时候做视频的人可能口误了，或者说的话有前提条件（比如只对算气体的情形成立，对于凝聚相不成立），如果机械性地理解，就容易把口误理解成真理，或者忽略前提条件。必须自己从原理上验证视频所说的是否有明显的不合理之处，所谓尽信书不如无书。如果不懂原理，宁可暂时不要做计算，把基础理论打扎实了，有判断视频说的话对不对的能力，再去看视频。

整不会了

wzkchem5 发表于 2022-7-16 16:53
这两个能量显然不一样，因为前者没有考虑cohesive energy而后者考虑了。如果你要算自由原子，显然只能用第 ...

我也感觉方案一比较合理，但不知道格子边长取多少才对

Mac0817

您好，我想请教一下，方案一是如何建模呢。晶格还是选择该元素单质的晶格类型（空间群）吗，建好模型后真空层的添加有具体要求没有呀

Mac0817

wzkchem5 发表于 2022-7-16 16:53
这两个能量显然不一样，因为前者没有考虑cohesive energy而后者考虑了。如果你要算自由原子，显然只能用第 ...

您好，我想请教一下，方案一是如何建模呢。晶格还是选择该元素单质的晶格类型（空间群）吗，建好模型后真空层的添加有具体要求没有呀

wzkchem5

Mac0817 发表于 2025-5-21 10:26
您好，我想请教一下，方案一是如何建模呢。晶格还是选择该元素单质的晶格类型（空间群）吗，建好模型后真 ...

对空间群没有要求
真空层要求在三个方向上都足够大，所谓足够大就是把真空层进一步加大也几乎不会改变算出来的能量。如果拿不准，就参考文献的设置