24小时热门版块排行榜    

查看: 2907  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 材料专硕326求调剂 +6 墨煜姒莘 2026-03-15 7/350 2026-03-17 17:10 by ruiyingmiao
[考研] 本人考085602 化学工程 专硕 +16 不知道叫什么! 2026-03-15 18/900 2026-03-17 17:05 by ruiyingmiao
[考研] 311求调剂 +8 冬十三 2026-03-15 8/400 2026-03-17 16:59 by ruiyingmiao
[考研] 求调剂,总分315,考的生物医药,一志愿湖南师范大学。调剂到任何专业都可以 +4 小丁想进步 2026-03-11 5/250 2026-03-17 16:05 by 外星文明
[考研] 302求调剂 +4 小贾同学123 2026-03-15 8/400 2026-03-17 10:33 by 小贾同学123
[考研] 材料专硕274一志愿陕西师范大学求调剂 +5 薛云鹏 2026-03-13 5/250 2026-03-17 10:15 by Sammy2
[考研] 274求调剂 +5 时间点 2026-03-13 5/250 2026-03-17 07:34 by 热情沙漠
[考研] 304求调剂 +3 曼殊2266 2026-03-14 3/150 2026-03-16 16:39 by houyaoxu
[考研] 0703化学调剂,求各位老师收留 +8 秋有木北 2026-03-14 8/400 2026-03-16 15:21 by 哦哦123
[考研] 085600材料与化工 求调剂 +13 enenenhui 2026-03-13 14/700 2026-03-16 15:19 by 了了了了。。
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 材料专硕288分求调剂 一志愿211 +4 在家想你 2026-03-11 4/200 2026-03-13 22:49 by JourneyLucky
[考研] 材料与化工求调剂一志愿 985 总分 295 +8 dream…… 2026-03-12 8/400 2026-03-13 22:17 by 星空星月
[考研] 315求调剂 +9 小羊小羊_ 2026-03-11 10/500 2026-03-13 21:13 by SXNU李老师
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 290求调剂 +3 柯淮然 2026-03-10 8/400 2026-03-11 13:48 by 柯淮然
信息提示
请填处理意见