24小时热门版块排行榜    

查看: 2488  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

luochuiyi90

木虫 (正式写手)

[求助] MS勾选自旋优化失败

各位虫友,如图,我用CASTEP几何优化镍酸锂,考虑电子自旋,所以勾选了spin polarized ,并且根据帖子里面说的把initial spin 设置为5,其他的设置均未改变。在未考虑自旋情况下优化没有什么问题,但是考虑自旋就会显示:Warning: electronic minimisation did not converge when finding ground state.然后优化失败。小弟刚学MS,导师也不懂,也没师兄姐可以请教,哪位高手可以帮帮我,小弟感激不尽!

MS勾选自旋优化失败
jietu2.jpg
回复此楼
努力提升自己ing
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冲天飞龙

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
luochuiyi90: 金币+5 2015-10-02 21:53:22
luochuiyi90: 金币+10 2015-10-03 08:59:18
引用回帖:
11楼: Originally posted by luochuiyi90 at 2015-10-02 19:01:54
谢谢你。我仔细看了计算好的文件,自旋设置为1和2时都会显示
BFGS : Geometry optimization completed successfully.这应该表明我的优化成功了吧?至于说自旋向上与向下的电子填充在那个能级上,我好像没有看到呢 ...

大概看了下,你的两个文件算出来的结果是一样的,你可以自己比较下。
我说的那个看电子填充在哪个能级上的事情,我是看的*.bands这个文件,你在文件夹里找找应该有。这个文件会列出这种东西
Fermi energy (in atomic units)     0.101340
Unit cell vectors
   18.838680    0.000000    0.000000
    0.000000   16.314775    0.000000
    0.000000    0.000000   21.345212
K-point     1  0.00000000  0.25000000  0.00000000  1.00000000
Spin component 1
   -0.69947963
   -0.68700675
   -0.68338554
   -0.68239331
   -0.68237069
   -0.68236994
   -0.67949367
   -0.67837160
   -0.67837103
每个k点和每个spin方向的能级都会列出来,而且文件头会列出fermi level。在fermi level能量之下的就是占据轨道。
至于电子排布和initial spin的问题,不必太过纠结,我之前说了,Ni和O直接电负性差别没那么大,所以可能电子对是共用的,说不清的。价键理论本身就是用来简单分析的模型,本质上原子的外层轨道还是会互相杂化形成新轨道的。所以还是看castep算出来的优化好的net spin是多少就可以了,初始值设成1还是3还是其他值按道理都可以优化到相同结果的。
最后给个建议,你的bands设置的有点少,25个电子只用了13条bands。最好在MS里面不要勾选fix occupancy,extra bands设成20%或者30%,保证有几条空的能带比较好

» 本帖已获得的红花(最新10朵)

14楼2015-10-02 20:30:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

冲天飞龙

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
luochuiyi90: 金币+6 2015-10-01 20:05:57
首先我没接触过这个材料,所以说的不一定对
net spin真的是5吗?这么大!能不能介绍一下这种材料的性质,为什么spin会是5,以及你算的时候用的晶胞大小?
此外,由于总电子数是偶数个,如果spin是5的话,那岂不是自旋向上以及自旋向下的电子数都不是整数了
2楼2015-09-30 18:37:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jimsmart

捐助贵宾 (职业作家)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
luochuiyi90: 金币+4 2015-10-01 20:06:19
如果报错的话, 应该是参数设置还是有不对的地方, 改到它不报错为止.
3楼2015-10-01 18:30:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luochuiyi90

木虫 (正式写手)

引用回帖:
2楼: Originally posted by 冲天飞龙 at 2015-09-30 18:37:06
首先我没接触过这个材料,所以说的不一定对
net spin真的是5吗?这么大!能不能介绍一下这种材料的性质,为什么spin会是5,以及你算的时候用的晶胞大小?
此外,由于总电子数是偶数个,如果spin是5的话,那岂不是 ...

谢谢你的回复。之所以设置为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呢?期待您的指教。金币不够的话我会追加,谢谢!
努力提升自己ing
4楼2015-10-01 20:05:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 7lpolszZVXgi 2026-08-14 6/300 2026-08-15 09:28 by sunzitan
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 4/200 2026-08-15 04:17 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
信息提示
请填处理意见