24小时热门版块排行榜

返回列表

qphll

金虫 (正式写手)

应助: 18 (小学生)
金币: 2708.7
散金: 2294
红花: 8
帖子: 554
在线: 227.1小时
虫号: 89654
注册: 2005-08-29
性别: GG
专业: 碳素材料与超硬材料

[求助] 收敛测试,　LaFeO3, 能量变化异常

尝试用vasp计算LaFeO3,　在收敛测试的时候发现随着encut的增加,　总能反而升高.

(1) INCAR脚本

m WAVECAR
for  i in 200 250 300 350 400 450 500 550 600
do
cat > INCAR << !
SYSTEM = LF orthorhombic
ENCUT = $i
ISMEAR = 0 ! Gaussian smearing
SIGMA = 0.1
LREAL= Auto
PREC= Accurate
!
echo  "ENCUT= $i eV" ;
time;
./vasp5211 ;

E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'`
echo $i $E >> compare
done

(2) KPOINTS 文件

Automatic generation
0
Monhkorst-Pack
4    4    4
0.0    0.0    0.0

(3) POSCAR 文件

rthorhombic
1.0
      5.5349998474       0.0000000000       0.0000000000
      0.0000000000       7.8880000114       0.0000000000
      0.0000000000       0.0000000000       5.5989999771
12 4 4 ! O Fe La
Direct
   0.270900011       0.034699999       0.229399994
   0.729099989       0.965300024       0.770600021
   0.229099989       0.965300024       0.729399979
   0.770900011       0.034699999       0.270600021
   0.729099989       0.534699976       0.770600021
   0.270900011       0.465299994       0.229399994
   0.770900011       0.465299994       0.270600021
   0.229099989       0.534699976       0.729399979
   0.493999988       0.250000000       0.564400017
   0.506000042       0.750000000       0.435599983
   0.006000012       0.750000000       0.064400017
   0.993999958       0.250000000       0.935599983
   0.000000000       0.000000000       0.000000000
   0.500000000       0.000000000       0.500000000
   0.000000000       0.500000000       0.000000000
   0.500000000       0.500000000       0.500000000
   0.518400013       0.250000000       0.008900000
   0.481599987       0.750000000       0.991100013
   0.981599987       0.750000000       0.508899987
   0.018400013       0.250000000       0.491100013

最后得到的图是这样的:

http://edu.muchong.com/attachment/16/b3/89654_1333900508_628.jpg

同样, 如果修改脚本, 对sigma测试, 最后得到的图中能量是线性的:

SIGMA 脚本

rm WAVECAR
for  i in 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30
do
cat > INCAR << !
SYSTEM = LF orthorhombic
ENCUT = 500
ISMEAR = 0 ! Gaussian smearing
SIGMA = $i
LREAL= Auto
PREC= Accurate
ISTART = 0
ICHARGE = 2
!
echo  "SIGMA= $i eV" ;
time;
./vasp5211 ;

E=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'`
echo $i $E >> sigma.compare
done

sigma测试能量图:

http://edu.muchong.com/attachment/87/87/89654_1333900718_282.jpg

这是怎么回事呢? 可能的问题:
a. INCAR中的参数设置??
b.赝势的选择?? PAW-PBE的O, La, Fe有几个, 我倒没有进行选择
c. 先优化结构, 再单点能计算??  POSCAR直接从文献报导的CIF得到. 我的打算是做完测试以后, 再做晶体的结构优化.

大家有什么建议?

多谢.

ENCUT.test

SIGMA.test

[ Last edited by qphll on 2012-4-9 at 01:03 ]

回复此楼