24小时热门版块排行榜    

汕头大学海洋科学、生物学、生物与医药等3个专业接受调剂
查看: 6591  |  回复: 12
本帖产生 1 个 1ST强帖 ,点击这里进行查看

xiao72379

木虫 (正式写手)

宅男甲

[求助] 求助vasp做过电荷密度差分的文献以及答疑,谢谢!

Sample Text
大家好  以前我分析成键都用DOS+ELF+Bader  最近的一个工作 我师姐建议做下电荷差分密度(charge density difference) 我没做过 请大家帮忙 谢了先!

1)求用vasp做过电荷密度差分的文献 这个用ISI实在不好查  老板催工作我时间又不多- - 如果同学手上有现成的文献给我参考 感激不尽

2)查了下论坛原来的帖子 找到一个我认为比较好的参考贴[http://muchong.com/bbs/viewthread.php?tid=2129228](我不知道如何给他金币,版主看到了教下我,呵呵),渔火江枫 的 先引用下:


VASP与电荷密度差
   一般分为:差分电荷密度图(deformation charge density),二次差分电荷密度图(difference charge density)
两者的区别有很多种说法,其中一种定义认为:
     差分电荷定义为成键后的电荷密度与对应的点的原子电荷密度之差。通过差分电荷密度的计算和分析,可以清楚地得到在成键和成键电子耦合过程中的电荷移动以及成键极化方向等性质。
     “二次”是指同一个体系化学成分或者几何构型改变之后电荷的重新分布。
如何做对电荷密度做差,目前有两个版本的公式:

人认为公式(1)和差分电荷密度的定义比较吻合,但是这种方法会产生很多干扰信息,而文献中常用的是公式(2)。

如何用VASP做电荷密度差图?
①        对优化后的结构做静态自洽计算
要注意的是,AB、A和B要分别放在同样大小的空间格子中并保证A和B与AB中相应坐标不变,计算时也要保证三次自洽计算所采用的FFT mesh 一致(NGXF,NGYF,NGZF)。
       INCAR中几个注意的参数:
         IBRION = -1;NSW = 0;
         NGXF,NGYF,NGZF
②        做差
    将第一步得到的三个CHGCAR文件依次命名为CHGCAR1、CHGCAR2、CHGCAR3,采用脚本程序按照公式②对三个文件做差,得到CHGdiff文件。
③        视图
     3D显示可以将得到的CHGdiff文件添加vasp后缀(CHGdiff.vasp),用VESTA读取,另外VESTA也可以做2D切面,缺点是没有标度。
     2D的显示也可以用lev00去做切面,然后用origin作图,得到的2D面电荷密度差带有标度便于分析(具体做法可参阅lev00说明书)。

好了  我的问题就是  比如我的体系 基底是A材料 然后呢是多个B原子吸附或者是嵌入(这里你视为化学键合的作用就行),形成了A(B)n,对上面的公式2应用到我的体系如何处理?比如说步骤中的1,2步  B做静态自洽计算的时候 是n个B一起算(我理解是扣除A的所有坐标得n个B的坐标,不知对否)还是一个一个B分开算,如果是分开,那脚本做CHGdiff如何做呢?如果都不是,请教下我该如何做,谢谢!

[ Last edited by xiao72379 on 2011-5-25 at 10:32 ]
回复此楼

» 收录本帖的淘帖专辑推荐

材料计算模拟实用技巧 第一性原理 第一原理资料汇编 VASP
vasp计算资料 第一性原理 VASP相关 学习资料
第一性原理计算-tg

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

天行健君子以自强不息
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)


franch(金币+1): 谢谢回帖交流 2011-05-25 21:37:30
引用回帖:
Originally posted by xiao72379 at 2011-05-25 10:17:31:
Sample Text
大家好  以前我分析成键都用DOS+ELF+Bader  最近的一个工作 我师姐建议做下电荷差分密度(charge density difference) 我没做过 请大家帮忙 谢了先!

1)求用vasp做过电荷密度 ...

差无定法,怎么差要看你要什么
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
2楼2011-05-25 19:46:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)

宅男甲

引用回帖:
Originally posted by cavediger at 2011-05-25 19:46:56:
差无定法,怎么差要看你要什么

要这些原子也就是上面说的B与体系A的成键信息 或者说看电荷转移 目的还是挺明确的 我现在是不知道怎么操作  以前没做过 现在组里也没人教  具体疑问上面已经提出来了 如是说算一个B+A---AB我按帖子应该做的出来  不过N个B该如何算呢 我现在一直想是不是就把优化的CONTCAR里面的A坐标全部抠出来 把N个B一起算了来差  不知对不对
天行健君子以自强不息
3楼2011-05-25 19:58:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-05-26 19:27:45
引用回帖:
Originally posted by xiao72379 at 2011-05-25 19:58:56:
要这些原子也就是上面说的B与体系A的成键信息 或者说看电荷转移 目的还是挺明确的 我现在是不知道怎么操作  以前没做过 现在组里也没人教  具体疑问上面已经提出来了 如是说算一个B+A---AB我按帖子应该做的出来 ...

