24小时热门版块排行榜    

查看: 2590  |  回复: 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的回帖

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的回帖

tangosnow

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
4楼: Originally posted by souledge at 2012-07-19 22:49:55
呵呵,又是本科提前进实验室的啊~除了简历能光辉点适合保研外,实在培养不出什么科学品格来,因为基础还差了太多了……
不多说了,现在都流行这个~总之想弄明白建模问题就要多看固体物理,晶体学,尤其是后者。推 ...

这个特殊k点,是怎么生成的?
5楼2012-07-20 00:43:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

frank_zhan

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
sunyang1988: 呵呵,好学没错但基础也挺重要的。不建议做计算的童鞋把计算工具完全当做黑盒子来用 2012-07-20 18:23:11
引用回帖:
4楼: Originally posted by souledge at 2012-07-19 22:49:55
呵呵,又是本科提前进实验室的啊~除了简历能光辉点适合保研外,实在培养不出什么科学品格来,因为基础还差了太多了……
不多说了,现在都流行这个~总之想弄明白建模问题就要多看固体物理,晶体学,尤其是后者。推 ...

你要帮助人家就好好帮人家,总是装一个德高望重的老长者一样教训别人干什么?人家是虚心求教,虚心交流。我支持楼主的好学精神,想学习啥时候都不晚。
frank_zhan,zy2zhan@gmail.com
6楼2012-07-20 10:31:52
已阅   回复此楼   关注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的回帖

hanyanli0475

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你好LZ,这里有点不太明白的,graphene是六角立方,所以ibrav=4吧,我有点不懂,你为什么设置为0呢,手册上写,设置为0,就是自由结构了,我太懂,请解释一下吧,谢谢哈。
什么时候能算完!!!!
8楼2012-11-21 16:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 volcano93 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +16 Lanmanbaby 2026-08-09 26/1300 2026-08-11 17:16 by 天神眷顾
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +7 Tide man 2026-08-10 8/400 2026-08-11 15:49 by 铁帽子农民
[基金申请] 分享一下我之前已中青C的计划书的filecode +3 布布和一二 2026-08-11 4/200 2026-08-11 15:20 by aasahr
[论文投稿] 职称评审,求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于Filecode分析方法 +3 majunge000 2026-08-10 3/150 2026-08-10 15:46 by lch2012
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 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邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
信息提示
请填处理意见