| 查看: 2287 | 回复: 10 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
Silicane静态计算之后利用Bader分析发现两个Si的电荷数居然不等，还相差很大…… 已有1人参与
|
|||
|
前几天计算Silicane（单层Silicene两边用H饱和的结构，Silicene是一种二维材料），發現一個問題…… 结构优化，静态计算之後，进行Bader分析，结果相当出人意料，两个Si的电荷数居然不等！ 要知道，这两个Si是完全等价的…… 分析结果如下： # X Y Z CHARGE MIN DIST ATOMIC VOL -------------------------------------------------------------------------------- 1 2.2286 0.0000 9.5609 3.5469 0.7023 20.9265 2 4.4571 0.0000 10.2849 3.2780 0.7023 19.4710 3 2.2286 0.0000 8.0619 1.5876 0.7280 107.8755 4 4.4571 0.0000 11.7838 1.5876 0.7280 107.8053 -------------------------------------------------------------------------------- VACUUM CHARGE: 0.0000 VACUUM VOLUME: 0.0000 NUMBER OF ELECTRONS: 10.0000 一个Si是3.5469，另一个是3.2780，这两个数值相差的还挺大 开始以为是精度不够，但是加大K点到45X45X1，还是不等…… 尝试加大NGX, NGY, NGZ 还是不行 将原胞扩展成2X2超胞之后，也还是不行，Si原子之间的电子数就是不等…… 不知道有没有前辈遇到过这个问题，还请各位多多赐教…… |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有123人回复
-
今天基金会出结果吗？20260819
已经有13人回复
-
时间戳变了，能看出什么问题？
已经有16人回复
-
filecode，4个jtjc了
已经有13人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有5人回复
-
今天放榜没戏了吧
已经有5人回复
-
今天放榜吗？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
» 本主题相关商家推荐: (我也要在这里推广)
|
呵呵，我计算了硅稀与硅烷(silicene,silicane)异质结，即将硅稀层放在硅烷层的上面。后来我也用了那个bader charge 分析，结果在最后一步bader CHGRG -ref CHGCAR_sum却出现了错误，should be no new maxima in edge refinement。不知道楼上大神遇见过这种情况没有，拜求解惑 发自小木虫Android客户端 |
9楼2016-02-06 22:41:45
2楼2016-01-29 13:08:41
lfhuang
木虫 (著名写手)
- 应助: 39 (小学生)
- 金币: 3643.6
- 散金: 2321
- 红花: 37
- 帖子: 1432
- 在线: 681.5小时
- 虫号: 505824
- 注册: 2008-02-17
- 性别: GG
- 专业: 金属材料表面科学与工程
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
甲斐之虎: 金币+100, ★★★★★最佳答案 2016-02-03 11:27:41
甲斐之虎: 金币+100, ★★★★★最佳答案 2016-02-03 11:27:41
|
这也正常，因为Bader以电荷密度极小处做为原子边界分割点，这种方法其实是不适合共价体系。Si-Si键当中有较大电荷密度，如以下文章图1(b)： http://journals.aps.org/prb/abstract/10.1103/PhysRevB.91.205433 所以Bader方法在分割原子边界时，可定会避开Si-Si中间电荷极大点，而选择靠近其中一个原子的某个地方做为边界，所以会导致两个原子不均等了。 |
» 本帖已获得的红花（最新10朵）
luzihen3楼2016-01-30 06:34:00
4楼2016-02-03 11:35:59
