24小时热门版块排行榜    

查看: 2854  |  回复: 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的回帖

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的回帖

Tea_Candy

金虫 (小有名气)

引用回帖:
3楼: Originally posted by gougou531 at 2012-05-17 16:36:08:
能量设置太低了,可能收敛到不正确的基态,我觉得和这个有很大的关系,还有就是你设置u的方式不对,同意2楼的说法。如果可以的话可以考虑以力为收敛标准,设置EDIFFG这个参数,一般设置为-0.01

谢谢回复。
我试试看
Life is more than Chemistry.
5楼2012-05-17 20:21:17
已阅   回复此楼   关注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中有两种类型原子

按照这个参数优化后, 结果又偏大了一点8.5333。实验值是~8.39,文献值GGA+U优化时8.473。请问我这个结果能接受吗?
Life is more than Chemistry.
6楼2012-05-17 23:07:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu108

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
WDD880227: 金币+1, 感谢交流提示 2012-05-18 08:56:03
计算值和实验值只差不超过5%是可以接受的,另外,没见过91+U计算的,只见过PBE+U计算的。

» 本帖已获得的红花(最新10朵)

会当凌绝顶,一览众山小。
7楼2012-05-18 08:39:52
已阅   回复此楼   关注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的回帖

chenyuchen

木虫 (正式写手)

【答案】应助回帖


fzx2008: 金币+1, 鼓励交流 2012-05-18 12:08:14
有个疑问,文献说GGA用的是PW91,GGA+U没说GGA用的是哪一种,难道默认和前面一样吗?
选择科学就要耐得住寂寞!
9楼2012-05-18 10:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)

送鲜花一朵
引用回帖:
7楼: Originally posted by liu108 at 2012-05-18 08:39:52:
计算值和实验值只差不超过5%是可以接受的,另外,没见过91+U计算的,只见过PBE+U计算的。

谢谢您的回复,那我计算的误差为~1.6%就是可以的了。
关于那个交换关联势的选择,我看文献中就是用的GGA-PW91,没提到PBE,我把文献上传上来,一起研究研究吧

» 本帖附件资源列表

Life is more than Chemistry.
10楼2012-05-18 13:21:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Tea_Candy 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料调剂 +4 爱擦汗的可乐冰 2026-02-28 4/200 2026-03-01 00:38 by 猫猫球alter
[考研] 272求调剂 +3 材紫有化 2026-02-28 3/150 2026-02-28 22:52 by ms629
[考研] 化工专硕348,一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 292求调剂 +3 yhk_819 2026-02-28 3/150 2026-02-28 21:57 by gaoxiaoniuma
[考研] 290求调剂 +5 材料专硕调剂; 2026-02-28 6/300 2026-02-28 21:40 by gaoxiaoniuma
[考研] 295求调剂 +5 19171856320 2026-02-28 5/250 2026-02-28 21:39 by gaoxiaoniuma
[考博] 26申博 +4 想申博! 2026-02-26 4/200 2026-02-28 21:37 by limorning
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 284求调剂 +4 天下熯 2026-02-28 4/200 2026-02-28 21:13 by gaoxiaoniuma
[考研] 085600材料工程一志愿中科大总分312求调剂 +8 吃宵夜1 2026-02-28 10/500 2026-02-28 20:27 by L135790
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[考博] 博士自荐 +3 kkluvs 2026-02-28 3/150 2026-02-28 16:59 by StarAura
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[基金申请] 面上可以超过30页吧? +12 阿拉贡aragon 2026-02-22 13/650 2026-02-26 22:09 by Hahaxia
信息提示
请填处理意见