24小时热门版块排行榜    

查看: 4181  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

[求助] CH4在Ir(111)表面的解离吸附过渡态计算:Henkelman的PRL文章

用vtstool的人大概都很熟悉Graeme Henkelman了,这个是他2001年的文章。看看他是怎么计算和分析过渡态来发PRL的。

我有不少问题没有看懂,希望路人指教!

1.作者把CH4放在表面外4埃左右的高度作为起点,解离成H和CH3的表面吸附态作为终点,插入若干image做NEB计算,然后在NEB的基础上做dimer计算找鞍点。
请问,为什么不直接用CNEB方法找鞍点呢?如果先用NEB方法,一般收敛到多高的精度(EDIFFG)再来做dimer计算?

2.作者通过dimer找到了鞍点,那么如何确定MEP呢?直接把dimer方法找到的鞍点放到前面收敛的NEB路径中?

3.如下图,作者做了零点能修正和dispersion interaction修正。零点能修正为什么降低总能呢?long range dispersioni interaction是什么意思?如何修正?

4.这篇文章给出的反应路径,非常明确的给出了CH4先有物理吸附,再越过能垒,转变为化学吸附,然后解离的路径。物理吸附的能量很小,这种若作用vasp的结果是否可靠呢?下图粗黑线表明物理吸附能阱的出现,依赖于ZPE修正,如果不做修正的话是不会出现这样的结果的,图中细实线就是直接把总能列出的结果,从中丝毫没有看到物理吸附的迹象。那么ZPE对气相分子和固体材料都有图示所列的这么大的影响?我觉得只对质量小,自由度大的小分子才有大的影响的,而对金属氧化物固体材料的影响可以忽略。不知道为什么这里给出的结果,ZPE修正的粗实线整体低于细实线如此明显?

5. MEP的ZPE修正是如何操作的呢?E+ZPE(=0.5hv)?如果是这样的话,为什么下图ZPE的粗实线在细实线的下面?难道仅仅是整体平移的“幻相”。问题见另外一个帖子:http://muchong.com/bbs/viewthread.php?tid=3507528&fpage=1


[ Last edited by cenwanglai on 2011-8-21 at 15:40 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理 我的饕餮 纳米及第一性原理方向 他山之玉
计算原理

» 本帖@通知

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)

★ ★ ★
zzy870720z(金币+3): 谢谢指点 2011-09-05 08:36:27
1.一阶鞍点有明确物理意义,高阶鞍点没有明确物理意义;
2.先NEB再DIMER,我想是因为NEB收敛后能量高的的IMAGE不一定就是TS,在此基础上取能量高的两个image,然后再dimer更容易收敛到TS;
3.找TS,CNEB是个很好的选择,命中率是很高的。CNEB后,找能量最高的image,在做一下频率计算以验证是否TS;
4.要确定能垒,需先确定TS,即一阶鞍点。
4楼2011-08-20 19:03:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)


zzy870720z(金币+1): 谢谢补充 2011-09-05 08:36:39
另外,对于包含H的体系,zpe校正是很明显的
5楼2011-08-20 19:05:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hakuna

木虫 (知名作家)


zzy870720z(金币+1): 谢谢提示 2011-09-11 22:56:07
引用回帖:
10楼: Originally posted by piaoxue001 at 2011-09-10 11:00:33:
您好,请问如何做零点能修正?

http://theory.cm.utexas.edu/henkelman/pubs/henkelman06_044706.pdf
这篇文章里有详细介绍
11楼2011-09-10 17:10:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cenwanglai 的主题更新
信息提示
请填处理意见