| 查看: 2489 | 回复: 19 | |||
luochuiyi90木虫 (正式写手)
|
[求助]
MS勾选自旋优化失败
|
||
各位虫友，如图，我用CASTEP几何优化镍酸锂，考虑电子自旋，所以勾选了spin polarized ,并且根据帖子里面说的把initial spin 设置为5，其他的设置均未改变。在未考虑自旋情况下优化没有什么问题，但是考虑自旋就会显示：Warning: electronic minimisation did not converge when finding ground state.然后优化失败。小弟刚学MS，导师也不懂，也没师兄姐可以请教，哪位高手可以帮帮我，小弟感激不尽！ jietu2.jpg |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
我这边最近估计比较忙，你自己看看我给你回复的研究一下吧？重要的地方我都跟你讲过了。 另外你看下我给你的附件，我算了几个例子在里面，initial spin=1,2,5的，以及还有一个是3倍原胞大小的计算，所有计算的结论都是每个原胞的net spin是1。你比较下你的参数有没有设置的不对的地方 |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : LiNiO2.zip
2015-10-02 23:22:59, 58.15 K
» 本帖已获得的红花（最新10朵）
luochuiyi9017楼2015-10-02 23:27:13
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
castep的默认设置是，如果检测到体系总电子数是奇数，就会自动进行自旋优化。 你这个体系是25个电子所以得到的结果已经是考虑了自旋的结果了。 如果你设置initial spin，是参考了别人的结论，可以看下别人的cell是怎么设置的，至少目前看来这个单胞得到net spin=1还是比较正常的，毕竟总电子数是奇数。 建议做下大的晶胞，比如2个，4个单胞大小的，看看这样设置后总电子数变成偶数后，还会不会有net spin。 另外后面如果需要考虑磁学性质的话，可以加上spin orbital coupling看看，vasp有这个功能，castep好像最新的版本有 |
» 本帖已获得的红花（最新10朵）
7楼2015-10-01 22:46:16
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2015-09-30 18:37:06
jimsmart
捐助贵宾 (职业作家)
- 应助: 109 (高中生)
- 贵宾: 0.2
- 金币: 3583.7
- 散金: 100
- 红花: 5
- 帖子: 3110
- 在线: 568.1小时
- 虫号: 57301
- 注册: 2005-01-29
- 性别: GG
- 专业: 光谱分析
3楼2015-10-01 18:30:54
luochuiyi90
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5624.3
- 散金: 1255
- 红花: 2
- 帖子: 869
- 在线: 157.1小时
- 虫号: 3515475
- 注册: 2014-11-03
- 性别: GG
- 专业: 其他无机非金属材料
|
谢谢你的回复。之所以设置为5，是因为看到有的帖子说参数可以设置大点，计算的时候会自行调整。但是，我发现设置为2或者1，对应的 2*Integrated Spin Density 与2*Integrated |Spin Density| 值会发生变化，比如，设置为5时： 2*Integrated Spin Density = 5.00000 2*Integrated |Spin Density| = 5.00000 ，但是，设置为2时： 2*Integrated Spin Density = 2.00000 2*Integrated |Spin Density| = 2.10030 设置为1时：2*Integrated Spin Density = 1.00000 2*Integrated |Spin Density| = 1.14647，不同的值会有不同的结果。另外LiNiO2中Ni离子为+3价，核外电子排布应该是3d74s0，所以我个人觉得是不是应该在Modify里面将Ni Spin设置为3呢？期待您的指教。金币不够的话我会追加，谢谢！ |
4楼2015-10-01 20:05:40
luochuiyi90
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5624.3
- 散金: 1255
- 红花: 2
- 帖子: 869
- 在线: 157.1小时
- 虫号: 3515475
- 注册: 2014-11-03
- 性别: GG
- 专业: 其他无机非金属材料
5楼2015-10-01 20:07:28
luochuiyi90
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5624.3
- 散金: 1255
- 红花: 2
- 帖子: 869
- 在线: 157.1小时
- 虫号: 3515475
- 注册: 2014-11-03
- 性别: GG
- 专业: 其他无机非金属材料
6楼2015-10-01 20:50:51
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
按道理来说，确实会即使设置不同的initial spin也都能优化到相同的结果，你看一看算好的.castep文件的末尾有没有说优化成功。然后分别看下自旋向上和自旋向下的电子分别填充在哪个能级上，castep会分别列出来的，自旋向上和自旋向下的电子数之差就是优化好的net spin。 关于Ni的外层电子，+3价的我觉得排成3d5,4s2似乎更合理，d轨道半满，s轨道全满能量应该比较低。不过Ni的金属性不强，和O的电负性差别没那么大，所以Ni和O之间可能是共价性的吧，所以Ni和O之间会公用电子吧，这一点我也不太肯定。 另外推荐个网站 https://www.materialsproject.org/ 我刚才看了，上面有每种LiNiO2晶体构型人家算的磁性Final Magnetic Moment，不同的晶体结构算出来的磁性也不同，1到5的都有 |
8楼2015-10-01 23:10:52
luochuiyi90
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 5624.3
- 散金: 1255
- 红花: 2
- 帖子: 869
- 在线: 157.1小时
- 虫号: 3515475
- 注册: 2014-11-03
- 性别: GG
- 专业: 其他无机非金属材料
9楼2015-10-02 17:19:48
冲天飞龙
木虫 (小有名气)
- 1ST强帖: 2
- 应助: 70 (初中生)
- 金币: 1888.6
- 红花: 9
- 帖子: 239
- 在线: 61.7小时
- 虫号: 1944534
- 注册: 2012-08-18
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★
luochuiyi90: 金币+5, ★★★很有帮助 2015-10-02 19:04:02
luochuiyi90: 金币+5, ★★★很有帮助 2015-10-02 19:04:02
|
对，如果是总电子数是奇数，即使没勾选spin polarized，castep算的时候也会自动考虑自旋，这是castep的自动纠错功能，你要证实的话就检测一下你的*.castep文件，看看里面会有这么一段 **************************** Electronic Parameters **************************** number of electrons : 348.0 net charge of system : 0.000 net spin of system : 0.000 number of up spins : 174.0 number of down spins : 174.0 treating system as non-spin-polarized 比如我这里是non-spin-polarized，就是算的时候没自旋。 只要这里写的是 treating system as spin-polarized就是算的时候考虑了自旋。 |
10楼2015-10-02 18:21:25
