| 查看: 9271 | 回复: 42 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
fentuo1203新虫 (初入文坛)
|
[交流]
MS吸附能计算的问题
|
||||
|
各位大神，本人用MS的CASTEP模块做吸附能计算，大概过程是： 1、根据FINDIT的数据构建的晶胞模型，没有做结构优化，因为优化后的参数跟实验值相差太大，所以没做优化 2、切面，选取的是（003）面，没有建立超晶胞，第一层和最低层原子是一样的。随后对该Slab进行了结构优化。计算了单点能。 3、计算了被吸附小分子的能量，直接建立分子结构优化后进行的单点能计算。 4、将优化后的小分子加到第二步优化后的Slab上，固定下面一半，上面一半弛豫。 定义吸附能=系统能量-之前单个的能量，根据这个计算得到的吸附能是正值，相当于该反应是吸热反应，这应该是不可能的。 备注：我计算的是NiCl2晶体表面的吸附。 请各位大神不吝赐教，本人感激不尽！ |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有298人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
1. 然后回到该贴的主题，还是计算NiCl2，优化完晶胞，我就切面优化表面，但是为什么优化表面的时候，表面SLAB的空间群会发生改变，比如我优化（003）面，切出来的时候是P1，优化完之后却是PM1（a=b），还有（101）面，原来是P1，优化之后都不知道是什么空间群了（a=b,alpha=beta）请问一下版主，这样正常吗。如果不正常应该怎么改正，如何设置计算参数呢。我的计算参数都是按照优化晶胞的参数设置的，唯一的就是没有勾选optimize cell。 回复：按我的理解，切出晶面，再加真空层得到的结构，与原来的晶体一般不会有相同的空间群，一般都是P1。NiCl2晶胞没问题，用MS切出来的晶面就应该没问题。 2. 切面的时候为了保证NiCl2个数是晶胞的整数倍，上下表面原子并不一致，优化时固定了下层原子。 回复：通常情况下，切出来的晶面与原晶面能保持整数的比例关系。 “上下表面原子并不一致”，这里的上下是什么意思？ |
27楼2015-01-21 10:19:42
DRP7
新虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 537.9
- 散金: 2
- 帖子: 116
- 在线: 14.6小时
- 虫号: 3443558
- 注册: 2014-09-26
- 性别: GG
- 专业: 自然语言理解与机器翻译
2楼2015-01-11 23:44:24
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68123.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
fentuo1203: 金币+10 2015-01-12 10:41:18
小木虫: 金币+0.5, 给个红包，谢谢回帖
fentuo1203: 金币+10 2015-01-12 10:41:18
|
第一步存在问题，模拟计算一个很关键的问题就是：你如何证明，你选用的计算参数设置是合理的？如果不能解决这个问题，后续的吸附能计算结果的准确性也会受质疑。 一般说来，为了保证计算的可信，会采用做几何优化时的参数设置来做后续的计算。所以建议，几何优化这一步时调整计算参数（实际上这一块要做的工作很多，比如交换能相关函数、赝势、E-cutoff、k-points等参数的设置，这些都是CASTEP计算的最基础步骤），使得计算值和实验值相差<2%。 |
3楼2015-01-12 08:59:00
fentuo1203
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 235.5
- 散金: 16
- 帖子: 44
- 在线: 18小时
- 虫号: 2832932
- 注册: 2013-11-28
- 专业: 无机非金属类电介质与电解
|
版主您好，非常感谢您的指点，我对这些参数的设置也进行了一些尝试，但是得到的结果都有一个倾向，就是能量值上下波动，没有规律，所以也不收敛，我想知道产生这样的现象有哪些特定的原因。感谢！下面是其中一个优化的能量变化图。http://image.keyan.cc/data/bcs/2 ... _1421038685_320.jpg  NiCl2 cell Energies.jpg |
4楼2015-01-12 12:58:27