计算化学公社

标题: 计算差分电荷密度后如何得到异质结界面处转移的电子数？ [打印本页]

作者
Author: jiandan 时间: 2024-4-21 14:07
标题: 计算差分电荷密度后如何得到异质结界面处转移的电子数？
请问各位老师，在计算异质结的差分电荷密度后，是通过CHGCAR中哪一个结果做差得到界面转移电子数呢？或者是通过其他计算得到？

作者
Author: sobereva 时间: 2024-4-22 06:28
绘制电荷位移曲线
使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线
http://sobereva.com/638（http://bbs.keinsci.com/thread-28225-1-1.html）
Multiwfn基于CHGCAR也可以按上文绘制，上文说了这点

作者
Author: jiandan 时间: 2024-4-22 10:21

sobereva 发表于 2024-4-22 06:28
绘制电荷位移曲线
使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线
http://soberev ...

老师您好，我用Multiwfn导入CHGCAR文件后输入13报错：
Error: Grid data has not been loaded or generated! If you want to load a grid data file now, input its path, e.g. C:\nico.cub, else input 0 to return
请问是INCAR中参数设置的问题吗？
INCAR文件如下：
ISTART = 1
ISPIN = 2

LREAL = .FALSE.
LWAVE = .TRUE.
LCHARG = .TRUE.
ADDGRID = .TRUE.

ISMEAR = 0
SIGMA = 0.05
NELM = 100
NELMIN = 6
EDIFF = 1E-5

NGXF = 120
NGYF = 120
NGZF = 320

LCHGCAR = T
LELF = T
LAECHG = T

作者
Author: 卡开发发 时间: 2024-4-22 11:20
本帖最后由卡开发发于 2024-4-22 11:30 编辑

转移电子数需要通过布居分析得到原子电荷，根据片段原子电荷积累进行计算，当然不同的原子电荷种类在划分策略有所不同，导致的结果也有所不同。目前VASP能够借助Henkelman组的Bader程序做Bader电荷（一种AIM电荷），VASP自身可以进行Hirsfeld系列的电荷计算，也有不太严格通过WS半径划分得到的电荷，虽然Bader电荷用的广泛，但确实有其缺陷。

作者
Author: sobereva 时间: 2024-4-23 06:44

jiandan 发表于 2024-4-22 10:21
老师您好，我用Multiwfn导入CHGCAR文件后输入13报错：
Error: Grid data has not been loaded or gener ...

把在Multiwfn里输入的所有命令（包括一开始载入文件写的文件路径）完整贴出来。最好把CHGCAR压缩后传网盘便于我检查

我不用VASP。这里只说CHGCAR本身的事

作者
Author: jiandan 时间: 2024-4-23 11:01

卡开发发发表于 2024-4-22 11:20
转移电子数需要通过布居分析得到原子电荷，根据片段原子电荷积累进行计算，当然不同的原子电荷种类在划分策 ...

谢谢老师，请问有没有可以参考的博文或者是教程？

作者
Author: jiandan 时间: 2024-4-23 11:12

sobereva 发表于 2024-4-23 06:44
把在Multiwfn里输入的所有命令（包括一开始载入文件写的文件路径）完整贴出来。最好把CHGCAR压缩后传网盘 ...

(, 下载次数 Times of downloads: 7)
将CHGCAR拖入框中，然后输入13出现错误。

CHGCAR文件百度网盘地址：
链接：https://pan.baidu.com/s/1S8bUX-5E4RqS6fs6_ZJnlg?pwd=7ml0
提取码：7ml0

感谢老师

作者
Author: 卡开发发 时间: 2024-4-23 14:26
本帖最后由卡开发发于 2024-4-23 14:27 编辑

jiandan 发表于 2024-4-23 11:01
谢谢老师，请问有没有可以参考的博文或者是教程？

VASP怎么去做Bader电荷你可以搜索一下，资料不太难找。剩下的问题是，假定两个片段相互作用，即A+B=AB，那么你需要单独计算A的Bader电荷累加（注意扣除核电荷）。当然你也可以用Hirshfeld，那个需要开着IVDW=2X独立做一个自洽场计算（也就是单点或者说静态计算）。不过这两种电荷其实都不怎么样，往往就是看个趋势。此外你也可以使用DDEC6电荷，可以参考这个链接，从文献表现来说可靠性还可以，使用方法你可以好好看他的Readme，门槛不太高。

作者
Author: jiandan 时间: 2024-4-23 22:28

卡开发发发表于 2024-4-23 14:26
VASP怎么去做Bader电荷你可以搜索一下，资料不太难找。剩下的问题是，假定两个片段相互作用，即A+B=AB， ...

老师您好，请问您说的第一种方法是分别将A组分、B组分、AB组分的ACF.dat文件中CHARGE相加（在加之前将氧原子减6，金属元素减10等处理）。分别形成A、B、AB的电荷累计，然后用AB-A-B得到界面电子转移数吗？

作者
Author: 卡开发发 时间: 2024-4-24 04:35

jiandan 发表于 2024-4-23 22:28
老师您好，请问您说的第一种方法是分别将A组分、B组分、AB组分的ACF.dat文件中CHARGE相加（在加之前将氧 ...

如果你做了和核电荷相加的话，其实单独求和A的就可以了，因为AB最终几乎是0（实际经过Bader划分和盆积分之后会有一点数值误差，但网格足够密这个误差就很小）。

作者
Author: sobereva 时间: 2024-4-24 06:21

jiandan 发表于 2024-4-23 11:12
将CHGCAR拖入框中，然后输入13出现错误。

CHGCAR文件百度网盘地址：

别用老掉牙的Multiwfn版本

作者
Author: jiandan 时间: 2024-4-24 10:03

卡开发发发表于 2024-4-24 04:35
如果你做了和核电荷相加的话，其实单独求和A的就可以了，因为AB最终几乎是0（实际经过Bader划分和盆积分 ...

您的意思是将AB异质结所产生的ACF.dat文件中的CHARGE扣除核电荷后，然后对A组分的原子求和得到转移数吗？
（比如说AB异质结的ACF.dat文件中A原子序数为1-7，b的为8-10。那么根据如下步骤计算：①首先将ACF.dat文件中的CHARGE列分别减去原子的核电荷（比如氧原子减去6，镍原子减去10等等）。②将A处理后的CHARGE相加，也就是对1-7个原子处理后的CHARGE求和。最终得到异质结界面电子转移数。）

欢迎光临计算化学公社 (http://bbs.keinsci.com/)