24小时热门版块排行榜    

查看: 2302  |  回复: 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

木虫 (著名写手)

爱你就等于爱自己

我今天看了看帮助文件,觉得第一个问题可能跟我没有设置initial spin有关系,我计算的时候,initial spin 设置的是0.
不知道这个参数该如何设置?请指教!
与超版斗,其乐无穷;与斑竹斗,其乐无穷;与虫子斗,其乐无穷。
4楼2007-11-28 20:42:54
已阅   回复此楼   关注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的回帖

zzgyb

荣誉版主 (文坛精英)

小木虫警察局局长

优秀版主

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
小红豆(金币+4,VIP+0):3Q!辛苦了!
cometring(金币+50,VIP+0):谢谢你多次帮忙!!
第一个文标题: 从他提示的关键词来看
NEXTRA_BANDS
This keywords controls the number of extra bands in addition to the number of occupied bands. These extra bands are necessary for metals or finite temperature insulators. The default value for this parameter is 0.

SMEARING_WIDTH
This keyword determines the width of the Fermi-surface smearing if the system is being treated as a metal.

因为不知道你的具体物质,所以我猜测有两种情况,一是你的结构有问题。二是因为你的结构的特殊性,所以你默认的设置计算是不可信的,就如提示而言。如果你的结构没问题,你可以按照提示调整一下参数算一算。

.param为后缀名的文件里,有以上两个参数,你打开调整一下,然后计算你调整后的文件,就是点Files,然后Run Files。如果你调整后没有提示,就说明你的结构需要调整参数。

第二个问题

原子在表面吸附就那么几种位置,你都要计算试一试,找到能量最小的吸附位置。至于吸附原子放在什么位置,你可以通过原子半径计算,这样计算能加快你的收敛速度。如果不是计算,而是大致放个位置,收敛就和你放的位置有关,如果放的位置很接近收敛位置,就会很快收敛。因为这这两种方法最终的计算收敛后,他们的位置应该是一样的(误差范围之内)。至于计算表面吸附total energy的时候,是选择"optimized cell"还是固定表面,我看文献上这两种都有,因为这要和你计算的体系有关,如果体系不是很大,应该选择optimized cell,这样更近现实验情况。

第三个问题,initial Spin设置,MS提供两种方法,你可以参考。
There are two ways of defining the initial magnetic configuration: either specify the total magnetic moment per unit cell, which gets uniformly distributed over the space, or provide detailed information on the absolute values and direction (up or down) of the spins for each atom in the unit cell. The former method can be used for relatively simple systems where only two solutions are expected (magnetic and non-magnetic). The latter method, which specifies the spin state of the atoms in the system, is more general and gives much more flexibility. It is possible to set up ferromagnetic, ferrimagnetic, or antiferromagnetic calculations to get different starting spin arrangements.
不过我记得好像他们公司说过,一般formal spin都可以,除非是计算一些磁性物质,比如Fe、Co这种。目前我都是默认。至于charge就是系统的电荷,如果不带电就是0,实际是多少就是多少。

自己的一点拙见,仅供参考
江南七怪+马玉+洪七公+一灯大师+周伯通教一个弱智儿童=郭靖,天下第一王重阳教七个天资不错的小孩=全真七子,全真七子×3《郭靖,这就是盲目扩招的后果
5楼2007-11-28 21:18:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见