原子化学势的计算

本文转载自：http://muchong.com/t-11236385-1-pid-2

 

可以参考一个以前的帖子http://muchong.com/bbs/viewthread.php?tid=10823604 怎么计算“缺陷形成能”公式中的化学势


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，但是我认为在楼主关注的情况下，是考虑楼主总结的方案。