求助电荷差分密度计算详细步骤
答案:2 悬赏:30 手机版
解决时间 2021-03-20 14:07
- 提问者网友：遁入空寂
- 2021-03-19 20:22
求助电荷差分密度计算详细步骤
最佳答案
- 五星知识达人网友：往事隔山水
- 2021-03-19 21:13
差分电荷密度是什么？
用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://emuch.net/bbs/viewthread.php?tid=4287522&fpage=1 这个帖子也不错，很详细 @沙风）
===============
转自：http://emuch.net/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很容易做出电荷密度差。
用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://emuch.net/bbs/viewthread.php?tid=4287522&fpage=1 这个帖子也不错，很详细 @沙风）
===============
转自：http://emuch.net/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很容易做出电荷密度差。
全部回答
- 1楼网友：刀戟声无边
- 2021-03-19 22:38
江上往来人,但爱鲈鱼美.
我要举报
如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
点此我要举报以上问答信息
大家都在看
推荐资讯