| 查看: 4562 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
slz0811木虫 (小有名气)
|
[求助]
小分子CH4吸附在金属Co100表面的参数问题
|
||
|
刚开始用vasp计算，先前进行了一些初步的计算，可是在算吸附能时，却得出负值。与以往castep计算结果，文献值（正的）不符。 我算的是小分子CH4在金属Co100表面的吸附解离。 INCAR: SYSTEM = name ENCUT = 400 PREC = Medium GGA = 91 VOSKWN = 1 ISYM = 0 ISPIN = 2 ISTART = 0 ICHARG = 2 ISMEAR = 2 (金属) SIGMA = 0.1 IBRION = 2 ISIF = 2 NSW = 1000 POTIM = 0.05 NELMIN = 5 NELM = 200 LWAVE = .FALSE. LCHARG = .FALSE. KPOINTS： surface 0 M 5 5 2 （castep计算的fine精度） 0 0 0 atom or molcular 1 Rec 0 0 0 1 POSCAR： CH4 1.0 10.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 10.000000000000000 0.000000000000000 0.000000000000000 0.000000000000000 10.000000000000000 4 1 D -0.4604751244820571 -0.5955083450831910 0.3572567372930910 -0.6097261178601441 -0.5191680930578321 0.4382028422581180 -0.4459877478070221 -0.4307015616245711 0.4426085768177100 -0.4708332035946150 -0.5855807188904930 0.5428642882608110 -0.4967554857621820 -0.5327396496507220 0.4452330771902660 我的问题是: （1）在优化小分子时，参数是否应和优化金属时一样？ 比如ISMEAR = 2 (金属) ISMEAR = 0 （小分子）是否这样？ （2）k点的选取，我以前用的是castep，一般fine精度就可以满足。是否可以这样选取。 （3）计算结果如图：CH4吸附为物理吸附，吸附能应很小，20kJ/mol左右。  大家看看是哪儿出问题了？ |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理 |
» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- MS金属氧化物表面吸附有机物平衡问题
已经有4人回复
- DFT计算和X射线发射光谱研究金属表面小分子吸附电子结构
已经有26人回复
- 【分享】有机分子结构查询，表面吸附小分子构型不用愁
已经有15人回复
- CH4在Ir（111）表面的解离吸附过渡态计算：Henkelman的PRL文章
已经有10人回复
- 【求助】关于金属表面吸附
已经有10人回复
- 【求助完毕】CO极性分子在金属氧化物表面吸附为什么是O上C下的构型？
已经有6人回复
- 【求助】重金属吸附问题
已经有17人回复
- 【讨论】原子/分子在金属表面的吸附过程中外层轨道初始宽化过程
已经有3人回复
- 【讨论】co做探针分子吸附问题，请懂该方面的进来解答下
已经有11人回复
- 【求助】CO化学吸附法测定贵金属催化剂分散度的问题
已经有10人回复
kmw.8668
金虫 (正式写手)
归期守望者
- 应助: 9 (幼儿园)
- 金币: 1731.6
- 散金: 217
- 红花: 2
- 帖子: 668
- 在线: 235.9小时
- 虫号: 1008177
- 注册: 2010-04-29
- 专业: 凝聚态物性 II ：电子结构
6楼2011-05-30 08:33:32
slz0811
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1967.7
- 散金: 140
- 红花: 2
- 帖子: 103
- 在线: 84小时
- 虫号: 761947
- 注册: 2009-05-03
- 性别: GG
- 专业: 理论和计算化学
2楼2011-05-29 11:18:31
kmw.8668
金虫 (正式写手)
归期守望者
- 应助: 9 (幼儿园)
- 金币: 1731.6
- 散金: 217
- 红花: 2
- 帖子: 668
- 在线: 235.9小时
- 虫号: 1008177
- 注册: 2010-04-29
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
franch(金币+1): 谢谢回帖交流 2011-05-29 21:22:36
slz0811(金币+3): 这个我也差不多看过，就是理解不太透彻。 2011-05-29 21:31:29
franch(金币+1): 谢谢回帖交流 2011-05-29 21:22:36
slz0811(金币+3): 这个我也差不多看过，就是理解不太透彻。 2011-05-29 21:31:29
|
进行任何的静态计算或态密度计算，且K点数目（从IBZKPT文件中读取）大于4时，取ISMEAR＝－5；当由于原胞较大而K点数目较少（小于4个）时，取ISMEAR＝0，并设置一个合适的SIGMA值。另外对半导体或绝缘体的计算（不论是静态还是结构优化），取ISMEAR＝－5;当体系呈现金属性时，取ISMEAR＝1和2，以及设置一个合适的SIGMA值，这是因为取ISMEAR＝-5这种方法没有对电子的部分占据进行变分计算，所以对金属的力的计算误差较大。在进行能带结构计算时，ISMEAR 和SIGMA用默认值（ISMEAR＝1。SIGMA＝0.2）就好。 我一般优化计算时，ISMEAR采用的都是默认值，我觉的不同的ISMEAR值算出来的结果差距不是很大，你可以比较一下的 |
3楼2011-05-29 20:54:01
kmw.8668
金虫 (正式写手)
归期守望者
- 应助: 9 (幼儿园)
- 金币: 1731.6
- 散金: 217
- 红花: 2
- 帖子: 668
- 在线: 235.9小时
- 虫号: 1008177
- 注册: 2010-04-29
- 专业: 凝聚态物性 II ：电子结构
4楼2011-05-29 20:59:17
