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

登录注册

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

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

引用回帖:

10楼: Originally posted by newton3915 at 2013-04-17 16:30:06
就是你MS算band不是也要用GMKG这些布里渊区的高对称点吗，算nscf的话，KPOINTS就是要写出这些点的坐标，你还是先看看qe的example吧，里面有算能带的例子。...

MS里的SCF、NSCF是不需要自己动手，而是软件自动连接起来两个计算的。还有高对称点的选取，也都是软件有默认设定的，所以N多用MS的人根本就不知道这个世界上算能带和DOS是先做SCF再做NSCF……

赞一下

回复此楼

思想重于技巧，内涵重于表象

11楼2013-04-17 22:09:51

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

lfhuang

木虫 (著名写手)

应助: 39 (小学生)
金币: 3583.6
散金: 2141
红花: 37
帖子: 1432
在线: 681.4小时
虫号: 505824
注册: 2008-02-17
性别: GG
专业: 金属材料表面科学与工程

【答案】应助回帖

感谢参与，应助指数 +1

搂住可以参考以下链接的图2。
http://blog.sciencenet.cn/blog-345795-370129.html

赞一下

回复此楼

先进使役材料计算（https://aidme.nimte.ac.cn）

12楼2013-04-18 07:02:55

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

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

引用回帖:

11楼: Originally posted by souledge at 2013-04-17 22:09:51
MS里的SCF、NSCF是不需要自己动手，而是软件自动连接起来两个计算的。还有高对称点的选取，也都是软件有默认设定的，所以N多用MS的人根本就不知道这个世界上算能带和DOS是先做SCF再做NSCF……...

你好，我是做了SCF再做的NSCF，但是得到的能带图很奇怪，完全看不懂啊。求指导下，感谢。
以下是我的input，不知是哪里出问题了。

&control
calculation= ‘band‘
restart_mode= ‘from_scratch‘
pseudo_dir= ‘/home/wxzhang/soft/espresso-4.2.1/pseudo/
outdir= ‘/home/zyy/tmp/’
prefix= ‘graphene’
tprnfor = .true.
tstress = .true.
&system
ibrav= 4,
celldm(1)=4.01, celldm(3)= 10.00,
nat= 2,
ntyp= 1,
ecutwfc= 40.0,
ecutrho= 480,
occupations=‘smearing’ , smearing=‘gaussian’, degauss=0.05，
nbnd= 16
/
&electrons
diagonalization=‘david’
mixing_beta= 0.7
mixing_mode= ‘plain’
conv_thr= 1.0d-10
/
ATOMIC_SPECIES
C 12.0107          C.UPF
ATOMIC_POSITIONS（crystal)
C 2/3 1/3 1/2
C 1/3 2/3 1/2
K_POINTS（automatic)
60 60 1 0 0 0

try.bands.png

赞一下

回复此楼

13楼2013-04-18 15:47:15

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

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

感谢参与，应助指数 +1

引用回帖:

13楼: Originally posted by 异域神迹 at 2013-04-18 15:47:15
你好，我是做了SCF再做的NSCF，但是得到的能带图很奇怪，完全看不懂啊。求指导下，感谢。
以下是我的input，不知是哪里出问题了。

&control
calculation= ‘band‘
restart_mode= ‘from_ ...

当做calculations = 'band'的计算时，需要自己指定K点路径。对于石墨烯，可以用：

引用回帖:

