| 查看: 971 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 爪机党 |
» 猜你喜欢
274求调剂
已经有11人回复
-
材料求调剂
已经有3人回复
-
化学工程085602 305分求调剂
已经有4人回复
-
07化学303求调剂
已经有4人回复
-
一志愿河工大 081700 276求调剂
已经有3人回复
-
312求调剂
已经有6人回复
-
321求调剂
已经有3人回复
-
086000生物与医药292求调剂
已经有7人回复
-
085600 材料与化工 329分求调剂
已经有5人回复
-
材料与化工304求B区调剂
已经有3人回复
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
|
上面的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 ：电子结构
4楼2009-11-04 22:31:45
5楼2009-11-04 23:23:46
