24小时热门版块排行榜    

查看: 1278  |  回复: 16
当前主题已经存档。

enola

捐助贵宾 (正式写手)

[交流] 【求助】pwscf作能带图时特殊k点的选择 已有4人参与

刚刚学pwscf,我计算硅材料,本来有例子,它是以fcc的结构输入的,但我换成sc的结构输入,然后每个原胞有8个原子,坐标如下:
Si 0.000 0.000 0.000
Si 0.000 0.500 0.500
Si 0.500 0.000 0.500
Si 0.500 0.500 0.000
Si 0.750 0.250 0.750
Si 0.250 0.250 0.250
Si 0.250 0.750 0.750
Si 0.750 0.750 0.750
其到空间还是SC结构,特殊k点是:
G 0.0 0.0 0.0
M 0.5 0.5 0.0
R 0.5 0.5 0.5
X 0.0 0.5 0.0
但是我按照G-M-R-X-G把能带画出来,和用fcc结构画出来的是不一样的,我知道是k点方向选择的关系,所以我想问大侠们,这个sc结构的Si怎么选取k点可以和fc就够的能带一样呢?若我构建的是一个硅的超胞,是以sc为结构的,那又该怎么k点来看能带呢?
总的来说当我构建一个新的结构(超胞,layer),我怎么才能获得我所需要的能带图呢,就是k点怎么选择。
非常感谢
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

timeflies..
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

goldenfisher

金虫 (著名写手)

★ ★
enola(金币+1):5 2010-04-12 20:57
qasd(金币+2):谢谢 2010-04-12 22:11
enola(金币+3): 2010-04-13 08:36
用xcrysden软件选择特殊K点,不同的晶报在里面会有不同的布里渊区,在其中选择你的特殊K点,然后设定总的K点个数即可。这个时候设定K_POINT cyrstal。
4楼2010-04-12 20:45:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

enola

捐助贵宾 (正式写手)

我的scf输入文件:
&control
    calculation='scf'
    restart_mode='from_scratch',
    prefix='si'
    pseudo_dir = '/mirror/espresso
    outdir='/mirror/tmp',
    tstress=.true.
    tprnfor=.true.
    disk_io='high'
/
&system
    ibrav= 1, celldm(1)= 10.21, na
    ecutwfc =36.0
/
&electrons
    conv_thr =  1.0d-8
    mixing_beta = 0.7
/

ATOMIC_SPECIES
Si  28.086  Si.vbc.UPF
ATOMIC_POSITIONS {crystal}
Si 0.000 0.000 0.000
Si 0.000 0.500 0.500
Si 0.500 0.000 0.500
Si 0.500 0.500 0.000
Si 0.750 0.250 0.750
Si 0.250 0.250 0.250
Si 0.250 0.750 0.750
Si 0.750 0.750 0.750
K_POINTS
。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
bands输入文件:
&control
    calculation='scf'
    restart_mode='from_scratch',
    prefix='si'
    pseudo_dir = '/mirror/espresso-4.0.5/pseudo/',
    outdir='/home/enola/tmp/'
/
&system
    ibrav= 1, celldm(1)= 10.2, nat= 2, ntyp= 1,
    ecutwfc =18.0
/
&electrons
    conv_thr =  1.0d-8
    mixing_beta = 0.7
/
ATOMIC_SPECIES
Si  28.086  Si.vbc.UPF
ATOMIC_POSITIONS
Si 0.00 0.00 0.00
Si 0.25 0.25 0.25
K_POINTS
下面的k点怎么选才能获得能带呢
timeflies..
2楼2010-04-12 18:23:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xirainbow

木虫 (正式写手)

★ ★
enola(金币+1):2 2010-04-12 20:57
qasd(金币+2):谢谢 2010-04-12 22:11
enola(金币+2): 2010-04-13 08:36
找到你希望得到的k点的直角坐标
不是以2pi/a为单位
是以i,j,k为单位
3楼2010-04-12 19:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