K_POINTS
151
0.00000 0.00000 0.00000  1.00
  -0.00667 0.01155 0.00000  1.00
  -0.01333 0.02309 0.00000  1.00
  -0.02000 0.03464 0.00000  1.00
  -0.02667 0.04619 0.00000  1.00
  -0.03333 0.05774 0.00000  1.00
  -0.04000 0.06928 0.00000  1.00
  -0.04667 0.08083 0.00000  1.00
  -0.05333 0.09238 0.00000  1.00
  -0.06000 0.10392 0.00000  1.00
  -0.06667 0.11547 0.00000  1.00
  -0.07333 0.12702 0.00000  1.00
  -0.08000 0.13856 0.00000  1.00
  -0.08667 0.15011 0.00000  1.00
  -0.09333 0.16166 0.00000  1.00
  -0.10000 0.17321 0.00000  1.00
  -0.10667 0.18475 0.00000  1.00
  -0.11333 0.19630 0.00000  1.00
  -0.12000 0.20785 0.00000  1.00
  -0.12667 0.21939 0.00000  1.00
  -0.13333 0.23094 0.00000  1.00
  -0.14000 0.24249 0.00000  1.00
  -0.14667 0.25403 0.00000  1.00
  -0.15333 0.26558 0.00000  1.00
  -0.16000 0.27713 0.00000  1.00
  -0.16667 0.28868 0.00000  1.00
  -0.17333 0.30022 0.00000  1.00
  -0.18000 0.31177 0.00000  1.00
  -0.18667 0.32332 0.00000  1.00
  -0.19333 0.33486 0.00000  1.00
  -0.20000 0.34641 0.00000  1.00
  -0.20667 0.35796 0.00000  1.00
  -0.21333 0.36950 0.00000  1.00
  -0.22000 0.38105 0.00000  1.00
  -0.22667 0.39260 0.00000  1.00
  -0.23333 0.40415 0.00000  1.00
  -0.24000 0.41569 0.00000  1.00
  -0.24667 0.42724 0.00000  1.00
  -0.25333 0.43879 0.00000  1.00
  -0.26000 0.45033 0.00000  1.00
  -0.26667 0.46188 0.00000  1.00
  -0.27333 0.47343 0.00000  1.00
  -0.28000 0.48497 0.00000  1.00
  -0.28667 0.49652 0.00000  1.00
  -0.29333 0.50807 0.00000  1.00
  -0.30000 0.51962 0.00000  1.00
  -0.30667 0.53116 0.00000  1.00
  -0.31333 0.54271 0.00000  1.00
  -0.32000 0.55426 0.00000  1.00
  -0.32667 0.56580 0.00000  1.00
  -0.33333 0.57735 0.00000  1.00
  -0.32667 0.57735 0.00000  1.00
  -0.32000 0.57735 0.00000  1.00
  -0.31333 0.57735 0.00000  1.00
  -0.30667 0.57735 0.00000  1.00
  -0.30000 0.57735 0.00000  1.00
  -0.29333 0.57735 0.00000  1.00
  -0.28667 0.57735 0.00000  1.00
  -0.28000 0.57735 0.00000  1.00
  -0.27333 0.57735 0.00000  1.00
  -0.26667 0.57735 0.00000  1.00
  -0.26000 0.57735 0.00000  1.00
  -0.25333 0.57735 0.00000  1.00
  -0.24667 0.57735 0.00000  1.00
  -0.24000 0.57735 0.00000  1.00
  -0.23333 0.57735 0.00000  1.00
  -0.22667 0.57735 0.00000  1.00
  -0.22000 0.57735 0.00000  1.00
  -0.21333 0.57735 0.00000  1.00
  -0.20667 0.57735 0.00000  1.00
  -0.20000 0.57735 0.00000  1.00
  -0.19333 0.57735 0.00000  1.00
  -0.18667 0.57735 0.00000  1.00
  -0.18000 0.57735 0.00000  1.00
  -0.17333 0.57735 0.00000  1.00
  -0.16667 0.57735 0.00000  1.00
  -0.16000 0.57735 0.00000  1.00
  -0.15333 0.57735 0.00000  1.00
  -0.14667 0.57735 0.00000  1.00
  -0.14000 0.57735 0.00000  1.00
  -0.13333 0.57735 0.00000  1.00
  -0.12667 0.57735 0.00000  1.00
  -0.12000 0.57735 0.00000  1.00
  -0.11333 0.57735 0.00000  1.00
  -0.10667 0.57735 0.00000  1.00
  -0.10000 0.57735 0.00000  1.00
  -0.09333 0.57735 0.00000  1.00
  -0.08667 0.57735 0.00000  1.00
  -0.08000 0.57735 0.00000  1.00
  -0.07333 0.57735 0.00000  1.00
  -0.06667 0.57735 0.00000  1.00
  -0.06000 0.57735 0.00000  1.00
  -0.05333 0.57735 0.00000  1.00
  -0.04667 0.57735 0.00000  1.00
  -0.04000 0.57735 0.00000  1.00
  -0.03333 0.57735 0.00000  1.00
  -0.02667 0.57735 0.00000  1.00
  -0.02000 0.57735 0.00000  1.00
  -0.01333 0.57735 0.00000  1.00
  -0.00667 0.57735 0.00000  1.00
