应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » MS » 【表面模拟】castep计算中遇到的几个问题
北京石油化工学院2026年研究生招生接收调剂公告
91/1返回列表
查看: 1918  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价6次，作者csfn增加金币 6 个
当前主题已经存档。

csfn

荣誉版主 (知名作家)

[资源] 【表面模拟】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要固定的话，肯定在你计算的时候你选择的任务是只算能量，我估计你应该是选择了结构优化，所以我想你是没固定原子位置的。
第三个问题：
选择自旋极化，推荐采用默认设置。正常情况下，都会收敛到比较合理的值。除非自己了解很透彻，再自定义。
回复此楼

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

模拟知识汇总

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2007-12-03 20:12:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhanglang51

木虫 (正式写手)

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



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

fegg7502

荣誉版主 (知名作家)

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

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

buble

金虫 (著名写手)
太好了，我又多了解了一些
赞一下 回复此楼
4楼2007-12-08 16:56:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cometring

木虫 (著名写手)

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

对于问题1，我试了以下，对于减小smearing影响不大，基本上还是出这个提示，只有增加空轨道数量，很有用，按照他提示的设定，就可以消除这个提示
赞一下 回复此楼
5楼2007-12-08 19:03:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

378287892

金虫 (小有名气)
非常感谢！！！第一个问题和我的一摸一样，一直不知道怎么回事，这次终于知道了，本人新手，望多多关照！再次谢谢楼上高手！
赞一下 回复此楼
6楼2008-06-10 21:34:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dongping001

木虫 (正式写手)

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

多多益善
回复此楼
7楼2008-06-10 22:56:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aylayl08

荣誉版主 (文坛精英)

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

辛苦啦
回复此楼
8楼2008-06-13 18:40:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 csfn 的主题更新
91/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 263求调剂 +3 李nihao 2026-04-08 3/150 2026-04-08 05:40 by 挽忆欢
[考研] 288求调剂 +12 没有答案_ 2026-04-05 12/600 2026-04-08 00:17 by T可可西里T
[考研] 调剂 +4 mcbbc 2026-04-06 5/250 2026-04-07 12:33 by upczlm1989
[考研] 295求调剂 +18 xndjjj 2026-04-04 19/950 2026-04-07 11:02 by wangjy2002
[考研] 338求调剂 +4 我想上岸ii 2026-04-05 4/200 2026-04-06 21:04 by 木子君1218
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 315求调剂 +5 ＆123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 322求调剂 +3 嗯哼哼恒 2026-04-05 3/150 2026-04-05 19:52 by nepu_uu
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 324求调剂 +9 想上学求调 2026-04-03 9/450 2026-04-04 23:57 by 果冻大王
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
[考研] 334求调剂 +8 曾仰之 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 400分求调剂 +3 尴尬且挠头 2026-04-04 3/150 2026-04-04 08:41 by jp9609
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 数一英一285求调剂 +7 AZMK 2026-04-03 9/450 2026-04-03 13:03 by ms629
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 一志愿山东大学，085600，344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 283求调剂 +3 jiouuu 2026-04-02 4/200 2026-04-02 14:08 by 哒哒哒呱呱呱
信息提示
关闭
请填处理意见
关闭