24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 5386  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

贝芷苡

铜虫 (小有名气)

[求助] 掺杂化学势及形成能计算

本人在算氧化锌中掺杂形成能,但有些计算参数找不到参照,特发此贴,希望各位高手给指导一下。一些

氧化锌形成焓:Hf ZnO=(1/M)EtotZnO-(1/N)EtotZn-1/2EtotO2, where Hf ZnO is the enthalpy of formation。EtotZnO为含M个锌原子的超胞总能,EtotZn为N个锌原子的超胞总能,超胞均为六方密排结构,EtotO2为边长为10埃的立方体中放一氧分子的总能,此时得出的形成焓满足Uzn+Uo=Hf ZnO,本人计算各参数值为Hf ZnO=-3.443eV此值与文献中的值符合,(1/N)EtotZn=-1.262eV,1/2EtotO2=-4.388eV,那么锌空位在富氧Uo=0,Uzn=Hf ZnO=-3.443eV及缺氧条件Uzn=0,Uo=Hf ZnO=-3.443eV下的形成能(计算公式:Ef=Etot(q)-Etot(p)+Uzn,Etot(q)
:为含一个锌空位的超胞总能,Etot(p):为完美晶胞)为6.677eV,3.233eV.

比文献中值大一个电子伏左右。不知道是不是对公式的理解错误,还是计算有问题。
一下为我的个输入文件
INCAR
SYSTEM =ZnO
ISTART =0
ICHARG =2
ENCUT  =380
EDIFF  =1E-5
EDIFFG =-0.01
ISMEAR=0
SIGMA=0.1
NSW=100
NGX=62
NGY=62
NGZ=66
ISPIN=2
IBRION = 2
ISIF   =3
PREC=Accurte
LREAL= Auto

KPOINTS
SYSTEM =ZnO
0
M
3  3  2
0.0  0.0  0.0
静态计算时KPOINTS
SYSTEM =ZnO
0
M
5 5 3
0.0  0.0  0.0
另EtotZn为N个锌原子的超胞总能,超胞均为六方密排结构,EtotZn=-1.262(晶胞含2个锌原子),EtotZn=-0.944,哪个更准确?

回复此楼

» 收录本帖的淘帖专辑推荐

带电掺杂体系形成能 vasp计算资料 @相关

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

坚持每天进步一点点
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pengliusss

木虫 (小有名气)

引用回帖:
4楼: Originally posted by 贝芷苡 at 2013-04-24 16:47:21
UGa+UN=UGaN(bulk),中的UGaN(bulk)就是生成焓吧?
在PRB 63,075205这篇文章中有这样的介绍Under thermal equilibrium
UZn+UO=DHZnO , where DHZnO is the formation enthalpy
of ZnO,在Journal of Crystal Gro ...

晕了。。。。
5楼2013-04-25 09:34:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

yourne

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
贝芷苡: 金币+5, ★★★很有帮助, 超胞大小,和q值一样,那个形成能与费米能级、价带能以及 矫正位势有关时,是求的带电的缺陷吧。是我理解错了吗?我再看看这两篇文献。谢谢 2013-04-23 16:52:24
你的q值、超胞大小都和文献一样吗?尤其是超胞大小对结果影响较大。PRB, 61, 15019
公式也少了一大项,一般情况下Eformation都和费米能级、价带顶能量以及矫正位势有关。PRB, 54, 16676
“自由有许多困难,民主亦非完美,然而我们从未建造一堵墙把我们的人民关在里面,不准他们离开我们。” ——肯尼迪
2楼2013-04-23 09:21:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yourne

木虫 (小有名气)

不带电缺陷的话,第二篇文献给出了化学势的计算方法,没有用到生成焓,可以试一下。
“自由有许多困难,民主亦非完美,然而我们从未建造一堵墙把我们的人民关在里面,不准他们离开我们。” ——肯尼迪
3楼2013-04-24 11:13:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

贝芷苡

铜虫 (小有名气)

引用回帖:
3楼: Originally posted by yourne at 2013-04-24 11:13:29
不带电缺陷的话,第二篇文献给出了化学势的计算方法,没有用到生成焓,可以试一下。

UGa+UN=UGaN(bulk),中的UGaN(bulk)就是生成焓吧?
在PRB 63,075205这篇文章中有这样的介绍Under thermal equilibrium
UZn+UO=DHZnO , where DHZnO is the formation enthalpy
of ZnO,在Journal of Crystal Growth 287 (2006) 58–65中则是这样描述的The upper bounds defined above also lead to lower
bounds by using the thermodynamic stability condition for
ZnO,UZn+UO = EtotZnO, (3)where EtotZnO is the total energy of a two-atom unit of bulk ZnO.
坚持每天进步一点点
4楼2013-04-24 16:47:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 277工科求调剂 +10 1915668 2026-04-04 11/550 2026-04-06 23:53 by 心媛意北
[考研] 一志愿吉大化学327求调剂 +6 王王白石 2026-04-06 7/350 2026-04-06 21:46 by dongzh2009
[考研] 085410人工智能 初试316分 求调剂 +7 残星拂曙 2026-03-31 7/350 2026-04-06 10:09 by 蓝云思雨
[考研] 332求调剂 +17 小小孟... 2026-04-05 18/900 2026-04-06 09:51 by 蓝云思雨
[考研] 一志愿南昌大学,085600,344分求调剂 +9 调剂上岸玘 2026-04-05 10/500 2026-04-06 09:30 by dongzh2009
[考研] 086000生物与医药298调剂求助 +9 元元青青 2026-03-31 12/600 2026-04-05 21:03 by 学员8dgXkO
[考研] 一志愿上海海洋大学083200食品学硕,求调剂,接受其他专业083200 +4 what张 2026-04-04 5/250 2026-04-05 14:07 by chw1980_0
[考研] 调剂 +3 好好读书。 2026-04-02 3/150 2026-04-05 13:02 by arrow8852
[考研] 数一英一274机械调剂 +5 星陨流霞 2026-04-04 6/300 2026-04-05 11:38 by arrow8852
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 调剂 +11 JLLLLLLLLLL 2026-04-03 11/550 2026-04-04 22:21 by hemengdong
[考研] 085600调剂 +4 1amJJ 2026-04-02 4/200 2026-04-04 21:53 by hemengdong
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考研] 273求调剂 +20 李芷新1 2026-03-31 20/1000 2026-04-03 09:58 by linyelide
[考研] 326求调剂 +3 9ahye 2026-04-02 4/200 2026-04-03 08:43 by Jaylen.
[考研] 314求调剂 +11 1xiaojun23 2026-03-31 12/600 2026-04-02 12:31 by 1xiaojun23
[考研] 一志愿346上海大学生物学 +3 上海大学346调剂 2026-04-01 3/150 2026-04-02 08:36 by w虫虫123
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 00:21 by 百秒光年
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
信息提示
请填处理意见