24小时热门版块排行榜    

查看: 3148  |  回复: 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的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见