Q1:朱老师,想看圈起来那三个原子的差分电荷二维图,应该切哪个面呢?

A:只能先算整体3D差分,再截面获取2D差分图;
Edit → Plane → New → By 3 points,在结构窗口点这 3 个原子;这会生成唯一经过三点的平面;
Properties → Isosurface/Slice 勾选 Slice,选择刚建的平面,Thickness 设很小(如 0.1–0.2 Å);
Contour 选择等高线;色条设置为对称范围(如 ±0.02 e/ų)
Q2:朱老师,单个原子的电荷密度数值利用vasp在哪个文件里?

A:VASP无法直接得到各原子电荷转移数目;
需要计算电荷,做bader电荷布居分析;
通过bader后处理得到ACF.dat。
Q3:朱老师,如果想比较两个模型表面的静电势,除了F(min)~F(max)范围保持一致,这里的颜色设置参数Saturation level (Max, Min)是否要保持相同取值范围?
A:要“可比”,不仅要统一 F(min)–F(max) 的显示范围,Saturation level(Max/Min)也必须设成同一组绝对数值;

VESTA 的 Saturation level 就是颜色条的上下限;超出这个范围的值会被截断为端点颜色。
如果两个模型用不同的饱和上下限,同样颜色对应的数值就不同,视觉对比会失真。
Q4:朱老师,如何导出两个原子间的电荷转移积分?
A:保证三份电荷密度在同一晶格和同一 FFT 网格下计算,然后做差得到Bader 原子盆内积分 Δρ;
用总体系的密度做 Bader 分区(建议用全电子密度 AECCAR0+2 作为 -ref);
在这些原子盆内对 CHGCAR_diff 积分:critic2或将权重乘到 差分电荷文件上并积分。
Q5:朱老师,这个六元环散了是什么原因呀?

A:上下原子作用太强了,离得远一些,设置在3埃左右;
自由度太多、一次性全放,不要一下子让环的所有原子与邻近无机位全部自由;
超胞太小,面内 1×1 或 2×1 让相邻影像强相互作用,分子会相互牵扯;
晶格放松过早,一上来 ISIF=3 放开晶格,软框架会被拉伸变形;
下一篇:没有了