| 查看: 1350 | 回复: 9 | |||
[交流]
【求助】做状态方程计算时候得到的EtVo.dat是空的 已有4人参与
|
|
下面是我按照侯博士教程写的脚本： #!/bin/sh for i in 7.846 7.946 8.046 8.146 8.246 8.346 8.446 8.546 8.646 8.746 8.846 8.946 9.046 do cat > POSCAR < aaa $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 2 4 8 Direct 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.2500000000000000 0.2500000000000000 0.2500000000000000 0.6205000281333946 0.6205000281333943 0.6205000281333943 0.6205000281333944 0.6205000281333944 -1.8615000844001834 0.6205000281333947 -1.8615000844001834 0.6205000281333946 -1.8615000844001834 0.6205000281333947 0.6205000281333944 0.3849999904632571 0.3849999904632571 0.3849999904632571 0.8499000072479258 0.8499000072479258 0.8499000072479256 0.3849999904632571 0.3849999904632571 -1.1549999713897712 0.8499000072479258 0.8499000072479258 -2.5497000217437775 0.3849999904632572 -1.1549999713897714 0.3849999904632571 0.8499000072479259 -2.5497000217437775 0.8499000072479259 -1.1549999713897714 0.3849999904632572 0.3849999904632571 -2.5497000217437775 0.8499000072479260 0.8499000072479258 ! 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 结果输出文件EtVo.dat是空的，看了OUTCAR，最终结果只有最后一个晶格参数的计算。请问脚本那里出了问题？ 还有如果不想得到体模量，只想优化，什么时候才有做状态方程拟合晶格参数必要，一般直接做全优化可以吗？还有做状态方程的时候，没有优化自由度，OUTCAR最后FORCE比较大，这个能量也不能说明问题吧 |
回复此楼» 猜你喜欢
材料专硕322
已经有7人回复
-
求调剂
已经有20人回复
-
一志愿郑州大学085600求调剂
已经有20人回复
-
化学调剂
已经有17人回复
-
3-溴吡啶-4-甲醛合成
已经有9人回复
-
一志愿南京航空航天大学 材料与化工329分求调剂
已经有3人回复
-
285求调剂
已经有9人回复
-
材料求调剂
已经有11人回复
-
调剂
已经有8人回复
-
材料调剂
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何求解描述振动的二阶微分方程
已经有21人回复
- 椭圆形偏微分方程标准化后该如何解？
已经有6人回复
- 有关磁流体方程计算的询问
已经有8人回复
- 怎么才能做出三个9的线性方程啊？？！！愁死我了
已经有11人回复
- 求助Ta2O5的色散方程
已经有4人回复
- PRSV方程中的k
已经有4人回复
- 大家看看我的excell拟合方程，怎么能合理的解释相关性呢？
已经有8人回复
- matlab拟合方程参数时初值的选择
已经有15人回复
- R语言图书-求解差分方程和并行计算！
已经有43人回复
- 求助这个非线性方程改怎么解啊？？
已经有4人回复
- 求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧
已经有20人回复
- 【求助】如何计算实验数据与已知方程的相关系数？
已经有6人回复
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
zhang668(金币+3): 2010-04-24 10:16
|
别人的脚本，不要拿过来就用，最好试一试。 你可以找一个算完的OUTCAR试一试，确定脚本没问题了，再进行计算。 EtVo.dat是空的，就说明是E和V的变量没有赋进去值。 我看你用的符号都是' '，应该是一个`,一个是' E='grep ''TOTEN'' OUTCAR | tail -1 | awk '{printf ''%12.6f \n'' , $5}' ' 应该是E=`grep ''TOTEN'' OUTCAR | tail -1 | awk '{printf ''%12.6f \n'' , $5}' ` 拟合物态方程绝对不是为了获得B0， 大部分的目的还是认为压力算的不准，不如拟合物态方程算出的压力准。至于拟合物态方程也紧紧限于立方结构，非立方结构由于加压之后，很可能晶格常数的比例变了，必须全优化后，才能拟合物态方程，也有不拟合的。 至于你说的FORCE很大， 那最大的原因就是你没有relax离子，否则不可能很大，或者你的收敛参数不够高。对应的就是isif=2， nsw=60 立方体系也可以全优化，但是那样比较浪费时间，没有这样的快！ |
2楼2010-04-23 23:35:28
wangzg9292
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3933
- 散金: 56
- 帖子: 348
- 在线: 89.6小时
- 虫号: 723844
- 注册: 2009-03-16
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
3楼2010-04-24 09:05:33
4楼2010-04-24 11:08:41
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
5楼2010-04-24 11:23:40
zhang668(金币+1): 2010-04-24 11:27
|
本帖内容被屏蔽 |
6楼2010-04-24 11:24:40
7楼2010-04-24 11:27:42
8楼2010-04-24 11:29:34
9楼2010-04-24 11:45:33
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
本帖内容被屏蔽 |
10楼2010-04-24 13:11:28
10