24小时热门版块排行榜    

查看: 8935  |  回复: 41
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

paopaotu326

金虫 (正式写手)

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

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

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


差分电荷密度图的疑惑-1
CHGCAR.jpg
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Joannaouc

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
franch: 金币+4, 谢谢回帖交流。。。 2014-06-10 09:10:58
总结了下你的问题,回答如下:
1. 电荷差分需要是总的减去每个分的。也就是说ABC-A-B-C。其中A、B、C分别的需要做一下自洽,不需要进行结构优化。这个是针对你刚开始提出的ABC-A-C,我觉得在你的处理过程没有减去B,不知道是不是这样。
2. 如果要分析原子之间成键或者因为掺杂而导致的电荷转移不同,我不推荐用切片方法。用create isosurface可能会更直观一些。我用VESTA的时候直接给出的就是等高面,如果你的没有可以试试objects->volumetric data->show isosurfaces
18楼2014-05-29 15:17:56
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:34 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:18 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-28 23:19 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:45 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:28 by ez6fGg9abYaj
[博后之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-28 22:21 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 17:30 by ZYXdhzDAy9ZX
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +7 zju2000 2026-09-26 8/400 2026-09-28 16:42 by beefly
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 15:32 by ZYXdhzDAy9ZX
信息提示
请填处理意见