24小时热门版块排行榜    

查看: 2243  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿山东大学药学学硕求调剂 +3 开开心心没烦恼 2026-03-23 4/200 2026-03-24 00:06 by 开开心心没烦恼
[考研] 求材料,环境专业调剂 +3 18567500178 2026-03-18 3/150 2026-03-23 23:50 by 热情沙漠
[考研] 求调剂 +7 十三加油 2026-03-21 7/350 2026-03-23 23:48 by 热情沙漠
[考研] 306求0703调剂一志愿华中师范 +7 纸鱼ly 2026-03-21 8/400 2026-03-23 23:31 by chixmc
[考研] 279分求调剂 一志愿211 +17 chaojifeixia 2026-03-19 19/950 2026-03-23 23:26 by 呆呆师姐
[考研] 335分 | 材料与化工专硕 | GPA 4.07 | 有科研经历 +4 cccchenso 2026-03-23 4/200 2026-03-23 23:00 by 徐ckkk
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 305分求调剂(食品工程) +4 Sxy112 2026-03-21 6/300 2026-03-22 15:26 by 无懈可击111
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 材料学硕301分求调剂 +7 Liyouyumairs 2026-03-21 7/350 2026-03-21 22:31 by peike
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 278求调剂 +9 烟火先于春 2026-03-17 9/450 2026-03-21 17:47 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 321求调剂 +9 何润采123 2026-03-18 11/550 2026-03-20 23:19 by JourneyLucky
[考研] 261求B区调剂,科研经历丰富 +3 牛奶很忙 2026-03-20 4/200 2026-03-20 19:34 by JourneyLucky
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
信息提示
请填处理意见