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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 按文献值优化Fe3O4晶胞后偏小（附有详细参数）
101/1返回列表
查看: 2242  |  回复: 20
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者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
回复此楼

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

DFT计算 建模

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
Life is more than Chemistry.
1楼2012-05-17 10:16:40
已阅   回复此楼   关注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的回帖

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

金虫 (小有名气)
引用回帖:
9楼: Originally posted by chenyuchen at 2012-05-18 10:24:31:
有个疑问，文献说GGA用的是PW91，GGA+U没说GGA用的是哪一种，难道默认和前面一样吗？

我把文献上传了，在9楼，一起研究研究吧
赞一下 回复此楼
Life is more than Chemistry.
11楼2012-05-18 13:22:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
引用回帖:
9楼: Originally posted by chenyuchen at 2012-05-18 10:24:31:
有个疑问，文献说GGA用的是PW91，GGA+U没说GGA用的是哪一种，难道默认和前面一样吗？

是在10楼  
回复此楼
Life is more than Chemistry.
12楼2012-05-18 13:23:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
引用回帖:
14楼: Originally posted by liu108 at 2012-05-18 14:24:38:
杩欎釜璇樊鏄?互鐨勶紝LDA+U鐨勭畻娉曚紭浜嶨GA+U锛屼负浠�涔堜�?鐢↙DA+U璁＄畻鍛€��

怎么是乱码呢？？
赞一下 回复此楼
Life is more than Chemistry.
15楼2012-05-18 18:14:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)
引用回帖:
16楼: Originally posted by 344526558 at 2012-06-08 09:56:19
您好！我想问一下 您有关于Fe3O4第一性原理方面的文献吗？我怎么没找到啊？我想计算一下。但是找不到相关参数的设置 谢谢您啦！由于刚刚接触这个方面 所以不是很了解 麻烦您啦！

我在10楼传了一篇哈
赞一下 回复此楼
Life is more than Chemistry.
17楼2012-06-11 15:51:56
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Tea_Candy 的主题更新
101/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 自然基金变化？？ +4 今晚推荐22 2024-06-21 4/200 2024-06-21 11:40 by jurkat.1640
[基金申请] 国自然青年基金，1A4B能上会吗？青年和面上的上会标准是一样的吗？ +11 今晚推荐22 2024-06-20 15/750 2024-06-21 11:36 by jurkat.1640
[基金申请] 听大佬说今年信息口本子数量大幅增加？ +5 wutzxt 2024-06-21 5/250 2024-06-21 11:12 by 1501311
[硕博家园] 难道我真的要放弃吗？ +22 133456 2024-06-20 22/1100 2024-06-21 11:10 by 记忆碎片2018
[基金申请] 青年和面上，哪个上会难度更大 +7 今晚推荐22 2024-06-21 7/350 2024-06-21 11:00 by gilonacv
[论文投稿] Scientific Reports二审90天了，催稿也没反应怎么办？ +4 智晨君 2024-06-18 7/350 2024-06-21 09:57 by 改日再写
[有机交流] 锌粉还原反应 20+3 明宝栋 2024-06-17 5/250 2024-06-21 08:26 by hptianyan
[考博] 有机化学迷茫学生 +4 佛系摸鱼5 2024-06-18 7/350 2024-06-20 23:56 by 佛系摸鱼5
[基金申请] 我标书代码变了 +63 学员NHuqdk 2024-06-16 86/4300 2024-06-20 22:03 by 841790061
[找工作] 高校两个offer选择 +15 cowox2021 2024-06-18 16/800 2024-06-20 15:24 by ase123456
[基金申请] 面青地会评时间？？？ +7 Axvdvbfs 2024-06-19 8/400 2024-06-20 11:16 by 路遥还有谁
[基金申请] F03青年基金函评结果 +3 暨阳一只柴 2024-06-19 4/200 2024-06-20 11:15 by 暨阳一只柴
[基金申请] 太卷了 +14 laoyuefubio 2024-06-17 27/1350 2024-06-20 09:52 by htjwqy
[催化] 镍负载氧化铝的保存问题 8+3 lwn0130 2024-06-15 6/300 2024-06-20 09:00 by lwn0130
[基金申请] 希望今年自己国自然面上项目和老婆青年项目能中！ +10 恐龙爸爸 2024-06-14 10/500 2024-06-20 08:42 by jsguo4127
[高分子] 寻找聚酯反应釜 +3 茕茕恭煮 2024-06-15 6/300 2024-06-18 14:15 by 茕茕恭煮
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[有机交流] 求教最后一步问题出在哪 20+3 苏州浩恒 2024-06-15 4/200 2024-06-18 08:33 by like_2005
[论文投稿] 二审返修送审10天了，原来一审的3个审稿人只有2个接受了审稿，会邀请新审稿人么？ 50+3 huanpo116 2024-06-15 5/250 2024-06-16 10:27 by bobvan
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
信息提示
关闭
请填处理意见
关闭