| 查看: 303 | 回复: 0 | ||
Gina88木虫 (正式写手)
|
[求助]
总能---长度 测试不是抛物线形状怎么办呢？
|
|
下面贴出了计算Si纳米线的INCAR和产生POSCAR的文件。 计算得到了一个能量-长度图。我有如下问题： 1. 按照网上所说，这个图应该是抛物线的，这个二次曲线系数则和杨氏模量相关。 可是我计算出来不是抛物线的，这是怎么回事啊？ 2. 还有，测试的时候，我的长度方向的伸长和压缩量应该选多少啊。 我选的是15%。是不是这个百分比太大，所以导致不是抛物线形的啊？ 3. 我这里将ISTART设置成了1，不知道这个设置有没有问题？ 看手册上说，要是开始计算时没有WAVECAR，计算时，程序会自动将ISTART当成0啊，后面的计算有WAVECAR，则ISTART=1. INCAR内容： SYSTEM=Si111_relax ISTART=1;ICHARG=2 PREC=Accurate ENCUT=300 GGA=91;VOSKOWN=1 ISMEAR=0;SIGMA=0.08 EDIFF=1E-4 EDIFFG=-0.08 NELMIN=6 LREAL=A ALGO=Fast ISYM=2 ISIF=2 NSW=200 IBRION=2 POTIM=0.5 !LWAVE=.FALSE. LCHARG=.FALSE. 产生POSCAR的文件： rm WAVECAR >& /dev/null for i in 0.85 0.88 0.91 0.94 0.97 1.00 1.03 1.06 1.09 1.12 1.15 do c=`echo "9.405986241642619*$i" |bc` cat >POSCAR< Si1112_1cell 1 22.680000000000000 -0.000000000000000 -0.000000000000000 -11.340000000000000 19.641456157831101 0.000000000000000 -0.000000000000019 0.000000000000010 $c 38 Direct 0.387122727907410 0.274245455814824 0.546291382498555 0.330684091861117 0.330684091861119 0.212948761497599 0.387122727907409 0.274245455814822 0.796298348249276 0.330684091861118 0.330684091861118 0.462955727248317 0.499999999999997 0.330684091861118 0.212948761497597 0.499999999999998 0.330684091861117 0.462955727248317 0.274245455814827 0.387122727907411 0.879634003499515 0.387122727907413 0.443561363953705 0.546291382498562 0.330684091861121 0.499999999999999 0.212948761497605 0.274245455814826 0.387122727907412 0.129613106247365 0.387122727907413 0.443561363953704 0.796298348249276 0.330684091861122 0.499999999999999 0.462955727248324 0.443561363953704 0.387122727907410 0.879634003499515 0.556438636046290 0.443561363953704 0.546291382498556 0.500000000000000 0.500000000000000 0.212948761497604 0.443561363953704 0.387122727907410 0.129613106247357 0.556438636046290 0.443561363953703 0.796298348249276 0.500000000000001 0.500000000000000 0.462955727248324 0.612877272092583 0.387122727907411 0.879634003499515 0.669315908138876 0.499999999999998 0.212948761497597 0.612877272092582 0.387122727907411 0.129613106247357 0.669315908138877 0.499999999999998 0.462955727248316 0.387122727907413 0.612877272092585 0.546291382498564 0.387122727907413 0.612877272092583 0.796298348249276 0.443561363953708 0.556438636046291 0.879634003499515 0.556438636046294 0.612877272092585 0.546291382498562 0.500000000000001 0.669315908138879 0.212948761497604 0.443561363953708 0.556438636046292 0.129613106247364 0.556438636046293 0.612877272092584 0.796298348249276 0.500000000000001 0.669315908138880 0.462955727248323 0.612877272092587 0.556438636046290 0.879634003499515 0.725754544185174 0.612877272092586 0.546291382498563 0.669315908138879 0.669315908138879 0.212948761497604 0.612877272092586 0.556438636046292 0.129613106247364 0.725754544185174 0.612877272092585 0.796298348249276 0.669315908138880 0.669315908138879 0.462955727248323 0.612877272092587 0.725754544185170 0.879634003499515 0.612877272092586 0.725754544185172 0.129613106247364 ! echo "c=9.405986241642619*$i angstrom" time mpirun -np $NP -machinefile $CURDIR/nodelist $HOME/vasp.5.2.11-intel9.1/vasp.5.2.11 en=`grep "entropy=" OUTCAR|tail -1|awk '{printf "%12.6f \n",$7 }'` 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 "$c $en" >>commentVE cp CONTCAR CONTCAR$i done     [ Last edited by Gina88 on 2011-11-10 at 15:48 ] |
回复此楼» 猜你喜欢
有时候，自然基金真的不能太认真 （我的申报经验）
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有68人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
关于弹性常数分析问题，C44特别大可能是什么原因造成的呢
已经有3人回复
- 在药企里面有生产型的企业，外包性质的企业，药品报批的企业，请教哪一个比较有前途呢
已经有4人回复
- 求助 这种情况我该怎么办呢
已经有5人回复
- 硕士毕业了，去事业单位好呢还是企业单位好呢
已经有25人回复
- 【求助】做聚氨酯时，用来滴加异氰酸酯的恒压漏斗活塞拔不出来了，该怎么办啊？？
已经有12人回复
- 【交流】做洗涤产品为什么一般是阴离子和非离子复配，阳离子主要是干什么的呢
已经有10人回复
- 急问：点“检查保护”后让输密码要输吗？输了别人怎么打开看呢
已经有4人回复
- 【求助/交流】PBS缓冲溶液高温灭菌后怎么变浑浊了呢？
已经有19人回复
- 【原创】发现教授论文有问题怎么办？
已经有11人回复
- 搞催化的在有机反应路线上创新投什么杂志好呢?Communiations之类的
已经有5人回复
- 请问文章并列第一作者是不是很不好？
已经有55人回复
- 【请教】培养晶体得到的是无机盐，该怎么办呢？
已经有9人回复
- 《寻求帮助》在985高校工作不到半个月的我，要不要开始异地恋呢？~~
已经有32人回复
- 急！！土壤中的硝酸盐的来源都有那些个呢？？？尿素施入后怎么变成硝酸盐的呢？？
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
