应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【求助】如何解决自旋极化计算时的难以收敛问题？
71/1返回列表
查看: 4148  |  回复: 16
查看全部回帖@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

forkn

银虫 (小有名气)

[交流] 【求助】如何解决自旋极化计算时的难以收敛问题？已有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），就可以收敛。但这样的结构不可信的，本来就是磁性体系，况且加自旋极化计算的能量都比不加的低。

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

» 收录本帖的淘帖专辑推荐

第一原理资料汇编 siesta问题 心理互助

» 本帖已获得的红花（最新10朵）

歇斯底琳

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
山阻石拦，大江毕竟东流去；雪辱霜欺，梅花依旧向阳开！
1楼2010-06-04 18:14:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)
引用回帖:
Originally posted by acridine at 2010-06-05 00:19:47:
你提到的energy cut off , k point set等参数，确实可以影响收敛，但他们最主要还是影响电子scf的收敛。从你的优化过程看，已经走了20多步了，所以说scf自洽的收敛情况还是相当漂亮的，这些参数的设置上不会有太大 ...

谢谢你的建议，原以为加上优化晶胞会使结构松弛的更好些，看来不一定是这样了。
还有，3层好像是太少了，那么原子层数的选取原则是什么啊？如何知道选择3层还是5层比较合理呢？
赞一下 回复此楼
山阻石拦，大江毕竟东流去；雪辱霜欺，梅花依旧向阳开！
4楼2010-06-06 11:31:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)
★ ★
cenwanglai(金币+2):谢谢参与讨论! 2010-06-06 15:07:57
引用回帖:
Originally posted by 163.com at 2010-06-05 06:29:45:

第一性原理计算指出，加不加自旋对对体系的结构变化不大。

是这样吗？我之前算过一些结构，好像加自旋和不加好像结构差的还不小呢。最外层的原子向真空层弛豫相差了0.4埃左右呢。
赞一下(1人) 回复此楼
山阻石拦，大江毕竟东流去；雪辱霜欺，梅花依旧向阳开！
5楼2010-06-06 11:48:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)
★ ★ ★
xiaohunhun(金币+3):谢谢 2010-06-08 23:49:18
引用cenwanglai的测试结果：  
“对于含有未配对电子的结构,计算能量时是否考虑自旋对结果影响不可忽略,而且最好在结构优化时就考虑自旋最好;但是作为计算时间成本的考虑,在结构优化时不考虑自旋,而在后续能量计算时考虑自旋可能对能量精度影响不大;但是是否在结构优化时考虑自旋对结构有一定影响.”
详见：http://muchong.com/bbs/viewthread.php?tid=1532270&fpage=0&view=&highlight=&page=2
本人觉得还是结构优化时就加上自旋比较好。至于先不加自旋优化，再此结构上再加自旋优化的方法，我正在试，希望能收敛吧,呵呵。
赞一下(3人) 回复此楼
山阻石拦，大江毕竟东流去；雪辱霜欺，梅花依旧向阳开！
8楼2010-06-07 09:57:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)
引用回帖:
Originally posted by zhenfeili at 2010-06-07 21:03:52:
我采用的pc机，有时有些复杂的胞或者原子数比较大的胞，
采用自旋极化优化超级胞时计算机算不动，所以就采用先结构优化再做自旋计算
如果主要为了得到定性的结果，节约时间，这种方法还是可以的

我觉得很多半 ...

呵呵，我主要关心的是能量问题，如果加自旋优化真的很难算的话，也只能先结构优化再做自旋计算能量了。其实第一原理计算的结果相对值才有意义，绝对值是很难准确定量描述问题了。
赞一下 回复此楼
山阻石拦，大江毕竟东流去；雪辱霜欺，梅花依旧向阳开！
10楼2010-06-08 10:03:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

forkn

银虫 (小有名气)
我已尝试在不加自旋优化出的结构上加自旋计算，交换关联能 GGA PBE ,E cut =350ev ，k point=7*7*5 , smearing=0.2 , emptybonds =13 .  density mixing 作电子最小化，还是SCF不收敛， 改用all bonds/ EDFT 作电子最小化，算了一天，能量只走了一步，我给kill掉了。
现在都不知该怎么算了，郁闷啊！
赞一下 回复此楼
山阻石拦，大江毕竟东流去；雪辱霜欺，梅花依旧向阳开！
12楼2010-06-12 10:46:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 forkn 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