24小时热门版块排行榜

返回列表

zhang668

木虫 (正式写手)

1ST强帖: 3
应助: 11 (小学生)
贵宾: 0.061
金币: 1735
散金: 1294
红花: 11
帖子: 843
在线: 324.9小时
虫号: 733242
注册: 2009-03-28
性别: GG
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】做状态方程计算时候得到的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比较大，这个能量也不能说明问题吧

回复此楼

» 猜你喜欢

自然科学基金委宣布启动申请书“瘦身提质”行动已经有4人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有6人回复
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）已经有4人回复
求个博导看看已经有18人回复
自荐读博已经有6人回复
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？已经有3人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复

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

如何求解描述振动的二阶微分方程已经有21人回复
椭圆形偏微分方程标准化后该如何解？已经有6人回复
有关磁流体方程计算的询问已经有8人回复
怎么才能做出三个9的线性方程啊？？！！愁死我了已经有11人回复
求助Ta2O5的色散方程已经有4人回复
PRSV方程中的k 已经有4人回复
大家看看我的excell拟合方程，怎么能合理的解释相关性呢？已经有8人回复
matlab拟合方程参数时初值的选择已经有15人回复
R语言图书-求解差分方程和并行计算！已经有43人回复
求助这个非线性方程改怎么解啊？？已经有4人回复
求教origin拟合简单的渗流方程（又称逾渗方程），做复合材料的应该会遇到吧已经有20人回复
【求助】如何计算实验数据与已知方程的相关系数？已经有6人回复

1楼 2010-04-23 23:10:12

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

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
立方体系也可以全优化，但是那样比较浪费时间，没有这样的快！