版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

hn19870519

木虫 (正式写手)

1ST强帖: 1
应助: 45 (小学生)
金币: 1225.2
散金: 1861
红花: 11
帖子: 379
在线: 203.1小时
虫号: 867945
注册: 2009-10-10
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 关于silicene能带计算的问题已有1人参与

最近刚开始学习pwscf，计算silicene的能带，碰到一个很奇怪的问题。我计算了高对称路径G-K-M和K-M，但是发现这两次计算，第三条和第四条能带的本征值(K-M)反转了。下面贴一下计算的输入文件，附件是两次计算的bands.dat文件。希望高手能够解惑。
&control
calculation = 'bands',
prefix = 'Si',
pseudo_dir='./'
outdir='./'
tstress=.true.
tprnfor=.true.
wf_collect=.true.
/
&system
ibrav =  4,
celldm(1) = 7.32477,celldm(3)=4.09569
nat =  2,
ntyp = 1,
ecutwfc =40
occupations= 'smearing'
smearing='gaussian'
degauss=0.020D0
nbnd=12
/
&electrons
mixing_beta = 0.7
conv_thr=1.D-8
/

ATOMIC_SPECIES
Si 28.086  Si.pbe-rrkj.UPF

ATOMIC_POSITIONS (crystal)
Si    0.666666667 0.333333333 0.486223239
Si    0.333333333 0.666666667  -0.486223239
K_POINTS crystal
40
  0.000000  0.000000  0.000000 1.0
-0.017526  0.035105  0.000000 1.0
-0.035053  0.070211  0.000000 1.0
-0.052579  0.105316  0.000000 1.0
-0.070105  0.140421  0.000000 1.0
-0.087632  0.175526  0.000000 1.0
-0.105158  0.210632  0.000000 1.0
-0.122684  0.245737  0.000000 1.0
-0.140211  0.280842  0.000000 1.0
-0.157737  0.315947  0.000000 1.0
-0.175263  0.351053  0.000000 1.0
-0.192789  0.386158  0.000000 1.0
-0.210316  0.421263  0.000000 1.0
-0.227842  0.456368  0.000000 1.0
-0.245368  0.491474  0.000000 1.0
-0.262895  0.526579  0.000000 1.0
-0.280421  0.561684  0.000000 1.0
-0.297947  0.596789  0.000000 1.0
-0.315474  0.631895  0.000000 1.0
-0.333000  0.667000  0.000000 1.0
-0.333000  0.667000  0.000000 1.0
-0.315474  0.658211  0.000000 1.0
-0.297947  0.649421  0.000000 1.0
-0.280421  0.640632  0.000000 1.0
-0.262895  0.631842  0.000000 1.0
-0.245368  0.623053  0.000000 1.0
-0.227842  0.614263  0.000000 1.0
-0.210316  0.605474  0.000000 1.0
-0.192789  0.596684  0.000000 1.0
-0.175263  0.587895  0.000000 1.0
-0.157737  0.579105  0.000000 1.0
-0.140211  0.570316  0.000000 1.0
-0.122684  0.561526  0.000000 1.0
-0.105158  0.552737  0.000000 1.0
-0.087632  0.543947  0.000000 1.0
-0.070105  0.535158  0.000000 1.0
-0.052579  0.526368  0.000000 1.0
-0.035053  0.517579  0.000000 1.0
-0.017526  0.508789  0.000000 1.0
  0.000000  0.500000  0.000000 1.0

回复此楼

» 猜你喜欢

