版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

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 ]

回复此楼

» 猜你喜欢

售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有4人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有5人回复
评审感受-评审感受-评审感受已经有20人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有4人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有4人回复
售SCI一区T0P文章，我:8.O.5.5.1.O.5.4,科目齐全,可+急已经有3人回复

Life, Love, Laugh.

1楼 2012-04-09 01:02:12

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

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

求助关于荧光共振能量转移的问题？已经有8人回复
用Dmol优化收敛的结构，用vasp不收敛？已经有9人回复
3.22文献解读：三元材料在高能量密度锂离子电池中的应用已经有236人回复
干涉光能量哪里去了？已经有5人回复
求助：castep计算参数设置没有收敛已经有9人回复
用CHEM3D软件PM3理论水平计算分子轨道却无法显示轨道能量已经有5人回复
关于dmol3的收敛性测试已经有12人回复
天然气燃烧的能量利用效率应该怎么计算？已经有4人回复
请教扣电测试异常已经有4人回复
两个原子靠近成键的过程中，系统总能量会不会发生变化？已经有3人回复
XPS中能量如何校准已经有7人回复
计算表面时：Vacuum真空层和slab层厚度收敛性测试，保持超晶胞体积不变吗？？已经有10人回复
【其他】关于收敛测试的问题已经有22人回复
【求助】收敛测试怎么做？？？已经有17人回复
【求助】收敛性测试的问题已经有40人回复

点击这里搜索更多相关资源

科研从小木虫开始，人人为我，我为人人

相关版块跳转我要订阅楼主 qphll 的主题更新