24小时热门版块排行榜    

查看: 1521  |  回复: 3

cuicui2397

新虫 (小有名气)

[求助] vasp中,结构参数优化时,出现错误,求助!谢谢!

侯老师的vasp教程中,针对复杂情况的结构参数优化出现问题,求助!

复杂的情况,以六角结构Mg的晶格常数为例
1.先进行一次体积保持不变的粒子驰豫的计算(ISIF=5)
2再把优化得到的CONTCAR拷贝成POSCAR,进行一次静态的计算
3对a的值取10个左右的点,每个点重复上面两步,得到静态计算下volume-Etot的关系。
这三部的脚步可以通过以下脚本run_cell来进行,内容如下:
#!/bin/sh
rm WAVECAR
for i in   5.20  5.23  5.26  5.29 5.32 5.35 5.38  5.41 5.44  5.47
do
cat > POSCAR << !
ureadopef
$i
0                      -1                0
0.8660254        0.5              0
0                      0               1.6230530
2
Direct
0.666666667  0.333333333    0.75
0.333333333  0.666666667    0.25
!
cp INCAR.relax INCAR
echo "a=$i angestrom"; mpiexec -np 8 vasp
cp CONTCAR POSCAR
cp INCAR.static INCAR
echo "a= $i angstrom ";mpiexec -np 8 vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf  "%12.6f \n", $5}'`
V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f \n", $5}'`
echo $V $E >>EtVo.dat
done
4.采用状态方程拟合得到的平衡态下的体积,体弹性模量。
5.在该体积下,重复上面第2步和第3步,得到平衡状态下的其他晶胞参数。这一步也就是:在的到了E(V)曲线后,通过状态方程拟合得到平衡状态下的体积,计算出上面脚本中的变量$i的值,并改变$i的循环值,再运行run_cell计算一次,得到其他的结构参数c和位置u。


我的问题是:在执行sh run.cell命令时,出现错误提示是这样的:
run_cell: line 15: syntax error near unexpected token `newline'
run_cell: line 15: `!'
请问这是什么错误啊?谢谢!我改了改!附近的内容,还是不行。
回复此楼

» 猜你喜欢

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

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

cuicui2397

新虫 (小有名气)

是我自己大意了,已经知道什么错误了。
2楼2013-11-15 20:57:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

不要烦恼

银虫 (小有名气)

引用回帖:
2楼: Originally posted by cuicui2397 at 2013-11-15 20:57:12
是我自己大意了,已经知道什么错误了。

您好楼主,我初涉vasp,我在运行这里时候通过状态方程拟合得到平衡状态下的体积,然后如何实现“计算出上面脚本中的变量$i的值,并改变$i的循环值,再运行run_cell计算一次,得到其他的结构参数c和位置u。”呢?
3楼2016-08-30 19:22:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
3楼: Originally posted by 不要烦恼 at 2016-08-30 19:22:02
您好楼主,我初涉vasp,我在运行这里时候通过状态方程拟合得到平衡状态下的体积,然后如何实现“计算出上面脚本中的变量$i的值,并改变$i的循环值,再运行run_cell计算一次,得到其他的结构参数c和位置u。”呢?...

你好,我也是刚学VASP,不知道你的问题解决了吗?如何实现“计算出上面脚本中的变量$i的值,并改变$i的循环值,再运行run_cell计算一次,得到其他的结构参数c和位置u。”呢?
4楼2016-10-31 18:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cuicui2397 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 大佬帮着分析今年的评审意见?该课题连续申请了三年了,明年要不要换个方向。 100+3 cai-11138 2026-08-31 10/500 2026-09-06 19:40 by babero
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 19:21 by eObJPa4gWmp6
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:21 by eObJPa4gWmp6
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 18:09 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 18:09 by eObJPa4gWmp6
[文学芳草园] 两块石头 +7 阿美_Lml888 2026-09-03 9/450 2026-09-06 17:48 by myrtle
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 6/300 2026-09-06 17:21 by eObJPa4gWmp6
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 9/450 2026-09-06 17:09 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 17:09 by eObJPa4gWmp6
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 13:09 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[基金申请] 面上没中,邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[硕博家园] 哈尔滨工业大学韩晓军教授课题组招收2027年硕士推免生及博士研究生 +3 灯灯灯灯DDDD 2026-09-03 3/150 2026-09-03 21:22 by 科研皮皮猪
[基金申请] 要骂人了,新模版改版就是要淡化问题凝练这种虚的东西,结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
请填处理意见