| 查看: 321 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
找到一些相关的精华帖子，希望有用哦~
关于弹性常数分析问题，C44特别大可能是什么原因造成的呢
已经有3人回复
- 在药企里面有生产型的企业，外包性质的企业，药品报批的企业，请教哪一个比较有前途呢
已经有4人回复
- 求助 这种情况我该怎么办呢
已经有5人回复
- 硕士毕业了，去事业单位好呢还是企业单位好呢
已经有25人回复
- 【求助】做聚氨酯时，用来滴加异氰酸酯的恒压漏斗活塞拔不出来了，该怎么办啊？？
已经有12人回复
- 【交流】做洗涤产品为什么一般是阴离子和非离子复配，阳离子主要是干什么的呢
已经有10人回复
- 急问：点“检查保护”后让输密码要输吗？输了别人怎么打开看呢
已经有4人回复
- 【求助/交流】PBS缓冲溶液高温灭菌后怎么变浑浊了呢？
已经有19人回复
- 【原创】发现教授论文有问题怎么办？
已经有11人回复
- 搞催化的在有机反应路线上创新投什么杂志好呢?Communiations之类的
已经有5人回复
- 请问文章并列第一作者是不是很不好？
已经有55人回复
- 【请教】培养晶体得到的是无机盐，该怎么办呢？
已经有9人回复
- 《寻求帮助》在985高校工作不到半个月的我，要不要开始异地恋呢？~~
已经有32人回复
- 急！！土壤中的硝酸盐的来源都有那些个呢？？？尿素施入后怎么变成硝酸盐的呢？？
已经有4人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
