| 查看: 978 | 回复: 2 | ||||
空谷新苗木虫 (小有名气)
|
[求助]
用vasp做结构参数的优化时输出文件中只有一组体积-能量数据，求解答。
|
|
我是参照侯老师《入门指南》中的一个脚本来算的，但输出的EtVo.dat文件中只有一组体积-能量数据，麻烦大家看看是不是我的脚本有问题。 脚本： #!bin/sh rm WAVECAR for i in 4.09 4.19 4.29 4.39 4.49 4.59 4.69 4.79 4.89 4.99 5.09do cat > POSCAR <<! TiO2 1 4.594000000000000 0.000000000000000 0.000000000000000 0.000000000000000 4.594000000000000 0.000000000000000 0.000000000000000 0.000000000000000 2.959000000000000 O Ti 4 2 Direct 0.305000 0.305000 0.000000 -0.305000 -0.305000 0.000000 0.195000 0.805000 0.500000 0.805000 0.195000 0.500000 0.000000 0.000000 0.000000 0.500000 0.500000 0.500000 ! cp INCAR.relax INCAR echo "a=$i angstrom ";time vasp cp CONTCAR POSCAR cp INCAR.static INCAR echo "a=$i angstrom ";time 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 INCAR.relax文件： SYSTEM = bulk rutile TiO2 unit cell ENCUT = 450 ISTART = 0 ICHARG = 2 ISMEAR = 1 SIGMA = 0.05 ISIF=5 IBRION = 2 NSW = 300 POTIM = 0.1 PREC = Accurate EDIFF = 1E-6 EDIFFG = 1E-4 INCAR.static文件： SYSTEM = bulk rutile TiO2 unit cell ENCUT = 450 PREC = Accurate ISTART = 0 ICHARG = 2 ISMEAR = -5 KPOINTS文件： Automatic 0 Monkhorst–Pack 5 5 8 0.0 0.0 0.0 输出的EtVo.dat: 62.4500 -53.101836 62.4500 -53.101836 62.4500 -53.101836 62.4500 -53.101836 62.4500 -53.101836 ······ |
回复此楼» 猜你喜欢
求有机合成大神指点三硫酸乙烯酯（CAS：2793408-99-6）的合成路线
已经有10人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
如何实现卤原子转化
已经有5人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有5人回复
-
考博自荐
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于vasp中，结构优化时所使用的参数POTIM
已经有3人回复
- ELK简单计算失败，求高手解答！
已经有4人回复
- Vasp相比经典分子动力学的优缺点
已经有11人回复
- 用vasp优化CO结构，机器不算，不知道哪出错了？高手来帮帮我！
已经有29人回复
- P4VASP做pdos交流，及后续origin计算d-band center
已经有34人回复
- in文件运行失败（lost atoms）求助！！
已经有11人回复
- 结构优化完结果老出现ZBRENT: fatal error in bracketing
已经有30人回复
- 用VASP优化结构收敛不了，求指教
已经有10人回复
- VASP结构优化时出现的POTIM问题
已经有4人回复
- vasp优化结构中出现的问题
已经有10人回复
- 重金请教采用准谐近似求热膨胀系数的方法
已经有16人回复
- 【求助完结】vasp优化结构运行一段时间后没有出现错误，但是停止了！
已经有20人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】vasp中如何做结构优化？
已经有10人回复
- 【求助成功】VASP中，如何固定磁矩做结构优化
已经有14人回复
- 【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能
已经有19人回复
卡开发发
专家顾问 (著名写手)
Ab Initio Amateur
-
专家经验: +224 - 1ST强帖: 2
- 应助: 47 (小学生)
- 金币: 3539.8
- 散金: 14
- 红花: 223
- 帖子: 2477
- 在线: 1254小时
- 虫号: 1369423
- 注册: 2011-08-16
- 性别: GG
- 专业: 金属材料的磨损与磨蚀
- 管辖: 第一性原理
★
ljw4010: 金币+1, 谢谢交流！ 2015-03-27 20:34:25
ljw4010: 金币+1, 谢谢交流！ 2015-03-27 20:34:25
|
在POSCAR需要让晶格常数成为变量，这里我只好固定一下c/a的值。 。。。。。。 rm WAVECAR for i in 4.09 4.19 4.29 4.39 4.49 4.59 4.69 4.79 4.89 4.99 5.09do cat > POSCAR <<! TiO2 1 $i 0.000000000000000 0.000000000000000 0.000000000000000 $i 0.000000000000000 0.000000000000000 0.000000000000000 $[$i*0.6441] 。。。。。。。 不过 $[$i*0.6441]这样的表达我不确定行不行，如果不行的话你可以找找linux脚本相关的资料进行修改。 |
» 本帖已获得的红花（最新10朵）
空谷新苗2楼2015-03-23 11:10:50
空谷新苗
木虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 3797.8
- 散金: 182
- 红花: 5
- 帖子: 238
- 在线: 158.8小时
- 虫号: 2218173
- 注册: 2013-01-02
- 性别: GG
- 专业: 光学
3楼2015-03-23 15:25:38
精华III: