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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 按文献值优化Fe3O4晶胞后偏小（附有详细参数）
211/3返回列表123下一页
查看: 2896  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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.
1楼 2012-05-17 10:16:40
已阅   回复此楼   关注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朵）

Tea_Candy
会当凌绝顶，一览众山小。
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 的主题更新
211/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 生物学308求调剂（一志愿华东师大） +5 相信必会光芒万� 2026-04-10 5/250 2026-04-10 21:36 by zhouxiaoyu
[考研] 22408 366分，本科211，一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 一志愿鲁东大学071000生物学学硕初试分数276求调剂 +20 慕绝cc 2026-04-09 24/1200 2026-04-10 10:22 by xujun0624
[考研] 调剂 +8 只叙离别辞 2026-04-09 10/500 2026-04-09 20:25 by zl8213662
[考研] 一志愿电子科技大学085600材料与化工 329分求调剂 +14 Naiko 2026-04-04 14/700 2026-04-09 16:56 by luoyongfeng
[考研] 085501机械专硕 302分 不挑专业求调剂 +5 汪某. 2026-04-09 5/250 2026-04-09 15:38 by 蒋皓禹
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 材料调剂 +14 一样YWY 2026-04-05 15/750 2026-04-09 13:36 by 故人??
[考研] 086000生物与医药调剂 +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by 北极159263
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
[考研] 293分求调剂，外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 调剂求助（生物与医药） +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
[考研] 机械专硕274求调剂，不挑专业学校 +6 泛泛2333 2026-04-05 8/400 2026-04-06 18:06 by 泛泛2333
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
信息提示
关闭
请填处理意见
关闭