| 查看: 17531 | 回复: 18 | ||||||||
[求助]
怎么计算“缺陷形成能”公式中的化学势 已有4人参与
|
||||||||
|
我推出来的计算氧化锌的缺陷形成能的公式里面需要用到氧原子、锌原子、氧化锌分子的化学势，想请问一下各位这个应该怎么计算，或者是能直接从VASP或者material studio中读取吗？有没有能够参考的书籍或资料？ 非常感谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | 带电掺杂体系形成能 | vasp计算资料 | 差分电荷密度 |
| 缺陷计算 | 计算模拟 |
» 猜你喜欢
可靠消息，8月21号(明天)放榜
已经有10人回复
-
只有每年这种时候来逛逛小木虫
已经有20人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有125人回复
-
今天基金会出结果吗？20260819
已经有15人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
|
最后做个结尾吧，以用VASP计算氧化锌为例子，推出缺陷形成能公式后需要用到Zn_bulk、O_bulk和ZnO_bulk的化学势，那么大多数情况下，若要算得O_bulk的化学势，通常直接建立氧气分子，然后计算该氧气分子的能量 读取OUTCAT中free energy TOTEN除以2就得到了 O_bulk的化学势，同理计算Zn_bulk 和ZnO_bulk也是建立结构后算的free energy TOTEN然后除以相应原子/分子数量即可； 在进行该计算的INCAR中 K点1X1X1即可，encut与你计算的supercell保持一致即可。 我目前就是这样算的，如若不对，欢迎指出~非常感谢！ |
» 本帖已获得的红花（最新10朵）
z8418880778楼2016-11-21 19:39:41
9楼2016-11-23 20:32:59
2楼2016-11-17 20:02:38
冯亚娟
金虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2694.5
- 红花: 7
- 帖子: 328
- 在线: 121.4小时
- 虫号: 2721395
- 注册: 2013-10-13
- 性别: MM
- 专业: 原子和分子物理
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
lison4(蛋蛋小童鞋代发): 金币+5, 谢谢参与，小木虫鼓励积极回答问题者 2016-11-18 10:32:18
感谢参与，应助指数 +1
lison4(蛋蛋小童鞋代发): 金币+5, 谢谢参与，小木虫鼓励积极回答问题者 2016-11-18 10:32:18
| 个人认为化学势是与你计算体系的化学环境相关联的。其中经常采用的包括：（1）单质所对应的每个原子的能量。比如Fe是fcc结构，可以算单个原胞的能量，然后除以每个原胞包含的原子数目4；（2）考虑元素化合物形式（感觉里面是包含成键特性的），比如如果你的目的是生成Fe3O4，但是体系中会有Fe2O3的杂质相产生，那你Fe3O4的化学势取值必然要小于Fe2O3的化学势，以此来表现你对与杂质相的考虑（这个是举例说明，至于Fe3O4与Fe2O3哪个化学势大我并不清楚，可千万不要产生误导）。 |
5楼2016-11-18 09:15:01
7楼2016-11-18 15:55:17
|
本帖内容被屏蔽 |
3楼2016-11-17 20:06:15
4楼2016-11-17 23:05:51
wuli8
荣誉版主 (知名作家)
…………
- 1ST强帖: 2
- 应助: 35 (小学生)
- 贵宾: 12.924
- 金币: 20192.4
- 散金: 15888
- 红花: 88
- 沙发: 4
- 帖子: 7840
- 在线: 1114.6小时
- 虫号: 465889
- 注册: 2007-11-23
- 专业: 物理学I
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
lison4(蛋蛋小童鞋代发): 金币+5, 谢谢交流 2016-11-21 19:40:58
感谢参与，应助指数 +1
lison4(蛋蛋小童鞋代发): 金币+5, 谢谢交流 2016-11-21 19:40:58
6楼2016-11-18 11:14:11
zhangfrank
金虫 (正式写手)
- 应助: 178 (高中生)
- 金币: 660.2
- 散金: 226
- 红花: 53
- 帖子: 891
- 在线: 344.5小时
- 虫号: 2679423
- 注册: 2013-09-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
lison4: 金币+15, ★★★很有帮助 2016-12-02 09:29:06
lison4: 金币+15, ★★★很有帮助 2016-12-02 09:29:06
|
楼主总结基本上是可以的，但是想做几个补充，不对的地方欢迎交流： 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
