24小时热门版块排行榜    

查看: 965  |  回复: 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的回帖

zdd2008

金虫 (职业作家)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by tjyl at 2009-11-4 23:30:
建议修改为使用$( command )结构来代替 ` command `。


呵呵!!也是一种好方法哈!!!呵呵!
7楼2009-11-04 23:31:59
已阅   回复此楼   关注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

铁杆木虫 (著名写手)

★ ★
小木虫(金币+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的回帖

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 359求调剂 +3 王了个楠 2026-03-25 3/150 2026-03-25 12:50 by Dyhoer
[考研] 303求调剂 +6 元夕元 2026-03-20 7/350 2026-03-25 12:00 by edmund7
[考研] 281求调剂 +4 Koxui 2026-03-24 5/250 2026-03-25 11:38 by userper
[考研] 0703化学求调剂 +6 奶油草莓. 2026-03-22 7/350 2026-03-25 10:00 by shangxh
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 一志愿211 初试270分 求调剂 +5 谷雨上岸 2026-03-23 6/300 2026-03-24 16:32 by laoshidan
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +4 曼111 2026-03-24 5/250 2026-03-24 15:44 by 星空星月
[考研] 277分求调剂,跨调材料 +3 考研调剂lxh 2026-03-24 3/150 2026-03-24 13:52 by JourneyLucky
[考研] 【双一流院校新能源、环境材料,材料加工与模拟招收大量调剂】 +4 Higraduate 2026-03-22 7/350 2026-03-24 11:23 by 种大树
[基金申请] 请教下大家 2026年国家基金申请是双盲审吗? +3 lishucheng1 2026-03-22 5/250 2026-03-24 08:22 by gltch
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 9/450 2026-03-23 12:57 by yqdszhdap-
[考研] 291求调剂 +5 孅華 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[考研] 求调剂一志愿海大,0703化学学硕304分,有大创项目,四级已过 +6 幸运哩哩 2026-03-22 10/500 2026-03-22 20:10 by edmund7
[考研] 一志愿华中农业071010,总分320求调剂 +5 困困困困坤坤 2026-03-20 6/300 2026-03-22 17:41 by hxsm
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
[基金申请] 山东省面上项目限额评审 +4 石瑞0426 2026-03-19 4/200 2026-03-22 08:50 by Wei_ren
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
信息提示
请填处理意见