24小时热门版块排行榜    

查看: 2074  |  回复: 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的回帖
信息提示
请填处理意见