24小时热门版块排行榜    

查看: 2096  |  回复: 5

zhangguangping

木虫 (著名写手)

[交流] 用VESTA做VASP的差分电荷密度图已有5人参与

转自:http://muchong.com/bbs/viewthread.php?tid=9046086&fpage=1
本帖子的内容来自于上面帖子中的附件。为了方便大家比较直接的看到,而不用下载附件,故重复发帖。@Jiansh

用VESTA 做VASP差分电荷(C-A+B)图方法:
1、导入 C 的CHGCAR(拖拽到vesta)面板框中
2、Edit——>Edit data——>Volumetric Data————>Import (选择A 的CHGCAR)——>Subtract
from current data——>OK
3、Edit——>Edit data——>Volumetric Data————>Import (选择B 的CHGCAR)——>Subtract
from current data——>OK

注意:A,B,C三个文件的格点要一一对应。
(后来发现,http://muchong.com/bbs/viewthread.php?tid=4287522&fpage=1 这个帖子也不错,很详细 @沙风)
===============
转自:http://muchong.com/bbs/viewthread.php?tid=7392808&fpage=1中的 “电荷密度差的计算.doc” @gulangyu

本文中提到电荷密度差,有时也称为差分电荷密度。
这里主要是讨论在VASP中如何得到用来分析成键前后电荷转移的电荷密度差。

此时电荷密度差定义为:delta_RHO = RHO_sc - RHO_atom
其中 RHO_sc 为自洽的面电荷密度,而 RHO_atom 为相应的非自洽的面电荷密度,是由理想的原子周围电荷分布堆彻得到的,即为原子电荷密度的叠加(a superposition of atomic charge densities)。

不同晶面的 RHO_sc 可由自洽计算的CHG或CHGCAR得到;
而计算 RHO_atom 所需的 CHG 或 CHGCAR 可由下述非自洽计算得到:
仍使用原来自洽计算时的四个输入文件,但INCAR中需要设置 ICHARG=12 和 NELM=0,其他设置不变。

需要特别注意的,应保持前后两次计算(自洽和非自洽)中的 FFT-mesh 一致。因为,只有维数一样,我们才能对两个RHO作相应的矩阵相减。不过,只要按上一段提到的设置方法做就行了,无须特别增加别的设置。

获得了RHO_sc和RHO_atom的CHG或者CHGCAR文件之后,可以利用开始提到的办法,用VESTA很容易做出电荷密度差。
回复此楼

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

综述专辑 wly的好资源专辑

» 猜你喜欢

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

弘德明志博学笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gowoaini131

新虫 (小有名气)

2楼2015-12-26 23:58:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
3楼2015-12-27 13:39:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jimsmart

捐助贵宾 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
非常感谢你的奉献精神
4楼2015-12-28 09:10:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

obaica

禁言 (著名写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

5楼2015-12-28 21:03:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无知!无畏?

木虫 (正式写手)

我怕哪天,现实饿了,就把梦想给吃了。
6楼2022-02-12 14:56:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhangguangping 的主题更新
普通表情 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 【45岁以上博士】柔性人才引进项目,有补贴 +6 Dreamsummit 2024-06-24 9/450 2024-06-25 14:17 by qsx04152006
[基金申请] 专科院校申请国自然 +12 孤单12站 2024-06-24 17/850 2024-06-25 13:14 by wangting_nju
[硕博家园] 联培博士文章第一单位署名问题交流~ +7 橙成成c 2024-06-23 18/900 2024-06-25 12:04 by 橙成成c
[考博] 对象没有,还非常想读博,难以抉择 +20 pvrw0224 2024-06-23 42/2100 2024-06-25 10:06 by huixiong0627
[考博] 关于读博感觉自己很抓马 +13 小九月 2024-06-19 19/950 2024-06-24 23:32 by 投必得科研顾问
[教师之家] 有没有今年的影响因子? +6 jurkat.1640 2024-06-22 10/500 2024-06-24 22:48 by 潇湘之迷
[有机交流] 求助析晶问题 20+4 dengdawang 2024-06-24 5/250 2024-06-24 21:22 by cc116
[硕博家园] 怎么带研究生? +21 豆豆小小2 2024-06-20 23/1150 2024-06-24 13:39 by bluewindmill
[基金申请] b口yq上会情况 +5 XUAAAAA 2024-06-19 6/300 2024-06-24 01:13 by sunrise3148
[基金申请] 国自然青年基金,1A4B能上会吗?青年和面上的上会标准是一样的吗? +19 今晚推荐22 2024-06-20 32/1600 2024-06-23 23:17 by andywei1028
[基金申请] F03青年基金函评结果 +5 暨阳一只柴 2024-06-19 6/300 2024-06-23 14:30 by adsqsj
[公派出国] CSC博士联培对将来就业有用吗 +3 也就这样 2024-06-22 3/150 2024-06-22 22:04 by 326lhpqk
[博后之家] 在国内某高校做全职博士后2年,现在找到新的单位,出站或退站对新工作有什么影响? +10 nxplfcc 2024-06-20 10/500 2024-06-22 07:52 by 徐长安
[论文投稿] 水果保鲜投稿 5+4 zhengjiandong 2024-06-19 6/300 2024-06-21 22:27 by 宋小爷
[基金申请] 听大佬说今年信息口本子数量大幅增加? +8 wutzxt 2024-06-21 9/450 2024-06-21 19:58 by wutzxt
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
[考博] 2025年博士申请,可先做科研助理 +5 limit888 2024-06-18 8/400 2024-06-21 12:45 by 半简体
[有机交流] 想要用氢化钠拔掉吲哚N上的氢取代酰氯 50+3 光敏剂 2024-06-19 4/200 2024-06-20 18:41 by HF111001
[论文投稿] 投稿求助 6+3 小短腿安琪拉 2024-06-19 4/200 2024-06-20 17:40 by 不一样烟火12345
[论文投稿] 投稿求助 +4 平凡的日子 2024-06-19 5/250 2024-06-20 16:24 by yueyueyue@
信息提示
请填处理意见