24小时热门版块排行榜    

查看: 2104  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

异域神迹

铁虫 (初入文坛)

[求助] 刚接触QE。做了石墨烯的能带和态密度,不是是否正确,求各位鉴定下




以前用MS算过能带,但感觉QE算出来的跟MS差别不是一般大啊。。。求各位帮助。
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 20 个回答

异域神迹

铁虫 (初入文坛)

能看到能带和态密度图么,由于第一次发,不知成功否
2楼2013-04-17 14:57:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

newton3915

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by 异域神迹 at 2013-04-17 14:57:53
能看到能带和态密度图么,由于第一次发,不知成功否

三个图片都不到哦
3楼2013-04-17 14:58:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

异域神迹

铁虫 (初入文坛)

不好意思,现在来上传试试

graphene.bands.png



Graphene(1).png



Graphene(2).png

4楼2013-04-17 15:07:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见