24小时热门版块排行榜    

查看: 1245  |  回复: 13
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jianchaoyv

金虫 (小有名气)

[交流] 【求助】初学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 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

y1ding

铁杆木虫 (著名写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
wuli8(金币+1,VIP+0):3q 11-4 22:41
E= `grep "TOTEN"  OUTCAR | tail -1 | awk '{printf   "%12.6f \n ", $5 } '`
3楼2009-11-04 22:30:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

jianchaoyv

金虫 (小有名气)

上面的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的回帖

y1ding

铁杆木虫 (著名写手)

★ ★
wuli8(金币+2,VIP+0):3ks 11-4 22:40
应该是`,不是'
这和linux有关,与vasp无关
初学vasp,不要用脚本
究竟是脚本问题,还是vasp的问题,初学者往往不会判断
4楼2009-11-04 22:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)

★
小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by y1ding at 2009-11-4 22:31:
应该是`,不是'
这和linux有关,与vasp无关
初学vasp,不要用脚本
究竟是脚本问题,还是vasp的问题,初学者往往不会判断

同意
是的!!!当时我也碰到了这样的情况!!呵呵!不是单引号!!而是`呵呵!!
5楼2009-11-04 23:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:37 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见