24小时热门版块排行榜    

查看: 2777  |  回复: 20
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Tea_Candy

金虫 (小有名气)

[求助] 按文献值优化Fe3O4晶胞后偏小(附有详细参数)

按文献上的参数设置后,优化结果差别比较大:
文献上优化后的晶胞参数是 8.473,我的是8.166
我先后尝试过直接用ISIF=3和ISIF=2(然后用状态方程拟合)优化,以及不加U的优化,但结果都在8.16左右。
不知道是哪儿设置错了,恳请大家指教!
文献参数:


这是我设置的参数
INCAR:
ISTART = 0
   ICHARG = 2
   ISPIN = 2

   IBRION = 2
   ISIF   = 3
   NSW   = 100

   PREC   = Medium
   ENCUT  = 500
   NELM   = 100
   NELMIN = 2
   EDIFF  = 1.0e-3
   ISMEAR =  0
   SIGMA  = 0.1
   LWAVE = .FALSE.
   LCHARG = .FALSE.
   IALGO  = 48

   LREAL = AUTO
   GGA = 91
   VOSKOWN = 1
   LMAXMIX = 4
  
   LDAU = .TRUE.
   LDAUL = 2
   LDAUU = 4.5
   LDAUJ = 0.89

POSCAR:
Title
1.0
   8.3941000    0.0000000    0.0000000
   0.0000000    8.3941000    0.0000000
   0.0000000    0.0000000    8.3941000
32 24
Direct
0.37990 0.37990 0.37990
0.37990 0.87990 0.87990
0.87990 0.37990 0.87990
0.87990 0.87990 0.37990
0.62010 0.12010 0.87990
0.62010 0.62010 0.37990
0.12010 0.12010 0.37990
0.12010 0.87990 0.62010
0.12010 0.37990 0.12010
......

KPOINTS:
Auto
0
M
6 6 6
0 0 0

赝势文件用的是PAW_GGA
回复此楼

» 收录本帖的淘帖专辑推荐

DFT计算 建模

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

Life is more than Chemistry.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenyuchen

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
fzx2008: 金币+1, 鼓励交流 2012-05-18 12:08:20
引用回帖:
7楼: Originally posted by liu108 at 2012-05-18 08:39:52:
计算值和实验值只差不超过5%是可以接受的,另外,没见过91+U计算的,只见过PBE+U计算的。

同意这个。
lz计算的和文献报道的差别应该是交换关联势不同造成的。
选择科学就要耐得住寂寞!
8楼2012-05-18 10:20:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 21 个回答

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢指教 2012-05-17 17:16:40
1.EDIFF  = 1.0e-3有点小,至少应该-4次方
2.GGA+U的设置,这样写应该会报错,因为POSCAR中有两种类型原子
2楼2012-05-17 15:39:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)


感谢参与,应助指数 +1
liliangfang: 金币+1, 谢谢指教 2012-05-17 17:16:48
本帖内容被屏蔽

3楼2012-05-17 16:36:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)

引用回帖:
2楼: Originally posted by fzx2008 at 2012-05-17 15:39:24:
1.EDIFF  = 1.0e-3有点小,至少应该-4次方
2.GGA+U的设置,这样写应该会报错,因为POSCAR中有两种类型原子

谢谢回复。
请问我想把U加在第二个原子上,LDA+U这样设置对么:

   LDAU = .TRUE.
   LDAUL = -1 2
   LDAUU = 0.00 4.5
   LDAUJ = 0.00 0.89
Life is more than Chemistry.
4楼2012-05-17 20:20:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见