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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 差分电荷密度图的疑惑
53/1返回列表
查看: 8988  |  回复: 41
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

paopaotu326

金虫 (正式写手)

[求助] 差分电荷密度图的疑惑 已有6人参与

大家好：
     我是新手，刚学着做差分电荷密度，作出来的图如下所示：
    我做的是掺杂体系，要看掺杂前后的电荷密度变化，我是分别计算出三个CHGCAR后直接在VESTA中相减的到的图片，不知道这个怎么分析啊？只知道红色区域表示电荷密度增大，蓝色减少，形成的键就不知道了，没学过量字化学啊，求高手帮忙下；
    第二个问题就大了，同样的三个CHGCCAR文件，在linix系统中用脚本CHGdiff.pl来依次相减，就是ABC－A，得到一个结果，再用它－C＝chgdiff，放在vesta中显示，才好玩，和直接在vesta里相减得出的图完全不一样，如第二张图，这是怎么回事？
    我的体系是在SnO2中用Pd原子取代Sn原子形成的。

差分电荷密度图的疑惑
CHGCAR_vesta.jpg


差分电荷密度图的疑惑-1
CHGCAR.jpg
回复此楼

» 猜你喜欢
1楼 2014-05-23 08:55:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZHC_Y

金虫 (小有名气)

Dr.
★ ★
franch: 金币+2, 鼓励交流。。。 2014-06-10 09:11:29
引用回帖:
24楼: Originally posted by paopaotu326 at 2014-06-09 08:27:06
嗯嗯，谢谢楼上，我再试试看。对了，您有没有用vesta做差分？奇怪的是，我用vesta做的差分图比较合理，但用chgdiff脚本去做两次减法后，得到的图却很奇怪啊，嘿嘿：）我在想是不是做错了？我是用掺杂后总的chgcar分 ...

我也用过chgdiff做过差分，得到的结果也很不理想，就是整个图看上去和vesta的差不多，但是图中有些锯齿，很不漂亮。用chgdiff脚本去差分的话，我也是按照你的步骤做的。不过，我也用过lev00去做过差分，得到的图就很好了，和vesta基本一样，你可以试试看。我也想问你一个问题，就是我用vesta去差分之后得到的结果，再去切片，但是始终得不到很好的结果。你能告诉我一下，你用的vesta版本和你去切片的方法吗？谢谢了！
赞一下 回复此楼
25楼2014-06-09 10:07:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 42 个回答

paopaotu326

金虫 (正式写手)
补充下第一张图的标尺，最上面是：0.1,最下面是－0。25
赞一下 回复此楼
2楼2014-05-23 08:56:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paopaotu326

金虫 (正式写手)
第三个疑惑，如果想把所获得的CHGCAR文件放在lev00里面处理，最后放在origin中作图，有个很大的问题，就是用lev00程序的时候，一开始就需要输入：outcar文件和contcar文件，之后才输入chgcar文件，那么最后输入的chgcar文件就用三个CHGCAR相减得到的差，可是前面地 outcar 和contcar该输入哪个？是ABC的，还是B的，或者C的？ 我的理解是计算的时候三个东东的空间格子都建立的是一样的，NGX、NGY和NGZ也一样，那么是不是输入ABC的outcar和contcar就行了？
赞一下 回复此楼
3楼2014-05-23 09:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paopaotu326

金虫 (正式写手)
求帮忙啦，不然难毕业啦，吼吼              
赞一下 回复此楼
4楼2014-05-23 09:02:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 42 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯，如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[基金申请] 申请科研项目纯靠实力能成功吗？ 200+5 David1982 2026-10-13 8/400 2026-10-14 21:38 by lhw_GPA4.0
[职场人生] 40+中年人失业后，还能找到好工作吗？ +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了，离AI做有机合成还有多远？ +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么？ +4 好运来～9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷！
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬，路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些？ +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
关闭
请填处理意见
关闭