24小时热门版块排行榜    

查看: 2494  |  回复: 4

金天赐

木虫 (小有名气)

[求助] 差分电荷密度图问题

我是照着张仁琴教授的教程做的,但做出来的图和原文献里的不一样,有懂的或做过的前辈帮帮忙,看看是哪里出错了,多谢指教。
结构基本上和原文献一样,差分电荷图就是用自洽的总电荷密度减去分的电荷密度,分的电荷密度我incar还是用的原来总的incar,里面只加了NGX,NGY,NGZ,POSCAR和POTCAR我也改了,还有什么遗漏吗?

差分电荷密度图问题
原文献.png


差分电荷密度图问题-1
CHGCAR1.jpg


差分电荷密度图问题-2
CHGCAR2.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金天赐

木虫 (小有名气)

后来我在别的帖子里看到分的电荷密度Incar里要将NELM = 0; ICHARG=12,但我这样改了后,分的电荷密度都跑到盒子上下表面了
2楼2017-03-24 09:58:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金天赐

木虫 (小有名气)

引用回帖:
2楼: Originally posted by 金天赐 at 2017-03-24 09:58:36
后来我在别的帖子里看到分的电荷密度Incar里要将NELM = 0; ICHARG=12,但我这样改了后,分的电荷密度都跑到盒子上下表面了

KPOINT我没改,和总的自洽时的k点一样,INCAI里也加了NGX,NGY,NGZ。哎,感觉INCAR里改了ICHARG和NELM还不如没改之前的图好。
有知道原因的前辈吗?求教
差分电荷密度图问题-3
A.jpg


差分电荷密度图问题-4
B.jpg

3楼2017-03-24 10:07:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heyong0104

金虫 (正式写手)

硕士研究生

引用回帖:
2楼: Originally posted by 金天赐 at 2017-03-24 09:58:36
后来我在别的帖子里看到分的电荷密度Incar里要将NELM = 0; ICHARG=12,但我这样改了后,分的电荷密度都跑到盒子上下表面了

楼主,我也在算差分电荷密度,按照张仁琴教授算出来的感觉不太对,你计算的是掺杂的还是没有掺杂的,计算时固定原子位置了吗?
待人要好,做事要专心,少说话,多做事。
4楼2017-03-24 20:53:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

金天赐

木虫 (小有名气)

引用回帖:
4楼: Originally posted by heyong0104 at 2017-03-24 20:53:48
楼主,我也在算差分电荷密度,按照张仁琴教授算出来的感觉不太对,你计算的是掺杂的还是没有掺杂的,计算时固定原子位置了吗?...

没参杂,现在我的问题解决了,我算B原子时是所有的B原子同时算的,这样不同位置的B原子有的会成键,电荷发生转移,而原文献是要求单独孤立的B原子,我又把所有B原子拆成四部分,分别算了一次,结果就和文献上的一致了。
位置是固定的,做自恰不是吗,而且还要保证FFT一致。

发自小木虫Android客户端
5楼2017-03-25 00:29:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 金天赐 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 找导师 +3 yuanjiabao 2026-08-29 4/200 2026-08-30 08:09 by zjb_nudt
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 06:57 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 06:13 by ZPa0EcMwuECS
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 06:02 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 基金不中,共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见