24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

forkn

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 232.5
红花: 1
帖子: 77
在线: 9.6小时
虫号: 626283
注册: 2008-10-14
性别: GG
专业: 材料

[交流] 【求助】如何解决自旋极化计算时的难以收敛问题？已有9人参与

小弟想用castep计算两种异质界面(Fe/Ni)结合的电子结构和能量学信息，
构建模型的过程如下：
1 .导入程序自带的Fe的晶胞，几何优化,切出(100)面，构建真空层后几何优化，得出Fe(100)的构型。晶胞参数a=b=2.775埃
2 .导入程序自带的Ni的晶胞，几何优化,切出(100)面，构建真空层后几何优化，得出Ni(100)的构型。晶胞参数c=d=2.513埃
3. 根据优化后的Ni(100)为底，向真空层添加Fe原子，界面间距选择为优化后的两种结构的层间距的平均值。Fe的层间距采用步骤1的数据。以此搭建Fe(100)/Ni(100)结构,晶胞边距选择为（a+c)/2，入下图：

现在的问题是：
（1）以上的几何优化都是在加spin polarized下进行的，对这个自己建的结构加自旋极化（自旋设为low和hight都试过）进行几何优化时，一直不收敛，我已经尝试了增加energy cut off (350ev), 增加k point set (7*7*2），empty bonds, 降低SCF 精度，降低smearing等，都不收敛，总是提示：
*Warning* max. SCF cycles performed but system has not reached the groundstate.
Warning: electronic minimisation did not converge when finding ground state
收敛曲线：

能量曲线：

（2）我计算过不加自旋的优化，很容易收敛，energy cut off (300ev), k point set (3*3*2），就可以收敛。但这样的结构不可信的，本来就是磁性体系，况且加自旋极化计算的能量都比不加的低。

是不是我的初始结构不太合理啊？或是计算参数的问题？希望大家帮帮忙吧，最好能结合我的情况给些加速收敛具体的措施。急等！

回复此楼