24小时热门版块排行榜

返回列表

Gina88

木虫 (正式写手)

1ST强帖: 1
应助: 9 (幼儿园)
金币: 3328.8
散金: 124
红花: 8
帖子: 395
在线: 113小时
虫号: 1047684
注册: 2010-06-26
专业: 理论和计算化学

[求助] 总能---长度测试不是抛物线形状怎么办呢？

下面贴出了计算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 ]

回复此楼