应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于vasp计算单个原子的能量的交流
北京石油化工学院2026年研究生招生接收调剂公告
211/3返回列表123下一页
查看: 6620  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

从心而觅

新虫 (小有名气)

[交流] 关于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）得到，不知这样是否合理？
请各位来交流一下，多谢！

关于vasp计算单个原子的能量的交流
incar.png


关于vasp计算单个原子的能量的交流-1
kopints.png


关于vasp计算单个原子的能量的交流-2
poscar.png


关于vasp计算单个原子的能量的交流-3
能量统计.png
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理计算-tg soft study

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
一切有为法，如梦幻泡影！如露亦如电，应作如是观！
1楼 2016-08-15 20:08:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

从心而觅

新虫 (小有名气)
引用回帖:
15楼: Originally posted by xxxxxxxfc at 2016-08-18 14:07:14
啥叫消除简并，麻烦解释一下。...

摘自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.
赞一下(1人) 回复此楼
一切有为法，如梦幻泡影！如露亦如电，应作如是观！
16楼2016-08-18 15:16:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

obaica

禁言 (著名写手)
本帖内容被屏蔽
2楼2016-08-15 23:03:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
结合能的E(A) E(B)是孤立原子的能量也叫气态原子能量吧，建模型时候用真空大盒子里面放原子来得到。
你这个单个原子能量计算得到的是形成能。formation energy
赞一下 回复此楼
路漫漫其修远兮
3楼2016-08-16 15:31:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)
引用回帖:
3楼: Originally posted by hs0218 at 2016-08-16 15:31:43
结合能的E(A) E(B)是孤立原子的能量也叫气态原子能量吧，建模型时候用真空大盒子里面放原子来得到。
你这个单个原子能量计算得到的是形成能。formation energy

我这里A、B确实是孤立原子放在大盒子呀，怎么不是结合能呢。
形成能的A、B是稳定单质，而不是孤立原子。
赞一下 回复此楼
一切有为法，如梦幻泡影！如露亦如电，应作如是观！
4楼2016-08-16 15:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
算单质能量，求平均就行了
赞一下 回复此楼
5楼2016-08-16 17:12:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)
引用回帖:
5楼: Originally posted by xxxxxxxfc at 2016-08-16 17:12:18
算单质能量，求平均就行了

不对吧，两个原子结合到一块形成单质，不得放出能量？
赞一下 回复此楼
一切有为法，如梦幻泡影！如露亦如电，应作如是观！
6楼2016-08-16 18:04:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
6楼: Originally posted by 从心而觅 at 2016-08-16 18:04:45
不对吧，两个原子结合到一块形成单质，不得放出能量？...

我也不太懂，但确实是这么算的
赞一下 回复此楼
7楼2016-08-16 23:24:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)
引用回帖:
7楼: Originally posted by xxxxxxxfc at 2016-08-16 23:24:37
我也不太懂，但确实是这么算的...

看图
关于vasp计算单个原子的能量的交流-4
结合能.png
赞一下 回复此楼
一切有为法，如梦幻泡影！如露亦如电，应作如是观！
8楼2016-08-17 08:48:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxxxxxxfc

金虫 (正式写手)
引用回帖:
8楼: Originally posted by 从心而觅 at 2016-08-17 08:48:19
看图

结合能.png
...

你这个可靠不
回复此楼
9楼2016-08-17 08:55:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)
引用回帖:
9楼: Originally posted by xxxxxxxfc at 2016-08-17 08:55:26
你这个可靠不...

这个是侯玉柱老师编写的《vasp软件包使用入门指南》里面的，我在网上搜的别人也是这样说的。
赞一下 回复此楼
一切有为法，如梦幻泡影！如露亦如电，应作如是观！
10楼2016-08-17 16:09:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 从心而觅 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +11 一样YWY 2026-04-02 11/550 2026-04-03 23:46 by hemengdong
[考研] 338求调剂 +7 晟功? 2026-04-03 7/350 2026-04-03 16:46 by wxiongid
[考研] 复试调剂 +5 春日来信- 2026-04-03 5/250 2026-04-03 15:01 by buqi613
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[考研] 一志愿中国科学院大学265求调剂 +9 恬淡ye 2026-03-31 10/500 2026-04-03 11:10 by txp1986
[考研] 材料调剂 +8 懒羊羊轻置玉臀 2026-04-02 8/400 2026-04-02 22:03 by liu823948201
[考研] 求生物学调剂 +10 15172915737 2026-04-01 10/500 2026-04-02 18:53 by 哦哦嗯哈
[考研] 301求调剂 +4 细胞相关蛋白 2026-04-02 8/400 2026-04-02 16:30 by 271179835
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 0710生物学，325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
[考研] 一志愿北交大材料工程，总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 085601材料工程找调剂 +20 oatmealR 2026-03-29 21/1050 2026-04-01 21:00 by lijunpoly
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 一志愿南昌大学324求调剂 +7 hanamiko 2026-03-30 7/350 2026-04-01 13:22 by JourneyLucky
[考研] 【调剂】一志愿厦大生物与医药调剂 +3 Echo虾米 2026-03-31 3/150 2026-04-01 08:40 by JourneyLucky
[考研] 085404 22408 315分 +5 zhuangyan123 2026-03-31 6/300 2026-03-31 13:48 by limeifeng
[考研] 哈尔滨工业大学材料与化工专硕378求调剂 +3 塔比乌斯 2026-03-30 3/150 2026-03-30 22:55 by 无际的草原
[考研] 11408总分309，一志愿东南大学求调剂，不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 药学105500求调剂 +3 Ssun。。 2026-03-28 3/150 2026-03-28 11:24 by lxf170613
信息提示
关闭
请填处理意见
关闭