24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

gleerat

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3926.9
帖子: 322
在线: 102.6小时
虫号: 860161

[交流] 【求助完成】如何使用pbnd.f画能带图

我计算能带时参考了http://blogold.chinaunix.net/u2/85704/showart_2075157.html，处理结果时参考了http://old.blog.edu.cn/user1/11542/archives/2006/1232097.shtml。
用pbnd.f处理的结果在附件中，但是结果却显然不是那么一回事。
这是我用来计算能带的文件：

引用回帖:

#!/bin/bash
echo "AFM1:static"
date
cat >INCAR < System = CaFe2As2 afm1:dos
PREC = Accurate
ISPIN = 2
MAGMOM = 2*0 2*-3 2*3 4*0
#ISMEAR = -5
ISTART = 0
ICHARG = 2
SIGMA = 0.18
EDIFF = 1E-5
EDIFFG = -1E-4
GGA = PE
ENCUT = 400
POTIM = 0.2
NSW = 0
ISIF = 2
IBRION = 2

LORBIT = 11
LREAL = Auto
NPAR = 1
LPLANE = .TRUE.
LSCALU = .FALSE.
NSIM = 4
!
cat > KPOINTS < Monhkorst-Pack
0
Monhkorst-Pack
8 8 8
0.0 0.0 0.0
!
/usr/local/bin/mpirun -np 4 /usr/bin/vasp
echo "AFM1:band"
date
cat >INCAR < System = CaFe2As2 afm1:dos
PREC = Accurate
ISPIN = 2
MAGMOM = 2*0 2*-3 2*3 4*0
#ISMEAR = -5
ISTART = 1
ICHARG = 11
SIGMA = 0.18
EDIFF = 1E-5
EDIFFG = -1E-4
GGA = PE
ENCUT = 400
POTIM = 0.2
NSW = 0
ISIF = 2
IBRION = 2

LORBIT = 11
LREAL = Auto
NPAR = 1
LPLANE = .TRUE.
LSCALU = .FALSE.
NSIM = 4
!
cat > KPOINTS < k-points along high symmetry lines
50 !50 intersections
Line mode
rec
0.0000000000    0.0000000000    0.0000000000 !G
0.5000000000    0.5000000000    0.0000000000 !Y
0.5000000000    0.5000000000    0.5000000000 !T
0.0000000000    0.0000000000    0.5000000000 !Z
0.0000000000    0.5000000000    0.0000000000 !S
0.0000000000    0.0000000000    0.0000000000 !G
!
/usr/local/bin/mpirun -np 4 /usr/bin/vasp

因为输入的能带文件中只有两列数据，所以搞出来的图中就只有一条“看似震荡”的曲线。杂乱无章，根本不是能带图形。
是结果有错误还是我处理的问题。

用pbnd.f处理的结果在附件中

[ Last edited by gleerat on 2011-3-15 at 09:42 ]

回复此楼

» 猜你喜欢

依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有13人回复
AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有11人回复
表哥与省会女结婚，父母去帮带孩子被省会女气回家生重病了已经有9人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有10人回复
依托企业入选了国家启明计划青年人才。有无高校可以引进的。已经有11人回复

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

【求助必读】文献版求助者必读（2013版) 已经有351人回复
“【有奖活动】说说大家都从哪里出发开始查文献”总结已经有21人回复
在求助完成时，你知道你帖子CEI是如何发放的吗已经有14人回复
【求助】求助关于GAR2PED指认苯环已经有3人回复
【求助】如何确定微区（直径约为3nm区域）成分已经有4人回复
【求助】吉非替尼已经有2人回复
【求助】Linux 下MS计算已经有7人回复
【求助】NMR基础问题求解已经有4人回复
【求助】Distillation design and control using ASPEN simulation 已经有3人回复
【求助】关于Dmol3的吸附问题已经有6人回复
【原创】分享考博的一些心得(更新在75楼，字数限制) ！已经有488人回复
【求助】求助真空干燥箱使用的相关注意事项已经有2人回复
【求助】两个MS在linux系统下安装和使用的问题已经有3人回复
【求助】有关使用Abinit能带计算的一些疑问已经有0人回复
【外文书籍求助版】版规（2010版）已经有34人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2011-03-14 10:58:56

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

gleerat

木虫 (正式写手)

应助: 1 (幼儿园)
金币: 3926.9
帖子: 322
在线: 102.6小时
虫号: 860161

引用回帖:

Originally posted by bingmou at 2011-03-14 15:44:12:
你应该说你想得到什么结果，而实际给出了什么结果，这两个之间的差异在哪里？

但是为什么给出这样的结果，具体原因则是摸不着头脑。

还有在bnd.dat中给出的两列数据是什么意思？第二列是能量，那么第一列是什么东西？

引用回帖:

  0.000000    -22.43620
  0.001751    -22.43620
  0.003502    -22.43610
  0.005252    -22.43610
  0.007003    -22.43610
  0.008754    -22.43610
  0.010505    -22.43610
  0.012255    -22.43610
  0.014006    -22.43610
  0.015757    -22.43610
  0.017508    -22.43610
  0.019258    -22.43610
  0.021009    -22.43610
  0.022760    -22.43610
  0.024511    -22.43600
  0.026262    -22.43600
  0.028012    -22.43600
  0.029763    -22.43600
  0.031514    -22.43600
  0.033265    -22.43600
  0.035015    -22.43590
  0.036766    -22.43590
  0.038517    -22.43590
  0.040268    -22.43590
  0.042019    -22.43590
  0.043769    -22.43580
  0.045520    -22.43580
  0.047271    -22.43580
  0.049022    -22.43580
  0.050772    -22.43580
  0.052523    -22.43580
  0.054274    -22.43570
  0.056025    -22.43570
  0.057775    -22.43570
  0.059526    -22.43570
  0.061277    -22.43570
  0.063028    -22.43570
  0.064779    -22.43560
  0.066529    -22.43560
  0.068280    -22.43560
  0.070031    -22.43560
……………………
……………………

[ Last edited by gleerat on 2011-3-14 at 17:32 ]