24小时热门版块排行榜    

查看: 4021  |  回复: 27
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

kattygirl

铁虫 (初入文坛)

[求助] 计算Cl-在Fe表面吸附问题时体系加电荷方法

计算Cl-离子在Fe表面吸附时,我的算法是对cl原子加一个负电荷后,对整个体系赋予-1价,几何优化、能量计算后、结果计算出的吸附能为正值。
想请大侠们帮忙解决一下,这是哪方面的问题,会导致最后的吸附能为正值。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★
sunyang1988(金币+2): 谢谢交流 2011-12-09 23:28:45
我认为有两点值得注意吧
首先,非电中性体系计算各种能量差的时候要特别小心。最好能详细给出你所采用的公式,以及每一项是怎么计算的,才能判断计算过程中是否有错误。
另一方面这种带电体系模型的意义,就是说你这里的模型是否应该带电,也是需要慎重考虑的。
2楼2011-12-09 02:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

★ ★ ★
sunyang1988(金币+3): 谢谢指点 2011-12-11 09:26:08
引用回帖:
3楼: Originally posted by kattygirl at 2011-12-10 21:10:23:
非常感谢您的帮助,首先我用的计算公式是:吸附能=(能量计算后)吸附体系的总能量—清洁Fe(100)的能量—Cl-的能量,计算Cl-的能量时我是把cl原子加一个负电荷后放到10埃的立方盒子里面,几何优化、能量计算后得 ...

第一个,关于公式,最好能有相关的参考文献对比。因为比较带电体系能量的时候需要进行能级的alignment,以使得能量具有可比性。这个过程在计算带电缺陷形成能的时候经常出现,不知道研究带电吸附的时候有没有类似的过程,找下相关文章吧。
第二个问题。首先,因为有些时侯离子是一种比较化学的、定性的描述,而实际体系中并不一定是这样的带电原子形式存在的。所以要先明确你模拟的体系是否真的是这种带负电的Cl原子参与。
另外就是模型的意义和合理性。我们想要模拟的是无限固体Fe表面的吸附,但由于计算的限制,我们只能建立吸附位置附近的结构模型。如果此时如果你让模型带电,那么其实模拟的是:电荷局域在吸附位置附近的情况。而实际上,即使真的有额外电荷引入,可能也不会局限在吸附位置附近,尤其这里的Fe还是导体。那么你的模型的合理性就有待商榷了。
以上只是我的看法了,我也没有具体研究过吸附,只当做建议就好了。
4楼2011-12-11 00:48:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
uuv2010(金币+2): 多谢指导! 2011-12-15 23:49:57
首先说,MS是个大杂烩,虽然很方便但也容易出现误会。比如你们讨论的这个设置原子电荷的功能,并不一定是针对dmol或castep模块的设置。所以在使用这些设置之前,一定要看好帮助文件中的介绍。
另外关于体系中的离子,其实第一原理就是寻找体系基态的过程,至于体系中的到底是离子还是原子,则需要看计算出来的电子得失情况。带电体系也是一样,就比如楼主的离子,cl原子吸附到Fe,你对体系加了一个电子,那么软件就会自动寻找这个电子最稳定分布。有可能是分布在cl原子上,使之成为离子。亦有可能是分布到Fe-slab上。究竟是那种情况,自然是要听从软件的计算结果。
限制这个电子停留在Cl上,我不知道意义何在,因为在实际的实验中你也不能保证这个电子始终停留在Cl上而没有转移到别处吧~
14楼2011-12-15 23:34:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

acridine

木虫 (正式写手)

引用回帖:
15楼: Originally posted by kattygirl at 2011-12-15 23:57:34:
我的意思是刚开始建模的时候对cl原子加一个负电荷,那经过结构优化和能量计算以后,这个电荷要重新分布,可能是cl失去一些电子而Fe得到了一些电子,那么cl和Fe表面是有相互作用的,或成键之类的。但是,现在困扰 ...

你用的是什么模块?
如果是dmol和castep的话,在设置计算参数的时候可以设置体系净电荷
16楼2011-12-16 00:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kattygirl 的主题更新
信息提示
请填处理意见