应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】初学VASP者的一个问题
141/2返回列表12下一页
查看: 1170  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

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 ]
回复此楼

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

爪机党

» 猜你喜欢
1楼 2009-11-04 21:13:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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 } '`
赞一下(21人) 回复此楼
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的问题,初学者往往不会判断
赞一下(20人) 回复此楼
4楼2009-11-04 22:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

金虫 (职业作家)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by y1ding at 2009-11-4 22:31:
应该是`,不是'
这和linux有关,与vasp无关
初学vasp,不要用脚本
究竟是脚本问题,还是vasp的问题,初学者往往不会判断

同意
是的！！！当时我也碰到了这样的情况！！呵呵！不是单引号！！而是`呵呵！！
赞一下(1人) 回复此楼
5楼2009-11-04 23:23:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tjyl

金虫 (正式写手)
★
aylayl08(金币+1,VIP+0):感谢解答，欢迎常来 11-5 08:40
建议修改为使用$( command )结构来代替 ` command `。
引用回帖:
Originally posted by zdd2008 at 2009-11-4 23:23:


同意
是的！！！当时我也碰到了这样的情况！！呵呵！不是单引号！！而是`呵呵！！

赞一下(6人) 回复此楼
6楼2009-11-04 23:30:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zdd2008

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


呵呵！！也是一种好方法哈！！！呵呵！
赞一下(1人) 回复此楼
7楼2009-11-04 23:31:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfan2192

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
呵呵 引号有时候真会害人
赞一下(1人) 回复此楼
相逢是缘！
8楼2009-11-05 08:30:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jianchaoyv

金虫 (小有名气)
再请问不用脚本运行，是怎么运行的？我刚学，只看到有脚本的例子，直接运行脚本就行了，其他方法不会啊，请指教！谢谢各位好心人！！
赞一下 回复此楼
9楼2009-11-05 08:32:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

至尊天蝎

优秀版主
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by y1ding at 2009-11-4 22:31:
应该是`,不是'
这和linux有关,与vasp无关
初学vasp,不要用脚本
究竟是脚本问题,还是vasp的问题,初学者往往不会判断

同样请教y1ding专家，初学VASP，不用脚本怎样运行啊？谢谢
赞一下(1人) 回复此楼
珍惜拥有，快乐人生!欢迎光临计算模拟区(*^__^*)。。。。。。。。。。。。。。
10楼2009-11-05 08:42:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jianchaoyv 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +37 医学老男孩 2026-08-13 87/4350 2026-08-18 06:35 by 医学老男孩
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了，能看出什么问题？ +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 今天系统多次维护，明天很可能放榜！ +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
关闭
请填处理意见
关闭