24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2065  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

clickme

木虫 (小有名气)

[交流] 【讨论】DMol3和CASTEP计算过渡金属表面吸附比较

上次yzcluster说用DMol3计算金属表面很不准确,建议用CASTEP。我接下来建立相同模型仔细比较了CASTEP以及DMol3计算Ni(100)面吸附乙烯分子(C2H4)的吸附能和构型,得到的结论是DMol3的精度能够符合计算要求,并且可以大幅度节约
计算时间。

模型是3层原子2x1晶格的Ni(100) slab构型,上面放置了C2H4分子,底层Ni原子固定。分别用DMol3和CASTEP进行了一步步Coarse,Medium和Fine的构型优化。

在最后计算能量的时候,DMol3的所有模型强制采用Orbital Cutoff 4.0eV, K point mesh:6x3x1其它采用Medium精度,对金属原子进行全电子计算时候,吸附能为1.083eV,换作ecp时候,吸附能为1.169eV.

CASTEP计算参数全都采用默认的Fine,强制采用Energy Cutoff 300.00eV.计算得到的吸附能为0.9812eV

CASTEP和DMol3计算中的Density Mixing Charge都根据MS帮助文档,计算金属表面吸附时候设定为小值0.01

DMol3采用全电子计算方法,与CASTEP的差值只有0.1eV。而且这是在DMol3用Medium精度,CASTEP为fine精度。而两者的计算速度相差很大。

这是我的初步的对比结果,希望大家能讨论一下
回复此楼

» 猜你喜欢

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

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

yzcluster

金虫 (著名写手)

小木虫扫盲人

★ ★ ★
zdhlover(金币+3,VIP+0):欢迎常来这里交流讨论 9-24 22:00
引用回帖:
Originally posted by clickme at 2009-9-24 21:55:
上次yzcluster说用DMol3计算金属表面很不准确,建议用CASTEP。我接下来建立相同模型仔细比较了CASTEP以及DMol3计算Ni(100)面吸附乙烯分子(C2H4)的吸附能和构型,得到的结论是DMol3的精度能够符合计算要求,并且可 ...

第一,你的slab太薄了,无论模拟什么材料的表面,用3层模拟肯定是不行的。其结果未必可靠。

第二,我上次说DMOL计算周期性的结果不好,只是大体上说的(比如,曾经有一个含Mg的材料的计算例子,结果居然连Mg的3d电子都明显参与成键,显然不太符合实际,我试了试,用CASTEP就不会出现这种情况)。而且,很多固体的性质是DMOL无法计算的,比如NMR等等。但是,对于某个具体的体系,究竟哪个软件更合适,就不好说了。这个你要试一试。
3楼2009-09-25 11:20:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

clickme

木虫 (小有名气)


zdhlover(金币+1,VIP+0):欢迎常来这里交流讨论 9-24 22:00
DMol3的基集函数用的是DNP,3.5版的
2楼2009-09-25 08:50:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

clickme

木虫 (小有名气)

谢谢yzcluster老大的回复,非常感谢!

我这个模型只是用来比较用的,不是真正用拉计算吸附能的。
另外问一下,如何在DMol3种看到参与成键的电子?是进行轨道分析么?
4楼2009-09-25 16:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)

小木虫扫盲人

★ ★ ★
zdhlover(金币+3,VIP+0):欢迎常来这里交流讨论 9-24 22:01
引用回帖:
Originally posted by clickme at 2009-9-25 16:29:
谢谢yzcluster老大的回复,非常感谢!

我这个模型只是用来比较用的,不是真正用拉计算吸附能的。
另外问一下,如何在DMol3种看到参与成键的电子?是进行轨道分析么?

对于DMOL能否显示出周期性体系“分子轨道”的原子成分,我不清楚,因为我没有用DMOL算过固体方面的东西。不过这个不难,你可以试一试,用print关键字,看看能否得到“布洛赫波”是由那些原子轨道组成的。
对于键合电子,一般看看态密度和穆里肯布局分析就行了。不过穆里肯布局分析对金属体系没有物理意义。这点要小心。
5楼2009-09-25 17:01:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工学硕294分,求导师收留 +33 yzyzx 2026-04-12 37/1850 2026-04-17 23:00 by wunaiy88
[考研] 求调剂 +9 小聂爱学习 2026-04-16 11/550 2026-04-17 22:34 by chixmc
[考研] 304求调剂 +7 castLight 2026-04-16 7/350 2026-04-17 20:05 by 关一盏灯cd
[考博] 求博导|生物质基多孔碳/超级电容方向,已有相关成果,寻能源材料/碳材料方向老师 +3 猪猪人Zzz 2026-04-12 3/150 2026-04-17 19:10 by 阳阳阳^_^
[考研] 0854求调剂 +21 门路摸摸 2026-04-15 25/1250 2026-04-17 15:45 by qzxyhcsy
[考研] 322求调剂 +6 tekuzu 2026-04-17 6/300 2026-04-17 13:48 by Espannnnnol
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 322求调剂 +8 123安康 2026-04-12 15/750 2026-04-16 11:07 by Espannnnnol
[考研] 327求调剂 +26 Xxjc1107. 2026-04-13 29/1450 2026-04-16 10:52 by Espannnnnol
[考研] 药学求调剂 +14 喽哈加油 2026-04-14 16/800 2026-04-16 10:15 by beilsong20
[考研] 289 分105500药学专硕求调剂(找B区学校) +4 白云123456789 2026-04-13 4/200 2026-04-16 00:18 by 粉沁若尘
[考研] 085404 22408 309分求调剂 +9 lzmk 2026-04-14 10/500 2026-04-15 20:02 by 学员JpLReM
[考研] 药学求调剂 +11 RussHu 2026-04-12 13/650 2026-04-15 19:07 by zhuwenxu
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] 302求调剂 +10 易!? 2026-04-13 10/500 2026-04-13 19:04 by lbsjt
信息提示
请填处理意见