卡开发发 发表于 2024-4-24 04:35 您的意思是将AB异质结所产生的ACF.dat文件中的CHARGE扣除核电荷后，然后对A组分的原子求和得到转移数吗？ （比如说AB异质结的ACF.dat文件中A原子序数为1-7，b的为8-10。那么根据如下步骤计算：①首先将ACF.dat文件中的CHARGE列分别减去原子的核电荷（比如氧原子减去6，镍原子减去10等等）。②将A处理后的CHARGE相加，也就是对1-7个原子处理后的CHARGE求和。最终得到异质结界面电子转移数。） |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 | 对。 |
jiandan 发表于 2024-4-23 11:12 别用老掉牙的Multiwfn版本 |
jiandan 发表于 2024-4-23 22:28 如果你做了和核电荷相加的话，其实单独求和A的就可以了，因为AB最终几乎是0（实际经过Bader划分和盆积分之后会有一点数值误差，但网格足够密这个误差就很小）。 |
卡开发发 发表于 2024-4-23 14:26 老师您好，请问您说的第一种方法是分别将A组分、B组分、AB组分的ACF.dat文件中CHARGE相加（在加之前将氧原子减6，金属元素减10等处理）。分别形成A、B、AB的电荷累计，然后用AB-A-B得到界面电子转移数吗？ |
本帖最后由 卡开发发 于 2024-4-23 14:27 编辑 jiandan 发表于 2024-4-23 11:01 VASP怎么去做Bader电荷你可以搜索一下，资料不太难找。剩下的问题是，假定两个片段相互作用，即A+B=AB，那么你需要单独计算A的Bader电荷累加（注意扣除核电荷）。当然你也可以用Hirshfeld，那个需要开着IVDW=2X独立做一个自洽场计算（也就是单点或者说静态计算）。不过这两种电荷其实都不怎么样，往往就是看个趋势。此外你也可以使用DDEC6电荷，可以参考这个链接，从文献表现来说可靠性还可以，使用方法你可以好好看他的Readme，门槛不太高。 |
卡开发发 发表于 2024-4-22 11:20 谢谢老师，请问有没有可以参考的博文或者是教程？ |
jiandan 发表于 2024-4-22 10:21 把在Multiwfn里输入的所有命令（包括一开始载入文件写的文件路径）完整贴出来。最好把CHGCAR压缩后传网盘便于我检查 我不用VASP。这里只说CHGCAR本身的事 |
sobereva 发表于 2024-4-22 06:28 老师您好，我用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 |
|
绘制电荷位移曲线 使用CP2K结合Multiwfn绘制密度差图、平面平均密度差曲线和电荷位移曲线 http://sobereva.com/638（http://bbs.keinsci.com/thread-28225-1-1.html） Multiwfn基于CHGCAR也可以按上文绘制，上文说了这点 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-24 18:40 , Processed in 0.191415 second(s), 27 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites