计算化学公社

 找回密码 Forget password
 注册 Register

用CP2K做过渡金属001表面HER台阶图时遇到的数值问题

查看数: 3073 | 评论数: 6 | 收藏 Add to favorites 4
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-2-15 10:26

正文摘要:

得到的结果和文献报道的结果有非常大的不同,经过人为粗暴的数值校正后,看起来比较类似;是否计算设置相关的精度需要进行进一步的提升? 体系是Ag        Au        Cu&n ...

回复 Reply

jiangning198511 发表于 Post on 2023-2-15 16:50:05
Eudaimonia 发表于 2023-2-15 16:33
1.这个不是很确定频率和k点之间的关系了,个人感觉很容易带来误差

2.有很多文献报道过,但是一时间可 ...

1.这个不是很确定频率和k点之间的关系了,个人感觉很容易带来误差
因为频率计算不能使用K点,你可以都用超胞计算以保持计算参数稳定性,用一个体系分析对吸附能的影响

2.有很多文献报道过,但是一时间可能不太好找
如果是发表文章的话,还是要找一下,不然会有审稿人质疑计算合理性

3.纵轴的获得方式有好几种,a.实验值 b.经由计算结果结合动力学进行模拟 c.使用能斯特方程进行校正,此处使用的是方法c,只考察趋势的话应该是可以对比的,可以参考本文的相关结论 10.1039/c1cp22271f
HER的文章很多直接给出过电位的对比图来分析不同金属的催化活性

Eudaimonia 发表于 Post on 2023-2-15 16:38:32
kubingxia23 发表于 2023-2-15 15:40
为什么是 “固定与H原子距离>3埃的所有原子” 呢,看输入文件似乎是还放开了表面的两个金属原子,可能这里 ...

此处是根据 该帖子 来设置的,主要是cp2k计算频率如果全放开计算量太大实在承受不了所以才如此进行的
Eudaimonia 发表于 Post on 2023-2-15 16:33:36
jiangning198511 发表于 2023-2-15 11:28
1. 几何优化用的K点是331,频率计算采用的是伽马点计算,这个是否会带来误差
2. 这些金属的吸附能数据是否 ...

1.这个不是很确定频率和k点之间的关系了,个人感觉很容易带来误差

2.有很多文献报道过,但是一时间可能不太好找

3.纵轴的获得方式有好几种,a.实验值 b.经由计算结果结合动力学进行模拟 c.使用能斯特方程进行校正,此处使用的是方法c,只考察趋势的话应该是可以对比的,可以参考本文的相关结论 10.1039/c1cp22271f
jiangning198511 发表于 Post on 2023-2-15 16:06:19
kubingxia23 发表于 2023-2-15 15:40
为什么是 “固定与H原子距离>3埃的所有原子” 呢,看输入文件似乎是还放开了表面的两个金属原子,可能这里 ...

一般好像是固定slab,只放开吸附分子
kubingxia23 发表于 Post on 2023-2-15 15:40:23
为什么是 “固定与H原子距离>3埃的所有原子” 呢,看输入文件似乎是还放开了表面的两个金属原子,可能这里导致你的能量校正有点问题

全部选择(001)面有什么理由吗,建议先做调整使Pt(111)的计算结果和文献相符

jiangning198511 发表于 Post on 2023-2-15 11:28:55
本帖最后由 jiangning198511 于 2023-2-15 11:42 编辑

1. 几何优化用的K点是331,频率计算采用的是伽马点计算,这个是否会带来误差
2. 这些金属的吸附能数据是否有文献报道过?
2. 文献里面Y轴用的是电流,你用的是电位,这两个可以直接对比吗?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-20 23:59 , Processed in 0.175088 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list