24小时热门版块排行榜    

查看: 2635  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统查不到 +9 董八千 2026-08-26 9/450 2026-08-26 12:30 by yanmofang
[基金申请] 怎么查啊 +4 huang1991js 2026-08-26 4/200 2026-08-26 11:26 by panjy.cn
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 8/400 2026-08-26 11:03 by zjlp
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
请填处理意见