24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2945  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 工科 267求调剂 +6 wanwan00 2026-04-02 8/400 2026-04-08 22:09 by 朱云虎202
[考研] 293分求调剂,外语为俄语 +7 加一一九 2026-04-07 10/500 2026-04-08 20:14 by yutian743
[考研] 285求调剂 +18 AZMK 2026-04-02 19/950 2026-04-08 10:33 by xingguangj
[考研] 求调剂 +11 wwwwabcde 2026-04-07 11/550 2026-04-07 23:16 by JourneyLucky
[考研] 专硕085403,291分,有两篇专利,一国一奖 +3 哈吉咪哈吉咪 2026-04-07 3/150 2026-04-07 18:21 by 蓝云思雨
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
[考研] 312求调剂 +4 LR6 2026-04-06 4/200 2026-04-07 08:42 by jp9609
[考研] 297分083200求助 +9 aekx 2026-04-05 9/450 2026-04-06 20:57 by flysky1234
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 10/500 2026-04-06 19:47 by jean5056
[考研] 求调剂 +11 xzghyuj 2026-04-04 11/550 2026-04-06 11:49 by lijunpoly
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 材料与化工306分找调剂 +12 沧海轻舟e 2026-04-03 13/650 2026-04-04 23:45 by lqwchd
[考研] 一志愿南农090401,268,求调剂 +5 一木鸟然 2026-04-04 5/250 2026-04-04 17:07 by babysonlkd
[考研] 一志愿北交大材料工程总分358 +6 cs0106 2026-04-03 6/300 2026-04-04 11:20 by w_xuqing
[考研] 材料调剂 +11 吴棂颖! 2026-04-03 11/550 2026-04-04 09:56 by 小小树2024
[考研] 材料科学与工程考研 +10 拯救皮特托先生 2026-04-02 10/500 2026-04-03 23:57 by userper
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 求调剂 +3 usbdndj 2026-04-03 3/150 2026-04-03 14:10 by dxiaoxin
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
信息提示
请填处理意见