24小时热门版块排行榜    

北京石油化工学院2025年硕士研究生招生接受调剂公告
查看: 6602  |  回复: 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

金虫 (正式写手)

★ ★
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的回帖
查看全部 13 个回答

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +3 wzn1234 2025-04-10 6/300 2025-04-11 00:05 by wzn1234
[基金申请] 广东省科技厅TZ计划出公示结果了没? +4 lcmdlut 2025-04-09 6/300 2025-04-10 21:40 by zp519
[考研] 085500数一英一283求调剂 +6 忧郁的番茄! 2025-04-05 7/350 2025-04-10 20:00 by 诗与自由
[考研] 296求调剂 +6 jiyu_111 2025-04-06 6/300 2025-04-10 18:23 by zengyanbo
[考研] 372求调剂 +3 yuiyiy 2025-04-09 3/150 2025-04-10 13:21 by orlando
[考研] 284求调剂 +5 变态小阿姨 2025-04-05 5/250 2025-04-10 12:57 by ajp861027
[考研] 材料调剂 +14 玄觞26 2025-04-09 14/700 2025-04-10 11:04 by a132412
[考研] 材料332求调剂 +5 核动力驴pro版 2025-04-08 5/250 2025-04-09 21:59 by 敲那个钟
[考研] 扬州大学材料专硕(085601)接收调剂 +3 夜斗VGA 2025-04-06 4/200 2025-04-07 23:09 by forever.C
[考研] 0854求调剂 +4 楚阳川木 2025-04-07 5/250 2025-04-07 21:47 by sunqijun
[考研] 课题组包进,实力雄厚,抓紧时间啦! +4 大黄鸭HQR 2025-04-07 5/250 2025-04-07 19:55 by yfxgugg
[考研] 材料科学与工程290求调剂 +27 双马尾痞老板2 2025-04-04 28/1400 2025-04-07 19:17 by 啥南丹你的梦想�
[考研] 一志愿南昌大学293材料工程085601 +11 蜘蛛侠呢 2025-04-04 11/550 2025-04-07 18:40 by liu85623724
[考研] 化工专硕360求调剂,本211,一志愿985 +5 kqatk 2025-04-06 6/300 2025-04-07 16:53 by 17313308663
[考研] 轻工技术与工程 082200 考了数学 263分 单科均已过A区线 +7 Pdz666 2025-04-05 8/400 2025-04-07 16:01 by 杨杰123
[考研] 中国矿业大学 总分:262 专业:工学->土木工程(专硕) +3 墨染hy 2025-04-06 3/150 2025-04-07 12:16 by zllcz
[考研] 求助085602化学工程调剂 304分 +8 妮小汁 2025-04-05 10/500 2025-04-07 11:31 by lijunpoly
[考研] 初始302,一志愿湖南大学专硕,强烈意愿读博 +7 高歌奋进 2025-04-04 7/350 2025-04-07 06:54 by houyaoxu
[考研] 材料与化工085600求收留 +13 cs0106 2025-04-04 14/700 2025-04-05 17:39 by 手札礼记
[考研] 河北大学生物与医药专业招收调剂!! +4 孙小北hea 2025-04-05 7/350 2025-04-05 10:14 by 一只嗷嗷嗷
信息提示
请填处理意见
澶嶈瘯璋冨墏澶嶈瘯璋冨墏
澶嶈瘯璋冨墏浜岀淮鐮�鑾峰彇澶嶈瘯 & 璋冨墏鎸囧崡锛�寰俊鎵竴鎵�
鍟嗗姟鍚堜綔鍟嗗姟鍚堜綔
鍟嗗姟鍚堜綔浜岀淮鐮�鍟嗗姟鍚堜綔鍜ㄨ寰俊鎵竴鎵�(澶囨敞锛氬晢鍔″悎浣�)
瀛︽湳骞茶揣瀛︽湳骞茶揣
瀛︽湳骞茶揣浜岀淮鐮�寰俊鎵竴鎵�浜嗚В鏈€鏂板鏈共璐�
璁烘枃杈呭璁烘枃杈呭
璁烘枃杈呭浜岀淮鐮�寰俊鎵竴鎵�娣诲姞鑰佸笀锛屽厤璐瑰挩璇�
甯堣祫鎷涘嫙甯堣祫鎷涘嫙