24小时热门版块排行榜    

查看: 2090  |  回复: 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很容易做出电荷密度差。
回复此楼
弘德明志博学笃行
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2024国社科通讯评审 +9 qsd10086 2024-06-13 15/750 2024-06-16 17:30 by zhtflysky19
[基金申请] 面青地会评时间 +6 tanjydd 2024-06-15 6/300 2024-06-16 15:42 by 我是王小帅
[找工作] 杭电、天津科技、青农和宁波工程学院如何选? +11 味道很好啊 2024-06-13 21/1050 2024-06-16 13:57 by wanglongzh
[博后之家] 山东大学(青岛)“天然药物生物智造”课题组 招聘“博士后”(年薪20.4-55.6万元) +5 第二种态度 2024-06-11 8/400 2024-06-16 10:14 by 午睡未进行
[考博] 2025考博 +5 自强不息a?a 2024-06-15 7/350 2024-06-16 10:07 by wang333666
[找工作] 江西双非一本和四川双一流高校如何选择? 5+7 寒山敲钟 2024-06-12 23/1150 2024-06-16 06:52 by puterde
[基金申请] 博后面上今天有bug可以看到是否资助? +20 lyfbangong 2024-06-12 31/1550 2024-06-15 21:18 by since—2010
[考博] 上海交大招收材料化学方向科研助理/“申请考核”博士(请勿回复帖子或站内投条) +3 灵梦and紫 2024-06-12 4/200 2024-06-15 20:58 by 1822836277
[基金申请] 为什么我的博后基金还在流动站审核中?不会是学院给我卡了吧? +14 王凯12 2024-06-13 26/1300 2024-06-15 15:22 by 好人与坏人
[教师之家] 饶议:什么制度能保障大学普通教师不用为领导拎包,不用看领导脸色 +8 zju2000 2024-06-12 14/700 2024-06-15 13:59 by chemhua
[基金申请] 有没有机械的前辈分享一下评上海优都是什么成果啊 +7 wulala800 2024-06-10 7/350 2024-06-15 09:33 by 晓目崇
[基金申请] 博后基金,以往的结果点不开,怎么回事呢?最后一次机会了,两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[考博] 申博找导师 +4 疏影横斜水清浅3 2024-06-13 6/300 2024-06-14 14:31 by zxl_1105
[基金申请] 国自然基金公布的时候基金号有吗 +8 潇洒怡惜 2024-06-13 11/550 2024-06-14 11:24 by JRfei
[基金申请] 75批博后基金 +10 kyukitu 2024-06-13 13/650 2024-06-14 10:31 by kyukitu
[论文投稿] 文章proof要求使用机构的邮箱 5+3 不可不信缘 2024-06-11 11/550 2024-06-14 07:00 by 3001160025
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[有机交流] 原料反应完了,怎么知道是产物还是中间体 +6 小胡在努力 2024-06-11 8/400 2024-06-13 13:33 by 091602
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
[论文投稿] water research状态咨询 5+3 Flyyawa 2024-06-10 6/300 2024-06-11 09:45 by bobvan
信息提示
请填处理意见