24小时热门版块排行榜    

查看: 2308  |  回复: 21
当前主题已经存档。
【有奖交流】积极回复本帖子,参与交流,就有机会分得作者 cometring 的 25 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cometring

木虫 (著名写手)

爱你就等于爱自己

[交流] 请教关于MS castep的几个问题

希望各位大侠帮忙,
1 计算表面过程中,计算成功,但是castep结果中提示

------------------------------------------------------------------------ <-- SCF
SCF loop      Energy           Fermi           Energy gain       Timer   <-- SCF
                               energy          per atom          (sec)   <-- SCF
------------------------------------------------------------------------ <-- SCF
Initial  -1.24983546E+004  5.76278315E+001                        55.44  <-- SCF
      1  -1.59703259E+004 -1.69750335E+000   2.89330938E+002     197.97  <-- SCF
Warning: There are no empty bands for at least one kpoint and spin; this may slow the convergence and/or lead to an inaccurate groundstate. If this warning persists, you should consider increasing nextra_bands and/or reducing smearing_width in the param file. Recommend using nextra_bands of 11 to 24.                             
  
      2  -1.62762609E+004 -6.23532889E+000   2.54945870E+001     356.62  <-- SCF
Warning: There are no empty bands for at least one kpoint and spin; this may slow the convergence and/or lead to an inaccurate groundstate. If this warning persists, you should consider increasing nextra_bands and/or reducing smearing_width in the param file.Recommend using nextra_bands of 11 to 24.   



2  要在表面上吸附小分子,如H2 ,或者CO,可以选择不同的位置,比如top, bridge, 等,但是如何确定具体每个原子的position? 直接根据原子半径按照几何知识来计算吗?另外在计算表面吸附total energy的时候,是否选择上"optimized cell" ?
如果不选择,就应该是每个原子的位置保持固定,那么如何知道所确定的吸附原子的位置是对的呢?

3 进行计算时选择上"spin polarized", 对于不同元素,其参数如何设置?"initial pin"和"charge"。
            

请各位高手不吝赐教!金币不够可以再加

[ Last edited by zzgyb on 2007-12-3 at 14:12 ]
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

与超版斗,其乐无穷;与斑竹斗,其乐无穷;与虫子斗,其乐无穷。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cometring

木虫 (著名写手)

爱你就等于爱自己

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
vzhang77(金币+10,VIP+0):多谢您的好帖子及好总结!:)
总结一下:
问题1
基本上是找到了问题,而且必须解决,通过软件的提示和大家的讨论,可以增加bands或者减小smearing,我会先试试减小smearing,看看结果如何

问题 2
看来大家大多数认为在吸附的relax过程中不选择optimized cell,原子的位置靠几何计算得到。我先这样试试,但是我还有一点疑虑,就是如果不选择optimized cell,吸附分子或者原子的位置是不是就保持固定了呢?csfn认为位置也不是保持固定,大家有经验的可以谈谈。我算完之后其实测量一下就知道到底固定不固定了:)

问题3
争议还比较大,
1 不设置,保持默认,但是默认的initial spin是0,这样不知道对不对?
2 设置spin=unpair electron+1,这个不知道是什么意思?有没有具体文献。
3 按照帮助文件设置,但是要考虑计算中超胞中很多原子是共用。

请大家继续讨论
与超版斗,其乐无穷;与斑竹斗,其乐无穷;与虫子斗,其乐无穷。
14楼2007-11-29 10:58:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 22 个回答

hufeng

新虫 (正式写手)


小红豆(金币+1,VIP+0):3Q!
你把empty band里的空带数目增加试试或者减少smearing 的数值,这两个办法都可以增加精度。
2楼2007-11-28 20:33:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hufeng

新虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小红豆(金币+1,VIP+0):3Q!
cometring(金币+10,VIP+0):谢谢~
你先把晶胞优化完成后再算吸附,只是注意两次计算所用方法要一致。确定原子的position是要根据能量高低来确定。
3楼2007-11-28 20:37:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cometring

木虫 (著名写手)

爱你就等于爱自己

我今天看了看帮助文件,觉得第一个问题可能跟我没有设置initial spin有关系,我计算的时候,initial spin 设置的是0.
不知道这个参数该如何设置?请指教!
与超版斗,其乐无穷;与斑竹斗,其乐无穷;与虫子斗,其乐无穷。
4楼2007-11-28 20:42:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见