24小时热门版块排行榜    

查看: 1112  |  回复: 6
当前主题已经存档。

xwyan

新虫 (正式写手)

[交流] 【求助】pwscf的能带计算中,occupations=' '还用设置吗?

pwscf的能带计算中,occupations=' '还用设置吗?
我看到scf的计算中是需要设置occupation=' '
在例子中,能带、费米面的‘nscf' 计算时,都没有设置occupation=' '

请问,nscf的计算中还用设置occupation=' '?
       在nscf的计算中设置了occupation=' '?, 结果还正确吗? 有什么影响?
        occupations选项的意义我不是太清楚?
十分感谢你的评论和指导
回复此楼

» 猜你喜欢

事了拂衣去,千里不留行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★ ★ ★ ★ ★
zxzj05(金币+2):谢谢,欢迎常来交流讨论 1-23 14:33
xwyan(金币+3): 1-23 22:11
xwyan(金币+1): 1-26 08:46
occupation是为了加速金属的scf收敛而设置的一个参数
并且最后要进行收敛测试(有关收敛测试内容可以参看这个http://blog.sina.com.cn/s/blog_5f15ead20100d1ki.html

nscf不是scf,我觉得应该不需要这个参数了
前提必须是你已经对这个参数做过收敛测试
2楼2010-01-23 14:04:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lfhuang

木虫 (著名写手)

★ ★ ★
zxzj05(金币+1):谢谢,欢迎常来交流讨论 1-23 14:33
xwyan(金币+1): 1-23 22:12
xwyan(金币+1): 1-26 08:45
很多都是有默认选项的。请楼主计算之前应该好好地读一下程序所给的一些资料,尤其是:(1)user-guide; (2) Quantum-Espresso/Doc/INPUT_*。
先进使役材料计算(https://aidme.nimte.ac.cn)
3楼2010-01-23 14:10:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwyan

新虫 (正式写手)

我对两种情况下做了计算,用diff 命令对比了输出的文件,相应k格点的处的能量值是完全一致的(degauss=0.002).
感谢上面两位的评论和指点。
正如xirainbow所说:
"occupation是为了加速金属的scf收敛而设置的一个参数"
nscf 的计算中是没有必要设置这个参数的。
即使设置也不会起什么作用,(因为是nscf)
所以设置occupations 和不设置occupations的计算结果是相同的。

[ Last edited by xwyan on 2010-1-27 at 08:45 ]
事了拂衣去,千里不留行
4楼2010-01-26 08:45:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Sunnyin

木虫 (初入文坛)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
qasd(金币+3):xiexie 2010-02-09 14:43
一般在计算半导体材料的能带时, 可以不设置smearing, 但是在计算金属材料时, 由于费米能级横跨某个能带, 为了准确描述费米面附件的电子态及其占据的情况, 需要用smearing或者tetrahedra技术. 特别是计算density of state,应该用更为精确的tetrahedra. smearing本身也分simple gaussian, MP, Vanderbilt-Marzari 等scheme, 这些都是gaussian smearing的变化, 主要是收敛的速度不一样.

NSCF 也是要设置smearing的, 比如你用smearing=0.02Ry 计算了SCF 的电子密度, 在此基础上,你可以用不同的smearing来算nscf (特别是在电声耦合计算里, 不同的smearing用于电声耦合的nscf计算,有不同的意义). smearing 本身可以理解为电子的温度,虽然不是绝对的意义上的温度.

在简单的能带结构计算种,smearing的意义体现的不是很多.
5楼2010-02-06 03:43:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwyan

新虫 (正式写手)

十分感谢sunnyin的回复!
事了拂衣去,千里不留行
6楼2010-02-06 10:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★ ★
qasd(金币+3):xiexie 2010-02-09 14:43
在pwscf中,计算能带应该首先设置calculation='bands'.然后设置好你要计算的k点方向。nscf是计算态密度的时候用的。在计算DOS时,occupation设置为tetrahedra为好,在计算能带时,可以与scf计算设置的一致。个人建议还是加上这一选项,因为只有加上这一选项,你才能得到费米能级的值。实际上,对于半导体而言,这一选项并不影响最终结果。对于金属而言,如果你不加这一选项,则你的结果是错误的甚至不能计算收敛。它的意义并不止于一个参数,而是存在实际的物理意义的。
7楼2010-02-08 20:44:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xwyan 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见