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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于形成能的计算
51/1返回列表
查看: 4794  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huichli

木虫 (正式写手)

[求助] 关于形成能的计算

小虫是vasp新手，对这方面都是一知半解，周围一个懂计算的人都没有，一路跌跌撞撞摸索过来，都不知道自己有没有入门。最近遇到一个关于形成能计算的问题，我研究的体系是掺杂其它金属元素的氧化铁，想把形成能计算一下，我查到形成能的计算公式是：
Ef=E(doped_Fe2O3)-E(pure_Fe2O3)+mE(dopant)-nE(Fe)，
其中E(doped_Fe2O3)和E(pure_Fe2O3)分别为掺杂Fe2O3和纯Fe2O3的总能量，E(dopant)和E(Fe)分别为掺杂原子和被取代的Fe原子的能量（也有化学势这一说），m和n分别为掺杂原子和被取代铁原子的数目（不知道是不是这样？）。我的问题是这里面掺杂原子和Fe原子的能量是怎么算的，我在相关的帖子和论文里查到这种单原子能量是通过计算纯体相金属（pure bulk metal）的总能而后平均到每个原子算得，不知道真的是不是这样，请大家指点一下，计算过程都要注意哪些啊？
还有侯博士的那个VASP入门指南里的关于单原子的计算跟这个有关系吗？
求指点、、求帮助啊
回复此楼

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

VASP 第一性原理计算-tg CASTEP

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-06-03 23:17:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huichli

木虫 (正式写手)
引用回帖:
2楼: Originally posted by 171713294 at 2013-06-04 08:16:44
掺杂原子和Fe原子的能量就是算他们的单点能，也就是着这两个原子分别放在立方超饱内，用单k点计算，你可以用15埃的超饱，k点选用1*1*1

麻烦大神再帮忙看一下我的INCAR和KPOINTS文件吧，不知道是不是这样写，谢谢啦！
INCAR文件
SYSTEM=Fe
ENCUT=270
ISYM=0
ISPIN=2
ISMEAR=0
SIGMA=0.1
AMIX=0.2
BMIX=0.0001
NELM=200
ICHARG=1
LMAXMIX=4
KPOINTS文件
KPOINTS文件
Automatic
0
Gama
1 1 1
0 0 0
赞一下 回复此楼
3楼2013-06-04 09:16:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

171713294

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2013-06-04 08:39:36
huichli: 金币+10, ★★★很有帮助, 非常感谢！ 2013-06-04 09:12:32
huichli: 金币+20, ★★★很有帮助 2013-06-04 09:16:53
掺杂原子和Fe原子的能量就是算他们的单点能，也就是着这两个原子分别放在立方超饱内，用单k点计算，你可以用15埃的超饱，k点选用1*1*1
赞一下 回复此楼
坚持，不放弃
2楼2013-06-04 08:16:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gkf高

木虫 (著名写手)

世界和平

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
huichli: 金币+10, ★★★很有帮助, 谢谢啊，大牛眼光果真犀利，昨天完善一不小心输错了！ 2013-06-04 10:26:02
sunyang1988: 金币+1, 谢谢帮主 2013-06-04 16:42:25
如果要算形成能的话 ，应该按相同的计算方法和设置计算所有能量才有可比较性。
你的形成能计算公式确定正确？要是取代Fe的话，应该为
Ef=E(doped_Fe2O3)-E(pure_Fe2O3)-mE(dopant)+nE(Fe)吧。
赞一下 回复此楼
天空不留下鸟的痕迹，但是我已飞过。
4楼2013-06-04 09:53:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huichli

木虫 (正式写手)
引用回帖:
4楼: Originally posted by gkf高 at 2013-06-04 09:53:03
如果要算形成能的话 ，应该按相同的计算方法和设置计算所有能量才有可比较性。
你的形成能计算公式确定正确？要是取代Fe的话，应该为
Ef=E(doped_Fe2O3)-E(pure_Fe2O3)-mE(dopant)+nE(Fe)吧。

您能否再详细的说明一下是那些地方，您的意思是指INCAR和KPOINTS文件都要和Fe2O3和掺杂Fe2O3体系的文件相同吗？小虫新手，很多地方不懂，还望能说的通俗易懂些，呵呵。
赞一下 回复此楼
5楼2013-06-04 10:28:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 调剂310 +3 温柔的晚安 2026-03-25 4/200 2026-03-25 23:16 by peike
[考研] 08工学调剂 +13 用户573181 2026-03-20 18/900 2026-03-25 22:00 by zbssa
[考研] 机械学硕总分317求调剂！！！！ +4 Acaciad 2026-03-25 4/200 2026-03-25 19:59 by hanserlol
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 278求调剂 +4 我可以上岸的对� 2026-03-19 4/200 2026-03-25 11:01 by userper
[考研] 286求调剂 +11 Faune 2026-03-21 11/550 2026-03-25 10:11 by 雾散后相遇lc
[考研] 07化学280分求调剂 +7 722865 2026-03-23 7/350 2026-03-25 09:29 by aa331100
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +6 @taotao 2026-03-19 6/300 2026-03-25 08:37 by 木托莫露露
[考研] 食品专硕 一志愿双一流 328 +3 xiaom99 2026-03-21 4/200 2026-03-24 21:20 by lailaisimei
[考研] 0703化学调剂，求导师收 +7 天天好运来上岸� 2026-03-24 7/350 2026-03-24 20:26 by peike
[考研] 070300化学求调剂 +9 苑豆豆 2026-03-20 9/450 2026-03-24 17:15 by licg0208
[考研] 求调剂一志愿武汉理工大学材料工程（085601） +5 WW.' 2026-03-23 7/350 2026-03-24 14:50 by sprinining
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 352求调剂 +3 大米饭！ 2026-03-22 3/150 2026-03-22 23:28 by king123！
[考研] 求调剂一志愿海大，0703化学学硕304分，有大创项目，四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 324求调剂 +6 lucky呀呀呀鸭 2026-03-20 6/300 2026-03-22 16:01 by ColorlessPI
[考研] 296求调剂 +4 www_q 2026-03-20 4/200 2026-03-21 17:26 by 学员8dgXkO
[考研] 一志愿 南京航空航天大学大学 ，080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 0856调剂，是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
信息提示
关闭
请填处理意见
关闭