版块导航
正在加载中...
客户端APP下载
论文辅导
申博辅导
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求,自2017年10月1日起,未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用,请尽快对帐号进行手机号验证,感谢您的理解与支持!
24小时热门版块排行榜
>
论坛更新日志
(1331)
>
考博
(47)
>
论文投稿
(35)
>
导师招生
(28)
>
公派出国
(28)
>
硕博家园
(26)
>
论文道贺祈福
(20)
>
教师之家
(18)
>
虫友互识
(18)
>
考研
(16)
>
博后之家
(14)
>
找工作
(13)
>
休闲灌水
(12)
>
基金申请
(11)
>
文献求助
(8)
>
招聘信息布告栏
(7)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
第一性原理
»
MS
»
【表面模拟】castep计算中遇到的几个问题
9
1/1
返回列表
查看: 1694 | 回复: 8
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
【奖励】
本帖被评价6次,作者csfn增加金币
6
个
当前主题已经存档。
csfn
荣誉版主
(知名作家)
1ST强帖: 1
应助: 10
(幼儿园)
贵宾: 28.97
金币: 20056.9
帖子: 9522
在线: 1047.1小时
虫号: 433110
[
资源
]
【表面模拟】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要固定的话,肯定在你计算的时候你选择的任务是只算能量,我估计你应该是选择了结构优化,所以我想你是没固定原子位置的。
第三个问题:
选择自旋极化,推荐采用默认设置。正常情况下,都会收敛到比较合理的值。除非自己了解很透彻,再自定义。
回复此楼
» 收录本帖的淘帖专辑推荐
模拟知识汇总
» 猜你喜欢
真诚求助:手里的省社科项目结项要求主持人一篇中文核心,有什么渠道能发核心吗
已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
论文投稿,期刊推荐
已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
孩子确诊有中度注意力缺陷
已经有14人回复
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
请问有评职称,把科研教学业绩算分排序的高校吗
已经有5人回复
2025冷门绝学什么时候出结果
已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关商家推荐:
(我也要在这里推广)
高级回复
1楼
2007-12-03 20:12:43
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
zhanglang51
木虫
(正式写手)
应助: 1
(幼儿园)
金币: 3104.8
帖子: 404
在线: 378.3小时
虫号: 352332
★★★★★ 五星级,优秀推荐
不错,很清楚。我也是新手,望多多关照,不吝赐教。
赞
一下
回复此楼
2楼
2007-12-04 14:59:34
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
fegg7502
荣誉版主
(知名作家)
应助: 20
(小学生)
贵宾: 13.133
金币: 18796.3
帖子: 7483
在线: 590小时
虫号: 352334
★★★★★ 五星级,优秀推荐
楼主,我支持你,呵呵!有空多奖励给我几个金币啊,哈哈!!
赞
一下
回复此楼
3楼
2007-12-04 15:56:08
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
buble
金虫
(著名写手)
应助: 1
(幼儿园)
金币: 3126.6
帖子: 1301
在线: 943.6小时
虫号: 467890
太好了,我又多了解了一些
赞
一下
回复此楼
4楼
2007-12-08 16:56:58
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
cometring
木虫
(著名写手)
应助: 0
(幼儿园)
贵宾: 1.03
金币: 4315.7
帖子: 2911
在线: 36.1小时
虫号: 195508
★★★★★ 五星级,优秀推荐
对于问题1,我试了以下,对于减小smearing影响不大,基本上还是出这个提示,只有增加空轨道数量,很有用,按照他提示的设定,就可以消除这个提示
赞
一下
回复此楼
5楼
2007-12-08 19:03:06
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
378287892
金虫
(小有名气)
应助: 0
(幼儿园)
金币: 217.4
帖子: 296
在线: 136.3小时
虫号: 477317
非常感谢!!!第一个问题和我的一摸一样,一直不知道怎么回事,这次终于知道了,本人新手,望多多关照!再次谢谢楼上高手!
赞
一下
回复此楼
6楼
2008-06-10 21:34:06
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
dongping001
木虫
(正式写手)
应助: 0
(幼儿园)
金币: 2018.4
帖子: 519
在线: 3.6小时
虫号: 140724
★★★★★ 五星级,优秀推荐
多多益善
回复此楼
7楼
2008-06-10 22:56:55
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
aylayl08
荣誉版主
(文坛精英)
1ST强帖: 1
应助: 3
(幼儿园)
贵宾: 4.825
金币: 18823.7
帖子: 16357
在线: 824小时
虫号: 533169
★★★★★ 五星级,优秀推荐
辛苦啦
回复此楼
8楼
2008-06-13 18:40:24
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
csfn
的主题更新
9
1/1
返回列表
☆ 无星级
★ 一星级
★★★ 三星级
★★★★★ 五星级
如果回帖内容含有宣传信息,请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
高级回复
(可上传附件)
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站,将文件上传后,然后将下载链接复制到帖子内容中就可以了。
信息提示
关闭
请填处理意见
关闭
确定