| 查看: 9448 | 回复: 42 | ||||
fentuo1203新虫 (初入文坛)
|
[交流]
MS吸附能计算的问题
|
|||
|
各位大神，本人用MS的CASTEP模块做吸附能计算，大概过程是： 1、根据FINDIT的数据构建的晶胞模型，没有做结构优化，因为优化后的参数跟实验值相差太大，所以没做优化 2、切面，选取的是（003）面，没有建立超晶胞，第一层和最低层原子是一样的。随后对该Slab进行了结构优化。计算了单点能。 3、计算了被吸附小分子的能量，直接建立分子结构优化后进行的单点能计算。 4、将优化后的小分子加到第二步优化后的Slab上，固定下面一半，上面一半弛豫。 定义吸附能=系统能量-之前单个的能量，根据这个计算得到的吸附能是正值，相当于该反应是吸热反应，这应该是不可能的。 备注：我计算的是NiCl2晶体表面的吸附。 请各位大神不吝赐教，本人感激不尽！ |
回复此楼» 收录本帖的淘帖专辑推荐
 掺杂计算 |
» 猜你喜欢
帮我的英语口语老师找学生
已经有4人回复
-
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
物理学I论文润色/翻译怎么收费?
已经有225人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有4人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 如何用MS计算分子晶体的晶格能
已经有5人回复
- 怎样用MS模拟分子在有机聚合物表面的吸附过程。
已经有3人回复
- MS中sorption模块得到吸附量的单位问题
已经有11人回复
- GCMC模拟气体吸附时各种吸附量该怎么求？
已经有5人回复
- MS 吸附能的计算
已经有3人回复
- MS吸附能如何取平均值
已经有6人回复
- MS计算光学吸收谱
已经有4人回复
- ms中吸附能的计算
已经有40人回复
- 用MS算的吸附能大概在什么数量级？
已经有5人回复
- 【求助】MS怎么计算吸附能
已经有12人回复
- 【求助】脱附能计算问题
已经有8人回复
- 【求助】MS模拟5A介孔微孔吸附CO2时吸附热的问题
已经有14人回复
- 【求助】ms怎么计算键能
已经有14人回复
- 【求助】求使用MS中sorption模块用于分子筛吸附小分子的问题
已经有5人回复
- 【求助成功】MS模拟沸石的吸附性能时的一些问题
已经有9人回复
RicKillMercy
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1541.1
- 散金: 709
- 红花: 1
- 帖子: 615
- 在线: 179.2小时
- 虫号: 5415362
- 注册: 2016-12-28
- 性别: GG
- 专业: 流体力学
41楼2017-04-05 10:46:52
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
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
- 金币: 68127.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
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
5楼2015-01-12 13:30:33
Jhone
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 333.6
- 帖子: 88
- 在线: 22.2小时
- 虫号: 3524968
- 注册: 2014-11-07
- 性别: MM
- 专业: 可再生与替代能源利用中的
6楼2015-01-12 15:30:33
fentuo1203
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 235.5
- 散金: 16
- 帖子: 44
- 在线: 18小时
- 虫号: 2832932
- 注册: 2013-11-28
- 专业: 无机非金属类电介质与电解
|
版主您好，按照您的方法增大了E-cutoff和k-points的精度，得到了一个结果如下图http://image.keyan.cc/data/bcs/2 ... _1421075280_255.jpg 该计算结果收敛，但是怎么看都觉得有点不靠谱的样子。 还有一点就是，我无意间发现优化晶胞的时候没有勾选optimize cell，计算很快就收敛了，大概在三四步就收敛算完了，试了好多次。不知道这样的算法正确不正确。实在是学识尚浅，恳请版主不吝赐教，谢谢！  nicl2 Energies.jpg |
7楼2015-01-12 23:12:40
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 计算强帖: 8
- 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
8楼2015-01-13 08:32:37
fentuo1203
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 235.5
- 散金: 16
- 帖子: 44
- 在线: 18小时
- 虫号: 2832932
- 注册: 2013-11-28
- 专业: 无机非金属类电介质与电解
|
版主您好，还是麻烦您一下。在MS里构建的晶胞原胞参数与文献上的实验值是一致的。这两天按照您的方法，我用GGA和LDA的交换关联函数都试过了，而且E cut-off也从小到大设置了一遍，都提高到了750eV，同时对k-point各个档位也一并试过了一遍。结果如下： 1、用GGA计算，不管是提高E cut-off和k-point精度，得到的结果趋势都差不多，能量都是上下波动，不能够收敛，大致如下图：http://image.keyan.cc/data/bcs/2 ... _1421258135_232.jpg http://image.keyan.cc/data/bcs/2 ... _1421258154_210.jpg http://image.keyan.cc/data/bcs/2 ... _1421258185_807.jpg http://image.keyan.cc/data/bcs/2 ... _1421258204_653.jpg 2、用LDA计算，结果好一点结果都收敛了，测试了一下，E cut-off的值大概是400eV是比较合适的，但是不管我怎么调节k-point精度和E cut-off的值都没有办法让优化得到的晶格参数与实验值误差达到2%以内，至多a参数在3%，c参数在4.6%。是在是无法再接近了。http://image.keyan.cc/data/bcs/2 ... _1421258243_597.jpg 我还把use LDA+U这个选项勾上了，因为看有人说含镍的体系属磁性体系。不过还是不能把误差拉近。 已经没有办法了，请版主再指点一下，还有没有其他可调节的。  nicl2 Energies.jpg  GGA 400ev.jpg  GGA 500ev k-point corse.jpg  GGA 500ev k-point fine.jpg  LDA 400ev kpoint corse.jpg |
9楼2015-01-15 01:57:54
阿Q~~
至尊木虫 (文坛精英)
- 应助: 77 (初中生)
- 金币: 33343.6
- 红花: 107
- 帖子: 36186
- 在线: 1028.1小时
- 虫号: 3237212
- 注册: 2014-05-27
- 性别: MM
- 专业: 宇宙学
10楼2015-01-15 08:44:24
