24小时热门版块排行榜    

查看: 1984  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

identation

金虫 (著名写手)

[交流] 【求助】PWSCF计算介电极化时为何出现如下提示? 已有7人参与

==================================================
                            POLARIZATION CALCULATION
                         !!! NOT THOROUGHLY TESTED !!!
               --------------------------------------------------


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
     from c_phase : error #         1
     Polarization only for insulators and no empty bands
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     stopping ...
计算的是铁磁性体系,自洽计算正常,非自洽计算中出现这样的提示,请问是怎么回事呢?谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
identation(金币+1):谢谢参与
zzy870720z(金币+1):谢谢指导 2010-06-17 18:26:07
identation(金币+5): 2010-06-17 23:27:35
electron polarization的计算一般针对有energy band gap(带隙)的体系。先做好基态的计算,确定它是否是有带隙的体系。
4楼2010-06-17 13:03:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
zzy870720z(金币+2):谢谢侯老师指导 2010-06-17 20:21:55
hedaors(金币+1):thanks 2010-06-17 23:32:10
引用回帖:
Originally posted by identation at 2010-06-17 13:18:56:

老师您好,那nbnd要指定成价电子数的一半么?我计算的体系如果用nscf的费米能级,费米能级正好在导带,但费米能级和价带最高处有空带呀,也就是DOS为零的带,这事怎么回事呢?

对有带隙的体系,nbnd可以设置为总的价电子数目的一半(假设为偶数的话)或总电子数(总电子数为奇数)+1的一半。注意,电子是按Pauli原理来从低到高来占据能带的。

如果是金属性体系,则nbnd一定要比总电子数的一半大(比如大到10以上),这样确保Fermi能级能被计算地正确。金属性的体系,Fermi能级会穿过某些未占据的能带。

empty bands指的未占据的能带。
9楼2010-06-17 19:03:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
qasd(金币+2):谢谢 2010-06-19 09:29:34
不一定,不论半导体还是金属性的体系,一般要设置的nbnd比总价电子数的一半要大一些,这是因为对半导体性的体系进行计算时,比如要计算band gap,这样需计算一些未占据的能带。

完全未占据态的能带或未占据的本征值对体系的总能是没有贡献的。

[ Last edited by valenhou001 on 2010-6-18 at 10:19 ]
12楼2010-06-18 10:17:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★
identation(金币+15): 2010-06-18 15:51:54
qasd(金币+2):谢谢 2010-06-19 09:29:43
引用回帖:
Originally posted by identation at 2010-06-18 12:59:08:
那介电极化为什么无法计算呢?我改过截断能和K网格,都会出现介电极化计算的错误提示;然后由于是FM自旋,我把occupation指定成smearing了;

前面有几个回帖提到了,你的体系可能是金属性的,或者gap非常小。 这套处理介电极化的,主要用来处理半导体或绝缘体的。换句话说,它们不能处理金属性的体系。金属性的体系或金属在外加电场下,体系内部的极化很小,电子会集中于体系的表面。这些知识,可以看看固体物理以及电磁学方面的书。

在前面的回帖中也提到了,先搞清楚该材料的基态性质,看它究竟是半导体性的还是金属性的,再看是否合适做介电极化的计算。

[ Last edited by valenhou001 on 2010-6-18 at 13:23 ]
14楼2010-06-18 13:21:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2):谢谢侯老师指点 2010-08-21 22:10:36
"另外,计算得到的能带和DOS图对不上呀,能带图的带隙比DOS的宽很多呀,这是怎么回事呢"

这个较正常,态密度中,由于采用展宽的方法来确定态密度值,在价带边或导带边的能量处而出现非0的值,使得带隙看起来变小了。

带隙的值,最好从能带图中确定。
16楼2010-08-21 20:48:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
cenwanglai(金币+2):谢谢侯老师答复 2010-08-22 09:56:07
那你在计算出来的out文件搜一下"of electrons",找出体系总的电子数,然后算出总的价电子的能带数,然后设置根据这个值设置nbnd。
18楼2010-08-22 09:40:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 identation 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见