| 查看: 4118 | 回复: 27 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
计算Cl-在Fe表面吸附问题时体系加电荷方法
|
|||
|
计算Cl-离子在Fe表面吸附时，我的算法是对cl原子加一个负电荷后，对整个体系赋予-1价，几何优化、能量计算后、结果计算出的吸附能为正值。 想请大侠们帮忙解决一下，这是哪方面的问题，会导致最后的吸附能为正值。 |
回复此楼» 猜你喜欢
青年基金C终止
已经有4人回复
-
青椒八年已不青，大家都被折磨成啥样了？
已经有7人回复
-
为什么nbs上溴 没有产物点出现呢
已经有10人回复
-
救命帖
已经有11人回复
-
招博士
已经有5人回复
-
26申博求博导推荐-遥感图像处理方向
已经有4人回复
-
限项规定
已经有7人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有5人回复
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
uuv2010(金币+2): 多谢指导！ 2011-12-15 23:49:57
感谢参与，应助指数 +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
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2011-12-09 02:04:33
3楼2011-12-10 21:10:23
acridine
木虫 (正式写手)
- 1ST强帖: 11
- 应助: 13 (小学生)
- 金币: 4901.6
- 红花: 13
- 帖子: 706
- 在线: 389.5小时
- 虫号: 305478
- 注册: 2006-12-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★
sunyang1988(金币+3): 谢谢指点 2011-12-11 09:26:08
sunyang1988(金币+3): 谢谢指点 2011-12-11 09:26:08
|
第一个，关于公式，最好能有相关的参考文献对比。因为比较带电体系能量的时候需要进行能级的alignment，以使得能量具有可比性。这个过程在计算带电缺陷形成能的时候经常出现，不知道研究带电吸附的时候有没有类似的过程，找下相关文章吧。 第二个问题。首先，因为有些时侯离子是一种比较化学的、定性的描述，而实际体系中并不一定是这样的带电原子形式存在的。所以要先明确你模拟的体系是否真的是这种带负电的Cl原子参与。 另外就是模型的意义和合理性。我们想要模拟的是无限固体Fe表面的吸附，但由于计算的限制，我们只能建立吸附位置附近的结构模型。如果此时如果你让模型带电，那么其实模拟的是：电荷局域在吸附位置附近的情况。而实际上，即使真的有额外电荷引入，可能也不会局限在吸附位置附近，尤其这里的Fe还是导体。那么你的模型的合理性就有待商榷了。 以上只是我的看法了，我也没有具体研究过吸附，只当做建议就好了。  |
4楼2011-12-11 00:48:00