24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

jianchaoyv

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.9
帖子: 198
在线: 94.8小时
虫号: 657809
注册: 2008-11-19
专业: 通信理论与系统

[交流] 【求助】初学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 ]

回复此楼

» 收录本帖的淘帖专辑推荐

爪机党

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

1楼2009-11-04 21:13:06

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.9
帖子: 198
在线: 94.8小时
虫号: 657809
注册: 2008-11-19
专业: 通信理论与系统

再请问不用脚本运行，是怎么运行的？我刚学，只看到有脚本的例子，直接运行脚本就行了，其他方法不会啊，请指教！谢谢各位好心人！！

赞一下

回复此楼

9楼2009-11-05 08:32:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

查看全部 14 个回答

jianchaoyv

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 964.9
帖子: 198
在线: 94.8小时
虫号: 657809
注册: 2008-11-19
专业: 通信理论与系统

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