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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 按文献值优化Fe3O4晶胞后偏小（附有详细参数）
211/3返回列表123下一页
查看: 2238  |  回复: 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的回帖

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下一页
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 工材口青年基金大概什么样能上会？ +12 今晚推荐22 2024-06-20 17/850 2024-06-20 21:50 by 今晚推荐22
[基金申请] 国产期刊影响因子大于12的有20多个个了 +7 babu2015 2024-06-20 7/350 2024-06-20 21:39 by 潇湘之迷
[考博] 一名额难求啊 +6 pinbo拼搏 2024-06-19 6/300 2024-06-20 21:11 by 82年拉菲
[基金申请] 国自然青年基金，1A4B能上会吗？青年和面上的上会标准是一样的吗？ +5 今晚推荐22 2024-06-20 7/350 2024-06-20 19:48 by goodboyycb
[有机交流] 跑板能跑开，过柱过不纯怎么办 +6 小胡在努力 2024-06-18 8/400 2024-06-20 18:54 by HF111001
[论文投稿] 求机械类四区sci推荐 5+4 迷茫小旷 2024-06-14 7/350 2024-06-20 17:42 by 不一样烟火12345
[找工作] 初始合伙人来啦！（生物试剂耗材标准品） +15 欢快的小科研人 2024-06-15 26/1300 2024-06-20 17:26 by 20086126
[教师之家] 试用期辞职 +10 ZHONGWU_U 2024-06-18 17/850 2024-06-20 17:07 by ZHONGWU_U
[精细化工] 试剂价格 +5 旭必上岸 2024-06-17 8/400 2024-06-20 16:25 by 旭必上岸
[硕博家园] 豫北虫友互识 +11 xuhongli903 2024-06-18 12/600 2024-06-20 13:32 by 1429524838
[考博] 有机化学迷茫学生 +4 佛系摸鱼5 2024-06-18 6/300 2024-06-20 10:05 by 295143924
[基金申请] 我标书代码变了 +62 学员NHuqdk 2024-06-16 85/4250 2024-06-19 11:52 by zhangjxnu
[高分子] 寻找聚酯反应釜 +3 茕茕恭煮 2024-06-15 6/300 2024-06-18 14:15 by 茕茕恭煮
[论文投稿] 投稿求助 10+3 2022_灵魂工程师 2024-06-16 8/400 2024-06-18 12:12 by 投必得科研顾问
[高分子] 烧瓶内合成聚酯 +3 大帝国乐 2024-06-17 7/350 2024-06-18 11:10 by 大帝国乐
[论文投稿] 审稿问题：为什么荧光激发波长和紫外吸收波长差的大？ 10+5 sdawege 2024-06-14 10/500 2024-06-17 18:54 by HH-探针
[基金申请] 面青地会评时间 +8 tanjydd 2024-06-15 8/400 2024-06-17 17:08 by 小龙虾2008
[论文投稿] 二审返修送审10天了，原来一审的3个审稿人只有2个接受了审稿，会邀请新审稿人么？ 50+3 huanpo116 2024-06-15 5/250 2024-06-16 10:27 by bobvan
[基金申请] 博士后创新人才支持计划公示 +9 aishida144 2024-06-14 15/750 2024-06-16 09:52 by msjy
[基金申请] 博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
信息提示
关闭
请填处理意见
关闭