24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2274  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿085802 323分求调剂 +9 drizzle_9 2026-04-12 10/500 2026-04-12 14:51 by drizzle_9
[考研] 299求调剂 +8 ZVVZ13 2026-04-08 8/400 2026-04-12 00:40 by 蓝云思雨
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 085400 328分 求调剂 +10 喂你一个大橙子 2026-04-09 14/700 2026-04-11 19:53 by lqspecial
[考研] 085501机械专硕 302分 不挑专业求调剂 +7 汪某. 2026-04-09 7/350 2026-04-11 14:37 by luhong1990
[考研] 还有化工二轮调剂的学校吗 5+14 化工人999 2026-04-09 48/2400 2026-04-11 10:27 by 89436494
[考研] 288求调剂 +15 代fish 2026-04-09 16/800 2026-04-11 10:26 by wwj2530616
[考研] 初试261 +3 Asht少 2026-04-10 6/300 2026-04-10 16:38 by Asht少
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 一志愿211 0703化学 346分求调剂 +22 土豆er? 2026-04-09 23/1150 2026-04-10 10:58 by 高维春
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 江苏大学 工科调剂 捡漏 +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[考研] 348求调剂 +3 candyyyi 2026-04-09 3/150 2026-04-09 17:20 by 段伟艳
[考研] 软件工程求调剂22软工296分求调剂,接受跨调 +4 yangchen2017 2026-04-08 5/250 2026-04-08 21:56 by 土木硕士招生
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 一志愿哈工大,初试329,求环境科学与工程调剂! +11 余未辛 2026-04-06 11/550 2026-04-08 15:21 by screening
[考研] 化工学硕 285求调剂 +26 Wisjxn 2026-04-07 26/1300 2026-04-08 14:42 by screening
[考研] 生物医药调剂|SCI中科院三区一作+多项科研成果 +8 likangxing 2026-04-07 11/550 2026-04-08 00:02 by lys0704
[考研] 生物学学硕求调剂:351分一志愿南京师范大学生物学专业 +6 …~、王…~ 2026-04-06 7/350 2026-04-06 18:54 by macy2011
信息提示
请填处理意见