2025冷门绝学什么时候出结果已经有5人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
孩子确诊有中度注意力缺陷已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用vasp计算能带时，出来的能带条数比较多的原因是？已经有11人回复
计算能带结构，高对称点的选择已经有14人回复
关于HSE06静态计算和HSE06能带计算的几个讨论已经有24人回复
关于态密度和能带计算的费米能级不统一的问题已经有6人回复
能带计算的问题已经有5人回复
手册中计算能带时k点设置的问题已经有12人回复
计算二维材料的能带，高对称点应该怎么选？已经有8人回复
求助我的有机半导体小分子计算出来的能带没有带隙已经有10人回复
Materials studio计算态密度、能带结构等等已经有12人回复
MS计算能带后怎样数据导出已经有3人回复
求助VASP能带计算的中断原因已经有6人回复
wien2k在终端下计算能带，改怎么做？已经有6人回复
[原创]一种基于能带计算结果的载流子有效质量计算方法已经有174人回复
HSE能带计算中ISYM参数的问题已经有7人回复
VASP+SOC后的能带计算已经有21人回复
关于使用MS中的CASTEP如何计算能带的问题已经有9人回复
HSE（GW类似）计算能带总结已经有51人回复
能带计算 vasp 已经有11人回复
【求助】VASP计算能带时KPOINTS设置出错已经有14人回复
【其他】不同版本软件计算Si的能带的结果比较，大家一起来做一下！已经有14人回复
【求助】如何通过能带计算电子和空穴的有效质量已经有7人回复
【求助】pwscf 能带计算已经有16人回复
【求助】PBE0 计算Si、Ge能带结构已经有9人回复

1楼2014-10-06 17:00:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

1ST强帖: 1
应助: 45 (小学生)
金币: 1225.2
散金: 1861
红花: 11
帖子: 379
在线: 203.1小时
虫号: 867945
注册: 2009-10-10
性别: GG
专业: 凝聚态物性 II ：电子结构

附件上传不了。发文本吧。截了两个数据。
G-K-M：
-0.333000  0.577928  0.000000
-14.538  -2.143 -14.525 -12.776  10.759  -2.124  10.766 8.251  12.752  11.215
  12.773  12.779
         -0.315474  0.577898  0.000000
-14.703  -2.376 -14.396 -12.730  10.661  -1.893  10.772 8.265  11.192  12.506
  12.560  12.506
K-M：
-0.333000  0.577928  0.000000
-14.538 -14.525 -12.776  -2.143  -2.124 8.251  10.759  10.766  11.215  12.752
  12.773  12.779
         -0.315474  0.577898  0.000000
-14.703 -14.396 -12.730  -2.376  -1.893 8.265  10.661  10.772  11.192  12.506
  12.560  12.506

赞一下

回复此楼

2楼2014-10-06 19:14:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2014-10-07 17:37:47
hn19870519: 金币+20, ★★★★★最佳答案, 太给力了！ 2014-10-08 10:34:33

算完scf -> bands之后，需要用bands.x提取能带信息。这个bands.x会根据当前k点对称性调整能带顺序。

可以在输入文件中加入no_onverlap = .true.禁用这一特性。

赞一下

回复此楼

3楼2014-10-07 10:13:03

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

1ST强帖: 1
应助: 45 (小学生)
金币: 1225.2
散金: 1861
红花: 11
帖子: 379
在线: 203.1小时
虫号: 867945
注册: 2009-10-10
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

3楼: Originally posted by liqizuiyang at 2014-10-07 10:13:03
算完scf -> bands之后，需要用bands.x提取能带信息。这个bands.x会根据当前k点对称性调整能带顺序。

可以在输入文件中加入no_onverlap = .true.禁用这一特性。

你好！再问一下，ph.x计算的时候，有没有对称性的开关。因为我计算silicene的声子谱，4×4×1的q点有虚频，但是3×3×1的就没有，所以我想测试一下是不是考虑了对称性之后K点数目变少的原因。

赞一下

回复此楼

4楼2014-10-08 21:10:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

liqizuiyang

木虫 (著名写手)

应助: 59 (初中生)
金币: 6859.5
散金: 261
红花: 94
帖子: 2093
在线: 708.9小时
虫号: 2366686
注册: 2013-03-21
性别: GG
专业: 凝聚态物性 II ：电子结构

引用回帖:

4楼: Originally posted by hn19870519 at 2014-10-08 21:10:37
你好！再问一下，ph.x计算的时候，有没有对称性的开关。因为我计算silicene的声子谱，4×4×1的q点有虚频，但是3×3×1的就没有，所以我想测试一下是不是考虑了对称性之后K点数目变少的原因。...

这个我不清楚，没用过ph模块。

赞一下

回复此楼

5楼2014-10-08 21:20:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 hn19870519 的主题更新

返回列表