0.00000 0.57735 0.00000  1.00
0.00000 0.56580 0.00000  1.00
0.00000 0.55426 0.00000  1.00
0.00000 0.54271 0.00000  1.00
0.00000 0.53116 0.00000  1.00
0.00000 0.51962 0.00000  1.00
0.00000 0.50807 0.00000  1.00
0.00000 0.49652 0.00000  1.00
0.00000 0.48497 0.00000  1.00
0.00000 0.47343 0.00000  1.00
0.00000 0.46188 0.00000  1.00
0.00000 0.45033 0.00000  1.00
0.00000 0.43879 0.00000  1.00
0.00000 0.42724 0.00000  1.00
0.00000 0.41569 0.00000  1.00
0.00000 0.40414 0.00000  1.00
0.00000 0.39260 0.00000  1.00
0.00000 0.38105 0.00000  1.00
0.00000 0.36950 0.00000  1.00
0.00000 0.35796 0.00000  1.00
0.00000 0.34641 0.00000  1.00
0.00000 0.33486 0.00000  1.00
0.00000 0.32332 0.00000  1.00
0.00000 0.31177 0.00000  1.00
0.00000 0.30022 0.00000  1.00
0.00000 0.28868 0.00000  1.00
0.00000 0.27713 0.00000  1.00
0.00000 0.26558 0.00000  1.00
0.00000 0.25403 0.00000  1.00
0.00000 0.24249 0.00000  1.00
0.00000 0.23094 0.00000  1.00
0.00000 0.21939 0.00000  1.00
0.00000 0.20785 0.00000  1.00
0.00000 0.19630 0.00000  1.00
0.00000 0.18475 0.00000  1.00
0.00000 0.17320 0.00000  1.00
0.00000 0.16166 0.00000  1.00
0.00000 0.15011 0.00000  1.00
0.00000 0.13856 0.00000  1.00
0.00000 0.12702 0.00000  1.00
0.00000 0.11547 0.00000  1.00
0.00000 0.10392 0.00000  1.00
0.00000 0.09238 0.00000  1.00
0.00000 0.08083 0.00000  1.00
0.00000 0.06928 0.00000  1.00
0.00000 0.05773 0.00000  1.00
0.00000 0.04619 0.00000  1.00
0.00000 0.03464 0.00000  1.00
0.00000 0.02309 0.00000  1.00
0.00000 0.01155 0.00000  1.00
0.00000 0.00000 0.00000  1.00

如果用MS时稍微认真点，多仔细看看每个button，也不至于不知道需要自己设置特殊K点路径～

赞一下

回复此楼

思想重于技巧，内涵重于表象

14楼2013-04-18 20:55:29

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

liqizuiyang

木虫 (著名写手)

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

【答案】应助回帖

感谢参与，应助指数 +1

石墨烯是两套六角格子套起来的，同一套格子的基矢长度是2.46 埃左右，大约是4.6487263 Bohr。

QE用的是呈120度夹角的基矢，把celldm(1)改成4.6487263 ，碳原子坐标改成 (0, 0, 1/2)和(2/3,1/3,1/2)试试。

另外，从QE5.0开始，计算bands有专门的关键字 calculation='bands'，而不再是calculation = 'nscf‘。K路径需要在bands计算时就制定。算完bands后再用band.x进行处理。

