应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 【求助】pwscf作能带图时特殊k点的选择
51/1返回列表
查看: 1282  |  回复: 16
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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..
1楼 2010-04-12 18:21:06
已阅   回复此楼   关注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为单位
赞一下(1人) 回复此楼
3楼2010-04-12 19:46:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答

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的回帖

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。
赞一下(2人) 回复此楼
4楼2010-04-12 20:45:46
已阅   回复此楼   关注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 ]
赞一下(1人) 回复此楼
timeflies..
5楼2010-04-13 09:29:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 17 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