计算化学公社

标题: VASP计算的Ni3P表面能与实验不符，求教！ [打印本页]

作者
Author: caoliqun 时间: 2018-6-1 16:17
标题: VASP计算的Ni3P表面能与实验不符，求教！
本帖最后由 caoliqun 于 2018-6-1 18:32 编辑

　　大家好！我最近在计算表面能与实验XRD图对比上遇到难以解决的问题！跪求各位大神赐教！十分感谢！感激涕零！
这是我用vasp计算表面能所用公式：

　　体系是Ni3P,分别取(3 2 1) (1 4 1)( 1 1 2)晶面，MS建模，为了对称，切面时取四层厚度，固定中间两层，放开上下两层做优化得Eslab,Ebulk通过计算原胞能量/原胞中分子数得到，即Ebulk是一个分子Ni3P的能量.A是面积．
　　实验结果XRD 图中峰高(3 2 1)> (1 4 1)> (1 1 2),结论是(3 2 1)面所占比例较大，我们认为表面能越低晶面越稳定，但是计算得到的(3 2 1)面的表面能最大．想求教大家几个问题：(1)用MS切面时不一定切到的原子是按Ni:P=3:1,求教应该怎么解决？(2)求大神指导我们用vasp计算得方法以及设置参数是否准确（附带INCAR 　KPOINTS文件）(3)晶面所占比例是否仅与表面能成反比，请问还有什么影响因素?

PS:我们还使用MS 粗略计算了各晶面表面能，使用相同公式，所得结果仍然是(3 2 1)晶面表面能最大．附加实验XRD图，计算表面能结果数据．
谢谢各位老师同学们不吝赐教！

作者
Author: 卡开发发 时间: 2018-6-1 23:21
本帖最后由卡开发发于 2018-6-1 23:22 编辑

（1）应当先保证化学计量比，不对称的情形可以用偶极修正来进行处理，否则表面能不对不说甚至SCF会很难收敛。
（2）参数问题并不严重，不过EDIFFG=-0.01这样的精度最好PREC稍微高点。
（3）严格来说不能说是反比，只能说是负相关。表面能应该是热力学方面的，如果从晶体劈裂成两个面或是从一种表面重构成另一种表面也许还有动力学的因素，但如何考虑我暂时也想不到。

PS：不知道你用的MS哪个模块，至少说明这个模块还是比较可靠而已。另外，Ni3P的Ni恐怕要小心，一般的DFT处理的结果有可能会不好。

作者
Author: caoliqun 时间: 2018-6-3 16:27
　　谢谢您的悉心指教，我用的MS模块Forcite算的，采用分子力场比较快，用CASTEP算不动．
　　在PREC设置比较高的时候会出现这个问题：
　　WARNING: aliasing errors must be expected set NGX to 94 to avoid them
　　WARNING: aliasing errors must be expected set NGY to 98 to avoid them
　　WARNING: aliasing errors must be expected set NGZ to 110 to avoid them
查的结果是把精度设置小一点可以避免这个问题，所以精度没有设置太高，这个问题应该也可以在INCAR中设置具体NGX等解决，您的回复我会认真学习，谢谢您的帮助！

作者
Author: 卡开发发 时间: 2018-6-3 16:55

caoliqun 发表于 2018-6-3 16:27
　　谢谢您的悉心指教，我用的MS模块Forcite算的，采用分子力场比较快，用CASTEP算不动．
　　在PREC设置 ...

这个warning不用管。

作者
Author: caoliqun 时间: 2018-6-3 18:34

卡开发发发表于 2018-6-3 16:55
这个warning不用管。

好的，我设置低精度后其实也没完全消除这个warning

欢迎光临计算化学公社 (http://ccc.keinsci.com/)