24小时热门版块排行榜    

查看: 2071  |  回复: 19

luochuiyi90

木虫 (正式写手)

引用回帖:
8楼: Originally posted by 冲天飞龙 at 2015-10-01 23:10:52
按道理来说,确实会即使设置不同的initial spin也都能优化到相同的结果,你看一看算好的.castep文件的末尾有没有说优化成功。然后分别看下自旋向上和自旋向下的电子分别填充在哪个能级上,castep会分别列出来的,自 ...

谢谢你。我仔细看了计算好的文件,自旋设置为1和2时都会显示
BFGS : Geometry optimization completed successfully.这应该表明我的优化成功了吧?至于说自旋向上与向下的电子填充在那个能级上,我好像没有看到呢?可能是我不会看结果。自旋向上与向下电子数之差是指2*Integrated Spin Density  与2*Integrated |Spin Density| 之差吗?
对于Ni的核外电子排布,如果3d轨道半满的话是不是应该设置Ni自旋spin为5,spin state 为high ,Direction 为up 呢?请指教。另外谢谢你推荐的网站。我把我算的.castep文件传上来,请你帮忙看一下,好吗?

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : spin=2.txt
  • 2015-10-02 19:01:01, 106.49 K
  • 附件 2 : 新建文本文档.txt
  • 2015-10-02 19:01:47, 106.49 K
努力提升自己ing
11楼2015-10-02 19:01:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luochuiyi90

木虫 (正式写手)

引用回帖:
10楼: Originally posted by 冲天飞龙 at 2015-10-02 18:21:25
对,如果是总电子数是奇数,即使没勾选spin polarized,castep算的时候也会自动考虑自旋,这是castep的自动纠错功能,你要证实的话就检测一下你的*.castep文件,看看里面会有这么一段
************************* ...

谢谢你,那请问我在Energy计算的时候应该要不要设置自旋呢?
努力提升自己ing
12楼2015-10-02 19:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冲天飞龙

木虫 (小有名气)

引用回帖:
12楼: Originally posted by luochuiyi90 at 2015-10-02 19:05:57
谢谢你,那请问我在Energy计算的时候应该要不要设置自旋呢?...

勾选上spin比较保险
13楼2015-10-02 19:11:35
已阅   回复此楼   关注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的回帖

luochuiyi90

木虫 (正式写手)

送红花一朵
引用回帖:
14楼: Originally posted by 冲天飞龙 at 2015-10-02 20:30:30
大概看了下,你的两个文件算出来的结果是一样的,你可以自己比较下。
我说的那个看电子填充在哪个能级上的事情,我是看的*.bands这个文件,你在文件夹里找找应该有。这个文件会列出这种东西
Fermi energy (in at ...

真的谢谢你!不好意思,没有及时回。我找全了所有结果文件夹,还是没有找到呢,会不会是版本的问题,我用的是MS7.0,还是优化的时候有的选择项没有勾上呢。另外那两个结果好像也不同,spin设置为1
时:**************************** Electronic Parameters
  
number of  electrons                           :  25.00   
net charge of system                           :  0.000   
net spin   of system                           :  1.000   
number of  up  spins                           :  13.00   
number of down spins                           :  12.00   
treating system as spin-polarized
number of bands                                :         13

且2*Integrated Spin Density   =     1.00000   
2*Integrated |Spin Density| =     1.14649
spin设置为2时:
**************************** Electronic Parameters
  
number of  electrons                           :  25.00   
net charge of system                           :  0.000   
net spin   of system                           :  2.000   
number of  up  spins                           :  13.50   
number of down spins                           :  11.50   
treating system as spin-polarized
number of bands                                :         14
且2*Integrated Spin Density   =     2.00000   
2*Integrated |Spin Density| =     2.10031
求指教。
现在在考虑你说的bands设置得有点少的问题,准备再算下。
努力提升自己ing
15楼2015-10-02 22:31:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冲天飞龙

木虫 (小有名气)

引用回帖:
15楼: Originally posted by luochuiyi90 at 2015-10-02 22:31:04
真的谢谢你!不好意思,没有及时回。我找全了所有结果文件夹,还是没有找到呢,会不会是版本的问题,我用的是MS7.0,还是优化的时候有的选择项没有勾上呢。另外那两个结果好像也不同,spin设置为1
时:********** ...

正常的Geometry Optimization task就会生成band文件的,你看我给你截屏的我的文件夹
你上传的那个spin=2的文件和你刚粘贴的那几行不一样
MS勾选自旋优化失败
Capture.JPG

16楼2015-10-02 22:41:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冲天飞龙

木虫 (小有名气)

引用回帖:
15楼: Originally posted by luochuiyi90 at 2015-10-02 22:31:04
真的谢谢你!不好意思,没有及时回。我找全了所有结果文件夹,还是没有找到呢,会不会是版本的问题,我用的是MS7.0,还是优化的时候有的选择项没有勾上呢。另外那两个结果好像也不同,spin设置为1
时:********** ...

我这边最近估计比较忙,你自己看看我给你回复的研究一下吧?重要的地方我都跟你讲过了。
另外你看下我给你的附件,我算了几个例子在里面,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朵)

17楼2015-10-02 23:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luochuiyi90

木虫 (正式写手)

引用回帖:
16楼: Originally posted by 冲天飞龙 at 2015-10-02 22:41:41
正常的Geometry Optimization task就会生成band文件的,你看我给你截屏的我的文件夹
你上传的那个spin=2的文件和你刚粘贴的那几行不一样

Capture.JPG
...

真的不好意思,我粘贴的时候出错了,,我这也太马虎了,,
努力提升自己ing
18楼2015-10-03 08:55:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luochuiyi90

木虫 (正式写手)

送红花一朵
引用回帖:
17楼: Originally posted by 冲天飞龙 at 2015-10-02 23:27:13
我这边最近估计比较忙,你自己看看我给你回复的研究一下吧?重要的地方我都跟你讲过了。
另外你看下我给你的附件,我算了几个例子在里面,initial spin=1,2,5的,以及还有一个是3倍原胞大小的计算,所有计算的结论 ...

真的谢谢你,除了说谢谢都不知道说啥了
努力提升自己ing
19楼2015-10-03 08:57:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

冲天飞龙

木虫 (小有名气)

引用回帖:
19楼: Originally posted by luochuiyi90 at 2015-10-03 08:57:13
真的谢谢你,除了说谢谢都不知道说啥了...

不客气~
大家都是从不懂慢慢学习的

发自小木虫IOS客户端
20楼2015-10-04 04:53:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luochuiyi90 的主题更新
信息提示
请填处理意见