如果B成团簇,讨论整个团簇与A相互作用,那就把B一起算;如果不是,B比较弥散,且周围环境等同,彼此分开较远,也可以一起算;如果B周围环境不等同,B、A作用也不同,那就一个一个地算;如果考虑B吸附或嵌入后对其他B的影响,先吸附/嵌入上去的B要同A一起作基底,其它B单独算......
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
4楼2011-05-26 05:43:58
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

漂泊四方

金虫 (小有名气)

【答案】应助回帖

★ ★
xiao72379(金币+10): 感谢回帖 2011-05-26 09:11:34
youzhizhe(金币+2, 1ST强帖+1): 谢谢交流。 2011-05-26 19:27:15
引用回帖:
Originally posted by xiao72379 at 2011-05-25 10:17:31:
Sample Text
大家好  以前我分析成键都用DOS+ELF+Bader  最近的一个工作 我师姐建议做下电荷差分密度(charge density difference) 我没做过 请大家帮忙 谢了先!

1)求用vasp做过电荷密度 ...

好了  我的问题就是  比如我的体系 基底是A材料 然后呢是多个B原子吸附或者是嵌入(这里你视为化学键合的作用就行),形成了A(B)n,对上面的公式2应用到我的体系如何处理?比如说步骤中的1,2 步  B做静态自洽计算的时候 是n个B一起算(我理解是扣除A的所有坐标得n个B的坐标,不知对否)还是一个一个B分开算,如果是分开,那脚本做CHGdiff如何做呢?如果都不是,请教下我该如何做,谢谢!

方案一里面:把B当做一个整体,主要考虑的是A整体和B的作用。
方案二里面:实际上考虑了A和B的作用,还考虑了B与B之间的作用。

所以重点看你关注什么,衬底与吸附原子之间的作用还是吸附原子之间的作用。两种作用之间的强弱关系?看你的意思。如果B形成了团簇,方案一可能是你关注的,如果B比较离散,作用力较弱,方案一也可以。如果BB之间作用力介于他们之间,和衬底的作用类似,分析起来就比较吃力了。

操作起来其实还不是太难,记得p4vasp好像可以做这个。
我与我纠缠久
5楼2011-05-26 08:43:26
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)

宅男甲

引用回帖:
Originally posted by 漂泊四方 at 2011-05-26 08:43:26:
好了  我的问题就是  比如我的体系 基底是A材料 然后呢是多个B原子吸附或者是嵌入(这里你视为化学键合的作用就行),形成了A(B)n,对上面的公式2应用到我的体系如何处理?比如说步骤中的1,2 步  B做静态自洽 ...

我的情况是:B在体系A表面不成团簇,B与A体系的作用从bader计算看是有一个|e|左右的电子从B转移到体系A(B是2价主族元素) DOS上看B的S电子与体系A的峰有杂化现象,呈现一定的共价成分,不过B的S电子在费米能级附近有明显裂分现象,一个峰出现在价带,一个峰出现在导带,还是很好的印证了bader电荷分析的结果。所以我初步判断了B与A体系的键合作用是主离子性辅共价性的混合键性。做电荷差分是进一步支持下我的论点。我在计算过程中除了考虑一个B与A的作用,也做了多个B与A的作用。当然一个B与A的作用做电荷差分现在看来我在步骤上已经是清楚了,但当多个的时候,情况应该是A与B是较强的化学作用,BB之间作用较弱,从ELF做的图看,没有化学键,顶多是静电排斥,我主要关心的是B与A的键合情况是不是我上面提到的情况,这个差分做出来也是为了多一个证据说明。写的有点乱,总结下就是:B不成团簇,B较分散,但不同B的化学环境是不同的,因为他们嵌入的位置不一样,不过成键性质是一样的(比如dos+ELF的结果显示定性是类似的,不过是定量上的不一样)。我确定B与A之间有化学键合,我希望做差分进一步研究成键情况,当然结果能匹配我上面提到的 我已经做过的一些成键分析结果是最好的了。如果是B一起做个差分(考虑到相同键性),我还能分析到每个B与A的电荷转移情况?还是该一个一个来?如你做过,能不能给个细致点的操作步骤,如何处理vasp文件进行计算得到差分结果?
天行健君子以自强不息
6楼2011-05-26 09:51:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)

宅男甲

引用回帖:
Originally posted by cavediger at 2011-05-26 05:43:58:
如果B成团簇,讨论整个团簇与A相互作用,那就把B一起算;如果不是,B比较弥散,且周围环境等同,彼此分开较远,也可以一起算;如果B周围环境不等同,B、A作用也不同,那就一个一个地算;如果考虑B吸附或嵌入后 ...

感谢回复,请看下我在6楼的回帖,我说明了下我已经得到的一些结果,在这个基础上能建议下该怎么操作么?若能就问题教下我,感激不尽,呵呵
天行健君子以自强不息
7楼2011-05-26 09:56:52
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

cavediger

金虫 (正式写手)

