24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 1915  |  回复: 8
【奖励】 本帖被评价6次,作者csfn增加金币 6
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[资源] 【表面模拟】castep计算中遇到的几个问题

CODE:
3个问题:
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"。

CODE:
第一个问题:
大家回答的差不多。你的结果肯定是有问题的,必须要调试到没这个警告出现。增加bands的做法不是很可靠,最好是调整smearing,怎么调?最好的方法是测试,就是大小都调。不过你这里提示是要你减小这个数值,可以一试。
第二个问题:
1,每个位置都要算,通过能量高低判断最佳吸附位置。对于分子可能有多种组合的位置,都要计算,不过最好参考下文献,一般是在高对称点咯!推荐下,计算的时候固定xy坐标,保持高度变化就ok!
2,根绝原子半径判断,呵呵,在表面就不要考虑,如果你考虑他可能吸附到体内去,比如那里是fcc位置,有一个洞,那是可以考虑原子半径。因为半径大了,那个洞容不下他肯定不会进去咯。不过你把他放到洞里算下就知道哒。一般是用吸附能判断,吸附能为正就是可以吸附,为负不吸附。不过吸附能定义的计算方法要搞正确就是。吸附能定义为吸附前后物质总能量的变化,是结合能的反数,也就是吸附体系的能量减去清洁表面的能量,再减去被吸附物质的能量。
3,是否选择“optimized cell”
一般的体系都可以不选择,选择上了计算量增大太多。如果不选择,也不是位置固定。只是那个合子,或者说整体的体积和外形不变化而已。我是这样的理解的!
4,怎么模拟表面
一般来说计算表面吸附是不会固定下面的基地的,固定了可信度小,发出去自己也没底气。^_^一般采用的2种办法,一种是楼上说的,固定的最下面层,然后上面和吸附原子接触的就relax,这个是相当于只考虑一个表面。固定的层数对半开。第二种方法是不固定层数,选择很大的层数,或者选择小一点的层数,固定中间层数,两面的原子层数不固定,这个就是有两个表面的情况。还有少数是全部固定原子层数。对于楼主来说,castep要固定的话,肯定在你计算的时候你选择的任务是只算能量,我估计你应该是选择了结构优化,所以我想你是没固定原子位置的。
第三个问题:
选择自旋极化,推荐采用默认设置。正常情况下,都会收敛到比较合理的值。除非自己了解很透彻,再自定义。

回复此楼

» 收录本帖的淘帖专辑推荐

模拟知识汇总

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buble

金虫 (著名写手)


太好了,我又多了解了一些
4楼2007-12-08 16:56:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

zhanglang51

木虫 (正式写手)


★★★★★ 五星级,优秀推荐



不错,很清楚。我也是新手,望多多关照,不吝赐教。
2楼2007-12-04 14:59:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

楼主,我支持你,呵呵!有空多奖励给我几个金币啊,哈哈!!
3楼2007-12-04 15:56:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cometring

木虫 (著名写手)


★★★★★ 五星级,优秀推荐

对于问题1,我试了以下,对于减小smearing影响不大,基本上还是出这个提示,只有增加空轨道数量,很有用,按照他提示的设定,就可以消除这个提示
5楼2007-12-08 19:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 288求调剂 +12 没有答案_ 2026-04-05 12/600 2026-04-08 00:17 by T可可西里T
[考研] 一志愿211,化学学硕,310分,本科重点双非,求调剂 +19 努力奋斗112 2026-04-07 20/1000 2026-04-07 23:12 by JourneyLucky
[考研] 085600材料与化工,求调剂 +5 won_qii 2026-04-07 5/250 2026-04-07 17:10 by 啵啵啵0119
[考研] 085602调剂 初试总分335 +10 19123253302 2026-04-05 10/500 2026-04-07 15:23 by 小乔同学ya
[考研] 372求调剂 +4 jj涌77 2026-04-02 4/200 2026-04-07 09:31 by 白云朵朵飞
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸! 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 求调剂 +4 wos666 2026-04-03 5/250 2026-04-06 15:22 by wos666
[考研] 一志愿211生物学280分 求调剂 +5 李rien 2026-04-05 5/250 2026-04-06 10:30 by zhyzzh
[考研] 319求调剂 +3 handrui 2026-04-05 3/150 2026-04-06 09:33 by jp9609
[考研] 一志愿哈尔滨工业大学085600英一数二337分求调剂 +12 lyz0427 2026-04-03 12/600 2026-04-06 06:37 by houyaoxu
[考研] 化学0703-一志愿211-338分求调剂 +8 vants 2026-04-05 8/400 2026-04-06 06:17 by houyaoxu
[考研] 一志愿北京交通大学材料工程总分358求调剂 +4 cs0106 2026-04-04 4/200 2026-04-05 18:46 by imissbao
[考研] 271分求调剂学校 +12 zph158488! 2026-04-02 13/650 2026-04-05 10:13 by lqwchd
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 324求调剂 +14 想上学求调 2026-04-02 15/750 2026-04-04 20:31 by 无际的草原
[考研] 26调剂 086003 +6 失活的细胞 2026-04-04 6/300 2026-04-04 09:50 by zhangdingwa
[考研] 调剂 +5 asdasdassda 2026-04-03 6/300 2026-04-03 20:27 by 岸上的一条鱼
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 293求调剂 +4 珂珂乐 2026-04-02 4/200 2026-04-02 20:10 by 6781022
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
信息提示
请填处理意见