| 查看: 4466 | 回复: 16 | ||||||
[交流]
【求助】如何解决自旋极化计算时的难以收敛问题？ 已有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朵）
歇斯底琳» 猜你喜欢
281求调剂（0805）
已经有8人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有3人回复
-
268求调剂
已经有3人回复
-
311求调剂
已经有10人回复
-
被我言中：新模板不强调格式了，假专家开始管格式了
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- castep 自旋极化问题
已经有10人回复
- 为什么我加自旋极化算出来的态密度用gnuplot画出来只有自旋向上态而没有自旋向下的态
已经有18人回复
- 关于自旋极化的设置
已经有4人回复
- 掺杂TiO2 自旋极化
已经有4人回复
- 【求助】自旋极化出错
已经有9人回复
- 【求助】关于拟合得到体弹性模量&自旋极化&LDA+U等问题，望大家指点指点
已经有4人回复
- 【求助】VASP怎么考虑自旋极化？
已经有14人回复
- 【求助】自旋极化计算设置的问题
已经有3人回复
- 【求助】什么样的图形能说明原子发生了自旋极化...？
已经有3人回复
- 【求助】自旋极化计算
已经有4人回复
- 【前沿】电子自旋极化首次在室温下实现
已经有4人回复
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
zzy870720z(金币+3):谢谢指导，很不错的注意，学习了 2010-06-05 08:32:23
forkn(金币+6):谢谢你的建议，我试试看 2010-06-06 10:53:22
zzy870720z(金币+3):谢谢指导，很不错的注意，学习了 2010-06-05 08:32:23
forkn(金币+6):谢谢你的建议，我试试看 2010-06-06 10:53:22
|
你提到的energy cut off , k point set等参数，确实可以影响收敛，但他们最主要还是影响电子scf的收敛。从你的优化过程看，已经走了20多步了，所以说scf自洽的收敛情况还是相当漂亮的，这些参数的设置上不会有太大的问题。那么问题就应该是来自结构的优化上。 其实俺一看到真空层+优化晶格就很晕，个人认为最好不要这么干。这就是相当于对电脑“下套”，电脑又不象我们这么聪明，它也不知道啥叫真空层，它会费力不讨好的不停优化c轴长度，而你应该根本就不关心这个方向的晶格长度吧？所以最好不要优化，也许不优化能很快收敛呢！ 另外如果非极化计算很容易收敛，那么以非极化优化出来的结构作为初结构，再进一步的使用自旋极化来优化，也许是个不错的主意。 另外你模型的原子层数可能太小了，模型的意义可能是个问题，但这和你问的问题没有什么关系，属于题外话了。 |
2楼2010-06-05 00:19:47
163.com
金虫 (著名写手)
实习版主
- 应助: 2 (幼儿园)
- 贵宾: 0.016
- 金币: 1489.8
- 散金: 1347
- 红花: 32
- 帖子: 1215
- 在线: 87小时
- 虫号: 923817
- 注册: 2009-12-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2010-06-05 06:29:45
★ ★ ★
xiaohunhun(金币+3):谢谢 2010-06-08 23:49:18
xiaohunhun(金币+3):谢谢 2010-06-08 23:49:18
|
引用cenwanglai的测试结果： “对于含有未配对电子的结构,计算能量时是否考虑自旋对结果影响不可忽略,而且最好在结构优化时就考虑自旋最好;但是作为计算时间成本的考虑,在结构优化时不考虑自旋,而在后续能量计算时考虑自旋可能对能量精度影响不大;但是是否在结构优化时考虑自旋对结构有一定影响.” 详见：http://muchong.com/bbs/viewthread.php?tid=1532270&fpage=0&view=&highlight=&page=2 本人觉得还是结构优化时就加上自旋比较好。至于先不加自旋优化，再此结构上再加自旋优化的方法，我正在试，希望能收敛吧,呵呵。 |
8楼2010-06-07 09:57:47
zhenfeili
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 73.6
- 散金: 1
- 帖子: 144
- 在线: 57.8小时
- 虫号: 922736
- 注册: 2009-12-08
- 专业: 生物电子学与生物信息处理
9楼2010-06-07 21:03:52
4楼2010-06-06 11:31:51
5楼2010-06-06 11:48:59
6楼2010-06-06 17:51:07
zhenfeili
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 73.6
- 散金: 1
- 帖子: 144
- 在线: 57.8小时
- 虫号: 922736
- 注册: 2009-12-08
- 专业: 生物电子学与生物信息处理
7楼2010-06-06 21:33:06
10楼2010-06-08 10:03:47