这是我的一个QE算能带的bands计算的例子：

&CONTROL
  prefix          = '1H-MoS2'
  calculation    = 'bands'
  restart_mode    = 'from_scratch'
  wf_collect    = .TRUE.
  tstress       = .TRUE.
  tprnfor       = .TRUE.
  outdir          = './'
  wfcdir          = './'
  pseudo_dir    = './'
/
&SYSTEM
  ibrav          = 4
  celldm(1)       = 5.8435958
  celldm(3)       = 5.8209138
  ntyp          = 2
  nat             = 3
  nbnd          = 21
  ecutwfc       = 50.0
/
&ELECTRONS
  electron_maxstep = 250
  conv_thr       = 1.0d-10
  mixing_mode    = 'plain'
  mixing_beta    = 0.7
  mixing_ndim    = 8
  diagonalization  = 'david'
  diago_david_ndim = 4
  diago_full_acc = .TRUE.
/
ATOMIC_SPECIES
  S  32.065  S.pz-mt_fhi.UPF
  Mo 95.96 Mo.pz-mt_fhi.UPF
ATOMIC_POSITIONS {crystal}
  S 0.000000000  0.000000000  0.000847125
  S 0.000000000  0.000000000  0.172854079
  Mo  0.666666700  0.333333300  0.086854396
K_POINTS {crystal_b}
  4
  0.0    0.0    0.0  15
  0.333333  0.333333  0.0  15
  0.5    0.0    0.0  15
  0.0    0.0    0.0  15

赞一下

回复此楼

15楼2013-04-19 11:44:38

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

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

引用回帖:

14楼: Originally posted by souledge at 2013-04-18 20:55:29
当做calculations = 'band'的计算时，需要自己指定K点路径。对于石墨烯，可以用：

如果用MS时稍微认真点，多仔细看看每个button，也不至于不知道需要自己设置特殊K点路径～...

这位兄台，看了你的回复，有几个问题想请教下：
1.这些k点是从nscf计算中选取出来的么？如果是，那么这是怎么选出来的啊，如果不是的话那是怎么选取出来的喃？
2.我用你说得k点计算一次，以下是计算结果的能带图。其中有两个问题，首先左下角的两个能带应该是简并的吧。其次在K点和费米能级处应该有焦点吧？
以下是我找到的石墨烯的能带图。感谢

grapheneatlast.png

201010683418248.JPG

赞一下

回复此楼

16楼2013-04-19 11:55:35

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

异域神迹

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5
帖子: 29
在线: 10.4小时
虫号: 2386648
注册: 2013-03-29
性别: GG
专业: 半导体电子器件

引用回帖:

15楼: Originally posted by liqizuiyang at 2013-04-19 11:44:38
石墨烯是两套六角格子套起来的，同一套格子的基矢长度是2.46 埃左右，大约是4.6487263 Bohr。

QE用的是呈120度夹角的基矢，把celldm(1)改成4.6487263 ，碳原子坐标改成 (0, 0, 1/2)和(2/3,1/3,1/2)试试。

另外 ...

感谢，请问K_POINTS后面的{crystal_b}是什么意思啊？

赞一下

回复此楼

17楼2013-04-19 12:31:24

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

liqizuiyang

木虫 (著名写手)

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

引用回帖:

17楼: Originally posted by 异域神迹 at 2013-04-19 12:31:24
感谢，请问K_POINTS后面的{crystal_b}是什么意思啊？...

这是算能带时指定按特定方向取k点的关键字。

另外，石墨烯能带在K点价带和导带应该是简并的。你算的有一个带隙，这可能是优化晶格参数的时候没优化好。你指定的celldm(1)只有4.00 bohr，而石墨烯同一套格子的基矢长换算成bohr有4.6左右。按4.0计算相当于施加了应力，这个应力会打开一个数十毫电子伏的带隙。

