24小时热门版块排行榜    

查看: 488  |  回复: 2

hisap-nano

新虫 (正式写手)

[求助] 怎样用ABINIT计算获得金属的功函数

最近在用ABINIT在算金属材料的力学性能时发现金属的功函数也是一个重要的物理量,然后就构建了一个slab算了表面能,但是不知道怎么得到相应表面的功函数,看了论坛里的关于功函数的计算的帖子,没有一个是用ABINIT的。按照功函数等于真空能减掉Fermi energy. 现在好像是ABINIT的输出文件里可以读到Fermi energy,但是真空能不知道从哪里获得,恳请高手指点!
下面的图是从输出的势函数中提出来的,但是图中的最高能级是负的,而Fermi energy是0.06Ha. 是不是我哪里设错了?!!




[ Last edited by hisap-nano on 2011-12-13 at 13:03 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
uuv2010(金币+1): 多谢指导! 2011-12-13 18:31:54
hisap-nano(金币+5): 2011-12-14 06:35:26
计算功函数的方法都差不多,所以没必要强调ABINIT中的处理方法
真空能级就是真空处的电势,能从平均静电势中得到。在远离slab的真空处,电势趋于常数,这个值就是真空能级。
你这个图是计算得到的电势图吗?如果是的话,对于周期的格子,最左侧和最右侧的电势应该相等啊,最好检查下你的数据处理过程有没有问题。
2楼2011-12-13 17:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

stractor

金虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
贺仪(金币+1): 多谢指教! 2011-12-14 09:42:01
功函数就是真空能级与费米能级的差值,费米能级输出文件有,而真空能级可以处理静电势的分布得到。
计算垂直于slab方向的平均静电势,在slab区域,静电势是震荡的,而在真空区域,平均静电势会是一个常数,这个常数就是真空能级。
3楼2011-12-14 01:19:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hisap-nano 的主题更新
信息提示
请填处理意见