| 查看: 7476 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理计算-tg | soft study |
» 猜你喜欢
极限拉扯，如何抉择
已经有9人回复
-
入职后第一次管课题组的专利缴费，对账才发现官费比别人多掏了一倍
已经有6人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有8人回复
-
课题组会帮你通过考核么？
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教关于氧原子化学势和用vasp计算氧分子的化学势的问题
已经有12人回复
- 关于用vasp计算磁矩的问题
已经有10人回复
- 关于Gibbs 程序的小问题请教（单位换算，E-V曲线）
已经有3人回复
- 关于VASP的并行
已经有10人回复
- 请教诸位前辈，关于总能计算，吸附能计算，修正能的一些疑问求教……
已经有5人回复
- 新手求助pwscf的进行磁性运算和LDA+U运算
已经有35人回复
- 关于vasp计算表面能中的Ebulk
已经有7人回复
- 【求助】计算单个Fe原子，无法收敛，请教
已经有19人回复
|
摘自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
|
本帖内容被屏蔽 |
2楼2016-08-15 23:03:41
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
6楼2016-08-16 18:04:45
xxxxxxxfc
金虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 7641.8
- 红花: 7
- 帖子: 431
- 在线: 310.5小时
- 虫号: 1821098
- 注册: 2012-05-16
- 性别: GG
- 专业: 无机材料化学
7楼2016-08-16 23:24:37
8楼2016-08-17 08:48:19
xxxxxxxfc
金虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 7641.8
- 红花: 7
- 帖子: 431
- 在线: 310.5小时
- 虫号: 1821098
- 注册: 2012-05-16
- 性别: GG
- 专业: 无机材料化学
9楼2016-08-17 08:55:26
10楼2016-08-17 16:09:11
