24小时热门版块排行榜    

查看: 2638  |  回复: 7

egg19900805

新虫 (小有名气)

[求助] NEDOS的數量會改變DOS的peak已有4人参与

我是计算80颗原子的DOS
KPOINT=4*4*4
NEDOS由1000 2000 3000 10000
band gap出现的peak很不一致
低NEDOS时没有peak 但2000-3000时有一个peak
到10000时有2个peak 请问哪一个图形才是正确的

谢谢!!


擷取1.PNG


-1
擷取2.PNG


-2
擷取3.PNG


-3
擷取4.PNG
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Joannaouc

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
egg19900805: 金币+10, ★★★很有帮助 2014-05-18 17:10:42
我的理解应该是NEDOS越高越准确。但是2000时候已经能够说明问题,在费米能级处有peak出现。
2楼2014-05-16 14:15:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egg19900805

新虫 (小有名气)

引用回帖:
2楼: Originally posted by Joannaouc at 2014-05-16 14:15:18
我的理解应该是NEDOS越高越准确。但是2000时候已经能够说明问题,在费米能级处有peak出现。

谢谢你的回覆 想请问
2000的时候只有一个peak 在费米能阶以上
10000的时候出现两个peak 在费米能阶上跟下各一个

所以是10000的时候比较准吗?
3楼2014-05-16 14:42:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

苦行客

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
egg19900805: 金币+10, ★★★很有帮助 2014-05-18 17:10:33
应该是1000比较准确的,但是2000也满足分析需求了,具体的正确还是不正确你还得看能带,毕竟dos是高斯积分出来的,能带图才是最准确的!
行有不得,反求诸己~
4楼2014-05-17 11:40:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

redhaier

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
egg19900805: 金币+20, ★★★很有帮助 2014-05-18 17:10:18
franch: 金币+2, 谢谢回帖交流。。。 2014-05-19 22:56:34
如上面几楼所述,2000反映大概的问题,而且可以看到DOS随能量点密度的增加而收敛,说明楼主计算正确,具体是不是要表达出费米能级附近的DOS细节取决于楼主的计算要求,如果你计算这个体系本来就是要很详细的研究费米能级附近的态,是需要把所有峰的细节包括进来的,但是如果峰间距过于小(比如说0.025eV),就没太大必要,因为一般计算所设的电子温度的能量已经达到了这个水平。
5楼2014-05-17 15:57:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

egg19900805

新虫 (小有名气)

引用回帖:
5楼: Originally posted by redhaier at 2014-05-17 15:57:06
如上面几楼所述,2000反映大概的问题,而且可以看到DOS随能量点密度的增加而收敛,说明楼主计算正确,具体是不是要表达出费米能级附近的DOS细节取决于楼主的计算要求,如果你计算这个体系本来就是要很详细的研究费米 ...

想请问“因为一般计算所设的电子温度的能量已经达到了这个水平”
是什么意思看不太懂
有什么参考资料吗
谢谢你!
6楼2014-05-18 17:15:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pangrui1985

铜虫 (小有名气)

引用回帖:
6楼: Originally posted by egg19900805 at 2014-05-18 17:15:38
想请问“因为一般计算所设的电子温度的能量已经达到了这个水平”
是什么意思看不太懂
有什么参考资料吗
谢谢你!...

这取决于你用的展宽。这些都是用eginvalue文件的本征值做的,如果这个没改变的话不会带来新结论
7楼2014-05-18 17:44:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★
franch: 金币+2, 谢谢回帖交流。。。 2014-05-19 22:56:40
个人觉得应该是一个峰,10000那个可能是因为展宽方式造成的。
建议楼主用ISMEAR=-5试下,另外楼主看下eginvalue文件,那块应该有一条很窄的能带。
8楼2014-05-19 20:30:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 egg19900805 的主题更新
信息提示
请填处理意见