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

登录注册

应《网络安全法》要求，自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

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有7人回复
上海工程技术大学激光智能制造课题组招收博士研究生已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有3人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有4人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有4人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有5人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有4人回复

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

用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 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3	QUjhNVAcOSff 2026-09-13	5/250	2026-09-14 21:25 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4	QUjhNVAcOSff 2026-09-13	5/250	2026-09-14 21:24 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4	两三岁ss 2026-09-14	4/200	2026-09-14 20:29 by 班主任的同桌
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	s3fFTmArrBt6 2026-09-14	3/150	2026-09-14 18:01 by Aj1rhIDL5ixY
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4	s3fFTmArrBt6 2026-09-14	4/200	2026-09-14 17:25 by Aj1rhIDL5ixY
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3	LwdutQ8HoqWP 2026-09-13	4/200	2026-09-14 16:13 by Aj1rhIDL5ixY
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4	LwdutQ8HoqWP 2026-09-13	4/200	2026-09-14 16:12 by Aj1rhIDL5ixY
[找工作] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4	3n8v2C8RimXI 2026-09-13	4/200	2026-09-14 16:00 by Aj1rhIDL5ixY
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4	3n8v2C8RimXI 2026-09-13	4/200	2026-09-14 15:48 by Aj1rhIDL5ixY
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	3n8v2C8RimXI 2026-09-13	4/200	2026-09-14 15:45 by Aj1rhIDL5ixY
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +5	QUjhNVAcOSff 2026-09-13	6/300	2026-09-14 15:21 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3	QUjhNVAcOSff 2026-09-13	3/150	2026-09-14 15:09 by Aj1rhIDL5ixY
[基金申请] 求教各位大神：2026教育部人文社科青年基金项目何时公示呀？ +5	云宴山人 2026-09-10	7/350	2026-09-14 12:11 by ecnu2013
[找工作] 浙江师范大学是怎么坑我的 +6	recruit123 2026-09-09	6/300	2026-09-14 12:00 by 萨仁其其格
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3	QUjhNVAcOSff 2026-09-13	4/200	2026-09-14 10:01 by 23jxep3nCNZb
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3	马方园 2026-09-10	4/200	2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3	Cahal1 2026-09-11	4/200	2026-09-14 08:50 by 北京莱茵编辑
[教师之家] 浙江师范大学是怎么坑我的 +10	recruit123 2026-09-09	15/750	2026-09-13 10:45 by seaskyy
[基金申请] 面上有专家说实验设备不是我们单位的 +6	zj_muchong 2026-09-08	8/400	2026-09-11 06:56 by suifengdao
[有机交流] 科研求助 10+4	我吃一只鱼 2026-09-07	4/200	2026-09-10 08:32 by 0940628

24小时热门版块排行榜

hn19870519

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

» 猜你喜欢

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

hn19870519

liqizuiyang

【答案】应助回帖

hn19870519

liqizuiyang

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