| 查看: 6430 | 回复: 20 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
关于vasp计算单个原子的能量的交流 已有8人参与
|
|||
|
为了计算体系的结合能，按照公式，结合能：E=E(A)+E(B)-E(AB)（注：此处A,B为孤立原子总能，即两个孤立原子结合在一起所释放的能量），我首先计算了“AB”的能量，接着分别计算单个原子“A”“B”的能量，文件如下： 此外POTCAR选用标准的PBE。得到的结果如图。 有如下总结： 1：得到H的能量是-13.51ev，与-13.6ev比较接近,说明计算相对可靠。 2：得到C的能量是-148.49ev，小于Si的-103.87ev，这个我无法理解，按说C要大于Si。 3：为了验证2，我查资料发现有文献报道C的能量在-37ev左右，我重新检查了输入文件，并重算，还跟原来一样，还请大神帮忙解释一下。 4：晶体的总能量我也是通过（energy without entropy）-（atomic energy EATOM）得到，不知这样是否合理？ 请各位来交流一下，多谢！  incar.png  kopints.png  poscar.png  能量统计.png |
回复此楼» 猜你喜欢
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
|
摘自vasp手册 大概意思就是为了得到原子的基态能量，将原子放在一个非正方体的晶格内，比如说11Ax10Ax9A的晶格（这样做是为了消除简并）。 9.5 Determining the groundstate energ of atoms The POTCAR file contains information on the energy of the atom in the reference configuration (i.e. the configuration for which the PP was generated). Total energies calculated by vasp are with respect to this configuration. The reference calculation, however, did not allow for spin-polarisation or broken symmetry solutions, which usually lower the energy for gradient corrected of hybrid functionals. To include these effects properly, it is required to calculate the lowest energy magnetic groundstate using VASP. Unfortunately convergence to the symmetry broken spin polarized groundstate can be relatively slow in VASP. The following INCAR file worked reasonably well for most elements: ISYM = 0 # no symmetry ISPIN = 2 # allow for spin polarisation VOSKOWN = 1 # this is important, in particular for GGA # but not required for PBE potentials ISMEAR = 0 # Gaussian smearing, otherwise negative occupancies might come up SIGMA = 0.002 # tiny smearing width to safely break symmetry AMIX = 0.2 # mixing set manually BMIX = 0.0001 NELM = 100 # often many steps are required ICHARG = 1 Execute VASP twice to three times, consecutively with this input file until energies are converged. Furthermore, we recommend to use large slightly non-cubic cells, i.e. 12 °A ×13 °A ×14 °A. |
16楼2016-08-18 15:16:53
3楼2016-08-16 15:31:43
4楼2016-08-16 15:37:30
xxxxxxxfc
金虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 7641.8
- 红花: 7
- 帖子: 431
- 在线: 310.5小时
- 虫号: 1821098
- 注册: 2012-05-16
- 性别: GG
- 专业: 无机材料化学
5楼2016-08-16 17:12:18