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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 求解，在用vasp脚本计算晶格常数时得到的结果一样
121/2返回列表12下一页
查看: 2710  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

天迹浪涯

金虫 (小有名气)

[求助] 求解，在用vasp脚本计算晶格常数时得到的结果一样

我用的INCAR 文件general:
SYSTEM  =  fcc Si
ISTART  =  0 ; ICHARG = 2
ENCUT   =  240
ISMEAR  = 0 ; SIGMAR = 0.1
和POSCAR文件fcc:
   4.3
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
   1
cartesian
0 0 0
KPOINTS文件K-Points
  0
Monkhorst Pack
11  11  11
0   0   0
使用的脚本
#! /bin/bash
BIN=~vw/bin/vasp
rm WAVECAR
for i in  3.5 3.6 3.7 3.8 3.9 4.0 4.1 4.2 4.3 ; do
cat >POSCAR <<!
fcc:
   $i
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
   1
cartesian
0 0 0
!
echo "a= $i" ; $BIN
E=`tail -1 OSZICAR` ; echo $i $E  >>SUMMARY.fcc
done
cat SUMMARY.fcc
得到的结果是这样的
.5 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
3.6 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
3.7 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
3.8 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
3.9 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
4.0 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
4.1 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
4.2 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
4.3 1 F= -.48810251E+01 E0= -.48765216E+01 d E =-.900717E-02
SUMMARY.fcc lines 1-9/9 (END)
这样正常吗
回复此楼

» 猜你喜欢
1楼 2017-07-27 08:53:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

失落的小白鸽

铜虫 (初入文坛)
从输出的文件里面找到最优晶格参数，在根据这个晶格参数写POSCAR，在下一步vasp计算

发自小木虫Android客户端
赞一下(1人) 回复此楼
5楼2017-07-27 09:39:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

失落的小白鸽

铜虫 (初入文坛)
★ ★
天迹浪涯(漫天飘雪代发): 金币+2, 谢谢交流 2017-07-27 11:28:29
你不要在文件夹里面加上POSCAR，你的脚本就是用来产生POSCAR啊，你在放个POSCAR他不就只用4.3一个数值了，根本没循环

发自小木虫Android客户端
赞一下 回复此楼
2楼2017-07-27 09:08:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天迹浪涯

金虫 (小有名气)
引用回帖:
2楼: Originally posted by 失落的小白鸽 at 2017-07-27 09:08:23
你不要在文件夹里面加上POSCAR，你的脚本就是用来产生POSCAR啊，你在放个POSCAR他不就只用4.3一个数值了，根本没循环

就是准备INCAR和KPOINTS和POTCAR三个文件，是先运行这个脚本呢？还是先提交任务进行计算后再运行脚本呢？

发自小木虫Android客户端
赞一下 回复此楼
3楼2017-07-27 09:13:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

失落的小白鸽

铜虫 (初入文坛)
引用回帖:
3楼: Originally posted by 天迹浪涯 at 2017-07-27 09:13:45
就是准备INCAR和KPOINTS和POTCAR三个文件，是先运行这个脚本呢？还是先提交任务进行计算后再运行脚本呢？
...

把INCAR KPOINTS POTCAR 和脚本放在一个文件夹，然后运行脚本，

发自小木虫Android客户端
赞一下 回复此楼
4楼2017-07-27 09:38:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天迹浪涯

金虫 (小有名气)
引用回帖:
5楼: Originally posted by 失落的小白鸽 at 2017-07-27 09:39:26
从输出的文件里面找到最优晶格参数，在根据这个晶格参数写POSCAR，在下一步vasp计算

运行以后，产生了一个完整的POSCAR 和一个SUMMARY文件里面只有从5.1到5.7这几个数，没有其他结果了，这是不是哪里出问题了？

发自小木虫Android客户端
赞一下 回复此楼
6楼2017-07-27 10:22:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

失落的小白鸽

铜虫 (初入文坛)
你脚本都错了，si怎么会一个原包只有一个原子

发自小木虫Android客户端
赞一下 回复此楼
7楼2017-07-27 11:03:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

失落的小白鸽

铜虫 (初入文坛)
你运行的时候linux会有错误提示，你把提示给我看看

发自小木虫Android客户端
赞一下 回复此楼
8楼2017-07-27 11:05:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天迹浪涯

金虫 (小有名气)
引用回帖:
8楼: Originally posted by 失落的小白鸽 at 2017-07-27 11:05:08
你运行的时候linux会有错误提示，你把提示给我看看

真是麻烦您了，这是运行后下面显示的./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
a= 5.2
./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
a= 5.3
./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
a= 5.4
./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
a= 5.5
./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
a= 5.6
./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
a= 5.7
./loop.sh: line 16: ~vw/bin/vasp: No such file or directory
tail: cannot open `OSZICAR' for reading: No such file or directory
5.1
5.2
5.3
5.4
5.5
5.6
5.7
赞一下 回复此楼
9楼2017-07-27 14:50:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

失落的小白鸽

铜虫 (初入文坛)
你把vasp执行文件放到主bin目录下，脚本里不用定义环境变量，直接time vasp，你这样可能是所在shell读不到，你可以把sh换成bash，用bash执行试试，参考侯老师的脚本

发自小木虫Android客户端
赞一下 回复此楼
10楼2017-07-27 22:46:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 天迹浪涯 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +16 爱看书的可乐 2026-08-26 18/900 2026-08-29 10:18 by winsaint
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