24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2803  |  回复: 10

相丽温

新虫 (著名写手)

[求助] DFT-D2优化晶格常数和计算bader电荷分析,结果与文献不符 已有1人参与

如题!!!
       在晶格常数优化时所用的截断能,精度,范德华修正的参数及KPOINTS,POSCAR,POTCAR均和文献一致,可是优化得到的结果还是和文献值有一定差距,计算值比文献值偏小,计算得到的bader电荷比文献值偏小更多,这可能会是由什么原因造成的?

      晶格常数实验值: a=b=3.674             c=5.033
      晶格常数文献值: a=b=3.587             c=4.906
      晶格常数计算值: a=b=3.549             c=4.886
  
     bader电荷分析文献值为0.83,计算值为0.54。
     
     附上晶格常数中的INCAR,其中红色线内均为文献给定值,还请大家帮忙看看是什么问题造成的,是参数设置问题吗?应该怎么解决,谢谢大家!!!

DFT-D2优化晶格常数和计算bader电荷分析,结果与文献不符
2.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
相丽温: 金币+5, ★★★很有帮助 2017-04-12 17:00:19
优化得到的结果还是和文献值有一定差距,计算值比文献值偏小,计算得到的bader电荷比文献值偏小更多
>>>>你说的POSCAR也一样,你确定么?是吸附还是取代体系看bader,杂质原子的位置也完全一样?
>>>>我觉得你的晶格常数部分无影响,1.07%的误差完全没有关系的。
>>>>bader可能就是跟你原子的位置是有关的,比如你的吸附了一个原子 是否吸附到合适位置
>>>>即便bader不一样,但是我们用bader是分析得失电子的,所以看的是delta电荷。你要有俩体系做比较,如果电荷变化不大那就没有关系了。【另外,看趋势骂,如果是得失电子规律是一样的,那么0.3e也没多大区别了】

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

计算顺利点吧
2楼2017-04-12 16:46:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

送红花一朵
引用回帖:
2楼: Originally posted by zhangfrank at 2017-04-12 16:46:15
优化得到的结果还是和文献值有一定差距,计算值比文献值偏小,计算得到的bader电荷比文献值偏小更多
>>>>你说的POSCAR也一样,你确定么?是吸附还是取代体系看bader,杂质原子的位置也完全一样?
> ...

谢谢你的回复!
>>>>我计算和文献中所用的实验值都是: a=b=3.674             c=5.033;
>>>>虽然晶格常数误差不是很大,但是bader计算太离谱,只好回过头来看是不是优化的问题了;
>>>>我计算的bader是体结构的,因为差太多就没继续计算表面和吸附体系的;
>>>>bader不止差0.3,理论上每个原子是1e,文献是0.83,我计算的只有0.54;
我尝试计算下表面的bader电荷,比较下试试!
3楼2017-04-12 17:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
2楼: Originally posted by zhangfrank at 2017-04-12 16:46:15
优化得到的结果还是和文献值有一定差距,计算值比文献值偏小,计算得到的bader电荷比文献值偏小更多
>>>>你说的POSCAR也一样,你确定么?是吸附还是取代体系看bader,杂质原子的位置也完全一样?
> ...

看错了,bader是差0.3,理论是1,计算得到0.54是不是差太多,还是说只要后续计算的电荷转移趋势是正确的就OK。
4楼2017-04-12 17:23:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

引用回帖:
3楼: Originally posted by 相丽温 at 2017-04-12 17:12:43
谢谢你的回复!
>>>>我计算和文献中所用的实验值都是: a=b=3.674             c=5.033;
>>>>虽然晶格常数误差不是很大,但是bader计算太离谱,只好回过头来看是不是优化的问题了;
& ...

>>>>我计算和文献中所用的实验值都是: a=b=3.674             c=5.033;
>>>>虽然晶格常数误差不是很大,但是bader计算太离谱,只好回过头来看是不是优化的问题了;
>>>>我计算的bader是体结构的,因为差太多就没继续计算表面和吸附体系的;
>>>>bader不止差0.3,理论上每个原子是1e,文献是0.83,我计算的只有0.54;


>>>>你的意思是,你现在在算的是理想块体的(没有任何缺陷参杂的)bader?那按理说结构不会有什么大错误。我也不好解释为什么差了0.3e了。 你说的理论上每个原子都是1e这个也只是理论上。我记得我以前也看过理想块体的bader,理论上是1的也是仅仅0.几。毕竟任何化学键都是既有共价键成分也有离子键成分。这点很好理解。
所以我一直觉得bader的绝对值意义不大,用来看电荷变化推测电荷转移就ok了
计算顺利点吧
5楼2017-04-12 18:26:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 相丽温 at 2017-04-12 17:23:29
看错了,bader是差0.3,理论是1,计算得到0.54是不是差太多,还是说只要后续计算的电荷转移趋势是正确的就OK。...

我觉得后续趋势正确就可以。你先全算完再比较看吧 可以后续再交流
计算顺利点吧
6楼2017-04-12 18:27:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

引用回帖:
6楼: Originally posted by zhangfrank at 2017-04-12 18:27:35
我觉得后续趋势正确就可以。你先全算完再比较看吧 可以后续再交流...

