VASP | 理论计算常见问题解答(93):溶剂、d 带中心、扩胞、晶格匹配与偶极修正答疑

Q1:朱老师,计算的时候需要加隐式溶剂MD/PCM,是需要在哪里加啊?

A:incar  lsol=.t.;

在 VASP 里它不是叫 PCM,而通常是 VASPsol(implicit solvation),加的位置就在 INCAR;

做AIMD计算的时候也是一样,增加对应参数LSOL;

总结就是:在INCAR中写入:ISOL=.TRUE.,默认溶剂是水。

Q2:朱老师,计算d带中心的时候 是要选择一个原子进行计算还是要统计全部的原子进行计算 我看课程是Pt单质 我的体系是AlNiYCo Al没有d轨道的贡献也要统计吗?

A:一般不是全体系所有原子无脑平均,而是选定有d态且是活性位点的原子集合;

算 d 带中心时,核心是:要让这个指标代表参与成键,吸附的d电子金属位点;

Al不需要统计,它的d态基本接近于0。

Q3:朱老师,为啥有的扩胞根号3那样扩的,这样得咋弄?

A:不是沿原来的晶格矢量整倍数扩胞(2×2×1、3×3×1),而是换了一组新的面内基矢,让新晶胞面积变成原来的3倍;
做 1/3覆盖度:在根号3超胞里放 1 个吸附物 → 相当于原 1×1 表面的 1/3 覆盖度,经典的吸附覆盖度构型;
具体做法是旋转晶胞,设置变换矩阵即可。

Q4:老师,120度变90度之后a,b晶格常数不一样了,那么a方向晶格常数不失配,但b方向失配了,还要继续调b吗?

A:一般是需要继续调整,直到a、b和角度都适配,不过也分情况;

如果只是 120°→90° 的等价变换:b 变不匹配多半是表示法问题,不用调 b,先统一基矢再比较;

如果要做周期界面模型:b 方向真失配就必须处理,继续调b看失配百分比,优先用超胞匹配。

Q5老师,如果我结构优化时加偶级修正收敛很不好,那可以先不加偶级修正吗?

A:可以的,但是如果要计算功函数的时候,那就需要加上,如果要严格测试,那就这样做:

先不加偶极修正,把结构优化收敛下来;

在最终结构上打开偶极修正做静态单点,比较开关偶极矫正的能量数值;

如果能量差异很大,则需要加上偶极修正重新优化。

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

相关推荐