24小时热门版块排行榜    

查看: 8053  |  回复: 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的回帖

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的回帖
信息提示
请填处理意见