| 查看: 978 | 回复: 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 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 爪机党 |
» 猜你喜欢
285求调剂
已经有6人回复
-
085701环境工程求调剂
已经有6人回复
-
291求调剂
已经有6人回复
-
315分求调剂
已经有6人回复
-
277跪求调剂
已经有8人回复
-
295求调剂
已经有5人回复
-
295材料工程专硕求调剂
已经有3人回复
-
083000学硕274求调剂
已经有4人回复
-
275求调剂
已经有8人回复
-
一志愿211院校 344分 东北农业大学生物学学硕，求调剂
已经有8人回复
|
上面的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
5楼2009-11-04 23:23:46
tjyl
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3218.1
- 红花: 2
- 帖子: 576
- 在线: 156.6小时
- 虫号: 765184
- 注册: 2009-05-07
- 专业: 生物无机化学
- 管辖: 程序语言
6楼2009-11-04 23:30:04
7楼2009-11-04 23:31:59
zhangfan2192
金虫 (著名写手)
- 应助: 9 (幼儿园)
- 金币: 1007.3
- 散金: 62
- 红花: 3
- 帖子: 1170
- 在线: 171.3小时
- 虫号: 648132
- 注册: 2008-11-07
- 性别: GG
- 专业: 理论和计算化学
8楼2009-11-05 08:30:05
9楼2009-11-05 08:32:04
aylayl08
荣誉版主 (文坛精英)
至尊天蝎
- 1ST强帖: 1
- 应助: 3 (幼儿园)
- 贵宾: 4.825
- 金币: 19183.7
- 散金: 2180
- 红花: 34
- 沙发: 48
- 帖子: 16357
- 在线: 824.1小时
- 虫号: 533169
- 注册: 2008-03-26
- 专业: 无机非金属基复合材料
- 管辖: 量子化学
10楼2009-11-05 08:42:07
6