24小时热门版块排行榜    

查看: 4760  |  回复: 27

kattygirl

铁虫 (初入文坛)

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

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

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的回帖

hakuna

木虫 (知名作家)

liliangfang(金币+1): 多谢指教 2012-02-21 09:15:54
都是形式电荷惹得祸,一个体系的电荷从来都不会集中在一个或几个特别原子上的.......
18楼2012-02-21 07:16:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

acridine

木虫 (正式写手)

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

kattygirl

铁虫 (初入文坛)

引用回帖:
2楼: Originally posted by acridine at 2011-12-09 02:04:33:
我认为有两点值得注意吧
首先,非电中性体系计算各种能量差的时候要特别小心。最好能详细给出你所采用的公式,以及每一项是怎么计算的,才能判断计算过程中是否有错误。
另一方面这种带电体系模型的意义,就是说 ...

非常感谢您的帮助,首先我用的计算公式是:吸附能=(能量计算后)吸附体系的总能量—清洁Fe(100)的能量—Cl-的能量,计算Cl-的能量时我是把cl原子加一个负电荷后放到10埃的立方盒子里面,几何优化、能量计算后得出的。其次,对于这种带电体系,对于加电荷的方式我不知道我这样对不对,但是我觉得为了保证电荷守恒,加一个带负电荷的氯离子在铁表面上,那么体系不是也应给再赋予给它一个负电荷吗?
希望您能指点我一下,万分感谢!
3楼2011-12-10 21:10:23
已阅   回复此楼   关注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的回帖

nwx_ky

木虫 (小有名气)


感谢参与,应助指数 +1
sunyang1988: 求助帖,请选择正确的回帖类型 2011-12-15 04:02:44
zzy870720z(金币-1, 应助指数-1): 请选择正确回帖类型,扣回应助指数并罚1BB,谢谢合作 2011-12-16 19:08:54
引用回帖:
3楼: Originally posted by kattygirl at 2011-12-10 21:10:23:
非常感谢您的帮助,首先我用的计算公式是:吸附能=(能量计算后)吸附体系的总能量—清洁Fe(100)的能量—Cl-的能量,计算Cl-的能量时我是把cl原子加一个负电荷后放到10埃的立方盒子里面,几何优化、能量计算后得 ...

你好,请问你的这个电荷是怎么加上去的。我试了很多次,MS中加不上去,不知你用的什么程序,我也再算离子吸附?可以讨论下。
梦洁
5楼2011-12-13 15:30:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kattygirl

铁虫 (初入文坛)

WDD880227(金币+1): 感谢交流 2012-02-21 08:29:26
引用回帖:
: Originally posted by nwx_ky at 2011-12-13 15:30:44:
你好,请问你的这个电荷是怎么加上去的。我试了很多次,MS中加不上去,不知你用的什么程序,我也再算离子吸附?可以讨论下。

我用的也是MS,第一种方式是:选中你要加电荷的原子修改properties中的charges。第二种方式:工具栏Modify——charges——set all charges/set total charge,视情况选择。
6楼2011-12-14 10:22:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxj_yc

金虫 (初入文坛)

求助:我也是做离子吸附的,问一下SO42-这类离子,吸附体系怎么建模呢?SO42-能量怎么计算呢?非常感谢
7楼2011-12-14 11:58:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwx_ky

木虫 (小有名气)

【答案】应助回帖

★ ★
贺仪(金币+2): 多谢你的热心帮助 2011-12-14 14:26:55
引用回帖:
6楼: Originally posted by kattygirl at 2011-12-14 10:22:52:
我用的也是MS,第一种方式是:选中你要加电荷的原子修改properties中的charges。第二种方式:工具栏Modify——charges——set all charges/set total charge,视情况选择。

首先第一种方法我试过,不可行,算出来和没设置的一样,第二种方法我也试过,也不可行,和第一种方法算出来一样吧,不放你再试试没有设置时的能量,可定和设置了的一样。我试过,就这个问题我用了很多办法,都不行。
梦洁
8楼2011-12-14 14:10:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangth08

金虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
sunyang1988(金币+1): 谢谢交流 2011-12-15 04:03:15
MS里面加电荷好像是赋予整个系统的,不能赋予单个原子,虽然是那么加的,最后结果好像是整个系统的,我同学计算的经验。
9楼2011-12-14 20:44:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwx_ky

木虫 (小有名气)

【答案】应助回帖

liliangfang(金币+1): 多谢交流 2012-02-21 09:15:20
引用回帖:
9楼: Originally posted by wangth08 at 2011-12-14 20:44:28:
MS里面加电荷好像是赋予整个系统的,不能赋予单个原子,虽然是那么加的,最后结果好像是整个系统的,我同学计算的经验。

恩,我也算试了试,却是是整个体系的,那你同学有没有什么办法才能算离子吗?可以问问嘛,我这个问题研究了很久了,非常感谢。
梦洁
10楼2011-12-15 08:52:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kattygirl 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 5/250 2026-08-18 11:25 by 無關想念
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +43 医学老男孩 2026-08-13 97/4850 2026-08-18 11:07 by 医学老男孩
[基金申请] 时间戳变了,能看出什么问题? +5 基诺咪客 2026-08-17 6/300 2026-08-18 11:06 by hunter无悔
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +4 Tide man 2026-08-14 4/200 2026-08-18 10:57 by marinewhale
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
信息提示
请填处理意见