| 查看: 1426 | 回复: 1 | ||||
| 【悬赏金币】回答本帖问题，作者xiaomao_16将赠送您 5 个金币 | ||||
xiaomao_16银虫 (小有名气)
|
[求助]
带电缺陷形成能的计算已有1人参与
|
|||
|
如题，看了论坛里大家推荐的J. Appl. Phys. 95, 3851 (2004)这篇文章，如图1-2，但是还是有些问题还是没太搞明白，请教一下大家。（我算的是MgO的表面体系，就代入表面体系的能量了） 1. 我在另一篇PRL文章里看到了图3所示的计算形成能的方法，大体形式类似，但是后面的电荷项有点差别。我可以认为是他把图1里的Ev取为0了吗？（因为文章里作“形成能-E_F”图时横轴是E_F=0~6eV），可以直接这样取吗，还是说Ev应该取相对于真空能级的值。 2. μ(O)的问题，图3里给出了这篇文章算μ(O)的办法，还有看到直接用氧气分子能量的一半的。这个哪种算法相对合理呢？有前辈说氧的化学势不该看做定值，可是这样的话一个形成能里，μ(O)和E_F两个变量了，还怎么描述呢？ 3. 两种算法都有一个Δ项，前者是电势的修正，后者是振动能之类的，但文章里都没细说这个东西怎么算，这个影响大吗，可以忽略吗？看论坛里大家用ΔV的比较多，可以指点一下怎么得到这个修正项吗？看了文章还是不懂怎么用VASP实操。 大概就是以上这些问题，请各路大神指教！！！  形成能1.jpg  形成能2.jpg  形成能3.jpg |
回复此楼» 收录本帖的淘贴专辑推荐
 缺陷计算 |
» 猜你喜欢
七级岗，到手不到8000，带本科论文12教学分/生，额定教学分280
已经有4人回复
-
找工作千万不能来差学校
已经有4人回复
-
北京师范大学人工智能学院2024年第二批博士研究生招生报名
已经有3人回复
-
高分子聚合物有机材料方向的合作/招聘科研助理
已经有3人回复
-
国家重点实验室基金，通知去答辩。希望大不？
已经有11人回复
-
博士毕业找工作选择！
已经有5人回复
-
五一假期评审申请书
已经有12人回复
-
投稿RSC旗下杂志突然间看不到投稿状态，里面投稿的文章信息不见了，有遇到过的吗？
已经有5人回复
-
有没有想聊天的
已经有12人回复
-
博士后和科研助理好难找啊
已经有12人回复
» 本主题相关商家推荐: (我也要在这里推广)
【答案】应助回帖
感谢参与，应助指数 +1
|
1.想要完全搞明白这些问题，你还得花一些时间，我这里不方便讲得那么细 1. 我在另一篇PRL文章里看到了图3所示的计算形成能的方法，大体形式类似，但是后面的电荷项有点差别。我可以认为是他把图1里的Ev取为0了吗？（因为文章里作“形成能-E_F”图时横轴是E_F=0~6eV），可以直接这样取吗，还是说Ev应该取相对于真空能级的值。 Ev 取相对于真空能级的值 2. μ(O)的问题，图3里给出了这篇文章算μ(O)的办法，还有看到直接用氧气分子能量的一半的。这个哪种算法相对合理呢？有前辈说氧的化学势不该看做定值，可是这样的话一个形成能里，μ(O)和E_F两个变量了，还怎么描述呢？ μ(O)应该取MgO中O的化学式，这个值μ(O)我们不太确定。但是它有两个极限值。极大值是氧气分子能量的一半，也就是通常说的O-rich condiction；极小值是当μ(Mg)取极大值时对应的μ(O)，也就是通常所说的O-poor condiction。你画图可以E_F为变量，而同时画两幅图，O-rich condiction和O-poor condiction，真正的formation energy只可能在这两者之间。 3. 两种算法都有一个Δ项，前者是电势的修正，后者是振动能之类的，但文章里都没细说这个东西怎么算，这个影响大吗，可以忽略吗？看论坛里大家用ΔV的比较多，可以指点一下怎么得到这个修正项吗？ 这一项有时候还挺大的（～2 eV），一般不能忽略。 看了文章还是不懂怎么用VASP实操。 如果只想得到计算值，pycdt可以帮你搞定一切！PyCDT: A Python toolkit for modeling point defects in semiconductors and insulators |
2楼2020-03-31 09:14:33
10(金币+5)