查看: 2238  |  回复: 20
【悬赏金币】回答本帖问题,作者Tea_Candy将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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
回复此楼
Life is more than Chemistry.
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangguiyang

新虫 (小有名气)

【答案】应助回帖


御剑江湖: 金币+1, 谢谢回帖 2012-06-15 08:56:46
EDIFF  = 1.0e-6


因为你优化体积了,应力需要更高收敛精度
18楼2012-06-15 00:11:21
已阅   回复此楼   关注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的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 怎么带研究生? +11 豆豆小小2 2024-06-20 11/550 2024-06-21 00:18 by smilerobin
[精细化工] 试剂价格 +7 旭必上岸 2024-06-17 10/500 2024-06-20 23:35 by 相与处于陆的鱼
[基金申请] 江苏省333人才工程 出校后被刷的概率大不大? +4 maxbirdzhang 2024-06-19 4/200 2024-06-20 22:38 by sunss
[基金申请] 我标书代码变了 +63 学员NHuqdk 2024-06-16 86/4300 2024-06-20 22:03 by 841790061
[论文投稿] ACS AMI 返回审稿意见,一个大修,两个据稿,编辑给的修改重投 +7 智商已更新 2024-06-19 7/350 2024-06-20 19:54 by kinlin13
[博后之家] 在国内某高校做全职博士后2年,现在找到新的单位,出站或退站对新工作有什么影响? +3 nxplfcc 2024-06-20 3/150 2024-06-20 19:50 by zhangysbad
[论文投稿] 投稿求助 6+3 小短腿安琪拉 2024-06-19 4/200 2024-06-20 17:40 by 不一样烟火12345
[有机交流] 求助 10+5 脂质纳米粒 2024-06-20 6/300 2024-06-20 17:16 by czyzsu
[有机交流] 怎么萃取出锡盐内包裹的化合物 +3 硕六过 2024-06-19 4/200 2024-06-20 15:45 by 2461777831
[教师之家] 每次骚扰女学生的都是院系领导,而不是普通教师,小编们要注意措辞正确 +11 zju2000 2024-06-15 13/650 2024-06-20 15:19 by ase123456
[硕博家园] 豫北虫友互识 +11 xuhongli903 2024-06-18 12/600 2024-06-20 13:32 by 1429524838
[精细化工] 切削液交流群 +5 SZ吴彦祖 2024-06-17 5/250 2024-06-20 12:10 by 毛爱波BBB
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[海外博后] 墨尔本大学博后offer要不要接 +3 kyxblmm 2024-06-18 3/150 2024-06-19 22:39 by blake1111
[基金申请] F口信息学部拿面上,大概需要什么样的成果 +7 _奋黎_ 2024-06-16 18/900 2024-06-19 20:02 by _奋黎_
[论文投稿] 审稿 +5 香瓜木香 2024-06-19 6/300 2024-06-19 17:44 by xli1984
[教师之家] 请问事业编制和年薪制冲突吗? +12 ZHONGWU_U 2024-06-14 12/600 2024-06-18 19:31 by fangyl2005
[考博] 2025考博 +8 自强不息a?a 2024-06-15 13/650 2024-06-18 18:12 by 投必得科研顾问
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[基金申请] 面青地会评时间 +8 tanjydd 2024-06-15 8/400 2024-06-17 17:08 by 小龙虾2008
信息提示
请填处理意见