VASP计算与可视化综合答疑:从电荷差分图切面、静电势比对到结构弛豫问题诊断

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

VASP计算与可视化综合答疑:从电荷差分图切面、静电势比对到结构弛豫问题诊断

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在哪个文件里?

VASP计算与可视化综合答疑:从电荷差分图切面、静电势比对到结构弛豫问题诊断

A:VASP无法直接得到各原子电荷转移数目;

需要计算电荷,做bader电荷布居分析;

通过bader后处理得到ACF.dat。

Q3:朱老师如果想比较两个模型表面的静电势,除了F(min)~F(max)范围保持一致,这里的颜色设置参数Saturation level (Max, Min)是否要保持相同取值范围?

VASP计算与可视化综合答疑:从电荷差分图切面、静电势比对到结构弛豫问题诊断

A:要“可比”,不仅要统一 F(min)–F(max) 的显示范围,Saturation level(Max/Min)也必须设成同一组绝对数值;

VESTA 的 Saturation level 就是颜色条的上下限;超出这个范围的值会被截断为端点颜色。

如果两个模型用不同的饱和上下限,同样颜色对应的数值就不同,视觉对比会失真。

Q4:老师,如何导出两个原子间的电荷转移积分?

A:保证三份电荷密度在同一晶格和同一 FFT 网格下计算,然后做差得到Bader 原子盆内积分 Δρ;

用总体系的密度做 Bader 分区(建议用全电子密度 AECCAR0+2 作为 -ref);

在这些原子盆内对 CHGCAR_diff 积分:critic2或将权重乘到 差分电荷文件上并积分。

Q5:朱老师这个六元环散了是什么原因呀?

VASP计算与可视化综合答疑:从电荷差分图切面、静电势比对到结构弛豫问题诊断

A:上下原子作用太强了,离得远一些,设置在3埃左右;

自由度太多、一次性全放,不要一下子让环的所有原子与邻近无机位全部自由;

超胞太小,面内 1×1 或 2×1 让相邻影像强相互作用,分子会相互牵扯;

晶格放松过早,一上来 ISIF=3 放开晶格,软框架会被拉伸变形;

为确保大家的学习成果得以沉淀和复用,我们特将朱老师本次及往期的VASP理论计算答疑精华,系统汇编为【理论计算答疑精选汇编】。它将成为你科研路上随取随用的“速查手册”,助你高效攻克难题。
🏅 华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。 🎯500+博士团队护航,累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果,计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。 👏👏👏

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 12小时前
下一篇 2025年6月6日 上午11:18

相关推荐