24小时热门版块排行榜

返回列表

当前主题已经存档。

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

jianchaoyv

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.9
帖子: 198
在线: 94.8小时
虫号: 657809
注册: 2008-11-19
专业: 通信理论与系统

[交流] 【求助】初学VASP者的一个问题

请各位帮忙指教一个VASP的例子，一直找不到原因，自己写的第一个例子就出错郁闷的不行了！请高手指点一下我这菜鸟吧。

INCAR:
SYSTEM=Si-Diamond
ENCUT=400
ISTART=0;ICHARG=2
ISMEAR=-5
PREC=Accurate

POSCAR:
Si-Diamond
5.430
0.0  0.5  0.5
0.5  0.0  0.5
0.5  0.5  0.0
2Direct
0.0  0.0  0.0
0.25 0.25 0.25

KPOINTS:
auto
0
Monhkorst-pack
9 9 9
0.0  0.0  0.0

SCRIPT:
#!/bin/bash
rm WAVECAR
for i in 150 200 250 300 350 400
do
cat >INCAR < SYSTEM=Si-Diamond
ENCUT=$i
ISTART=0;ICHARG=2
ISMEAR=-5
PREC=Accurate
!
echo "ENCUT=$i ev";time vasp
E= ' grep "TOTEN"  OUTCAR | tail -1 | awk '{printf "%12.6f \n ", $5 } ''
echo $i $E >>comment
done

最后得到的结果文件comment：
150
200
250
300
350
400
没有每个EVNCUT对应的总能，十分感谢你能耐心指点！

[ Last edited by jianchaoyv on 2009-11-5 at 15:59 ]

回复此楼

» 收录本帖的淘帖专辑推荐

爪机党

» 猜你喜欢

课题组招2027级博士上海工程技术大学激光智能制造方向已经有4人回复
现代”学阀”该如何界定已经有8人回复
上海工程技术大学激光智能制造课题组 2027级博士研究生招生已经有3人回复
实名举报：华南理工大学物理与光电学院副院长李志远婚内两次出轨女学生（博士）已经有16人回复
上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告已经有11人回复
上海工程技术大学激光智能制造课题组招收博士研究生已经有10人回复
申博发邮件已经有11人回复
各位大神，目前国内有哪些比较好用的逆合成软件？已经有12人回复
我的奶奶已经有3人回复
国社科系统bug了，是不是要放榜了？已经有9人回复

1楼 2009-11-04 21:13:06

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

zdd2008

金虫 (职业作家)

应助: 1 (幼儿园)
金币: 5.4
散金: 2613
红花: 12
帖子: 3845
在线: 108.4小时
虫号: 708180
注册: 2009-02-25
专业: 金属结构材料

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

引用回帖:

Originally posted by y1ding at 2009-11-4 22:31:
应该是`,不是'
这和linux有关,与vasp无关
初学vasp,不要用脚本
究竟是脚本问题,还是vasp的问题,初学者往往不会判断

同意
是的！！！当时我也碰到了这样的情况！！呵呵！不是单引号！！而是`呵呵！！

赞一下(1人)

回复此楼

5楼2009-11-04 23:23:46

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

查看全部 14 个回答

jianchaoyv

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.9
帖子: 198
在线: 94.8小时
虫号: 657809
注册: 2008-11-19
专业: 通信理论与系统

上面的POTCAR：
PAW Si 02Apr1999
4.00000000000000000
parameters from PSCTR are:
VRHFIN =Si: s2p2
LEXCH  = CA
EATOM  = 103.4338 eV, 7.6022 Ry

TITEL  = PAW Si 02Apr1999
LULTRA =       F use ultrasoft PP ?
IUNSCR =       1 unscreen: 0-lin 1-nonlin 2-no
RPACOR = 1.500 partial core radius
POMASS = 28.085; ZVAL = 4.000 mass and valenz
RCORE  = 1.900 outmost cutoff radius
RWIGS  = 2.480; RWIGS  = 1.312 wigner-seitz radius (au A)
ENMAX  =  245.704; ENMIN  =  184.278 eV
ICORE  =       2 local potential
LCOR =       T correct aug charges
LPAW =       T paw PP
EAUG =  322.069
DEXC = -.006
RMAX = 2.941 core radius for proj-oper
RAUG = 1.300 factor for augmentation sphere
RDEP = 1.993 core radius for depl-charge
QCUT = -4.250; QGAM = 8.499 optimization parameters

Description
   l    E    TYP  RCUT TYP  RCUT
   0 .000    23  1.900
   0 .000    23  1.900
   1 .000    23  1.900
   1 .000    23  1.900
   2 .000    7  1.900
Error from kinetic energy argument (eV)
NDATA  =    100
STEP = 20.000 1.050
  10.3    9.30    8.81    7.90    7.48    6.69    5.98    5.65
  5.03    4.48    4.23    3.76    3.34    2.97    2.64    2.35
  2.08    1.85    1.64    1.38    1.22    1.09    .910    .808
  .718    .601    .502    .445    .371    .309    .256    .212
  .176    .145    .120    .989E-01  .768E-01  .637E-01  .501E-01  .422E-01
  .342E-01  .282E-01  .237E-01  .204E-01  .178E-01  .157E-01  .139E-01  .119E-01
  .104E-01  .877E-02  .723E-02  .610E-02  .485E-02  .380E-02  .281E-02  .221E-02
  .178E-02  .147E-02  .130E-02  .123E-02  .120E-02  .116E-02  .111E-02  .101E-02
  .905E-03  .765E-03  .630E-03  .514E-03  .416E-03  .362E-03  .326E-03  .307E-03
  .291E-03  .271E-03  .240E-03  .208E-03  .173E-03  .144E-03  .123E-03  .112E-03
  .104E-03  .968E-04  .876E-04  .752E-04  .630E-04  .536E-04  .476E-04  .437E-04
  .403E-04  .361E-04  .308E-04  .263E-04  .228E-04  .209E-04  .192E-04  .170E-04
  .145E-04  .126E-04  .112E-04  .103E-04
END of PSCTR-controll parameters
local part
98.2657514061040160
.59263674E+01 .59326757E+01 .59422539E+01 .59582113E+01 .59805393E+01
.60092267E+01 .60442607E+01 .60856276E+01 .61333132E+01 .61873035E+01

赞一下

回复此楼

2楼2009-11-04 21:14:52

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