| 查看: 977 | 回复: 13 | ||||
| 当前主题已经存档。 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【求助】初学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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 爪机党 |
» 猜你喜欢
材料调剂
已经有7人回复
-
化学工程085602 305分求调剂
已经有10人回复
-
289求调剂
已经有15人回复
-
291 求调剂
已经有7人回复
-
274求调剂
已经有14人回复
-
北京林业大学硕导招生广告
已经有3人回复
-
309求调剂
已经有5人回复
-
292求调剂
已经有8人回复
-
求调剂
已经有4人回复
-
一志愿 西北大学 总分282 英语一62 求调剂
已经有3人回复
5楼2009-11-04 23:23:46
|
上面的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
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
3楼2009-11-04 22:30:03
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
4楼2009-11-04 22:31:45
