24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3447  |  回复: 13
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

wzh上善若水

木虫 (著名写手)

[交流] 【求助】vasp优化晶格参数a 已有8人参与

vasp优化晶格参数a,侯柱峰教程中:
准备好一个poscar文件,得到Volume-Etotal的数据,运行脚本程序
#!/bin/sh
rm WAVECAR
for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30
do
cat > POSCAR < Al-fcc
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
1
Direct
0.0 0.0 0.0
!
echo / a = $i angstrom /; time vasp
E=‘grep /TOTEN0 OUTCAR | tail -1 | awk -{printf /%12.6f \n0, $5 }.‘
V=‘grep /volume0 OUTCAR | tail -1 | awk -{printf /%12.4f \n0 , $5}.‘
echo $V $E >>EtVo.dat
done
小弟菜鸟,请问如何运行脚本
回复此楼

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

VASP 分子动力学

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

知 止 而 后 有 定 ; 定 而 后 能 靜 ; 靜 而 后 能 安 ; 安 而 后 能 慮 ; 慮 而 后 能 得
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowu787

木虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
引用回帖:
Originally posted by bingmou at 2010-10-25 15:15:40:
你要注意V=`后面的引号是反引号`,也就是键盘上1左边的纳和键,而不是单引号’。
同样最后的‘也是反引号`

time是用来看这个计算要花多少时间,属于让自己多点耐心的设置

谢谢您的解答:
试着做这个例子时,发现LDA和GGA的差别很大,下图是一篇文献上用VASP-GGA-PW91算的fcc(最下面的那根红线)的结果,经过单位23.06换算后,我用GGA-PW91算的结果在85Kcal/mol,差别怎么这么大?请版主帮忙分析一下原因。谢谢!

LDA-的计算结果
3.80 13.720000 -4.071701
3.85 14.270000 -4.117942
3.90 14.830000 -4.147312
3.95 15.410000 -4.162867
4.00 16.000000 -4.165279
4.05 16.610000 -4.156421
4.10 17.230000 -4.137987
4.15 17.870000 -4.111258
4.20 18.520000 -4.077605
4.25 19.190000 -4.038274
4.30 19.880000 -3.993437

GGA-的计算结果
a= 3.80    -3.533891      13.720000
a= 3.85    -3.598180      14.270000
a= 3.90    -3.644692      14.830000
a= 3.95    -3.675400      15.410000
a= 4.00    -3.692204      16.000000
a= 4.05    -3.696934      16.610000
a= 4.10    -3.691228      17.230000
a= 4.15    -3.676476      17.870000
a= 4.20    -3.653888      18.520000
a= 4.25    -3.624569      19.190000
a= 4.30    -3.589560      19.880000



[ Last edited by xiaowu787 on 2010-10-25 at 17:25 ]
9楼2010-10-25 17:22:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

1003307944

银虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
bingmou:什么意思? 2010-10-21 19:44:43
除非你一个人用超算计算机
3楼2010-10-21 17:03:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
bingmou(金币+1):感谢交流 2010-10-21 19:44:07
先学习一下linux的的基本操作。运行脚本,首先要保证脚本有执行权限
$chmod 777 filename
$./filename
4楼2010-10-21 17:06:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1003307944

银虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1):谢谢指教 2010-10-22 08:22:05
你这个脚本执行的时候可能会算,但是有多人提交任务的时候,你的任务会被停止,,还不如把这个脚本中对参量的改变,采用手动的方式,算完一个,再变动下。。
5楼2010-10-21 21:09:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂 +7 嗯滴的基本都 2026-03-27 7/350 2026-03-28 04:19 by fmesaito
[考研] 265求调剂11408 +3 刘小鹿lu 2026-03-27 3/150 2026-03-27 20:53 by nihaoar
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 339求调剂 +4 烤麦芽 2026-03-27 5/250 2026-03-27 13:23 by 752105528
[考研] 348求调剂 +4 小懒虫不懒了 2026-03-27 5/250 2026-03-27 12:47 by 果果妈咪
[考研] 085601 材料工程 313分 求调剂 +5 Ong3 2026-03-27 5/250 2026-03-27 12:24 by goldfish51
[考研] 286求调剂 +4 lim0922 2026-03-26 4/200 2026-03-27 10:28 by guoweigw
[考研] 求调剂 +6 林之夕 2026-03-24 6/300 2026-03-27 08:38 by hypershenger
[硕博家园] 招收生物学/细胞生物学调剂 +3 IceGuo 2026-03-26 4/200 2026-03-27 05:35 by user003
[考研] 325求调剂 +5 李嘉图·S·路 2026-03-23 5/250 2026-03-27 00:42 by wxiongid
[考研] 336材料求调剂 +7 陈滢莹 2026-03-26 9/450 2026-03-27 00:20 by wxiongid
[考研] 071000生物学求调剂,初试成绩343 +6 小小甜面团 2026-03-25 6/300 2026-03-26 23:01 by 不吃魚的貓
[考研] 求调剂 +8 Auroracx 2026-03-22 8/400 2026-03-26 19:55 by 不吃魚的貓
[考研] 环境专硕324分求调剂推荐 +5 轩小宁—— 2026-03-26 5/250 2026-03-26 12:05 by i_cooler
[考研] 0854AI CV方向招收调剂 +4 章小鱼567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[考研] 282求调剂 +3 wcq131415 2026-03-24 3/150 2026-03-25 12:16 by userper
[考研] 086003食品工程求调剂 +6 淼淼111 2026-03-24 6/300 2026-03-25 10:29 by 3Strings
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 一志愿北京化工大学070300 学硕336求调剂 +5 vv迷 2026-03-21 8/400 2026-03-22 14:20 by ColorlessPI
信息提示
请填处理意见