24小时热门版块排行榜    

CyRhmU.jpeg
查看: 16443  |  回复: 18
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lison4

金虫 (小有名气)

[求助] 怎么计算“缺陷形成能”公式中的化学势已有4人参与

我推出来的计算氧化锌的缺陷形成能的公式里面需要用到氧原子、锌原子、氧化锌分子的化学势,想请问一下各位这个应该怎么计算,或者是能直接从VASP或者material studio中读取吗?有没有能够参考的书籍或资料?
非常感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lison4: 金币+15, ★★★很有帮助 2016-12-02 09:29:06
引用回帖:
8楼: Originally posted by lison4 at 2016-11-21 19:39:41
最后做个结尾吧,以用VASP计算氧化锌为例子,推出缺陷形成能公式后需要用到Zn_bulk、O_bulk和ZnO_bulk的化学势,那么大多数情况下,若要算得O_bulk的化学势,通常直接建立氧气分子,然后计算该氧气分子的能量 读取O ...

楼主总结基本上是可以的,但是想做几个补充,不对的地方欢迎交流:
1.要注意化学势是一个约束条件,不能是任意取值的。在很多计算形成能的论文中都可以看到一个公式,就是Uo+Uzn=Uzno,bulk,其中Uo, Uzn表示化学势,而Uzno,bulk表示块体的能量。也就是说是二者化学势之和是定值,根据材料生长环境,两个元素的化学势是彼此约束的,以Uo为例,其取值上限就是氧气原子的能量/2,能量取值下限应该是zno块体能量减去反应形成晗。所以在讨论块体缺陷形成能时候必须说明是在什么条件下取的化学势。
2. O_bulk的化学势,除了是用楼主总结的“直接建立氧气分子,然后计算该氧气分子的能量 读取OUTCAT中free energy TOTEN除以2就得到了 O_bulk的化学势”,也有一种计算方法是算单个氧原子的,参见:https://cms.mpi.univie.ac.at/wiki/index.php/O_atom,但是我认为在楼主关注的情况下,是考虑楼主总结的方案。
计算顺利点吧
10楼2016-12-01 21:53:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lison4 的主题更新
信息提示
请填处理意见