24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1298  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料科学与工程320求调剂,080500 +4 黄瓜味薯片 2026-04-06 4/200 2026-04-06 19:25 by 无际的草原
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +18 披星河 2026-04-03 18/900 2026-04-06 13:55 by BruceLiu320
[考研] 一志愿南航,数一英一学硕317求调剂!! +6 Acaciad 2026-04-04 6/300 2026-04-06 12:13 by 考研学校招点人
[考研] 272分求调剂 +4 wangyile2233 2026-04-02 4/200 2026-04-05 22:21 by 286640313
[考研] 280求调剂 +4 李rien 2026-04-04 4/200 2026-04-05 18:44 by imissbao
[考研] 工科08专硕机械275求调剂 +3 AaAa7420 2026-04-02 3/150 2026-04-05 13:26 by jp9609
[考研] 求调剂 +3 电气小神童 2026-04-04 3/150 2026-04-05 10:17 by barlinike
[考研] 11408,335分,本科211,求调剂,可转专业 +5 鳄梨大鳄鱼 2026-04-03 5/250 2026-04-04 22:49 by chongya
[考研] 材料专硕306英一数二 +8 z1z2z3879 2026-03-31 8/400 2026-04-04 22:08 by hemengdong
[考研] 277求调剂 +4 12A3 2026-04-02 5/250 2026-04-04 20:28 by 蓝云思雨
[考研] 求调剂 +3 ffyyu 2026-04-02 3/150 2026-04-04 19:03 by 蓝云思雨
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 322求调剂 +6 FZAC123 2026-04-03 6/300 2026-04-03 22:23 by 科研小专家
[考研] 一志愿重庆大学085404,总分314分,求调剂 +4 zf83hn 2026-04-03 4/200 2026-04-03 21:25 by 啵啵啵0119
[考研] 考研调剂 +3 Draa 2026-04-03 3/150 2026-04-03 17:37 by hgwz7468
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 26考研调剂 +4 Wnz.20030617 2026-04-01 5/250 2026-04-02 16:11 by 1939136013狗壮
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 301求调剂 +8 axibli 2026-04-01 8/400 2026-04-01 09:51 by 我的船我的海
[考研] 一志愿中海洋320化学工程与技术学硕求调剂 +8 披星河 2026-03-30 8/400 2026-03-31 08:53 by lbsjt
信息提示
请填处理意见