计算化学公社

标题: vaspkit进行热力学校正，校正完能垒变成负的咋整？ [打印本页]

作者
Author: 我是大帅逼 时间: 2021-3-6 17:02
标题: vaspkit进行热力学校正，校正完能垒变成负的咋整？
如题。就是很简单的算了个Ni（111）上CO2的分解，DFT数据能垒0.7左右吧，然后始态，过渡态算了频率，用vaspkit501功能进行自由能校正，矫正完成负的了！！？？？？？哈哈哈哈哈哈哈哈哈哈哈哈
这也太搞笑了吧
这应该是不对的吧，反正文献里自由能垒是正的

作者
Author: 卡开发发 时间: 2021-3-9 09:14
有可能的原因是不合理的固定导致的，尤其是分子和表面有显著作用的情况，我不认为约束整个基底只计算吸附原子的自由度对Hessian的贡献是合理的。

作者
Author: 我是大帅逼 时间: 2021-3-9 16:11

卡开发发发表于 2021-3-9 09:14
有可能的原因是不合理的固定导致的，尤其是分子和表面有显著作用的情况，我不认为约束整个基底只计算吸附原 ...

谢谢，我后来又做了几个反应发现化学吸附的反应过程经过校正得到的数据是合理的，比如Ni（111）上C+O=CO.而涉及物理吸附的过程经过校正得到的数据就很不合理，像前面说的CO2=CO+O(开始时CO2在表面的吸附是物理吸附),还有后来我又试了CH4的分解，也是这样，现在我猜想是跟物理吸附有关。不知道您有没有类似的经验，涉及物理吸附的反应过程使用vaspkit进行校正结果是否准确。

作者
Author: 卡开发发 时间: 2021-3-9 16:16

我是大帅逼发表于 2021-3-9 16:11
谢谢，我后来又做了几个反应发现化学吸附的反应过程经过校正得到的数据是合理的，比如Ni（111）上C+O=CO. ...

我不能确定是否是vaspkit本身的问题，你可以试试phonopy，如果也给出不合理的结果，那么还是声子计算部分的操作有问题。

作者
Author: 我是大帅逼 时间: 2021-3-9 16:21

卡开发发发表于 2021-3-9 16:16
我不能确定是否是vaspkit本身的问题，你可以试试phonopy，如果也给出不合理的结果，那么还是声子计算部分 ...

好的，多谢

欢迎光临计算化学公社 (http://bbs.keinsci.com/)