好的好的  我先算完看看趋势  有问题再向你请教
7楼2017-04-12 18:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

-0.01的收敛限似乎不够,恐怕还不如默认的能量收敛限
POTIM有点大了
ISIF怎么回事?要变胞吗?通常优化晶格参数可以设一系列变值,然后拟合;也可以去ISIF=7直接优化,条件是不改变晶胞对称性。
K点也会影响到优化出来的晶格参数

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

8楼2017-04-13 08:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

相丽温

新虫 (著名写手)

送红花一朵
引用回帖:
8楼: Originally posted by hakuna at 2017-04-13 08:24:33
-0.01的收敛限似乎不够,恐怕还不如默认的能量收敛限
POTIM有点大了
ISIF怎么回事?要变胞吗?通常优化晶格参数可以设一系列变值,然后拟合;也可以去ISIF=7直接优化,条件是不改变晶胞对称性。
K点也会影响到优 ...

谢谢你的回复!
首先收敛标准和K值的设定都是参考文献的;另外我所给出的INCAR是在对体结构晶格参数优化时的INCAR,所以设置ISIF=3;同时由于FeS并不是金属,因此修改ISMEAR=0,POTIM=0.2,可是结果依然如此,最后查看计算结束时的OUTCAR中也有相关DFT-D2的信息,说明也不是范德华修正参数的问题,再次尝试计算bader电荷仍旧在0.5-0.6范围内。
我想请教下原理上每个原子的bader电荷为1,但是我计算的为0.5几。同时我计算体结构和表面以及吸附体系的bader电荷变化趋势正确,是否可以说明问题???
9楼2017-04-13 08:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

引用回帖:
9楼: Originally posted by 相丽温 at 2017-04-13 08:57:44
谢谢你的回复!
首先收敛标准和K值的设定都是参考文献的;另外我所给出的INCAR是在对体结构晶格参数优化时的INCAR,所以设置ISIF=3;同时由于FeS并不是金属,因此修改ISMEAR=0,POTIM=0.2,可是结果依然如此,最后 ...

不一定对,仅供参考
1.Fe磁性,需要考虑
2.vDW不必写那么多参数,都默认的,查你的VASP版本,设一下IVDW即可
3.如果你那个INCAR是优化晶格参数,ISIF取值7试一下

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

10楼2017-04-14 08:49:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 相丽温 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26应届毕业生考博求助 +3 wo一定上岸 2026-05-13 3/150 2026-05-14 21:47 by 明海天涯
[基金申请] 重磅!青年科学基金项目(C类)资助增幅预计超过50% +5 水和泥不是水泥 2026-05-13 7/350 2026-05-14 20:57 by 水和泥不是水泥
[教师之家] 教学课件你会给同学吗 +6 硕士研究生吗 2026-05-13 6/300 2026-05-14 20:44 by pin4501
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 5/250 2026-05-14 20:15 by 一切都是空工
[高分子] 本人最近太闲了,谁有问题可以提,每天会统一回复 +8 一切都是空工 2026-05-12 19/950 2026-05-14 20:03 by 一切都是空工
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +7 dd921ww 2026-05-12 8/400 2026-05-14 14:22 by dd921ww
[考博] 材料类只有一篇综述能申博么 +4 乐逍遥谷 2026-05-13 4/200 2026-05-14 12:05 by zhyzzh
[文学芳草园] 风把牡丹吹跑了 +3 myrtle 2026-05-12 3/150 2026-05-14 10:09 by xingzhewozhi
[基金申请] 请问大佬b0816评完了吗 +3 市民华南虎 2026-05-12 7/350 2026-05-14 07:41 by 市民华南虎
[基金申请] 精华III评审感受-评审感受-评审感受 +12 ferrarichen 2026-05-11 16/800 2026-05-14 07:33 by 2000zf36392
[基金申请] 青C资助名额大幅增加! +8 西葫芦炒鸡蛋 2026-05-13 12/600 2026-05-13 21:42 by 虫眼000
[论文投稿] 有带发论文的吗 +3 山楂之术 2026-05-09 3/150 2026-05-13 17:56 by Cyhcl2629
[硕博家园] 导师各种操作恶心咋办 +11 苍白的小青天 2026-05-09 13/650 2026-05-13 17:11 by 六两废铜
[论文投稿] 护理论文 晋升 +5 Taylor1990, 2026-05-08 5/250 2026-05-13 14:40 by tegsgjy20
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
[考博] 西南大学考核制博士 +3 lijunjie84 2026-05-11 6/300 2026-05-12 18:09 by lijunjie84
[文学芳草园] 窗边初夏的小雨 +7 阿美_Lml888 2026-05-09 10/500 2026-05-12 15:27 by 阿美_Lml888
[考博] 现在不知道怎么办,感觉很痛苦 +4 qweww 2026-05-11 5/250 2026-05-11 20:23 by Oversize
[考博] 生物学博士 +3 17749024330 2026-05-08 6/300 2026-05-11 14:29 by 17749024330
[考博] 2026年申博-电池方向 +7 沃克尔维特 2026-05-08 11/550 2026-05-10 13:56 by 沃克尔维特
信息提示
请填处理意见