计算化学公社

标题: 关于二维材料真空层优化的问题 [打印本页]

作者
Author: 汪杰 时间: 2024-10-23 09:23
标题: 关于二维材料真空层优化的问题
我使用一层二维材料，真空层初始设为25埃，结构优化后，为什么它的真空层就只有9埃了。
我看很多人说，这个优化最多也就优化掉几埃，为什么我这会优化这么多？
有大佬能够给点建议吗？不用重新编译vasp的那种。

作者
Author: sobereva 时间: 2024-10-23 12:18
我不用VASP，对于CP2K的情况，此时完全没必要也不应该令垂直于界面方向的尺寸被优化

作者
Author: 卡开发发 时间: 2024-10-24 10:46
本帖最后由卡开发发于 2024-10-24 10:48 编辑

如果不考虑重新编译或者增加其他的代码，在6.4.3之后的版本有个关键字叫LATTICE_CONSTRAINTS，老版本可以考虑使用isif=4进行恒定体积的计算，但有可能数值上没那么稳定，还有别的办法就是用ase的优化器去调用VASP进行结构优化，在ase中晶胞信息通过掩码参数要求不对c方向晶胞参数进行优化。

二维体系如果不做额外约束，除非c周足够足够长，否则原子间在c方向仍然会存在一些作用导致真空层被压缩甚至结构优化难收敛。

作者
Author: 汪杰 时间: 2024-10-25 13:39

卡开发发发表于 2024-10-24 10:46
如果不考虑重新编译或者增加其他的代码，在6.4.3之后的版本有个关键字叫LATTICE_CONSTRAINTS，老版本可以考 ...

感谢大佬，我还是费点功夫给重新编译完了

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