VASP理论计算常见问题答疑(88):晶格匹配、文件解析、表面能、溶剂化效应与电荷密度单位

Q1:朱老师,a方向晶格常数不失配,但b方向失配了,还要继续调b吗?

A:①需要,直到都匹配;

②做supercell配对;

③施加应变;

④换取向或换切面方向

Q2:朱老师,我粘贴POSCAR相关原子坐标,前面信息总是不显示,尝试windows里面建好上传仍然不显示,这种情况应该如何解决?

VASP理论计算常见问题答疑(88):晶格匹配、文件解析、表面能、溶剂化效应与电荷密度单位

A:通常不是 VASP 的问题,而是文本格式或不可见字符把解析器搞崩了;

检查修改Windows的隐藏字符/编码/行尾;

还有可能是粘贴时聊天窗口把多空格或行首空格自动删除,需要进行检查。

Q3:朱老师,如果我想构建一个不同晶面在Li+暴露下稳定性的Wulff构造模型,我应该是计算吸附能还是表面形成能呢?

A:用于 Wulff 构造模型,要算表面形成能或表面能:

吸附能只能说明在某个晶面上放一个 Li(或一定覆盖度)是否容易吸附;

但 Wulff 构造需要的是:不同晶面在相同环境下谁的表面自由能更低,低的面长得更大。

Q4:老师,CP-VASP做SCF计算时,巨正则系综的能量有没有考虑溶剂化效应?

VASP理论计算常见问题答疑(88):晶格匹配、文件解析、表面能、溶剂化效应与电荷密度单位

A:一般都是考虑的,隐式溶剂;

最稳的办法是看 OUTCAR/OSZICAR 是否出现溶剂模型的能量分解输出;

仅切换 LSOL=.TRUE.  / .FALSE. ,如果开启后 E0/F 和 GCE 都系统性变化,就说明 GCE 确实包含了溶剂模型贡献。

Q5老师,请教一下CHGCAR里面的电荷密度,要变成e/Bohr3或者其它真实电荷密度,是要用CHGCAR的数值÷还是×晶胞实际体积?

A:VESTA可以转换单位,edit data  volumetric data

不需要乘/除晶胞体积,只做单位换算即可:从 e/ų 转 e/Bohr³:乘以 0.1481847;

如果想从密度变成真实电荷(电子数目):乘以ΔV即可。

 

华算科技提供专业的第一性原理、分子动力学、生物模拟、量子化学、机器学习、有限元仿真等代算服务。

500+博士团队护航,累计助力5️⃣0️⃣0️⃣0️⃣0️⃣➕篇科研成果,计算数据已发表在Nature & Science正刊及大子刊、JACS、Angew、PNAS、AM系列等国际顶刊。

声明:如需转载请注明出处(华算科技旗下资讯学习网站-学术资讯),并附有原文链接,谢谢!
(0)
上一篇 4天前
下一篇 2025年8月13日 上午10:20

相关推荐