计算化学公社

标题: 求助cp2k计算此反应OOH*+H→O*+H2O 吉布斯，这样算是否正确？ [打印本页]

作者
Author: qzm 时间: 2024-5-1 21:11
标题: 求助cp2k计算此反应OOH*+H→O*+H2O 吉布斯，这样算是否正确？
本帖最后由 qzm 于 2024-5-1 22:35 编辑

1式子： *+O2+H→OOH*
2式子：OOH*+H→O*+H2O
吉布斯自由能
步骤1：先计算O2分子能量 -31.881592786405143  （把优化好的一个O2分子放在一个盒子里，cp2k计算获得能量）
         先计算H2分子能量 -1.1613144308840  （把优化好的一个H2分子放在一个盒子里，cp2k计算获得能量）
      先计算H2O分子能量-17.219987932384019  （把优化好的H2O分子放在一个盒子里，cp2k计算获得能量）
      计算频率  求出他们的 G(O2)   G(H2)   G(H2O)
   不确定是否准确，
步骤2：计算OOH*的能量  *的能量，O*的能量分别计算他们的频率求出 G（OOH*） G（*） G（OH*）

计算  1式子吉布斯自由能= G（OOH*） - G（*）- G(O2) - 1/2G(H2)

      2式子吉布斯自由能= G（O*）+ G (H2O) - G（OOH*） - 1/2G(H2)

求助：1吉布斯自由能这样计算是否正确？
      2给出了O2 H2 H2O 计算的输入和输出文件，这样计算获得的结果是否可以写论文直接用

作者
Author: sobereva 时间: 2024-5-1 22:23
没具体OOH*的输入文件没法说
O2基态是三重态
孤立体系还用k点莫名其妙
孤立体系用非周期性明显是更好选择
注意对于算能量目的，DZVP-MOLOPT-SR-GTH太烂，这只适合opt freq目的

作者
Author: qzm 时间: 2024-5-1 22:38
本帖最后由 qzm 于 2024-5-1 22:40 编辑

sobereva 发表于 2024-5-1 22:23
没具体OOH*的输入文件没法说
O2基态是三重态
孤立体系还用k点莫名其妙

感谢社长，更新了OOH的输入，计算能量选TZ2P-MOLOPT 可以吗，请问计算步骤是否存在错误呢？

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