★ ★ ★
xiaohunhun(金币+3):很好的实例 2010-04-13 19:30
引用回帖:
Originally posted by goldenfisher at 2010-04-12 20:45:46:
用xcrysden软件选择特殊K点,不同的晶报在里面会有不同的布里渊区,在其中选择你的特殊K点,然后设定总的K点个数即可。这个时候设定K_POINT cyrstal。

你好阿,刚刚按照你的方法计算拉,但当我设定为计算8条能级的时候,似乎计算出来的都是费米能级以下的能级,直到计算32条能带,本征值才从(-5-》12),所以我想问一下,假如我们预定能带是计算8条,那从那个本征值开始到那里结束呢?是由软件自己决定的吗,这8条能带它是怎么分布的呢,等间距的还是其他的。。。。。要是我计算的能带刚刚不是导带底和价带顶,那不是画出来的能带部队拉。
这是我画出来的能带图,k点选择是X-Gamma-M-R-Gamma,倒格子是SC,就是和例子的不一样:



[ Last edited by enola on 2010-4-13 at 09:31 ]
timeflies..
5楼2010-04-13 09:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

引用回帖:
Originally posted by xirainbow at 2010-04-12 19:46:39:
找到你希望得到的k点的直角坐标
不是以2pi/a为单位
是以i,j,k为单位

谢谢阿,不过为什么计算出来还是和例子不一样呢?
timeflies..
6楼2010-04-13 09:39:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

watermall

金虫 (小有名气)


enola(金币+2): 2010-04-13 14:09
xiaohunhun(金币+1):xiexie 2010-04-13 19:29
因为你用的结构是sc
而例子里的结构是fcc
所以不一样
假设你用的Si的赝势里只考虑了Si的4个价电子
sc结构里有8个原子,每个原子有4个电子,所以最少需要32个能带fcc里面有2个原子,所以最少需要8个能带
如果你算能带,可以不设置nbnd,程序自动设置能带数目
7楼2010-04-13 09:48:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

引用回帖:
Originally posted by watermall at 2010-04-13 09:48:33:
因为你用的结构是sc
而例子里的结构是fcc
所以不一样
假设你用的Si的赝势里只考虑了Si的4个价电子
sc结构里有8个原子,每个原子有4个电子,所以最少需要32个能带fcc里面有2个原子,所以最少需要8个能带
如果 ...

如果结构选的不一样,出来的能带结构也不一样啊?
我用MS算了一下,如果是SC结构,其能带和我用pwscf的计算出来的差不多,但换成FC,就和另外一个差不多啦。那对于同一种材料,其能带结构怎么区分,或者由于选取模型的不同,得出了不同的能带?怎么解释呢?
关于能带数目,谢谢你拉,感觉你的说法是对的。
谢谢你的帮组
timeflies..
8楼2010-04-13 14:12:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
xiaohunhun(金币+2):xiexie 2010-04-13 21:42
算能带的时候一般都应该将能带数目选择多一些,可以看INPUT_PW。这是因为你需要计算导带电子的能带,而这部分在基态时是没有电子占据的。
另外,你说的几条能带这个说法很不准确,nbnd并不是画出来的能带的数目,而是电子占据的能带数目。你的超原胞越大,画出来的能带数目越多,这是由于Bloch定理决定的。当你选择的是最小单胞时,画出来的能带数目最少。nbnd决定能带能量范围,也就是电子占据的能量范围。
9楼2010-04-13 20:59:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

goldenfisher

金虫 (著名写手)

★ ★
xiaohunhun(金币+2):xiexie 2010-04-13 21:42
enola(金币+1): 2010-04-14 09:31
你给的图中价带的能带貌似和例子一样,但是导带就完全不一样,我怀疑你的计算有问题。你可以将所有的特殊k点都包含在bands计算中,然后将图画出来对照一下看看。
10楼2010-04-13 21:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 enola 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见