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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp中，结构参数优化时，出现错误，求助！谢谢！
41/1返回列表
查看: 1253  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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: `!'
请问这是什么错误啊？谢谢！我改了改！附近的内容，还是不行。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2013-11-15 20:38:39
已阅   回复此楼   关注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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