| 查看: 9265 | 回复: 42 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
fentuo1203新虫 (初入文坛)
|
[交流]
MS吸附能计算的问题
|
|||
|
各位大神，本人用MS的CASTEP模块做吸附能计算，大概过程是： 1、根据FINDIT的数据构建的晶胞模型，没有做结构优化，因为优化后的参数跟实验值相差太大，所以没做优化 2、切面，选取的是（003）面，没有建立超晶胞，第一层和最低层原子是一样的。随后对该Slab进行了结构优化。计算了单点能。 3、计算了被吸附小分子的能量，直接建立分子结构优化后进行的单点能计算。 4、将优化后的小分子加到第二步优化后的Slab上，固定下面一半，上面一半弛豫。 定义吸附能=系统能量-之前单个的能量，根据这个计算得到的吸附能是正值，相当于该反应是吸热反应，这应该是不可能的。 备注：我计算的是NiCl2晶体表面的吸附。 请各位大神不吝赐教，本人感激不尽！ |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有189人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
fentuo1203
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 235.5
- 散金: 16
- 帖子: 44
- 在线: 18小时
- 虫号: 2832932
- 注册: 2013-11-28
- 专业: 无机非金属类电介质与电解
|
版主您好，关于NiCl2晶体结构确实只有那一篇。最近计算用的LDA函数，大概能算到2.89%和3.99%。 我用了一下Discover模块优化，发现我如果按照添加原子的方法建的晶胞会有很多的键，但实际上NiCl2应该是离子晶体不会出现共价键的呀，所以Discover总是无法进行优化，之后我把键都给删掉，就可以进行了。不过优化后的参数相差实在太大。 我师兄给我发了个NiCl2的cif文件，直接导入到MS中，出来的晶胞模型中确实没有共价键的存在。如图。 我现在不敢肯定那一种才是正确的构图方法，还是说两种都是可以的，参数都是一样的。请版主示教！！ http://image.keyan.cc/data/bcs/2 ... _1421600117_374.jpg http://image.keyan.cc/data/bcs/2 ... _1421600124_232.jpg  3D Atomistic.jpg  NiCl2#14208.jpg |
20楼2015-01-19 00:57:16
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
- 金币: 68121.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