计算化学公社

标题: Fe2P的slab优化后结构变化很大，阴离子向外跑，阳离子向内跑 [打印本页]

作者
Author: wang5360848 时间: 2022-11-10 22:39
标题: Fe2P的slab优化后结构变化很大，阴离子向外跑，阳离子向内跑
构建了Fe2P的slab，如果不扩胞直接进行结构优化可以正常收敛并且结构也没问题，但是扩胞成2x2的超胞后再进行结构优化，slab结构变化都非常大，P向层外跑，Fe往层内压缩的趋势（尝试固定和不固定部分层数、更改暴露表面均是该结果）。
INCAR在附件中，扩胞后K点只算111，赝势是正常的赝势，用的vasp_gam版去运行的。
想请教一下各位问题出在哪里呢？

作者
Author: sobereva 时间: 2022-11-11 09:33
以当前POSCAR里的盒子尺寸来说，只考虑gamma点不合适，至少2*2*2

slab太薄
并且考虑实际可能出现的表面重构的可能性

作者
Author: wang5360848 时间: 2022-11-11 10:02

sobereva 发表于 2022-11-11 09:33
以当前POSCAR里的盒子尺寸来说，只考虑gamma点不合适，至少2*2*2

slab太薄

谢谢sob大佬的回复！如此简化确实是出于无奈，之前slab是4层，依旧有P逸出slab的情况。
我有一点比较好奇，扩胞前后真空层厚度都是一样的，为什么扩胞后才出现这么明显的slab变化呢？是因为自由度改变了吗？

作者
Author: ChemG 时间: 2022-11-11 10:33
优化slab模型的时候至少也要固定最下面三层的原子，哪能不固定就优化呢。主要是为了尽可能的模拟真实的晶体内部结构，而对表面的原子进行弛豫。还要在当前的结构下对k点、截断能进行测试，最好也设置一下磁性的初猜。表面重构也很正常，不过要你根据经验或者查找文献看看重构后的表面长啥样。

作者
Author: wang5360848 时间: 2022-11-11 11:01

ChemG 发表于 2022-11-11 10:33
优化slab模型的时候至少也要固定最下面三层的原子，哪能不固定就优化呢。主要是为了尽可能的模拟真实的晶体 ...

谢谢大佬的回复！固定3层放开4层的话，以我那个小破机器估计得跑一年 :-(，K点、截断能以及色散校正等对于表面重构影响大不大呢？我这个重构太夸张了，继续结构优化，一些P跑到距离slab7-8A的位置上了感觉实在不合理。

作者
Author: sobereva 时间: 2022-11-11 11:14

wang5360848 发表于 2022-11-11 10:02
谢谢sob大佬的回复！如此简化确实是出于无奈，之前slab是4层，依旧有P逸出slab的情况。
我有一点比较好 ...

没准之前的结果本来就是不合理的，过小的晶胞造成的周期性约束抑制了理应有的表面重构

作者
Author: sobereva 时间: 2022-11-11 11:16

wang5360848 发表于 2022-11-11 11:01
谢谢大佬的回复！固定3层放开4层的话，以我那个小破机器估计得跑一年 :-(，K点、截断能以及色散校正等对 ...

扩胞后，k点2*2*2够了
色散校正跟这个没直接联系
截断能做个收敛性测试便知
嫌VASP慢不如用又免费又极快的CP2K

作者
Author: wang5360848 时间: 2022-11-11 11:27

sobereva 发表于 2022-11-11 11:16
扩胞后，k点2*2*2够了
色散校正跟这个没直接联系
截断能做个收敛性测试便知

好的，谢谢sob大佬！最近是想转CP2K来着，只是刚学会vasp[捂脸]

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