★ ★
youzhizhe(金币+2): 谢谢交流。 2011-05-26 19:27:59
引用回帖:
Originally posted by xiao72379 at 2011-05-26 09:56:52:
感谢回复,请看下我在6楼的回帖,我说明了下我已经得到的一些结果,在这个基础上能建议下该怎么操作么?若能就问题教下我,感激不尽,呵呵

一个B的时候就不用说了,多个B的时候,如果彼此间没有显著影响的话,一起算也可以,这样省点事儿,当然你要是想逐个算,也未尝不可,反正就是要一个定性的证据......
我们都老得太快,却聪明得太迟。当我一路跌跌撞撞登上山巅,蓦然发现,本来就有一条笔直的路,直达山顶....
8楼2011-05-26 17:34:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

漂泊四方

金虫 (小有名气)

【答案】应助回帖


gzqdyouxia(金币+1): 谢谢提示 2011-05-27 10:41:50
xiao72379(金币+3): 谢谢,体系还是有几十个原子的 不过我后面会都试试 2011-05-31 19:05:14
xiao72379(金币+10): 2011-08-31 09:15:58
引用回帖:
Originally posted by xiao72379 at 2011-05-26 09:56:52:
感谢回复,请看下我在6楼的回帖,我说明了下我已经得到的一些结果,在这个基础上能建议下该怎么操作么?若能就问题教下我,感激不尽,呵呵

其实做计算的很多情况下都不是做之前就能了解结果的。对于你这个问题最简单的办法就是都去做一遍,然后看看那个结果更符合你的要求,因为感觉上判断这两个方法对你的体系来说没有太大的区别。只有几个原子的体系,计算应该不会花很多时间吧?
我与我纠缠久
9楼2011-05-26 22:01:44
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

gloomy2004

金虫 (著名写手)

很有收获, 谢谢!
10楼2011-07-20 03:37:07
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiao72379 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 315求调剂 +3 2777084671 2025-04-03 3/150 2025-04-03 10:39 by 褚迎松小小
[考研] 085600材料与化工310英一求调剂 +17 几顿饭 2025-03-27 17/850 2025-04-03 00:53 by 诗与自由
[考研] 本科河南农业大学,一志愿苏州大学071000生物学289分求调剂 +4 潮潮玫瑰少年 2025-04-02 6/300 2025-04-02 22:47 by 佩琪宝贝
[考研] 接受跨专业调剂 +4 没学上咋整 2025-04-02 5/250 2025-04-02 17:18 by 我爱我家娃
[考研] 一志愿东北大学 材料科学与工程 学硕 323分求调剂 30+3 殊不知02 2025-04-01 7/350 2025-04-02 16:53 by 1104338198
[考研] 271求调剂 +12 yy1232123 2025-03-27 12/600 2025-04-02 15:36 by Cloris-y
[考研] 288求调剂 +7 张yc11 2025-04-02 8/400 2025-04-02 15:34 by Cloris-y
[考研] 求调剂一志愿东华大学材料与化工268 +5 泫弟弟啊 2025-04-01 5/250 2025-04-02 15:29 by nanaliuyun
[考研] 085500数一英一求调剂283 +4 忧郁的番茄! 2025-04-01 4/200 2025-04-02 14:48 by ljs628
[考研] 286求调剂 +11 白橄榄树 2025-04-01 16/800 2025-04-01 20:44 by syjjj0321
[考研] 一志愿河北工业大学294,专业代码085801求调剂 +4 895821003 2025-04-01 11/550 2025-04-01 15:07 by zhshch
[考研] 一志愿大连理工大学求调剂 +3 xtzz2003 2025-04-01 3/150 2025-04-01 10:44 by yujianx
[考研] 招08调剂 +16 laoshidan 2025-03-30 21/1050 2025-04-01 09:42 by 及时行乐fan
[考研] 085404 293分 希望调剂到材料相关专业 +8 Fizzes 2025-03-28 10/500 2025-03-31 23:33 by 褚迎松小小
[考研] 286求调剂 +10 123页 2025-03-31 10/500 2025-03-31 20:44 by Iveryant
[考研] 276求调剂 +12 荔枝荔枝思思 2025-03-30 13/650 2025-03-31 14:09 by 摇曳百合7
[考研] 311求调剂 +7 小张真的想读研 2025-03-28 8/400 2025-03-30 19:25 by zhouxiaoyu
[考研] 材料与化工考研336分求调剂 +4 张天傲121 2025-03-29 9/450 2025-03-30 07:28 by microscopy
[考研] 307求调剂 +4 宋超宇 2025-03-28 4/200 2025-03-29 20:18 by 褚迎松小小
[考研] 调剂收 生态学 资源利用与植物保护 +3 Bebrave_11 2025-03-28 8/400 2025-03-28 22:06 by It’spossible
信息提示
请填处理意见
复试调剂复试调剂
复试调剂二维码获取复试 & 调剂指南!微信扫一扫
商务合作商务合作
商务合作二维码商务合作咨询微信扫一扫(备注:商务合作)
学术干货学术干货
学术干货二维码微信扫一扫了解最新学术干货
论文辅导论文辅导
论文辅导二维码微信扫一扫添加老师,免费咨询
师资招募师资招募