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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » 刚接触QE。做了石墨烯的能带和态密度，不是是否正确，求各位鉴定下
202/2返回列表上一页12
查看: 2096  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

souledge

专家顾问 (著名写手)
引用回帖:
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

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +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的回帖

异域神迹

铁虫 (初入文坛)
引用回帖:
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

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +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

木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +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的回帖

异域神迹

铁虫 (初入文坛)
引用回帖:
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的回帖

异域神迹

铁虫 (初入文坛)
引用回帖:
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

木虫 (著名写手)
引用回帖:
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

铜虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
请问楼上的大虾  如果加石墨带加氢钝化后 要先优化晶格常数吗  如果要 怎么优化呢 ？？
赞一下 回复此楼
努力
19楼2013-04-19 14:37:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gemucai

木虫 (正式写手)
木有狄拉克点啊，赝势C.UPF是从哪里得到的？
赞一下 回复此楼
20楼2013-04-24 16:23:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 异域神迹 的主题更新
202/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