求助：六方晶系Fe-S化物构建超胞表面后，结构优化不收敛，请大神指教

昆明十一 · 发表于 Post on 2022-6-15 14:30:02

计算一个六方晶系的铁硫化物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）

超胞表面（Fe27S27）

附件1:优化后的单胞：

04_Spin=High= 4_FormalSpin=On_Optimized.xsd (45.34 KB, 下载次数 Times of downloads: 1)

附件2：（001）面切出来的超胞表面

15_(001)_(U=1-0-0)_0.3_(V=0-1-0)_0.3_[3-3-1]_SuperCell.xsd (56.73 KB, 下载次数 Times of downloads: 0)

附件3：测试计算时不收敛的.castep文件

02_(001)_[2-2-1]SuperCell_Spin= 4[-4].castep (93.15 KB, 下载次数 Times of downloads: 0)

昆明十一 · 发表于 Post on 2022-6-15 14:34:37

还试过将所有Fe原子的自旋态改为低自旋态，Spin=Low=0，情况还是依旧出现……

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[CASTEP/Dmol3/MS] 求助：六方晶系Fe-S化物构建超胞表面后，结构优化不收敛，请大神指教

浏览过的版块