| 查看: 3039 | 回复: 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，其中红色线内均为文献给定值，还请大家帮忙看看是什么问题造成的，是参数设置问题吗？应该怎么解决，谢谢大家！！！  2.png |
回复此楼» 收录本帖的淘帖专辑推荐
 材料相关（第一性vasp等） |
» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
相丽温
新虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 8830.6
- 红花: 1
- 帖子: 1048
- 在线: 144.8小时
- 虫号: 3400328
- 注册: 2014-09-05
- 性别: MM
- 专业: 理论和计算化学
4楼2017-04-12 17:23:29
zhangfrank
金虫 (正式写手)
- 应助: 178 (高中生)
- 金币: 660.2
- 散金: 226
- 红花: 53
- 帖子: 891
- 在线: 344.5小时
- 虫号: 2679423
- 注册: 2013-09-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★
感谢参与，应助指数 +1
相丽温: 金币+5, ★★★很有帮助 2017-04-12 17:00:19
感谢参与，应助指数 +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
相丽温
新虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 8830.6
- 红花: 1
- 帖子: 1048
- 在线: 144.8小时
- 虫号: 3400328
- 注册: 2014-09-05
- 性别: MM
- 专业: 理论和计算化学
3楼2017-04-12 17:12:43
zhangfrank
金虫 (正式写手)
- 应助: 178 (高中生)
- 金币: 660.2
- 散金: 226
- 红花: 53
- 帖子: 891
- 在线: 344.5小时
- 虫号: 2679423
- 注册: 2013-09-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
>>>>我计算和文献中所用的实验值都是： 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
20