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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于形成能的计算
101/1返回列表
查看: 4528  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖

huichli

木虫 (正式写手)
嘿嘿，不好意思计算公式输错了，应该是：
Ef=E(doped_Fe2O3)-E(pure_Fe2O3)-mE(dopant)+nE(Fe)
再次感谢“gkf高”的纠正！
赞一下(1人) 回复此楼
6楼2013-06-04 10:36:17
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

171713294

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
huichli: 金币+20, ★★★★★最佳答案, 再次感谢！ 2013-06-04 14:09:57
sunyang1988: 金币+1, 谢谢帮助 2013-06-04 16:42:42
AMIX=0.2
BMIX=0.0001这两个采用默认值
ICHARG设置为2
赞一下(1人) 回复此楼
坚持，不放弃
7楼2013-06-04 11:02:46
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duantian0529

金虫 (小有名气)
想问下，你查的这个文献是哪篇？可以告诉我捏文献名字吗？我想看看，谢啦！
赞一下 回复此楼
加油哦！
8楼2015-12-15 14:54:32
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

duantian0529

金虫 (小有名气)
m和n分别为掺杂原子和被取代铁原子的数目，而掺杂原子和取代原子数难道不一样吗？实在是不懂，求指教
赞一下 回复此楼
加油哦！
9楼2015-12-15 15:37:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LYJJFL

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by huichli at 2013-06-04 10:36:17
嘿嘿，不好意思计算公式输错了，应该是：
Ef=E(doped_Fe2O3)-E(pure_Fe2O3)-mE(dopant)+nE(Fe)
再次感谢“gkf高”的纠正！

我怎么觉得他把你给忽悠了呢？
赞一下 回复此楼
10楼2018-12-04 21:00:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huichli 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