tpiba_b : Used for band-structure plots. k-points are in units of 2 pi/a. nks points specify nks-1 lines in reciprocal space. Every couple of points identifies the initial and final point of a line. pw.x generates N intermediate points of the line where N is the weight of the first point.

赞一下

回复此楼

18楼2013-04-19 14:15:51

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

chenkunkun

铜虫 (小有名气)

应助: 14 (小学生)
金币: 269.6
帖子: 145
在线: 286.7小时
虫号: 2370184
注册: 2013-03-23
性别: GG
专业: 凝聚态物性 II ：电子结构

【答案】应助回帖

感谢参与，应助指数 +1

请问楼上的大虾如果加石墨带加氢钝化后要先优化晶格常数吗如果要怎么优化呢？？

赞一下

回复此楼

努力

19楼2013-04-19 14:37:44

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

gemucai

木虫 (正式写手)

应助: 16 (小学生)
金币: 1774.6
散金: 157
红花: 4
帖子: 693
在线: 201.7小时
虫号: 816478
注册: 2009-07-27
专业: 凝聚态物性 II ：电子结构

木有狄拉克点啊，赝势C.UPF是从哪里得到的？

赞一下

回复此楼

20楼2013-04-24 16:23:04

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

相关版块跳转我要订阅楼主异域神迹的主题更新

返回列表

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

[考研] 332求调剂 +15	小小孟... 2026-04-05	16/800	2026-04-05 23:29 by 来看流星雨10
[考研] 复试调剂 +13	呼呼？~+123456 2026-04-05	13/650	2026-04-05 22:07 by 醉翁wl
[考研] 一志愿青科085500，初试295分，公共课213分 +3	遇到的人愿望都� 2026-04-05	3/150	2026-04-05 18:45 by 蓝云思雨
[考研] 326求调剂 +3	顾若浮生 2026-04-05	3/150	2026-04-05 18:32 by 蓝云思雨
[考研] 材料专硕调剂 +14	CXN123456 2026-04-03	14/700	2026-04-05 17:18 by Hdyxbekcb
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +8	哇呼哼呼哼 2026-04-02	9/450	2026-04-05 17:14 by lbsjt
[考研] 346分的生物与医药08600求调剂 +5	常雨阳上岸 2026-04-05	6/300	2026-04-05 13:42 by imissbao
[考研] 0703调剂 +11	拾玖壹 2026-04-04	12/600	2026-04-05 10:29 by 果冻大王
[考研] 材料调剂 +9	革微桂 2026-04-04	9/450	2026-04-05 08:27 by 544594351
[考研] 男生，一志愿沪9生物学071000，初试308求调剂 +3	刘墨墨 2026-04-04	3/150	2026-04-05 08:26 by barlinike
[考研] 311分 22408 求调剂 +3	bing_bot 2026-04-03	3/150	2026-04-05 00:43 by chongya
[考研] 材料求调剂 +10	呢呢妮妮 2026-04-01	10/500	2026-04-04 23:12 by 无际的草原
[考研] 材料科学与工程339求调剂 +12	hyz0119 2026-03-31	13/650	2026-04-03 18:33 by ls刘帅
[考研] 274求调剂 +9	顺理成张 2026-04-03	10/500	2026-04-03 15:10 by 啊俊！
[考研] 260求调剂 +3	朱芷琳 2026-04-02	3/150	2026-04-03 08:44 by yulian1987
[考研] 326求调剂 +3	9ahye 2026-04-02	4/200	2026-04-03 08:43 by Jaylen.
[考研] 材料专业求调剂 +10	月月鸟木 2026-04-01	10/500	2026-04-02 12:57 by wxiongid
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +9	双马尾痞老板2 2026-04-01	9/450	2026-04-02 12:14 by oooqiao
[考研] 0703一志愿南师大334求调剂 +4	seven7yu 2026-03-30	4/200	2026-04-01 16:10 by oooqiao
[考研] 本科211总分289，08工学真心求调剂 +3	utopiaE 2026-03-30	3/150	2026-03-30 23:42 by ms629