24小时热门版块排行榜    

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

volcano93

新虫 (初入文坛)

[交流] 算graphene能带,很乱看不懂 已有5人参与

头一次用pwscf算能带,算的是graphene的,得到能带图很乱看不懂

输入文件是这样的:
&control
    calculation='bands'
    pseudo_dir = '/home/zc/espresso-4.3.1/pseudo/',
    outdir='/home/zc/tmp'
    prefix='g'
/
&system

    ibrav= 0, celldm(1) = 2.460031762, nat=  2, ntyp= 1,
    nbnd = 16,
    ecutwfc = 60,
/
&electrons
/
ATOMIC_SPECIES
C 12.0107 C.pz-vbc.UPF
ATOMIC_POSITIONS {crystal}
C   0.666666667   0.333333333   0.000000000  
C   0.333333333   0.666666667   0.000000000  
CELL_PARAMETERS { cubic }
1            0.000000000  0.000000000  
-0.50000000  0.866025404  0.000000000  
0.000000000  0.000000000 10.000000000  
K_POINTS
61
         0.0  0.0            0.0  1
         0.0250    0.0250         0    1.0000
    0.0500    0.0500         0    1.0000
    0.0750    0.0750         0    1.0000
    0.1000    0.1000         0    1.0000
    0.1250    0.1250         0    1.0000
    0.1500    0.1500         0    1.0000
    0.1750    0.1750         0    1.0000
    0.2000    0.2000         0    1.0000
    0.2250    0.2250         0    1.0000
    0.2500    0.2500         0    1.0000
    0.2750    0.2750         0    1.0000
    0.3000    0.3000         0    1.0000
    0.3250    0.3250         0    1.0000
    0.3500    0.3500         0    1.0000
    0.3750    0.3750         0    1.0000
    0.4000    0.4000         0    1.0000
    0.4250    0.4250         0    1.0000
    0.4500    0.4500         0    1.0000
    0.4750    0.4750         0    1.0000
    0.5000    0.5000         0    1.0000
    0.4917    0.5083         0    1.0000
    0.4833    0.5167         0    1.0000
    0.4750    0.5250         0    1.0000
    0.4667    0.5333         0    1.0000
    0.4583    0.5417         0    1.0000
    0.4500    0.5500         0    1.0000
    0.4417    0.5583         0    1.0000
    0.4333    0.5667         0    1.0000
    0.4250    0.5750         0    1.0000
    0.4167    0.5833         0    1.0000
    0.4083    0.5917         0    1.0000
    0.4000    0.6000         0    1.0000
    0.3917    0.6083         0    1.0000
    0.3833    0.6167         0    1.0000
    0.3750    0.6250         0    1.0000
    0.3667    0.6333         0    1.0000
    0.3583    0.6417         0    1.0000
    0.3500    0.6500         0    1.0000
    0.3417    0.6583         0    1.0000
    0.3333    0.6667         0    1.0000
    0.3167    0.6333         0    1.0000
    0.3000    0.6000         0    1.0000
    0.2833    0.5667         0    1.0000
    0.2667    0.5333         0    1.0000
    0.2500    0.5000         0    1.0000
    0.2333    0.4667         0    1.0000
    0.2167    0.4333         0    1.0000
    0.2000    0.4000         0    1.0000
    0.1833    0.3667         0    1.0000
    0.1667    0.3333         0    1.0000
    0.1500    0.3000         0    1.0000
    0.1333    0.2667         0    1.0000
    0.1167    0.2333         0    1.0000
    0.1000    0.2000         0    1.0000
    0.0833    0.1667         0    1.0000
    0.0667    0.1333         0    1.0000
    0.0500    0.1000         0    1.0000
    0.0333    0.0667         0    1.0000
    0.0167    0.0333         0    1.0000
         0         0         0    1.0000

建模是平行四边形的,高对称点是从(0,0,0)--(1/2,1/2,0)--(1/3,2/3,0)--(0,0,0)
画图的那部分是这样:
gbands.dat
-6.0 10
gbands.xmgr
gbands.ps
4.6745
1.0 4.6745
于是,怎么破……[ Last edited by volcano93 on 2012-7-19 at 21:32 ]
回复此楼

» 猜你喜欢

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

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

gemucai

木虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
5楼: Originally posted by tangosnow at 2012-07-20 00:43:47
这个特殊k点,是怎么生成的?...

自己写一个小程序,应该不难。里面的20和100自己定义吧,只能作为一个参考。
dx=(Xt-X0)/n
dy=(Yt-Y0)/n
dz=(Zt-Z0)/n
do i=1,n,1
write(20,100)X0+(i-1)*dx, Y0+(i-1)*dy, Z0+(i-1)*dz
end do
7楼2012-08-09 10:56:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

souledge

专家顾问 (著名写手)

★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
sunyang1988: 金币+1, 谢谢交流 2012-07-19 23:10:43
石墨烯是六方,这些高对称点没个靠谱的。同时,对称性跟你建的模型是个虾米样子是俩概念,CELL_PARAMETER中明明给的就是六方结构的可好?
我个人觉得,没啥要破的,硬要破,就破破LZ重视计算本身而忽视基础理论的惯性吧……
思想重于技巧,内涵重于表象
2楼2012-07-19 18:41:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

volcano93

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by souledge at 2012-07-19 18:41:09
石墨烯是六方,这些高对称点没个靠谱的。同时,对称性跟你建的模型是个虾米样子是俩概念,CELL_PARAMETER中明明给的就是六方结构的可好?
我个人觉得,没啥要破的,硬要破,就破破LZ重视计算本身而忽视基础理论的惯 ...

呃,意思是,那个基矢打开来是平行四边形的模样……
表骂我……呜呜,我也想学好基础理论啊,可是大三刚进实验室,固体物理课都没上过一节……我用功去补了,但还是迷茫的样子……
求教看什么书靠谱?然后这个能带究竟要怎么看?
3楼2012-07-19 21:39:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

souledge

专家顾问 (著名写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
sunyang1988: 金币+2, 谢谢帮助 2012-07-19 23:11:22
引用回帖:
3楼: Originally posted by volcano93 at 2012-07-19 21:39:53
呃,意思是,那个基矢打开来是平行四边形的模样……
表骂我……呜呜,我也想学好基础理论啊,可是大三刚进实验室,固体物理课都没上过一节……我用功去补了,但还是迷茫的样子……
求教看什么书靠谱?然后这个能 ...

呵呵,又是本科提前进实验室的啊~除了简历能光辉点适合保研外,实在培养不出什么科学品格来,因为基础还差了太多了……
不多说了,现在都流行这个~总之想弄明白建模问题就要多看固体物理,晶体学,尤其是后者。推荐Walter Borchardt-Ott的Crystallography,Springer出版社出版的英文翻译版(原始版本是德文版)。
至于能带嘛~用下面这一系列的K点坐标吧:
   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
思想重于技巧,内涵重于表象
4楼2012-07-19 22:49:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +11 yufeiwaner 2026-08-09 12/600 2026-08-12 01:27 by 南京中西医
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 29/1450 2026-08-12 00:24 by zhanghaozhu
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +7 布布和一二 2026-08-10 9/450 2026-08-11 22:12 by 医学老男孩
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +8 Tide man 2026-08-10 9/450 2026-08-11 20:39 by beefly
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[硕博家园] 读博的好处 +3 lnee 2026-08-11 3/150 2026-08-11 18:10 by 希望我好好的
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 13/650 2026-08-11 09:49 by 家与远方
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见