计算化学公社

标题: 求助：六方晶系Fe-S化物构建超胞表面后，结构优化不收敛，请大神指教 [打印本页]

作者
Author: 昆明十一 时间: 2022-6-15 14:30
标题: 求助：六方晶系Fe-S化物构建超胞表面后，结构优化不收敛，请大神指教
计算一个六方晶系的铁硫化物Fe1S1，使用的晶体文件从Findit软件库中导出，编号为ICSD-053526，空间群号：No.194_P63-MMC。对其单胞进行几何优化，主要参数为：截断能450eV，PW91泛函，k点为7*7*4，设置自旋极化计算，Fe原子的自旋态为高自旋态+4；并且使用Hubbard U校正，Fe原子的U值设置为1.90eV，单胞几何优化能够收敛。但是切（001）面，扩展3*3*1超胞之后，再按上述参数（仅K点对应降低）进行几何优化计算，第一个离子步的SCF就一直在反复震荡，跑完设定的768个步骤，无法收敛，且试过以下方法均不收敛：

（1）增加空带数目：由默认20%增加至50%，不收敛；
（2）修改Density Mixing参数：Charge=0.3, Spin=1.50，不收敛；
（3）All Bands/EDFT：试过几个SCF，成本高到离谱，且没有明显的收敛迹象；
（4）切超胞时候更改Surface Mesh参数，不收敛；
（5）切超胞时候给底部原子加赝H原子，依旧不收敛……

单胞（Fe2S2）
(, 下载次数 Times of downloads: 22)

超胞表面（Fe27S27）
(, 下载次数 Times of downloads: 16)

附件1:优化后的单胞： (, 下载次数 Times of downloads: 1)

附件2：（001）面切出来的超胞表面 (, 下载次数 Times of downloads: 0)

附件3：测试计算时不收敛的.castep文件 (, 下载次数 Times of downloads: 0)

作者
Author: 昆明十一 时间: 2022-6-15 14:34
还试过将所有Fe原子的自旋态改为低自旋态，Spin=Low=0，情况还是依旧出现……

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