| 查看: 9270 | 回复: 42 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
fentuo1203新虫 (初入文坛)
|
[交流]
MS吸附能计算的问题
|
|||
|
各位大神，本人用MS的CASTEP模块做吸附能计算，大概过程是： 1、根据FINDIT的数据构建的晶胞模型，没有做结构优化，因为优化后的参数跟实验值相差太大，所以没做优化 2、切面，选取的是（003）面，没有建立超晶胞，第一层和最低层原子是一样的。随后对该Slab进行了结构优化。计算了单点能。 3、计算了被吸附小分子的能量，直接建立分子结构优化后进行的单点能计算。 4、将优化后的小分子加到第二步优化后的Slab上，固定下面一半，上面一半弛豫。 定义吸附能=系统能量-之前单个的能量，根据这个计算得到的吸附能是正值，相当于该反应是吸热反应，这应该是不可能的。 备注：我计算的是NiCl2晶体表面的吸附。 请各位大神不吝赐教，本人感激不尽！ |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有154人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵 |
大神您好，我也是做吸附的。我现在是研究几种溶剂在硫酸钡和硫酸锶的晶体表面的物理吸附和化学吸附作用，用第一性原理和分子动力学解释这几种溶剂对硫酸钡和硫酸锶晶体的晶面生长的影响。现在遇到的问题是，我看了类似文献里用forcite和morphology这两个模块里面的力场，其中博士论文里用compass算的。然而钡和锶并不在compass的力场参数文件包含的元素里，是自己编译力场文件么。硕士论文里用的 universal这个力场。但我感觉universal力场的话是不是相比compass力场准确性差不少。但compass力场也没有钡和锶单独的力场参数文件。请问这种情况，怎么处理呢。谢谢您。 |
43楼2021-03-18 18:17:13
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