24小时热门版块排行榜    

查看: 778  |  回复: 0

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[交流] 【求助完毕】超胞DOS计算与K-mesh及ISMEAR选取

结构优化
INCAR:

ENCUT=480
EDIFF = 1E-4
EDIFFG = -0.01
ISMEAR =0
IBRION =2
ISPIN = 2
MAGMON = 2*3 64*0 30*0
ALGO =F
LREAL =  A

KPOINTS
A
0
M-K
1 1 1
0 0 0

静态计算
INCAR:
#IBRION = 2
ENCUT = 480
EDIFF = 1E-4
#EDIFFG = -0.01
ISMEAR = -5
SIGMA =0.05
ISPIN = 2
MAGMOM = 2*3 64*0 30*0

KPOINTS:
A
0
M-K
3 3 3
0 0 0

******************************************
问题:
(1)结构优化后收敛,静态计算电子自洽慢很多,最后得到的力也远大于0.01
(2)结构优化后自旋mag=6.00;静态计算mag=1.x
为什么会这样呢?结构优化后,我修改了ISMEAR, SIGMA和k-mesh,其它基本没变呢。

后来发现是这样的:
ISMEAR的选取与K点取样是有关系的。
对于半导体超胞结构,K点取样数可以适当减少,与之对应的ISMEAR比较适合取0。

超胞Gamma点结构优化后,加密k-mesh(如3×3×3),ISMEAR=0做静态自洽计算,得到DOSCAR。要得到高精度的DOS,有两个办法。一个是取较大的k-mesh,如 8×8×8,然后ISMEAR=-5;第二个是选取特别的k点,然后全部列出这些k点,ISMEAR = -5.对于超胞,第二种方法计算量可以接受,第一种则可能吃不消。

[ Last edited by cenwanglai on 2010-8-26 at 12:56 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见